*KEYWORD $ strengths, including compressive strengths, should be input positive $ 15 layers $ cmpflg=1 $ 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 $ fiber tensile strain reaches DFAILT $ mat55: nonzero FBRT has dubious effect >> set to 0. $ IP fails only if sig-aa reaches XT. $ Other strengths will cap the stress except XC may not serve as cap (YC $ may come into play for capping fiber compression. $ mat58: peak fiber tension at XT and then softens in a controlled manner $ mat59: fiber tension capped at XT * SR $ *database_bndout .01 *control_shell ,,,1 $,,,,,1 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $*control_accuracy $1,2 *CONTROL_TERMINATION $ ENDTIM ENDCYC DTMIN ENDNEG ENDMAS .55000E-0 $---+----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,15,1 1 *DATABASE_BINARY_D3PLOT $ DT/CYCL LCDT/NR BEAM NPLTC PSETID ISTATS TSTART IAVG 0.1000E-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 1 2 15 -2020 1 .1,.1,.1,.1 0. 90. 45. -45 0. 90. 45. 0 45 90. 0. -45 45 90. 0. *INTEGRATION_SHELL 2020 15 1 $---+----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 global sig-zz = 1e6 whether XC is neg or pos 1.e6 1.e6 1.e6 1.e6 1.e6 54 0.0 *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 1.0 0. 1.0 0.2 $ 0.2 0.05 1.0 0.0123008 0.0119812 0.0121410 $ XC XT YC YT SC CRIT BETA 1.e6 1.e6 1.e6 1.e6 1.e6 55 0.0 *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 0.0 -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 $ XC XT YC YT SC 1.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) $ XXT YYT ZZT (solids) $ XC XT YC YT SC (shells) 1.e6 1.e6 1.e6 1.e6 1.e6 *END