鄧肯張模型FORTRAN子程序源代碼_第1頁
鄧肯張模型FORTRAN子程序源代碼_第2頁
鄧肯張模型FORTRAN子程序源代碼_第3頁
鄧肯張模型FORTRAN子程序源代碼_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

1、鄧肯張模型FORTRAN子程序源代碼SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,MATERL,NDI,NSHR,NTENS,NSTATVPROPS,NPROPS,COORDS,DROT,PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,KSLAYKSPT,KSTEP,KINC)CINCLUDE ABA_PARAM.INCCCHARACTER*80 MATERLDIMENSI

2、ON STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)CDIMENSION PS(3),DSTRESS(NTENS),TDSTRESS(NTENS),TSTRESS(NTENS)PARAMETER (ONE=1.0D0,TWO=2.0D0,THREE=3.0D0,SIX=

3、6.0D0)K=PROPS(1)N=PROPS(2)RF=PROPS(3)C=PROPS(4)FAI=PROPS(5)/180.0*3.1415926G=PROPS(6)D=PROPS(7)F=PROPS(8)KUR=PROPS(9)PA=PROPS(10)DFAI=PROPS(11)/180.0*3.1415926S1S3O=STATEV(1)S3O=STATEV(2)SSS=STATE V(3)CALL GETPS(STRESS,PS,NTENS)FAI=FAI-DFAI*LOG10(S3O/PA)CALL GETEMOD(PS,K,N,RF,C,FAI,ENU,PA,KUR,EMOD,S

4、,S3O,G,D,F1 ,SSS,S1S3O)EBULK3=EMOD/(ONE-TWO*ENU)EG2=EMOD/(ONE+ENU)EG=EG2/TWOEG3=THREE*EGELAM=(EBULK3-EG2)/THREECALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)TDSTRESS=0.0CALL GETSTRESS(DDSDDE,TDSTRESS,DSTRAN,NTENS)DO 701 I1=1,NTENSTSTRESS(I1)=STRESS(I1)+TDSTRESS(I1)*0.5CONTINUECALL GETPS(TSTRESS,PS,NTE

5、NS)CALL GETEMOD(PS,K,N,RF,C,FAI,ENU,PA,KUR,EMOD,S,S3O,G,D,F,1 SSS,S1S3O)EBULK3=EMOD/(ONE-TWO*ENU)EG2=EMOD/(ONE+ENU)EG=EG2/TWOEG3=THREE*EGELAM=(EBULK3-EG2)/THREECALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)DSTRESS=0.0CALL GETSTRESS(DDSDDE,DSTRESS,DSTRAN,NTENS)DO 702 I1=1,NTENSSTRESS(I1)=STRESS(I1)+DST

6、RESS(I1)CONTINUECALL GETPS(STRESS,PS,NTENS)CALL GETEMOD(PS,K,N,RF,C,FAI,ENU,PA,KUR,EMOD,S,S3O,G,D,F,1 SSS,S1S3O)EBULK3=EMOD/(ONE-TWO*ENU)EG2=EMOD/(ONE+ENU)EG=EG2/TWOEG3=THREE*EGELAM=(EBULK3-EG2)/THREECALL GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)IF(PS(3).GT.S3O)S3O=PS(3)IF(PS(1)-PS(3).GT.S1S3O)S1S3O=P

7、S(1)-PS(3)IF(S.GT.SSS)SSS=SSTATEV(1)=S1S3OSTATEV(2)=S3OSTATE V(3)=SSSENDSUBROUTINE GETPS(STRESS,PS,NTENS)INCLUDE ABA_PARAM.INCDIMENSION PS(3),STRESS(NTENS)CALL SPRINC(STRESS,PS,1,3,3)DO 310 I=1,2DO 320 J=I+1,3IF(PS(I).GT.PS(J)THENPPS=PS(I)PS(I)=PS(J)PS(J)=PPSEND IF320 CONTINUE310 CONTINUEDO 330 K1=1

8、,3PS(K1)=-PS(K1)330 CONTINUERETURNENDSUBROUTINE GETEMOD(PS,K,EN,RF,C,FAI,ENU,PA,KUR,EMOD,S,S3O1 ,G,D,F,SSS,S1S3O)INCLUDE ABA_PARAM.INCDIMENSION PS(3)S=(1-SIN(FAI)*(PS(1)-PS(3)IF(PS(3).LT.0) THENPSFEI=0.1ELSEPSFEI=PS (3)END IFS=S/(2*C*COS(FAI)+2*PSFEI*SIN(FAI)IF(S.GE.0.99) THENS=0.99END IFAA=D*(PS(1)

9、-PS(3)AA=AA/(K*PA*(S3O/PA)*N)AA=AA/(1-RF*S)ENU=G-F*LOG10(S3O/PA)ENU=ENU/(1-AA)/(1-AA)IF(ENU.GT.0.49)ENU=0.49EMOD=K*PA*(S3O/PA)*N)*(1-RF*S)*2)IF(S.LT.SSS.AND.(PS(1)-PS(3).LT.S1S3O)THENEMOD=KUR*PA*(S3O/PA)*N)END IFENDSUBROUTINE GETFEI(STRESS,SMISES,NTENS,NDI)INCLUDE ABA_PARAM.INCDIMENSION STRESS(NTENS

10、)SMISES=STRESS(1)*2+STRESS(2)*2+STRESS(3)*2+1STRESS(4)*2+STRESS(5)*2+STRESS(6)*2SMISES=SQRT(SMISES)RETURNENDSUBROUTINE GETDDSDDE(DDSDDE,NTENS,NDI,ELAM,EG2,EG)INCLUDE ABA_PARAM.INCDIMENSION DDSDDE(NTENS,NTENS)DO 20 K1=1,NTENSDO 10 K2=1,NTENSDDSDDE(K2,K1)=0.0CONTINUECONTINUEDO 40 K1=1,NDIDO 30 K2=1,NDIDDSDDE(K2,K1)=ELAMCONTINUEDDSDDE(K1,K1)=EG2+ELAMCONTINUEDO 50 K1=NDI+1,NTENSDDSDDE(K1,K1)=EGCONTINUERETURNENDSUBROUTINE GETSTRESS(DDSDDE,STRESS,DSTRAN,NTENS)INCLUDE ABA_PARAM.INCDIMENSI

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽，若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論