随动双线性硬化本构模型.docxVIP

2 2 2 2 C ===================================================== C清华大学研究生精品课程《钢筋混凝土有限元》教学程序 C随动双线性硬化本构模型 C Kin ematics hardi ng plastic subrout ine for MSC.MARC C陆新征江见鲸 C LU Xin zhe ng , Jia ng Jianj ing C清华大学土木工程系，北京， 100084 C Dept. Civil En grg. of Tsin ghua Uni versity C last revised: Sep. 2003. C主要参考文献： C 1. MARC V olumn D user subroutine and Special Routines C 2. ABAQUS Writ ing UMATs, VUMA Ts, and UELs C 3.江见鲸”钢筋混凝土结构非线性有限元分析 ” C ===================================================== SUBROUTINE HYPELA(D,G,E,DE,S,TEMP, 1 DTEMP,NGENS,N,NN,KC,MATS,NDI,NSHEAR) c implicit none c* * * * * * c c user subroutine to define youngs modulus and poissons ratio c as function of stress in non-li near elastic small stra in c material. c c d stress strain law to be formed by user c g cha nge in stress due to temperature effects c e total strain c de in creme nt of strain c s stress - should be updated by user c temp state variables c dtemp in creme nt of state variables c ngens size of stress - stra in law c n eleme nt nu mber c nn in tegrati on point nu mber c kc layer nu mber c mats material i.d. c ndi nu mber of direct comp onents c n shear nu mber of shear comp onents c c* * * * * * implicit real*8 (a-h,o-z) INCLUDE ../common/concom ! 通过concom模块得到当前的计算步数 in teger :: ngens,nn ,kc,mats ,n di, nshear real*8 :: e(1),de(1),temp(40),dtemp(40),g(1),d(ngens,ngens),s(1) ! temp(40)需要在前面initial condition 里面设置完成 in teger :: n(2) C* * * * * *C Local Arrays C C EELAS - Elastic Strai ns ! 弹性应变 C EPLAS - Plastic Strai ns ! 塑性应变 C ALPHA - Shift Tensor ! 硬化参数 C Flow - Plasitc Flow Directi ons ! 塑性流动方向 C SIG - Stress at start of in creme nt ! 增量步开始时的应力 C EPSPL - Plastic Strai ns at start of in creme nt ! 增量步开始时的塑性应 变 C real*8 EELAS( nge ns),EPLAS( nge ns),ALPHA( nge ns) 1 ,FLOW( nge ns),SIG( nge ns),EPSPL (n ge ns) c in teger :: inc, in csub, n cycle real*8 E0,ENU,EBULK,EG,ELAM初始弹性模量，泊松比，体积模量，剪切模量 in teger :: K1, K2 ! 循环变量 real*8 :: Smises, SYIELD, hard ! VM 应力，屈服应力，硬化模量 real*8 :: SIGMQEQPL !