abaqus材料子程序Word文档下载推荐.docx
- 文档编号:15494255
- 上传时间:2022-11-02
- 格式:DOCX
- 页数:5
- 大小:16.08KB
abaqus材料子程序Word文档下载推荐.docx
《abaqus材料子程序Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《abaqus材料子程序Word文档下载推荐.docx(5页珍藏版)》请在冰豆网上搜索。
DIMENSIONSTRANT(6),TSTRANT(4),PT
(1)
DIMENSIONOLD_STRESS(6)
DIMENSIONDOLD_STRESS(6),D_STRESS(6)
DIMENSIONC(6,6),CD(6,6),DSTRESS(6),BSTRESS(6),ROOT(3),
+DFMNDE(6),DDMDE(6),DCDDM(6,6),ATEMP1(6),ATEMP2(6)
PARAMETER(ZERO=0.D0,ONE=1.D0,TWO=2.D0,FOUR=4.D0,HALF=0.5D0)
Cstart
CIF(NPROPS.LT.2)THEN
CWRITE(7,*)'
**ERROR:
UMATREQUIRES*NPROPS=2'
CSTOP
CENDIF
E11=PROPS
(1)
V12=PROPS
(2)
G12=PROPS
(1)/TWO/(ONE+PROPS
(2))
CCriticalvaluesofstresses
XT=PROPS(3)
XC=PROPS(4)
XS=PROPS(5)
GX=PROPS(6)!
Fractureenergyinmatrix
ETA=0.001
CCurrentstrain
DOI=1,NTENS
STRANT(I)=STRAN(I)+DSTRAN(I)
ENDDO
CStiffness
DOI=1,6
DOJ=1,6
C(I,J)=ZERO
ENDDO
ATEMP=(1+V12)*(1-TWO*V12)
C(1,1)=E11*(1-V12)/ATEMP
C(2,2)=E11*(1-V12)/ATEMP
C(3,3)=E11*(1-V12)/ATEMP
C(1,2)=E11*V12/ATEMP
C(1,3)=E11*V12/ATEMP
C(2,3)=E11*V12/ATEMP
C(4,4)=G12
C(5,5)=G12
C(6,6)=G12
DOI=2,6
DOJ=1,I-1
C(I,J)=C(J,I)
CCriticalvaluesofstrains
XET=XT/(C(1,1)-2*V12*C(1,2))
XEC=XC/(C(1,1)-2*V12*C(1,2))
XES=XS/C(4,4)
DMOLD=STATEV
(1)
CStraininitiationcriterion
A11=STRANT
(1)**TWO+STRANT
(2)**TWO+STRANT(3)**TWO
A12=A11/XET/XEC
A21=STRANT
(1)+STRANT
(2)+STRANT(3)
A22=(XEC-XET)/XEC/XET*A21
A31=STRANT(4)**TWO+STRANT(5)**TWO+STRANT(6)**TWO
A32=A31/XES**TWO
A1=A12+A22+A32
CB11=STRANT
(2)**TWO
CB12=B11/XET/XEC
CB21=STRANT
(2)
CB22=(XEC-XET)/XEC/XET*B21
CB31=STRANT(5)**TWO
CB32=B31/XES**TWO
CB1=B12+B22+B32
CC11=STRANT(3)**TWO
CC12=C11/XET/XEC
CC21=STRANT(3)
CC22=(XEC-XET)/XEC/XET*C21
CC31=STRANT(6)**TWO
CC32=C31/XES**TWO
CC1=C12+C22+C32
STATEV
(2)=A1
CSTATEV(3)=B1
CSTATEV(4)=C1
FMN=ZERO
IF(A1.GT.ZERO)THEN
FMN=SQRT(A1)
CIF(B1.GT.ONE)THEN
CFMN=FMN+SQRT(B1)
CIF(C1.GT.ONE)THEN
CFMN=FMN+SQRT(C1)
CENDIF
CENDIF
ENDIF
STATEV(5)=FMN
Cwrite(*,*)FMN
DM=ZERO
DDMDFMN=ZERO
DFMNDE(I)=ZERO
DDMDE(I)=ZERO
IF(FMN.GT.ONE)THEN
CCALCULATEDM,DDMDFMN
CWRITE(6,*)FMN
T1=(C(1,1)-2*V12*C(1,2))*XET**2*CELENT/GX
T2=(ONE-FMN)*T1
DM=ONE-EXP(T2)/FMN
CWRITE(6,*)'
T1'
T1,'
T2'
T2,'
DM'
DM
Cwrite(*,*)DM
CCALCULATETHEDERIVATIVEOFDAMAGEVARIABLEWITHRESPECTTOFAILURE
CRITERION
DDMDFMN=(ONE/FMN+T1)*(ONE-DM)
CCALCULATEDFMNDE
IF(DM.GT.DMOLD)THEN
DFMNDE
(1)=HALF/FMN*(TWO*STRANT
(1)+XEC-XET)/XET/XEC
DFMNDE
(2)=HALF/FMN*(TWO*STRANT
(2)+XEC-XET)/XET/XEC
DFMNDE(3)=HALF/FMN*(TWO*STRANT(3)+XEC-XET)/XET/XEC
DFMNDE(4)=ONE/FMN*TWO*STRANT(4)/XES**TWO
DFMNDE(5)=ONE/FMN*TWO*STRANT(5)/XES**TWO
DFMNDE(6)=ONE/FMN*TWO*STRANT(6)/XES**TWO
DOI=1,6
DDMDE(I)=DFMNDE(I)*DDMDFMN
ENDDO
ENDIF
DM=MAX(DM,DMOLD)
Cwrite(6,*)DM
CSAVETHEOLDSTRESSTOOLD_STRESS
OLD_STRESS(I)=STRESS(I)
CEffectivestiffness
CD(I,J)=C(I,J)
IF(DM.NE.ZERO)THEN
CD(1,1)=(ONE-DM)*C(1,1)
CD(1,2)=(ONE-DM)*C(1,2)
CD(2,1)=CD(1,2)
CD(2,2)=(ONE-DM)*C(2,2)
CD(1,3)=(ONE-DM)*C(1,3)
CD(3,1)=CD(1,3)
CD(2,3)=(ONE-DM)*C(2,3)
CD(3,2)=CD(2,3)
CD(4,4)=(ONE-DM)*C(4,4)
CD(5,5)=
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- abaqus 材料 子程序