hashin失效vumat_hashin破坏准则_vumatfailure_vumathashin失效_hashin_vumat

c http://abaqus-users.1086179.n5.nabble.com/Hashin-damage-theory-td12032.html subroutine vumat( c Read only - 1 nblock, ndir, nshr, nstatev, nfieldv, nprops, lanneal, 2 stepTime, totalTime, dt, cmname, coordMp, charLength, 3 props, density, strainInc, relSpinInc, 4 tempOld, stretchOld, defgradOld, fieldOld, 5 stressOld, stateOld, enerInternOld, enerInelasOld, 6 tempNew, stretchNew, defgradNew, fieldNew, c Write only - 7 stressNew, stateNew, enerInternNew, enerInelasNew ) c include 'vaba_param.inc' c c 3D Orthotropic Elasticity with Hashin 3d Failure criterion c c The state variables are stored as: c state(*,1) = material point status c state(*,2:7) = damping stresses c c User defined material properties are stored as c * First line: c props(1) --> Young's modulus in 1-direction, E1 c props(2) --> Young's modulus in 2-direction, E2 c props(3) --> Young's modulus in 3-direction, E3 c props(4) --> Poisson's ratio, nu12 c props(5) --> Poisson's ratio, nu13 c props(6) --> Poisson's ratio, nu23 c props(7) --> Shear modulus, G12 c props(8) --> Shear modulus, G13 c c * Second line: c props(9) --> Shear modulus, G23 c props(10) --> beta damping parameter c props(11) --> "not used" c props(12) --> "not used" c props(13) --> "not used" c props(14) --> "not used" c props(15) --> "not used" c props(16) --> "not used" c c * Third line: c props(17) --> Ultimate tens stress in 1-direction, sigu1t c props(18) --> Ultimate comp stress in 1-direction, sigu1c c props(19) --> Ultimate tens stress in 2-direction, sigu2t c props(20) --> Ultimate comp stress in 2-direction, sigu2c c props(21) --> Ultimate tens stress in 2-direction, sigu3t c props(22) --> Ultimate comp stress in 2-direction, sigu3c c props(23) --> "not used" c props(24) --> "not used" c c * Fourth line: c props(25) --> Ultimate shear stress, sigu12 c props(26) --> Ultimate shear stress, sigu13 c props(27) --> Ultimate shear stress, sigu23 c props(28) --> "not used" c props(29) --> "not used" c props(30) --> "not used" c props(31) --> "not used" c props(32) --> "not used" c dimension props(nprops), density(nblock), 1 coordMp(nblock,*), 2 charLength(*), strainInc(nblock,ndir+nshr), 3 relSpinInc(nblock,nshr), tempOld(nblock), 4 stretchOld(nblock,ndir+nshr), defgradOld(nblock,ndir+nshr+nshr), 5 fieldOld(nblock,nfieldv), stressOld(nblock,ndir+nshr), 6 stateOld(nblock,nstatev), enerInternOld(nblock), 7 enerInelasOld(nblock), tempNew(*), 8 stretchNew(nblock,ndir+nshr), defgradNew(nblock,ndir+nshr+nshr), 9 fieldNew(nblock,nfieldv), stressNew(nblock,ndir+nshr), 1 stateNew(nblock,nstatev), 2 enerInternNew(nblock), enerInelasNew(nblock) * character*80 cmname * parameter( zero = 0.d0, one = 1.d0, two = 2.d0, half = .5d0 ) * parameter( * i_svd_DmgFiberT = 1, * i_svd_DmgFiberC = 2, * i_svd_DmgMatrixT = 3, * i_svd_DmgMatrixC = 4, * i_svd_statusMp = 5, * i_svd_dampStress = 6, c * i_svd_dampStressXx = 6, c * i_svd_dampStressYy = 7, c * i_svd_dampStressZz = 8, c * i_svd_dampStressXy = 9, c * i_svd_dampStressYz = 10, c * i_svd_dampStressZx = 11, * i_svd_Strain = 12, c * i_svd_StrainXx = 12, c * i_svd_StrainYy = 13, c * i_svd_StrainZz = 14, c * i_svd_StrainXy = 15, c * i_svd_StrainYz = 16, c * i_svd_StrainZx = 17, * n_svd_required = 17 ) * parameter( * i_s33_Xx = 1, * i_s33_Yy = 2, * i_s33_Zz = 3, * i_s33_Xy = 4, * i_s33_Yz = 5, * i_s33_Zx = 6 ) * * Structure of property array parameter ( * i_pro_E1 = 1, * i_pro_E2 = 2, * i_pro_E3 = 3, * i_pro_nu12 = 4, * i_pro_nu13 = 5, * i_pro_nu23 = 6, * i_pro_G12 = 7, * i_pro_G13 = 8, * i_pro_G23 = 9, * * i_pro_beta = 10, * * i_pro_sigu1t = 17, * i_pro_sigu1c = 18, * i_pro_sigu2t = 19, * i_pro_sigu2c = 20, * i_pro_sigu3t = 21, * i_pro_sigu3c = 22, * i_pro_sigu12 = 25, * i_pro_sigu13 = 26, * i_pro_sigu23 = 27 ) * Temporary arrays dimension eigen(maxblk*3) * * Read material properties * E1 = props(i_pro_E1) E2 = props(i_pro_E2) E3 = props(i_pro_E3) xnu12 = props(i_pro_nu12) xnu13 = props(i_pro_nu13) xnu23 = props(i_pro_nu23) G12 = props(i_pro_G12) G13 = props(i_pro_G13) G23 = props(i_pro_G23) * xnu21 = xnu12 * E2 / E1 xnu31 = xnu13 * E3 / E1 xnu32 = xnu23 * E3 / E2 * * * Compute terms of stiffness matrix gg = one / ( one - xnu12*xnu21 - xnu23*xnu32 - xnu31*xnu13 * - two*xnu21*xnu32*xnu13 ) C11 = E1 * ( one - xnu23*xnu32 ) * gg C22 = E2 * ( one - xnu13*xnu31 ) * gg C33 = E3 * ( one - xnu12*xnu21 ) * gg C12 = E1 * ( xnu21 + xnu31*xnu23 ) * gg C13 = E1 * ( xnu31 + xnu21*xnu32 ) * gg C23 = E2 * ( xnu32 + xnu12*xnu31 ) * gg * f1t = props(i_pro_sigu1t) f1c = props(i_pro_sigu1c) f2t = props(i_pro_sigu2t) f2c = props(i_pro_sigu2c) f3t = props(i_pro_sigu3t) f3c = props(i_pro_sigu3c) f12 = props(i_pro_sigu12) f13 = props(i_pro_sigu13) f23 = props(i_pro_sigu23) * beta = props(i_pro_beta) * * Assume purely elastic material at the beginning of the analysis * if ( totalTime .eq. zero ) then if (nstatev .lt. n_svd_Required) then call xplb_abqerr(-2,'Subroutine VUMAT requires the '// * 'specification of %I state variables. Check the '// * 'definition of *DEPVAR in the input file.', * n_svd_Required,zero,' ') call xplb_exit end if call OrthoEla3dExp ( nblock, * stateOld(1,i_svd_DmgFiberT), * stateOld(1,i_svd_DmgFiberC), * stateOld(1,i_svd_DmgMatrixT), * stateOld(1,i_svd_DmgMatrixC), * C11, C22, C33, C12, C23, C13, G12, G23, G13, * strainInc, * stressNew ) return end if * * Update total elastic strain call strainUpdate ( nblock, strainInc, * stateOld(1,i_svd_strain), stateNew(1,i_svd_strain) ) * * Stress update call OrthoEla3dExp ( nblock, * stateOld(1,i_svd_DmgFiberT), * stateOld(1,i_svd_DmgFiberC), * stateOld(1,i_svd_DmgMatrixT), * stateOld(1,i_svd_DmgMatrixC), * C11, C22, C33, C12, C23, C13, G12, G23, G13, * stateNew(1,i_svd_strain), * stressNew ) * * Failure evaluation * call copyr ( nblock, * stateOld(1,i_svd_DmgFiberT), stateNew(1,i_svd_DmgFiberT) ) call copyr ( nblock, * stateOld(1,i_svd_DmgFiberC), stateNew(1,i_svd_DmgFiberC) ) call copyr ( nblock, * stateOld(1,i_svd_DmgMatrixT), stateNew(1,i_svd_DmgMatrixT) ) call copyr ( nblock, * stateOld(1,i_svd_DmgMatrixC), stateNew(1,i_svd_DmgMatrixC) ) nDmg = 0 call eig33Anal ( nblock, stretchNew, eigen ) call Hashin3d ( nblock, nDmg, * f1t, f2t, f3t, f1c, f2c, f3c, f12, f23, f13, * stateNew(1,i_svd_DmgFiberT), * stateNew(1,i_svd_DmgFiberC),