add 1/2 GC

This commit is contained in:
cosmic_power 2023-10-31 23:10:39 +08:00
parent 89c0d9cd1e
commit 80d3ec831c
4 changed files with 29 additions and 18 deletions

View File

@ -5,7 +5,7 @@
0 # nlages
0 # n_add_eq
0 # n_roots
2.000000 # tend
1.130000 # tend
1 # integration method
9.9999997E-06 # atolg(1)
1.0000000E-07 # rtolg(1)
@ -112,17 +112,17 @@
0.0000000E+00 # vipar
14 # N-Values for IntegrationOutputInfos
20 # nrhs_total
19 # nrhs_nojac
132 # nrhs_total
131 # nrhs_nojac
1 # n jacobi
0 # n mass
19 # n steps
131 # n steps
0 # n roots
0 # n error test failors
0 # n convergence test failors
0 # n steps_rejected
0.6892850 # last step size
9.9999998E-03 # last step size
1 # last order
0 # rhs_error_flag
0 # integrator_error_flag
2.000000 # tout required
1.130000 # tout required

Binary file not shown.

Binary file not shown.

View File

@ -57,18 +57,18 @@ prop2d.proj.dir = 2
! Views
!**********************************************************************
view.type ( $V_WorkingView ) = 0 ! Projection type
view.pos ( 1 , $V_WorkingView ) = 1.1536542177200317E+00 ! Position
view.pos ( 2 , $V_WorkingView ) = -1.8443827629089355E+00 ! Position
view.pos ( 3 , $V_WorkingView ) = 6.6502451896667480E-01 ! Position
view.orient ( 1 , $V_WorkingView ) = 5.7451611757278442E-01 ! Orientation
view.orient ( 2 , $V_WorkingView ) = 3.0375838279724121E-02 ! Orientation
view.orient ( 3 , $V_WorkingView ) = 2.9174527153372765E-02 ! Orientation
view.orient ( 4 , $V_WorkingView ) = 8.1740903854370117E-01 ! Orientation
view.pos ( 1 , $V_WorkingView ) = -7.4012970924377441E-01 ! Position
view.pos ( 2 , $V_WorkingView ) = -6.2400698661804199E-01 ! Position
view.pos ( 3 , $V_WorkingView ) = 6.7095392942428589E-01 ! Position
view.orient ( 1 , $V_WorkingView ) = 4.5991823077201843E-01 ! Orientation
view.orient ( 2 , $V_WorkingView ) = -3.4318101406097412E-01 ! Orientation
view.orient ( 3 , $V_WorkingView ) = -4.4884416460990906E-01 ! Orientation
view.orient ( 4 , $V_WorkingView ) = 6.8501162528991699E-01 ! Orientation
view.angle ( $V_WorkingView ) = 7.8539818525314331E-01 ! Lens angle
view.fclip.auto ( $V_WorkingView ) = 1 ! Front clipping plane auto
view.bclip.auto ( $V_WorkingView ) = 1 ! Back clipping plane auto
view.fclip.value ( $V_WorkingView ) = 9.9999997764825821E-03 ! Front clipping plane value
view.bclip.value ( $V_WorkingView ) = 3.8589730262756348E+00 ! Back clipping plane value
view.bclip.value ( $V_WorkingView ) = 3.0534176826477051E+00 ! Back clipping plane value
view.rotcenter.type ( $V_WorkingView ) = 0 ! Rotation center type
view.rotcenter.adjust ( $V_WorkingView ) = 1 ! Adjust rotation center
view.motion.active ( $V_WorkingView ) = 0 ! Camera is moved
@ -253,7 +253,7 @@ prim.par ( 3 , $P_Isys ) = 1.00000000
!**********************************************************************
body.m ( $B_Body1 ) = 1.0000000000000000E+00 ! Mass of the Body
body.mp ( $B_Body1 ) = 0 ! 0=manual; 1=auto (based on geometry); 2=mass manual, CG & Inertia auto
body.cg.pos ( 1 , $B_Body1 ) = 0.0000000000000000E+00 ! Center of gravity
body.cg.pos ( 1 , $B_Body1 ) = { 1/2 } ! Center of gravity
body.cg.kind ( $B_Body1 ) = 1 ! Kind of CG specification: 0=wrt CG; 1=wrt Marker
body.cg.ref ( $B_Body1 ) = $M_Body1_BRF ! Reference Marker for center of gravity
body.I.tens ( 1 , 1 , $B_Body1 ) = 1.0000000000000000E+00 ! Moments of inertia
@ -280,6 +280,10 @@ marker.pos ( 1 , $M_Body1_Top ) = $G_DH1.$_a
marker.pos ( 2 , $M_Body1_Top ) = 0.0000000000000000E+00 ! Position
marker.ang ( 1 , $M_Body1_Top ) = 0.0000000000000000E+00 ! Angles
marker.type ( $M_Body1_GC ) = 2 ! Type
marker.parent ( $M_Body1_GC ) = $B_Body1 ! Body
marker.pos ( 1 , $M_Body1_GC ) = { $G_DH1.$_a/2 } ! Position
joint.from ( $J_Body1 ) = $M_Isys ! From Marker
joint.to ( $J_Body1 ) = $M_Body1_BRF ! To Marker
joint.type ( $J_Body1 ) = 29 ! Type
@ -353,7 +357,7 @@ prim.par ( 11 , $P_Body1_Link ) = 0.00000000
body.m ( $B_Body2 ) = 1.0000000000000000E+00 ! Mass of the Body
body.mp ( $B_Body2 ) = 0 ! 0=manual; 1=auto (based on geometry); 2=mass manual, CG & Inertia auto
body.cg.pos ( 1 , $B_Body2 ) = 0.0000000000000000E+00 ! Center of gravity
body.cg.pos ( 1 , $B_Body2 ) = { 1/2 } ! Center of gravity
body.cg.kind ( $B_Body2 ) = 1 ! Kind of CG specification: 0=wrt CG; 1=wrt Marker
body.cg.ref ( $B_Body2 ) = $M_Body2_BRF ! Reference Marker for center of gravity
body.I.tens ( 1 , 1 , $B_Body2 ) = 1.0000000000000000E+00 ! Moments of inertia
@ -376,6 +380,10 @@ marker.pos ( 1 , $M_Body2_Top ) = $G_DH2.$_a
marker.pos ( 2 , $M_Body2_Top ) = 0.0000000000000000E+00 ! Position
marker.ang ( 1 , $M_Body2_Top ) = 0.0000000000000000E+00 ! Angles
marker.type ( $M_Body2_GC ) = 2 ! Type
marker.parent ( $M_Body2_GC ) = $B_Body2 ! Body
marker.pos ( 1 , $M_Body2_GC ) = { $G_DH2.$_a/2 } ! Position
joint.from ( $J_Body2 ) = $M_Body1_Top ! From Marker
joint.to ( $J_Body2 ) = $M_Body2_BRF ! To Marker
joint.type ( $J_Body2 ) = 29 ! Type
@ -443,7 +451,10 @@ prim.par ( 11 , $P_Body2_Link ) = 0.00000000
!**********************************************************************
! Sensors
!**********************************************************************
sensor.from ( $S_end ) = $M_Body2_BRF ! From Marker
sensor.to ( $S_end ) = $M_Body1_BRF ! To Marker
sensor.from ( $S_1 ) = $M_Isys ! From Marker
sensor.to ( $S_1 ) = $M_Body1_GC ! To Marker
sensor.from ( $S_2 ) = $M_Isys ! From Marker
sensor.to ( $S_2 ) = $M_Body2_GC ! To Marker