Graphical Verification of Smoothness and Deployment

>

> IVP := subs( k=kk, param, EOM[quad] union IC );

>

> SOLN := motion( IVP, t );

>

> scale[X] := { x=1 }:

> scale[XV] := { x=1000, v=10 }:

> scale[VAJ] := { v=10, a=981/100, j=981/100 }:

> scale[XVA] := { x=1000, v=10, a=981/100 }:

> scale[VA] := { v=10, a=981/100 }:

> scale[all] := { x=1000, v=10, a=981/100, j=981/100 }:

> OPTS := style=POINT, symbol=CIRCLE, numpoints=200:

>

> interface( plotdevice=default );

> PLOT1 := display(

> {

> plot_motion( SOLN, [0..30], scale[VA], numpoints=200, xtickmarks=3 ),

> textplot( [ [ 14, -2.0, `v/10` ],

> [ 14, 2.5, `a/g` ] ] )

> },

> title="Initial Motion for T-10" ):

> display( PLOT1 );

>

> PLOT2 := display(

> {

> plot_motion( SOLN, [9.7..13.5], scale[VAJ], OPTS, xtickmarks=3 ),

> textplot( [ [ 10.6, -5.5, `v/10` ],

> [ 11.6, 5.0, `a/g` ],

> [ 10.4, 6.2, `j/g` ] ] )

> },

> title="T-10 Canopy Deployment" ):

> display( PLOT2 );

>

> PLOT3 := display(

> {

> plot_motion( SOLN, [10.25..11.75], scale[VAJ], OPTS, xtickmarks=3 ),

> textplot( [ [ 10.7, -5.9, `v/10` ],

> [ 10.8, 3.8, `a/g` ],

> [ 10.6, 7.7, `j/g` ] ] )

> },

> title="Snatch and Opening Forces" ):

> display( PLOT3 );

>