*KEYWORD *database_extent_binary ,1 $ center of rotation at x,y,z = -8,0,0 $ $ add rigid body to end of blade *constrained_nodal_rigid_body 60 *set_node 60 16,32,48,64 $ ******************** begin DR input *database_binary_d3drlf 1 *load_body_rz 1,1,2, -8,0,0 *define_curve $ body load for transient phase 1 0,0 1000,0 *define_curve $ body load for DR phase 2,1 0,0 1e-5,300. 100,300. *control_dynamic_relaxation 20,1.e-7 $ ************************ end DR input *control_structured *control_accuracy 1,2 $ stiffness damping does not damp RB motion *damping_part_stiffness 1,0.10 *CONTROL_TIMESTEP $ .7 is ok for 10% damping $ , .5 *part_inertia rigid hub 2,1,2 -8,0,0, .6e-4 1,,,1,,1 0,0,0, 0,0,300 $*initial_velocity_generation $1,2, 300. $$-8.,0.,0., 0.,0.,1., 0 $$ phase=1 ONLY if dr is invoked (otherwise, velocities wont' take) $-8.,0.,0., 0.,0.,1., 1 *initial_velocity_generation 1,2, 300. $ phase=0 is also required so that V0 of rigid bodies is correct -8.,0.,0., 0.,0.,1., 0 *database_glstat 0.0001 *DATABASE_BINARY_D3PLOT 0.00018 *CONTROL_TERMINATION 0.009 *database_history_node 1,2,3 $*database_nodout $1.e-9 *TITLE blade spin-up *MAT_ELASTIC 1 3.00000-4 1.00000+7 0.3000000 *mat_rigid 2,3.e-4, 1.e7,.3 1,7,4 *SECTION_SHELL 1 16 0.0000000 0.0000000 0.0000000 0.0000000 0 $ 1 6 0.0000000 0.0000000 0.0000000 0.0000000 0 0.1000000 0.1000000 0.1000000 0.1000000 0.0000000 *control_hourglass 8 $*control_shell $,,,-2 $ *PART 1 1 1 0 0 0 0 0 *NODE 1-5.000000000E+00-1.000000000E+00 0.000000000E+00 0 0 2-4.333333492E+00-1.000000000E+00 0.000000000E+00 0 0 3-3.666666746E+00-1.000000000E+00 0.000000000E+00 0 0 4-3.000000000E+00-1.000000000E+00 0.000000000E+00 0 0 5-2.333333254E+00-1.000000000E+00 0.000000000E+00 0 0 6-1.666666627E+00-1.000000000E+00 0.000000000E+00 0 0 7-1.000000000E+00-1.000000000E+00 0.000000000E+00 0 0 8-3.333333433E-01-1.000000000E+00 0.000000000E+00 0 0 9 3.333333433E-01-1.000000000E+00 0.000000000E+00 0 0 10 1.000000000E+00-1.000000000E+00 0.000000000E+00 0 0 11 1.666666627E+00-1.000000000E+00 0.000000000E+00 0 0 12 2.333333254E+00-1.000000000E+00 0.000000000E+00 0 0 13 3.000000000E+00-1.000000000E+00 0.000000000E+00 0 0 14 3.666666746E+00-1.000000000E+00 0.000000000E+00 0 0 15 4.333333492E+00-1.000000000E+00 0.000000000E+00 0 0 16 5.000000000E+00-1.000000000E+00 0.000000000E+00 0 0 17-5.000000000E+00-3.333333433E-01 0.000000000E+00 0 0 18-4.333333492E+00-3.333333433E-01 0.000000000E+00 0 0 19-3.666666746E+00-3.333333433E-01 0.000000000E+00 0 0 20-3.000000000E+00-3.333333433E-01 0.000000000E+00 0 0 21-2.333333254E+00-3.333333433E-01 0.000000000E+00 0 0 22-1.666666627E+00-3.333333433E-01 0.000000000E+00 0 0 23-1.000000000E+00-3.333333433E-01 0.000000000E+00 0 0 24-3.333333433E-01-3.333333433E-01 0.000000000E+00 0 0 25 3.333333433E-01-3.333333433E-01 0.000000000E+00 0 0 26 1.000000000E+00-3.333333433E-01 0.000000000E+00 0 0 27 1.666666627E+00-3.333333433E-01 0.000000000E+00 0 0 28 2.333333254E+00-3.333333433E-01 0.000000000E+00 0 0 29 3.000000000E+00-3.333333433E-01 0.000000000E+00 0 0 30 3.666666746E+00-3.333333433E-01 0.000000000E+00 0 0 31 4.333333492E+00-3.333333433E-01 0.000000000E+00 0 0 32 5.000000000E+00-3.333333433E-01 0.000000000E+00 0 0 33-5.000000000E+00 3.333333433E-01 0.000000000E+00 0 0 34-4.333333492E+00 3.333333433E-01 0.000000000E+00 0 0 35-3.666666746E+00 3.333333433E-01 0.000000000E+00 0 0 36-3.000000000E+00 3.333333433E-01 0.000000000E+00 0 0 37-2.333333254E+00 3.333333433E-01 0.000000000E+00 0 0 38-1.666666627E+00 3.333333433E-01 0.000000000E+00 0 0 39-1.000000000E+00 3.333333433E-01 0.000000000E+00 0 0 40-3.333333433E-01 3.333333433E-01 0.000000000E+00 0 0 41 3.333333433E-01 3.333333433E-01 0.000000000E+00 0 0 42 1.000000000E+00 3.333333433E-01 0.000000000E+00 0 0 43 1.666666627E+00 3.333333433E-01 0.000000000E+00 0 0 44 2.333333254E+00 3.333333433E-01 0.000000000E+00 0 0 45 3.000000000E+00 3.333333433E-01 0.000000000E+00 0 0 46 3.666666746E+00 3.333333433E-01 0.000000000E+00 0 0 47 4.333333492E+00 3.333333433E-01 0.000000000E+00 0 0 48 5.000000000E+00 3.333333433E-01 0.000000000E+00 0 0 49-5.000000000E+00 1.000000000E+00 0.000000000E+00 0 0 50-4.333333492E+00 1.000000000E+00 0.000000000E+00 0 0 51-3.666666746E+00 1.000000000E+00 0.000000000E+00 0 0 52-3.000000000E+00 1.000000000E+00 0.000000000E+00 0 0 53-2.333333254E+00 1.000000000E+00 0.000000000E+00 0 0 54-1.666666627E+00 1.000000000E+00 0.000000000E+00 0 0 55-1.000000000E+00 1.000000000E+00 0.000000000E+00 0 0 56-3.333333433E-01 1.000000000E+00 0.000000000E+00 0 0 57 3.333333433E-01 1.000000000E+00 0.000000000E+00 0 0 58 1.000000000E+00 1.000000000E+00 0.000000000E+00 0 0 59 1.666666627E+00 1.000000000E+00 0.000000000E+00 0 0 60 2.333333254E+00 1.000000000E+00 0.000000000E+00 0 0 61 3.000000000E+00 1.000000000E+00 0.000000000E+00 0 0 62 3.666666746E+00 1.000000000E+00 0.000000000E+00 0 0 63 4.333333492E+00 1.000000000E+00 0.000000000E+00 0 0 64 5.000000000E+00 1.000000000E+00 0.000000000E+00 0 0 *ELEMENT_SHELL 1 2 1 17 18 2 2 1 2 18 19 3 3 1 3 19 20 4 4 1 4 20 21 5 5 1 5 21 22 6 6 1 6 22 23 7 7 1 7 23 24 8 8 1 8 24 25 9 9 1 9 25 26 10 10 1 10 26 27 11 11 1 11 27 28 12 12 1 12 28 29 13 13 1 13 29 30 14 14 1 14 30 31 15 15 1 15 31 32 16 16 2 17 33 34 18 17 1 18 34 35 19 18 1 19 35 36 20 19 1 20 36 37 21 20 1 21 37 38 22 21 1 22 38 39 23 22 1 23 39 40 24 23 1 24 40 41 25 24 1 25 41 42 26 25 1 26 42 43 27 26 1 27 43 44 28 27 1 28 44 45 29 28 1 29 45 46 30 29 1 30 46 47 31 30 1 31 47 48 32 31 2 33 49 50 34 32 1 34 50 51 35 33 1 35 51 52 36 34 1 36 52 53 37 35 1 37 53 54 38 36 1 38 54 55 39 37 1 39 55 56 40 38 1 40 56 57 41 39 1 41 57 58 42 40 1 42 58 59 43 41 1 43 59 60 44 42 1 44 60 61 45 43 1 45 61 62 46 44 1 46 62 63 47 45 1 47 63 64 48 *END