$ plot hisv1, hisv5 from d3thdt and compare to x-stress from d3thdt $ Illustrates that stress release after reaching XT is more abrupt (though NOT instantaneous) than $ drop in efail (hisv5). The latter occurs over 100 time steps. *KEYWORD *database_history_shell 54 *database_binary_d3thdt 1e-10 $ use c=pcom.a-dir to postprocess $ NOTE: Use DP $ Set DFAILT in mat54 so that some plasticity is allowed after $ reaching XT. $ tension, xt set to 1e6 $ aopt=2 (a = global z, b = global x) $ results: $ mat22: fails at XT $ mat54: fails at XT unless DFAILT given in which case stress holds at XT until $ strain reaches DFAILT, then element is deleted. $ mat55: z-stress capped briefly at FBRT * XT and then stress drops $ to zero withOUT element deletion; nodes shoot to side and $ error results unless NFAIL1 is set to 1 so that element is $ deleted due to bad shape: remedy is to set FBRT to 0 (no reduction $ in tensile fiber strength) $ mat58: peak z-stress at XT and then softens in a controlled manner $ mat59: z-stress capped at XT * SR $ *database_bndout 0.01 *control_shell ,,,0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $*control_accuracy $1,2 *CONTROL_TERMINATION $ ENDTIM ENDCYC DTMIN ENDNEG ENDMAS $ .55000E-0 $0.3 0.1 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *CONTROL_TIMESTEP $ DTINIT TSSFAC ISDO TSLIMIT DT2MS LCTM ERODE MS1ST $ 1.0000E-7 .7 ,.7 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *database_extent_binary ,6,,1 1,,,,,,1 *DATABASE_BINARY_D3PLOT $ DT/CYCL LCDT/NR BEAM NPLTC PSETID ISTATS TSTART IAVG 0.5000E-2 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ Geometry $ *ELEMENT_SHELL $ EID PID NID1 NID2 NID3 NID4 2 2 1 2 3 4 22 22 11 12 13 14 54 54 21 22 23 24 55 55 31 32 33 34 58 58 41 42 43 44 59 59 51 52 53 54 $ *NODE $ NID X Y Z TC RC 1 0.00 0. 0.0 2 1.00 0. 0.0 3 1.00 0. 1.0 4 0.00 0. 1.0 11 2.00 0. 0.0 12 3.00 0. 0.0 13 3.00 0. 1.0 14 2.00 0. 1.0 21 4.00 0. 0.0 22 5.00 0. 0.0 23 5.00 0. 1.0 24 4.00 0. 1.0 31 6.00 0. 0.0 32 7.00 0. 0.0 33 7.00 0. 1.0 34 6.00 0. 1.0 41 8.00 0. 0.0 42 9.00 0. 0.0 43 9.00 0. 1.0 44 8.00 0. 1.0 51 10.00 0. 0.0 52 11.00 0. 0.0 53 11.00 0. 1.0 54 10.00 0. 1.0 *set_node 101 1,11,21,31,41,51 *set_node 102 2,12,22,32,42,52 *set_node 103 3,13,23,33,43,53 *set_node 104 4,14,24,34,44,54 *boundary_spc_set 101,0, 1,1,1,1,0,1 102,0, 0,1,1,1,0,1 103,0, 0,1,0,1,0,1 104,0, 1,1,0,1,0,1 $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ Boundary conditions including constant strain condition $ *BOUNDARY_PRESCRIBED_MOTION_SET $ SETID DOF VAD LCID SF VID DEATH BIRTH 2 3 0 2 1. $ Top nodes *SET_NODE_LIST $ SID DA1 DA2 DA3 DA4 2 $ NID1 NID2 NID3 NID4 NID5 NID6 NID7 NID8 3,4,13,14,23,24,33,34 43,44,53,54 $ *DEFINE_CURVE $ LCID SIDR SFA SFO OFFA OFFO DATTYP 1 $ ABCISSA ORDINATE 0. 0. 10. 0. $ *DEFINE_CURVE $ LCID SIDR SFA SFO OFFA OFFO DATTYP 2 $ ABCISSA ORDINATE 0. 1.00 10. 1.00 *define_curve 22 0,0 10,10 $ $ Material *PART 2,1,2 22,1,22 54,1,54 55,1,55 58,1,58 59,1,59 *SECTION_SHELL $ SECID ELFORM SHRF NIP PROPT QR/IRID ICOMP SETYP $ 1 2 0.8333 3.0 0.0 1 1 1 2 0.8333 1.0 0.0 1 1 $ T1 T2 T3 T4 NLOC MAREA 0.1000 0.1000 0.1000 0.1000 0.0000 0.0000 $ B1 B2 B3 B4 B5 B6 B7 B8 $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *MAT_orthotropic_elastic $ MID RO EA EB (EC) PRBA 2 1.4763E-4 11.400E+6 11.100E+6 1.6500E+6 0.4500 $ GAB GBC GCA aopt 0.7660E+6 0.7500E+6 0.7500E+6 2 $ XP YP ZP A1 A2 A3 0.0 0.0 1.0 $ V1 V2 V3 D1 D2 D3 1.0 0.0 0.0 *MAT_composite_damage $ MID RO EA EB (EC) PRBA 22 1.4763E-4 11.400E+6 11.100E+6 1.6500E+6 0.4500 $ GAB GBC GCA KFAIL aopt 0.7660E+6 0.7500E+6 0.7500E+6 2 $ XP YP ZP A1 A2 A3 0.0 0.0 1.0 $ V1 V2 V3 D1 D2 D3 1.0 0.0 0.0 $ SC XT YT YC ALPH SN SYZ SZX 1.e6 1.e6 1.e6 1.e6 *MAT_ENHANCED_COMPOSITE_DAMAGE $ MID RO EA EB EC PRBA PRCA PRCB 54 1.4763E-4 11.400E+6 11.100E+6 1.6500E+6 0.4500 $ GAB GBC GCA KF AOPT 0.7660E+6 0.7500E+6 0.7500E+6 2.0 $ A1 A2 A3 MANGLE 0.0 0.0 1.0 $ V1 V2 V3 D1 D2 D3 DFAILM DFAILS 1.0 0.0 0.0 $ 1.0 0.0 0.0 0.0121410 0.1600000 $ TFAIL ALPH SOFT FBRT YCFAC DFAILT DFAILC EFS 0.2 0.05 1.0 $ 0.2 0.05 1.0 0.2 $ XC XT YC YT SC CRIT BETA $ fails at sig-zz = 1e6 whether XC is neg or pos $ -2.e6 1.e6 -1.e6 1.e6 1.e6 54 0.05 2.e6 1.e6 1.e6 1.e9 1.e6 54 0.05 $130.00E+3 155.00E+3 124.00E+3 138.00E+3 18.100E+3 54 0.05 *MAT_ENHANCED_COMPOSITE_DAMAGE $ MID RO EA EB EC PRBA PRCA PRCB 55 1.4763E-4 11.400E+6 11.100E+6 1.6500E+6 0.4500 $ GAB GBC GCA KF AOPT 0.7660E+6 0.7500E+6 0.7500E+6 2.0 $ A1 A2 A3 MANGLE 0.0 0.0 1.0 $ V1 V2 V3 D1 D2 D3 DFAILM DFAILS 1.0 0.0 0.0 $ 1.0 0.0 0.0 0.0121410 0.1600000 $ TFAIL ALPH SOFT FBRT YCFAC DFAILT DFAILC EFS $ 0.2 0.05 1.0 $ FBRT limits tensile strength 0.2 1.0 0.2 $ XC XT YC YT SC CRIT BETA $ fails at sig-zz = 1e6 whether XC is neg or pos $ -2.e6 1.e6 -1.e6 1.e6 1.e6 55 0.05 2.e6 1.e6 1.e6 1.e9 1.e6 55 0.05 $130.00E+3 155.00E+3 124.00E+3 138.00E+3 18.100E+3 55 0.05 *MAT_LAMINATED_COMPOSITE_FABRIC $ MID RO EA EB (EC) PRBA TAU1 GAMMA1 58 1.4763E-4 11.400E+6 11.100E+6 1.6500E+6 0.4500 15.000E+3 0.0600000 $ GAB GBC GCA SLIMT1 SLIMC1 SLIMT2 SLIMC2 SLIMS 0.7660E+6 0.7500E+6 0.7500E+6 0.05 1.0 0.05 1.0 0.05 $ AOPT TSIZE ERODS SOFT FS 2.0 .25 -1.0 $ XP YP ZP A1 A2 A3 0.0 0.0 1.0 $ V1 V2 V3 D1 D2 D3 BETA 1.0 0.0 0.0 $ E11C E11T E22C E22T GMS $ .05 .05 0.0116981 0.0118966 0.1600000 .05 .10 0.0116981 0.0118966 0.1600000 $ XC XT YC YT SC $ 130.00E+3 155.00E+3 124.00E+3 138.00E+3 18.100E+3 2.e6 1.e6 1.e6 1.e6 1.e6 *MAT_COMPOSITE_FAILURE_SHELL_MODEL $ MID RO EA EB (EC) PRBA PRBC PRCB $ 59 1.4763E-4 11.400E+6 11.100E+6 1.6500E+6 0.4500 0.31 0.31 59 1.4763E-4 11.400E+6 11.100E+6 1.6500E+6 0.4500 $ GAB GBC GCA KF AOPT MAFLAG 0.7660E+6 0.7500E+6 0.7500E+6 2.0 $ XP YP ZP A1 A2 A3 0.0 0.0 1.0 $ V1 V2 V3 D1 D2 D3 BETA 1.0 0.0 0.0 $ SBA SCA SCB XXC YYC ZZC (solids) $ TSIZ ALP SOFT FBRT SR SF (shells) .8 $ XXT YYT ZZT (solids) $ XC XT YC YT SC (shells) 2.e6 1.e6 1.e6 1.e6 1.e6 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ *END