邓肯模型-中点增量法.docVIP

下载本文档

15
0
约3.69千字
约 4页
2017-06-08 发布于重庆
举报
版权申诉

邓肯模型-中点增量法.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

邓肯模型-中点增量法

SUBROUTINE UMAT STRESS,STATEV,DDSDDE,SSE,SPD,SCD, 1 RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN, 2 TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,MATERL,NDI,NSHR,NTENS, 3 NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT,CELENT, 4 DFGRD0,DFGRD1,NOEL,NPT,KSLAY,KSPT,KSTEP,KINC C INCLUDE ABA_PARAM.INC C CHARACTER*80 MATERL DIMENSION STRESS NTENS ,STATEV NSTATV , 1 DDSDDE NTENS,NTENS ,DDSDDT NTENS ,DRPLDE NTENS , 2 STRAN NTENS ,DSTRAN NTENS ,TIME 2 ,PREDEF 1 ,DPRED 1 , 3 PROPS NPROPS ,COORDS 3 ,DROT 3,3 , 4 DFGRD0 3,3 ,DFGRD1 3,3 C DIMENSION PS 3 ,DSTRESS NTENS ,TDSTRESS NTENS ,TSTRESS NTENS PARAMETER ONE 1.0D0,TWO 2.0D0,THREE 3.0D0,SIX 6.0D0 EK PROPS 1 EN PROPS 2 RF PROPS 3 C PROPS 4 FAI PROPS 5 /180.0*3.1415926 UG PROPS 6 UD PROPS 7 UF PROPS 8 EKUR PROPS 9 PA PROPS 10 DFAI PROPS 11 /180.0*3.1415926 S1S3O STATEV 1 S3O STATEV 2 SSS STATEV 3 CALL GETPS STRESS,PS,NTENS FAI FAI-DFAI*LOG10 S3O/PA CALL GETEMOD PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,UG,UD,UF 1 ,SSS,S1S3O EBULK3 EMOD/ ONE-TWO*ENU EG2 EMOD/ ONE+ENU EG EG2/TWO EG3 THREE*EG ELAM EBULK3-EG2 /THREE CALL GETDDSDDE DDSDDE,NTENS,NDI,ELAM,EG2,EG TDSTRESS 0.0 CALL GETSTRESS DDSDDE,TDSTRESS,DSTRAN,NTENS DO 701 I1 1,NTENS TSTRESS I1 STRESS I1 +TDSTRESS I1 *0.5 701 CONTINUE CALL GETPS TSTRESS,PS,NTENS CALL GETEMOD PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,UG,UD,UF, 1 SSS,S1S3O EBULK3 EMOD/ ONE-TWO*ENU EG2 EMOD/ ONE+ENU EG EG2/TWO EG3 THREE*EG ELAM EBULK3-EG2 /THREE CALL GETDDSDDE DDSDDE,NTENS,NDI,ELAM,EG2,EG DSTRESS 0.0 CALL GETSTRESS DDSDDE,DSTRESS,DSTRAN,NTENS DO 702 I1 1,NTENS STRESS I1 STRESS I1 +DSTRESS I1 702 CONTINUE CALL GETPS STRESS,PS,NTENS CALL GETEMOD PS,EK,EN,RF,C,FAI,ENU,PA,EKUR,EMOD,S,S3O,UG,UD,UF, 1 SSS,S1S3O EBULK3 EMOD/ ONE-TWO*ENU EG2 EMOD/ ONE+ENU EG EG2/TWO EG3 THREE*EG ELAM EBULK3-EG2 /THREE CALL GETDDSDDE DDSDDE,NTENS,NDI,ELAM,EG2,EG IF PS 3 .GT.S3O S3O PS 3 IF PS 1 -PS 3 .GT.S1S3O S1S3O PS 1 -PS 3 IF S.GT.SSS SSS S STATEV 1 S1S3O STATEV 2 S3O STAT