$ model is driven by displacement $ *KEYWORD *comment The attached one element model loaded to failure in fiber tension and then compressed shows the effect of KFAIL. Stresses are output in the material coordinate system: x-stress = fiber stress = material a-direction = global z-direction. Here's what I observe: KFAIL.gt.0: When change in volume first becomes negative subsequent to element failure in fiber tension, pressure immediately is reset to zero and then increases in proportion to further change in volume according to the bulk modulus = KFAIL. The physical basis for the compressive stress discontinuity is not obvious to me. KFAIL.eq.0: Compressive stresses do not undergo a discontinuity, i.e., pressure is not reset to zero when volume change becomes negative. Jim 11/29/11 *boundary_spc_set_birth_death $ constrain in transverse direction at t=0.01 so that change in volume can become negative 18,,1,1 0.01 *set_node_list 18 1,2,3,4,5,6,7,8 *MAT_COMPOSITE_DAMAGE $# mid ro ea eb ec prba prca prcb $ 1 1.4800E-4 3.7000E+6 1.2100E+7 1.7500E+6 0.300000 0.300000 0.300000 $$ isotropic 1 1.4800E-4 3.7000E+6 3.700E+06 3.700E+6 0.300000 0.300000 0.300000 $# gab gbc gca kfail aopt macf $ 3.0000E+6 5.7500E+5 4.8500E+5 0.000000 3.000000 $ 3.0000E+6 5.7500E+5 4.8500E+5 1.E7 0.000000 3.000000 $$ isotropic 1.42e6 1.42e6 1.42e6 1.E7 0.000000 3.000000 $ 1.42e6 1.42e6 1.42e6 0 0.000000 3.000000 $# xp yp zp a1 a2 a3 $# v1 v2 v3 d1 d2 d3 beta $# sc xt yt yc alph sn syz szx $ 92500.000 92500.000 1.8500E+5 1.8500E+5 0.000 0.000 0.000 0.000 1.e9 92500. 1.e9 1.e9 0.000 1.e9 1.e9 1.e9 *database_binary_d3thdt 10.e-5 *database_history_solid 1 *database_glstat 0.2 *TITLE mat22.pull.english.k *DEFINE_CURVE $ Bolt material curve $ i i f f f f i i $ lcid sidr sfa sfo offa offo dattyp 2 0 0.000, 48600 0.012, 80475 0.022, 96125 0.032, 107174 0.042, 115943 0.052, 123313 0.062, 129723 0.072, 135429 0.082, 140591 0.092, 145319 0.102, 149692 0.112, 153767 0.122, 157590 0.132, 161194 0.142, 164608 0.152, 167854 0.162, 170951 0.172, 173914 0.182, 176755 0.192, 179488 0.202, 182120 0.212, 184661 0.222, 187117 0.232, 189496 0.242, 191802 0.252, 194041 0.262, 196217 0.272, 198335 0.282, 200398 0.292, 202409 0.302, 204372 0.312, 206289 0.322, 208162 0.332, 209994 0.342, 211787 0.352, 213543 0.362, 215264 0.372, 216951 0.382, 218606 *define_curve 13 0,45600 0.006,47000 0.012,55000 0.132616, 89965.83 *DEFINE_CURVE $ Plate material curve $ i i f f f f i i$ lcid sidr sfa sfo offa offo dattyp 3 0 0, 45600 0.002837, 45603.43 0.002858, 45706.18 0.002892, 45800.74 0.002932, 45881.94 0.002965, 45933.79 0.002994, 45997.75 0.003023, 46067.36 0.003062, 46114.62 0.003090, 46155.22 0.003124, 46183.74 0.003152, 46243.01 0.003189, 46285.72 0.003215, 46310.08 0.003245, 46348.11 0.003280, 46372.66 0.003314, 46400.40 0.003330, 46358.57 0.003355, 46377.35 0.003386, 46405.76 0.003419, 46438.18 0.003443, 46473.96 0.003476, 46527.19 0.003514, 46600.99 0.003549, 46667.11 0.003589, 46737.11 0.003626, 46805.92 0.003679, 46879.67 0.003719, 46906.83 0.003746, 46930.09 0.003780, 46976.64 0.003812, 47026.31 0.003844, 47073.05 0.003877, 47125.57 0.003925, 47208.56 0.003947, 47240.95 0.003985, 47277.60 0.004011, 47332.01 0.004047, 47384.06 0.004080, 47441.37 0.004113, 47488.56 0.004152, 47538.66 0.004188, 47579.05 0.004222, 47617.58 0.004265, 47656.64 0.004297, 47691.77 0.004332, 47704.93 0.004368, 47739.16 0.004404, 47770.38 0.004449, 47809.89 0.004481, 47826.28 0.004495, 47847.52 0.004524, 47898.17 0.004568, 47956.96 0.004592, 47995.21 0.004619, 48011.24 0.004656, 48063.33 0.004694, 48114.97 0.004730, 48146.60 0.004770, 48178.65 0.004801, 48211.26 0.004846, 48223.33 0.004878, 48246.13 0.004903, 48253.61 0.004940, 48292.52 0.004972, 48310.73 0.005007, 48336.21 0.005036, 48374.03 0.005074, 48402.15 0.005103, 48437.11 0.005139, 48469.77 0.005153, 48497.49 0.005194, 48542.88 0.005228, 48601.88 0.005264, 48625.77 0.005299, 48649.46 0.005330, 48698.23 0.005363, 48733.08 0.005408, 48766.40 0.005447, 48792.78 0.005471, 48790.18 0.005494, 48818.39 0.005530, 48878.13 0.005566, 48908.75 0.005606, 48957.74 0.005644, 48988.30 0.005684, 49016.83 0.005706, 49042.69 0.005743, 49065.95 0.005780, 49082.42 0.005806, 49110.77 0.005840, 49134.03 0.005868, 49166.02 0.005892, 49199.63 0.005918, 49212.83 0.005941, 49269.66 0.005988, 49357.25 0.006019, 49384.23 0.006060, 49420.97 0.006099, 49444.59 0.006141, 49467.80 0.006167, 49472.41 0.006207, 49496.50 0.006235, 49512.54 0.006264, 49536.64 0.006293, 49575.10 0.006321, 49601.29 0.006354, 49661.99 0.006396, 49696.03 0.006433, 49733.01 0.006459, 49762.27 0.006498, 49792.10 0.006532, 49815.10 0.006557, 49828.69 0.006592, 49867.58 0.006627, 49902.77 0.006658, 49935.13 0.006693, 49960.09 0.006729, 49988.55 0.006758, 50016.03 0.006790, 50052.17 0.006837, 50086.11 0.006873, 50114.14 0.006906, 50143.17 0.006940, 50152.28 0.006967, 50170.07 0.006995, 50193.64 0.007026, 50224.54 0.007057, 50271.32 0.007095, 50295.39 0.007120, 50325.05 0.007147, 50354.54 0.007180, 50409.09 0.007202, 50429.39 0.007235, 50466.31 0.007280, 50515.66 0.007310, 50547.76 0.007354, 50574.11 0.007396, 50625.93 0.007425, 50625.14 0.007461, 50638.66 0.007494, 50650.87 0.007524, 50680.42 0.007553, 50714.86 0.007580, 50738.56 0.007614, 50778.57 0.007649, 50809.58 0.007679, 50851.21 0.007709, 50870.03 0.007747, 50921.28 0.007785, 50954.62 0.007831, 50966.67 0.007858, 50986.70 0.007884, 51005.52 0.007919, 51020.30 0.007945, 51044.69 0.007979, 51086.25 0.008005, 51114.16 0.008037, 51159.85 0.008063, 51204.92 0.008096, 51254.24 0.008140, 51291.74 0.008181, 51327.64 0.008220, 51342.29 0.008256, 51347.25 0.008286, 51365.09 0.008311, 51379.10 0.008347, 51408.88 0.008375, 51427.80 0.008415, 51468.36 0.008443, 51510.43 0.008474, 51540.93 0.008495, 51483.21 0.008527, 51551.73 0.008562, 51596.02 0.008594, 51636.98 0.008629, 51682.79 0.008670, 51711.14 0.008714, 51736.55 0.008745, 51741.70 0.008771, 51759.61 0.008803, 51786.80 0.008836, 51814.29 0.008872, 51830.61 0.008907, 51865.59 0.008939, 51910.19 0.008966, 51932.81 0.008997, 51948.70 0.009038, 51991.17 0.009069, 52018.61 0.009103, 52037.90 0.009130, 52064.52 0.009152, 52086.16 0.009193, 52122.30 0.009234, 52167.42 0.009266, 52205.44 0.009292, 52227.31 0.009323, 52263.14 0.009361, 52299.40 0.009396, 52322.44 0.009434, 52346.96 0.009470, 52380.28 0.009502, 52412.51 0.009533, 52432.37 0.009563, 52469.82 0.009594, 52489.05 0.009634, 52506.94 0.009660, 52532.38 0.009689, 52562.71 0.009729, 52591.66 0.009755, 52618.06 0.009779, 52634.94 0.009814, 52655.67 0.009847, 52683.76 0.009884, 52730.60 0.009911, 52765.53 0.009950, 52784.58 0.009981, 52820.50 0.010004, 52854.93 0.010044, 52891.74 0.010082, 52923.80 0.010123, 52963.81 0.010156, 52996.14 0.010186, 53013.40 0.010233, 53034.61 0.010263, 53055.43 0.010297, 53068.59 0.010332, 53084.94 0.010358, 53088.25 0.010394, 53122.47 0.010416, 53132.18 0.010450, 53166.79 0.010475, 53202.45 0.010511, 53249.49 0.010547, 53292.38 0.010584, 53314.23 0.010623, 53346.97 0.010651, 53364.99 0.010683, 53393.98 0.010727, 53425.63 0.010759, 53440.08 0.010787, 53470.41 0.010815, 53491.40 0.010837, 53515.35 0.010865, 53543.36 0.010905, 53598.69 0.010942, 53630.80 0.010963, 53654.38 0.010998, 53692.92 0.011041, 53723.02 0.011060, 53718.02 0.011110, 53768.37 0.011134, 53784.19 0.011164, 53811.36 0.011201, 53851.11 0.011239, 53873.57 0.011276, 53888.27 0.011297, 53918.82 0.011336, 53938.69 0.011362, 53970.54 0.011400, 54013.78 0.011437, 54029.11 0.011477, 54073.33 0.011515, 54096.99 0.011556, 54109.66 0.011583, 54120.97 0.011605, 54118.50 0.011634, 54140.81 0.011665, 54161.01 0.011700, 54203.72 0.011734, 54234.46 0.011769, 54269.86 0.011796, 54296.33 0.011830, 54332.53 0.011874, 54353.48 0.011909, 54381.63 0.011936, 54407.76 0.011964, 54431.35 0.011989, 54439.03 0.012026, 54482.55 0.012054, 54493.60 0.012087, 54533.35 0.012121, 54570.76 0.012150, 54586.46 0.012181, 54617.36 0.012221, 54645.87 0.012249, 54664.34 0.012288, 54716.76 0.012319, 54751.65 0.012343, 54762.65 0.012376, 54792.86 0.012417, 54806.98 0.012450, 54830.63 0.012488, 54863.03 0.012557, 54892.46 0.012571, 54915.21 0.012648, 54949.32 0.012653, 54964.45 0.012653, 54990.15 0.012653, 55020.45 0.012735, 55037.94 0.012736, 55069.39 0.012808, 55088.59 0.012818, 55109.91 $0.012818, 55128.38 starts duplicating strain values - Not allowed! 0.012897, 55165.07 0.012901, 55197.93 0.012982, 55223.31 0.012983, 55245.88 $0.012983, 55258.26 $0.012983, 55280.83 0.013066, 55312.63 0.013066, 55345.76 0.013148, 55380.85 0.013148, 55393.58 0.013148, 55424.36 0.013231, 55447.95 0.013274, 55472.37 0.013313, 55505.44 0.013313, 55561.03 0.013468, 55960.92 0.013641, 56395.47 0.013962, 56674.53 0.014303, 56946.89 0.014738, 57253.91 0.015194, 57576.15 0.015644, 57871.45 0.016088, 58223.57 0.016578, 58595.61 0.017108, 58991.99 0.017699, 59382.87 0.018268, 59740.09 0.018859, 60102.97 0.019393, 60448.23 0.019896, 60792.45 0.020416, 61160.07 0.020936, 61518.69 0.021492, 61859.89 0.022360, 62355.14 0.022929, 62695.58 0.023439, 62989.87 0.023945, 63328.99 0.024494, 63649.36 0.025097, 63948.25 0.025724, 64231.51 0.026305, 64514.13 0.026830, 64779.84 0.027288, 65082.14 0.027764, 65400.05 0.028283, 65705.53 0.028830, 65994.61 0.029372, 66277.70 0.029921, 66569.95 0.030454, 66853.12 0.031000, 67135.16 0.031556, 67401.03 0.032095, 67676.04 0.032620, 67956.93 0.033176, 68224.61 0.033700, 68486.64 0.034241, 68744.07 0.034788, 68987.64 0.035290, 69255.54 0.035798, 69513.96 0.036340, 69768.46 0.036891, 70014.49 0.037442, 70265.99 0.038008, 70497.72 0.038561, 70725.92 0.039137, 70964.76 0.039720, 71187.09 0.040289, 71394.18 0.040776, 71606.95 0.041265, 71826.59 0.041768, 72045.02 0.042285, 72277.47 0.042858, 72492.03 0.043410, 72676.43 0.043888, 72889.35 0.044416, 73104.80 0.044929, 73314.57 0.045502, 73516.01 0.046075, 73699.31 0.046649, 73890.15 0.047155, 74078.74 0.047653, 74271.17 0.048158, 74460.21 0.048632, 74668.88 0.049158, 74863.35 0.049707, 75056.83 0.050236, 75243.34 0.050801, 75429.02 0.051353, 75597.22 0.051907, 75764.23 0.052425, 75930.73 0.052959, 76095.09 0.053480, 76279.25 0.053978, 76450.48 0.054495, 76626.41 0.055014, 76791.49 0.055550, 76951.89 0.056077, 77124.88 0.056645, 77280.26 0.057177, 77437.70 0.057745, 77581.31 0.058263, 77722.47 0.058750, 77881.24 0.059225, 78056.89 0.059744, 78201.18 0.060282, 78351.70 0.060792, 78490.29 0.061322, 78646.12 0.061861, 78801.84 0.062390, 78938.21 0.062965, 79057.21 0.063493, 79192.29 0.063990, 79320.37 0.064479, 79640.79 0.064995, 79810.76 0.065500, 79937.57 0.066011, 80056.78 0.066497, 80180.58 0.066977, 80332.37 0.067520, 80484.56 0.068068, 80613.50 0.068648, 80736.17 0.069214, 80827.60 0.069725, 80950.85 0.070250, 81064.87 0.070767, 81183.72 0.071275, 81314.44 0.071810, 81424.90 0.072345, 81534.04 0.072902, 81627.73 0.073433, 81742.23 0.073886, 81873.26 0.074365, 81982.59 0.074901, 82094.82 0.075462, 82194.61 0.076020, 82293.47 0.076533, 82381.56 0.076999, 82503.95 0.077479, 82614.63 0.077984, 82731.22 0.078481, 82820.20 0.078982, 82928.15 0.079478, 83043.49 0.079988, 83159.06 0.080516, 83247.42 0.081028, 83351.56 0.081562, 83463.67 0.082091, 83552.04 0.082592, 83659.60 0.083102, 83741.67 0.083630, 83831.48 0.084133, 83923.87 0.084685, 84008.29 0.085207, 84082.74 0.085694, 84171.25 0.086201, 84266.03 0.086702, 84358.26 0.087187, 84478.63 0.087709, 84551.07 0.088227, 84666.85 0.088710, 84742.54 0.089278, 84846.21 0.089797, 84884.36 0.090300, 84979.14 0.090815, 85070.28 0.091309, 85136.62 0.091827, 85211.90 0.092328, 85300.88 0.092825, 85389.39 0.093331, 85456.20 0.093842, 85548.73 0.094353, 85624.81 0.094873, 85682.32 0.095366, 85756.01 0.095843, 85853.67 0.096314, 85943.55 0.096807, 86026.73 0.097348, 86093.67 0.097850, 86136.93 0.098395, 86239.46 0.098905, 86303.03 0.099405, 86352.13 0.099881, 86457.43 0.100388, 86488.49 0.100857, 86580.41 0.101365, 86666.71 0.101917, 86730.89 0.102460, 86779.40 0.102977, 86829.01 0.103471, 86884.98 0.103941, 86971.80 0.104404, 87063.00 0.104927, 87140.88 0.105438, 87176.88 0.105962, 87237.15 0.106454, 87299.69 0.106924, 87346.58 0.107382, 87449.77 0.107901, 87520.17 0.108425, 87572.07 0.108981, 87611.23 0.109508, 87654.57 0.110037, 87712.46 0.110528, 87744.76 0.110999, 87809.58 0.111463, 87871.27 0.111947, 87943.54 0.112433, 88017.84 0.112901, 88064.10 0.113393, 88126.59 0.113928, 88190.26 0.114454, 88250.83 0.114961, 88284.71 0.115426, 88334.27 0.115901, 88412.45 0.116373, 88470.83 0.116873, 88516.45 0.117357, 88566.07 0.117853, 88633.03 0.118365, 88679.66 0.118835, 88748.38 0.119321, 88801.64 0.119854, 88857.12 0.120378, 88894.34 0.120891, 88920.57 0.121369, 88947.88 0.121856, 89018.85 0.122333, 89068.26 0.122830, 89129.05 0.123306, 89170.53 0.123789, 89212.04 0.124289, 89250.96 0.124754, 89297.49 0.125210, 89370.47 0.125704, 89431.56 0.126218, 89457.60 0.126736, 89488.06 0.127262, 89519.72 0.127742, 89562.26 0.128258, 89611.94 0.128745, 89649.88 0.129254, 89690.23 0.129729, 89742.97 0.130214, 89767.07 0.130707, 89802.41 0.131183, 89859.78 0.131680, 89906.19 0.132155, 89904.39 0.132616, 89965.83 $ $ *CONTROL_TERMINATION $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8 $ endtim endcyc dtmin endneg endmas 0.02 $ *CONTROL_TIMESTEP $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8 $ dtinit scft isdo tslimt dtms lctm erode ms1st .10 $ $ *DEFINE_CURVE $ lcid sidr scla sclo offa offo 5 $ abscissa ordinate 0.0 0.001 100.0 0.001 $ $ *CONTROL_ENERGY $ hgen rwen slnten rylen 2 2 $ *CONTROL_OUTPUT $ npopt neecho nrefup iaccop opifs ipnint ikedit 1000 $ $ *DATABASE_BINARY_D3PLOT $ dt/cycl lcdt 10e-5 $ *DATABASE_EXTENT_BINARY $ neiph neips maxint strflg sigflg epsflg rltflg engflg 3 1 $ $ cmpflg ieverp beamip 1 $ $ $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $$$$ Define Cross Sections $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $ *DATABASE_SECFORC $ dt/cycl lcdt 0.1 $ $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8 $ *DATABASE_CROSS_SECTION_PLANE $ i $ psid xct yct zct xch ych zch 0 0.0 0.0 1.0 0.0 0.0 5.0 $ xhev yhev xhev lenl lenm 1.0 0.0 1.0 $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $$$$ Boundary Motion Conditions $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8 $ *BOUNDARY_PRESCRIBED_MOTION_SET $ nid dof vad lcid sf vid 1 3 2 1 1.0 0 $ *DEFINE_CURVE $ lcid sidr scla sclo offa offo 1 0,0 0.005,.3 0.01,0 0.02,-.4 1.00,-.4 *SET_NODE_LIST $ sid da1 da2 da3 da4 1 $ $ nid1 nid2 nid3 nid4 nid5 nid6 nid7 nid8 5 6 7 8 $ $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $$$$ Define Parts and Materials $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8 $ *PART Solid Element $ pid sid mid eosid hgid adpopt 1 1 1 0 0 0 $ $ $ *SECTION_SOLID $ sid elform 1 $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $$$$ Define Nodes and Elements $ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $ $ *NODE $ i8 f f f i i $ node x y z tc rc 1 0.0 0.0 0.0 3 7 2 10.0 0.0 0.0 3 7 3 0.0 10.0 0.0 3 7 4 10.0 10.0 0.0 3 7 5 0.0 0.0 10.0 0 7 6 10.0 0.0 10.0 0 7 7 0.0 10.0 10.0 0 7 8 10.0 10.0 10.0 0 7 *ELEMENT_SOLID $ i i i i i i i i i i $ eid pid n1 n2 n3 n4 n5 n6 n7 n8 1 1 1 2 4 3 5 6 8 7 $ *END