各向同性材料损伤本构模型SUBROUTINEUMATSTRESSSTATEVDDSDDESSESPDSCDRPLDDSDDTDRPLDEDRPLDTSTRANDSTRANTIMEDTIMETEMPDTEMPPREDEFDPREDCMNAMENDINSHRNTENSNSTATVPROPSNPROPSCOORDSDROTPNEWDTCELENTDFGRD0DFGRD1NOELNPTLAYERKSPTKSTEPKINC
INCLUDEABA_PARAMINCCHARACTER80CMNAME
DIMENSIONSTRESSNTENSSTATEVNSTATVDDSDDENTENSNTENSDDSDDTNTENSDRPLDENTENSSTRANNTENSDSTRANNTENSTIME2PREDEF1DPRED1PROPSNPROPSCOORDS3DROT33DFGRD033DFGRD133
DIMENSIONSTRANT6TSTRANT4PT1DIMENSIONOLD_STRESS6DIMENSIONDOLD_STRESS6D_STRESS6
DIMENSIONC66CD66DSTRESS6BSTRESS6ROOT3
DFMNDE6DDMDE6DCDDM66ATEMP16ATEMP26
PARAMETERZERO0D0ONE1D0TWO2D0FOUR4D0HALF05D0
Cstart
C
IFNPROPSLT2THEN
C
WRITE7ERRORUMATREQUIRESNPROPS2
C
STOP
C
ENDIF
E11PROPS1V12PROPS2G12PROPS1TWOONEPROPS2
CCriticalvaluesofstressesXTPROPS3XCPROPS4XSPROPS5GXPROPS6Fracturee
ergyi
matrixETA0001
CCurre
tstrai
DOI1NTENSSTRANTISTRANIDSTRANIENDDO
CStiff
essDOI16DOJ16CIJZEROENDDOENDDO
fATEMP1V121TWOV12C11E111V12ATEMPC22E111V12ATEMPC33E111V12ATEMPC12E11V12ATEMPC13E11V12ATEMPC23E11V12ATEMPC44G12C55G12C66G12DOI26
DOJ1I1CIJCJI
ENDDOENDDO
CCriticalvaluesofstrai
sXETXTC112V12C12XECXCC112V12C12XESXSC44
DMOLDSTATEV1
CStrai
i
itiatio
criterio
A11STRANT1TWOSTRANT2TWOSTRANT3TWOA12A11XETXECA21STRANT1STRANT2STRANT3A22XECXETXECXETA21A31STRANT4TWOSTRANT5TWOSTRANT6TWOA32A31XESTWO
A1A12A22A32
C
B11STRANT2TWO
C
B12B11XETXEC
C
B21STRANT2
C
B22XECXETXECXETB21
C
B31STRANT5TWO
C
B32B31XESTWO
C
B1B12B22B32
C
C11STRANT3TWO
C
C12C11XETXEC
C
C21STRANT3
C
C22XECXETXECXETC21
C
C31STRANT6TWO
C
C32C31XESTWO
fC
C1C12C22C32
STATEV2A1
C
STATEV3B1
C
STATEV4C1
FMNZERO
IFA1GTZEROTHEN
FMNSQRTA1
C
IFB1GTONETHEN
C
FMNFMNSQRTB1
C
IFC1GTONETHEN
C
FMNFMNSQRTC1
C
ENDIF
C
ENDIF
ENDIF
STATEV5FMN
C
writeFMN
DMZERO
DDMDFMNZERO
DOI16
DFMNDEIZERO
DDMDEIZERO
ENDDO
IFFMNGTONETHEN
CCALCULATEDMDDMDFMN
C
WRITE6FMN
T1C112V12C12XET2CELENTGX
T2ONEFMNT1
DMONEEXPT2FMN
C
WRITE6T1T1T2T2DMDM
C
writeDM
CCALCULATETHEDERIVATIVEOFDAMAGEVARIABLEWITHRESPECTTOFAILURE
CRITERION
DDMDFMNONEFMNT1ONEDM
CCALCULATEDFMNDEIFDMGTDMOLDTHENr