资源描述：
《随动双线性硬化本构模型》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、C=====================================================C清华大学研究生精品课程《钢筋混凝土有限元》教学程序C随动双线性硬化本构模型CKinematicshardingplasticsubroutineforMSC.MARCC陆新征江见鲸CLUXinzheng,JiangJianjingC清华大学土木工程系，北京，100084CDept.CivilEngrg.ofTsinghuaUniversityClastrevised:Sep.2003.C主要参考文献：C1.MARCVolumnDusersubroutineandSpecialRo
2、utines"C2.ABAQUS"WritingUMATs,VUMATs,andUELs"C3.江见鲸"钢筋混凝土结构非线性有限元分析"C=====================================================SUBROUTINEHYPELA(D,G,E,DE,S,TEMP,1DTEMP,NGENS,N,NN,KC,MATS,NDI,NSHEAR)cimplicitnonec******ccusersubroutinetodefineyoung'smodulusandpoisson'sratiocasfunctionofstressinnon-linea
3、relasticsmallstraincmaterial.ccdstressstrainlawtobeformedbyusercgchangeinstressduetotemperatureeffectscetotalstraincdeincrementofstraincsstress-shouldbeupdatedbyuserctempstatevariablescdtempincrementofstatevariablescngenssizeofstress-strainlawcnelementnumbercnnintegrationpointnumberckclayernumber
4、cmatsmateriali.d.cndinumberofdirectcomponentscnshearnumberofshearcomponentscc******implicitreal*8(a-h,o-z)INCLUDE'../common/concom'!通过concom模块得到当前的计算步数integer::ngens,nn,kc,mats,ndi,nshearreal*8::e(1),de(1),temp(40),dtemp(40),g(1),d(ngens,ngens),s(1)!temp(40)需要在前面initialcondition里面设置完成integer::n(2
5、)C******CLocalArraysC-------------------------------------CEELAS-ElasticStrains!弹性应变CEPLAS-PlasticStrains!塑性应变CALPHA-ShiftTensor!硬化参数CFlow-PlasitcFlowDirections!塑性流动方向CSIG-Stressatstartofincrement!增量步开始时的应力CEPSPL-PlasticStrainsatstartofincrement!增量步开始时的塑性应变Creal*8EELAS(ngens),EPLAS(ngens),ALPHA(n
6、gens)1,FLOW(ngens),SIG(ngens),EPSPL(ngens)cinteger::inc,incsub,ncyclereal*8E0,ENU,EBULK,EG,ELAM!初始弹性模量，泊松比，体积模量，剪切模量integer::K1,K2!循环变量real*8::Smises,SYIELD,hard!VM应力,屈服应力，硬化模量real*8::SIGM,DEQPL!平均正应力,等效应变增量real*8::effg,efflam,effhard!等效剪切模量，等效lame常数，等效硬化模量Creal*8,parameter::ENUMAX=.49999D0,TOLER
7、=1.D-6CCtemp(1)-E!弹性模量Ctemp(2)-NU!泊松比Ctemp(3)-SYIELD!屈服强度Ctemp(4)-HARD!硬化模量C------------------------------C初始状态赋值Cif(inc==0.and.incsub==0.and.ncycle==0)thenTEMP(1)=200e3;!弹性模量200GPaTEMP(2)=0.27;!泊松比为0.27TEMP(3)=210.0;!
