$ See email from Borrvall to jday dated 9/2/14 $ $ Analytic soln (when loaded only in z and erosion turned off): $ zdisp = pl^3/3EI = 8.88e-2 $ hand calc for longitudinal stress (sigx) at center of 1st element is Mc/I = +-73. $ (or 75.4 if moment is calculated at IP nearest support rather than center of element) $ sigz = 2.22 (same sign for all 8 IP) $ sigy= 0 $ $ While displacement is not bad, stresses are NOT what I expected!! $ *KEYWORD *comment *control_implicit_general 1,0.1 $___________________________________________ $ ignore weird stresses and proceed with mat_add_erosion example anyway $ erosion OK in dev; not working in R7.1.1. I see no mention of a fix in the SVN log. *comment *mat_add_erosion $ card 1, NUMFIP in 7th field 1,,,,,,1 $1,,,,,,4 $1,,,,,,8 $ card 2 (use eloutdet to plot ipt values) $$ max prin strain $,,,1.0e-4 $$ max tensorial shear strain = (max prin eps - min prin eps)/2 ,,,,0.05e-3 $___________________________________________ *control_hourglass 6,1.0 *section_solid $1,1 1,-1 $1,-2 $1,2 $$ Non-S/R 8-ipt solid (SPC out rotational DOF at fixed end) $$ - y-stress is ~zero (good); x-stress is a bit high; z-stress is much too low and varies element-to-element $1,3 *comment *constrained_global $7,,1, 0,0,0 6,,1, 0,0,0 *boundary_spc_node 1,,1,,1 12,,1,,1 23,,1 34,,1 $ store stress for all 8 IP *DATABASE_EXTENT_BINARY ,,,1 $ cmpflg $ nintsld 8,,,,,,ALL *PART 1,1,1 $1,1,22 *interface_springback_lsdyna 9 *set_part 9 1 *database_history_solid 1 *database_elout 3.e-3 *TITLE ELFORM -1 failure test *CONTROL_TERMINATION $ 0.3990000 0 0.0000000 0 0.0000000 3.0 *CONTROL_TIMESTEP 0.0000000 0.6700000 0 0.0000000 0.0000000 0 0 0 *CONTROL_ENERGY 1 2 1 1 *DATABASE_BINARY_D3PLOT 3.00000-2 0 *MAT_ELASTIC $ 1 7.00000-4 3.00000+5 0.3000000 1 7.00000-4 3.00000+5 0.0000000 *mat_022 22,7.e-4, 3e5,3e5,3e5,0,0,0 1.5e5,1.5e5,1.5e5, 1e5 1e9,1e9,1e9,1e9, ,1e9,1e9,1e9 *ELEMENT_solid 1 1 23 24 35 34 1 2 13 12 2 1 24 25 36 35 2 3 14 13 3 1 25 26 37 36 3 4 15 14 4 1 26 27 38 37 4 5 16 15 5 1 27 28 39 38 5 6 17 16 6 1 28 29 40 39 6 7 18 17 7 1 29 30 41 40 7 8 19 18 8 1 30 31 42 41 8 9 20 19 9 1 31 32 43 42 9 10 21 20 10 1 32 33 44 43 10 11 22 21 *DEFINE_CURVE 3 0 0.0000000 0.0000000 0.0000000 0.0000000 0.00000000E+00 0.00000000E+00 2. 1.00000000E+00 1.00000000E+01 1.00000000E+00 *SET_NODE 100 11,22,33,44 *SET_NODE 101 33,44 *LOAD_NODE_SET $ axial tension $ 100 1 3 $$ y load $ 100 3 3 $$ z load $ 100 3 3 -5.0 101 3 3 -10. *NODE 1 0.0000000e+00 0.0000000e+00 0.0000000e+00 2 3.0000000e+00 0.0000000e+00 0.0000000e+00 3 6.0000000e+00 0.0000000e+00 0.0000000e+00 4 9.0000000e+00 0.0000000e+00 0.0000000e+00 5 1.2000000e+01 0.0000000e+00 -5.5511151e-17 6 1.5000000e+01 0.0000000e+00 0.0000000e+00 7 1.8000000e+01 0.0000000e+00 0.0000000e+00 8 2.1000000e+01 0.0000000e+00 -2.7755576e-17 9 2.4000000e+01 0.0000000e+00 2.7755576e-17 10 2.7000000e+01 0.0000000e+00 -4.1633363e-17 11 3.0000000e+01 0.0000000e+00 0.0000000e+00 12 0.0000000e+00 -3.0000000e+00 0.0000000e+00 13 3.0000000e+00 -3.0000000e+00 0.0000000e+00 14 6.0000000e+00 -3.0000000e+00 0.0000000e+00 15 9.0000000e+00 -3.0000000e+00 0.0000000e+00 16 1.2000000e+01 -3.0000000e+00 0.0000000e+00 17 1.5000000e+01 -3.0000000e+00 0.0000000e+00 18 1.8000000e+01 -3.0000000e+00 0.0000000e+00 19 2.1000000e+01 -3.0000000e+00 0.0000000e+00 20 2.4000000e+01 -3.0000000e+00 0.0000000e+00 21 2.7000000e+01 -3.0000000e+00 0.0000000e+00 22 3.0000000e+01 -3.0000000e+00 0.0000000e+00 23 0.0000000e+00 0.0000000e+00 3.0000000e+00 24 3.0000000e+00 0.0000000e+00 3.0000000e+00 25 6.0000000e+00 0.0000000e+00 3.0000000e+00 26 9.0000000e+00 0.0000000e+00 3.0000000e+00 27 1.2000000e+01 0.0000000e+00 3.0000000e+00 28 1.5000000e+01 0.0000000e+00 3.0000000e+00 29 1.8000000e+01 0.0000000e+00 3.0000000e+00 30 2.1000000e+01 0.0000000e+00 3.0000000e+00 31 2.4000000e+01 0.0000000e+00 3.0000000e+00 32 2.7000000e+01 0.0000000e+00 3.0000000e+00 33 3.0000000e+01 0.0000000e+00 3.0000000e+00 34 0.0000000e+00 -3.0000000e+00 3.0000000e+00 35 3.0000000e+00 -3.0000000e+00 3.0000000e+00 36 6.0000000e+00 -3.0000000e+00 3.0000000e+00 37 9.0000000e+00 -3.0000000e+00 3.0000000e+00 38 1.2000000e+01 -3.0000000e+00 3.0000000e+00 39 1.5000000e+01 -3.0000000e+00 3.0000000e+00 40 1.8000000e+01 -3.0000000e+00 3.0000000e+00 41 2.1000000e+01 -3.0000000e+00 3.0000000e+00 42 2.4000000e+01 -3.0000000e+00 3.0000000e+00 43 2.7000000e+01 -3.0000000e+00 3.0000000e+00 44 3.0000000e+01 -3.0000000e+00 3.0000000e+00 *END