umat程序代码及注释.docx

  1. 1、本文档共11页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
umat程序代码及注释

SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD, 1 RPL,DDSDDT,DRPLDE,DRPLDT, 2 STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME, 3 NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT, 4 CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)C INCLUDE ABA_PARAM.INCC CHARACTER*80 CMNAME DIMENSION STRESS(NTENS),STATEV(NSTATV), 1 DDSDDE(NTENS,NTENS), 2 DDSDDT(NTENS),DRPLDE(NTENS), 3 STRAN(NTENS),DSTRAN(NTENS),TIME(2),PREDEF(1),DPRED(1), 4 PROPS(NPROPS),COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3) C DIMENSION STRANT(6),CFULL(6,6),CDFULL(6,6),DCDDFV(6,6), 1 DCDDMV(6,6),DCDDDV(6,6),DFFDE(6),DFMDE(6),DFDDE(6), 2 TEMP1(6),TEMP2(6),TEMP3(6),TEMP4(6),TEMP5(6),TEMP6(6), 3 TEMP7(6,6) PARAMETER (ZERO = 0.D0,ONE = 1.D0,TWO = 2.D0,THREE = 3.D0)C 输入材料弹性参数 e1=props(1) e2=props(2) e3=props(3) u12=props(4) u13=props(5) u23=props(6) g12=props(7) g13=props(8) g23=props(9) C 输入材料的强度参数,用于失效判断。 xt=props(10) xc=props(11) yt=props(12) yc=props(13) s12=props(14) s13=props(15) s23=props(16) zt=props(17) zc=props(18) eta=props(19)粘性系数 gc1=props(20)复合材料的1方向的临界能量释放率 gc2=props(21) gc3=props(22)!CFULL 无损伤下的刚度阵 u21=(u12/e1)*e2 u31=(u13/e1)*e3 u32=(u23/e2)*e3 fact=1-u12*u21-u23*u32-u13*u31- 2*u21*u32*u13 factor=fact/(e1*e2*e3) CFULL(1,1)=(1-u23*u32)/(e2*e3*factor) CFULL(2,2)=(1-u13*u31)/(e1*e3*factor) CFULL(3,3)=(1-u12*u21)/(e1*e2*factor) CFULL(1,2)=(u12+u32*u13)/(e1*e3*factor) CFULL(2,3)=(u32+u12*u31)/(e1*e3*factor) CFULL(1,3

文档评论(0)

xjj2017 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档