共轭双线性硬化本构模型.docVIP

C =====================================================C 清华大学研究生精品课程《钢筋混凝土有限元》教学程序C 随动双线性硬化本构模型 C Kinematics harding plastic subroutine for HYPERLINK MSC.MARCC HYPERLINK \t _blank 陆新征 HYPERLINK /TCHSchool/%7Ejiangjj_c.htm 江见鲸C HYPERLINK \t _blank LU Xinzheng, Jiang JianjingC HYPERLINK /TCHSchool/index.asp 清华大学土木工程系，北京，100084C HYPERLINK /TCHSchool/index.asp Dept. Civil Engrg. of Tsinghua UniversityC last revised: Sep. 2003.C 主要参考文献：C 1. MARC Volumn D user subroutine and Special RoutinesC 2. ABAQUS Writing UMATs, VUMATs, and UELsC 3. 江见鲸 钢筋混凝土结构非线性有限元分析C ===================================================== SUBROUTINE HYPELA(D,G,E,DE,S,TEMP,1 DTEMP,NGENS,N,NN,KC,MATS,NDI,NSHEAR)c implicit nonec* * * * * *cc user subroutine to define youngs modulus and poissons ratioc as function of stress in non-linear elastic small strainc material.cc d stress strain law to be formed by userc g change in stress due to temperature effectsc e total strainc de increment of strainc s stress - should be updated by userc temp state variablesc dtemp increment of state variablesc ngens size of stress - strain lawc n element numberc nn integration point numberc kc layer numberc mats material i.d.c ndi number of direct componentsc nshear number of shear componentscc* * * * * *implicit real*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)需要在前面initial condition 里面设置完成integer :: n(2) C* * * * * *C Local ArraysC-------------------------------------C EELAS - Elastic Strains ! 弹性应变C EPLAS - Plastic Strains ! 塑性应变C ALPHA - Shift Tensor ! 硬化参数C Flow - Plasitc Flow Directions ! 塑性流动方向C SIG - Stress at start of increment ! 增量步开始时的应力C EPSPL - Plastic Strains at start of increment ! 增量步开始时的塑性应变Creal*8 EELAS(ngens),EPLAS(ngens),ALPHA(ngens)1 ,FLOW(ngens),SIG(ngens),EPSPL(ngens)c integer :: inc, incsub, ncyclereal*8 E0,ENU