资源描述：
《MSC.MARC等强硬化弹塑性本构源程序》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、MSC.MARC等强硬化弹塑性本构源程序C=====================================================C清华大学研究生精品课程《钢筋混凝土有限元》教学程序C等强双线性硬化本构模型C各向同性哈丁塑料子程序模拟C陆新征江见鲸C清华大学土木工程系，北京，100084C路新郑江见鲸，土木工程系。清华大学c最后修订：2009年9月2003oC主要参考文献：CE马克体积D用户子程序和特殊程序”C2。江见鲸”钢筋混凝土结构非线性有限元分析”CC==========================
2、===========================c子程序说明：C：得到等效应力零C：得到当前屈服应力产量Cwkslp:得到强化模量C公司：得到流动方向子公司(弹起，辛格，SC,T,ngens,NDI,nshear,N,NNN,层)C***Cc提供通用可塑性选项的流规则。C流动法则CC随机应力阵列！应力矩阵C从流动方向阵列！流动向量静水压力！静水压力未使用Cngens应力分量的数目CNDI(3)直接应力分量数Cnshear(3)数量的剪应力分量n元素数CNNN积分点数层层数cc***隐房X8（A-H0-Z）DP维斯托特（*
3、）,辛格（*）26我二1,NDI26从（我）二0.5*（3*斯托特（I）-SC）如果（nshearoeq.0）去30F3=3o28我二1,nshearJ二NDI+我28止弦（j）二F3*斯托特（J）30继续返回结束子程序wkslp（M、NN、KC、垫、坡度、ebarp,eqrate,stryt,DT,笫一次）隐房X8（A-H0-Z）DPC***cc用户子程序定义加工硬化斜率。CC硬化了程序c是当前的用户元素号。神经网络是积分点数C是层数（非分层元素=1）。C垫子是当前的材质idc边坡工作硬化斜坡应在此例程中定义Cebarp是等
4、效塑性应变Ceqrate是等效塑性应变率Cstryt电流屈服应力，也可以在这个程序的定义是当前的总温度。C第一标志区分第十个周期特性CORNL选项C可以获得内部元素号薄荷。C薄荷二ielint(M)CC***20e3斜率！强化模量为ES二2gpastryt二210+ebarp*20e3!强化后的屈服应力为210mpa+eps_eq_pl*ES返回结束真止的8功能分子(M,NN,KC,你，第一次，DT,eplas,同时，该,*用来存储）C***C得到屈服应力C求屈服应力的电流值C碳元素指数C你初始屈服应力C笫一旗笫十周期产量CD
5、T温度Ceplas等效塑性应变C,等效塑性应变率C用来存储表携产量CC***隐房X8(A-H0-Z)c***cc效率二你ccC用户子程序wkslpC调用做功强化了程序得到强化模量和屈服应力Cstryt二0.dO电话wkslp(M、NN、KC、垫、坡度、eplas,,stryt,DT,第一次)效率二stryt返回结束真正的8函数零点(NDI,nshear,T,IORT,ianiso,yrdir,yrshr,AMM,AO)c***cc得到等效应力求应力张量的等效应力。C应力的直接成分CNDI数Cnshear数量的剪应力分量c(I)
6、应力分量IC中说明如果使用曲线坐标C二0放疗没有使用曲线坐标C二1T获得IORT混合组件C=2假定IORT已经混合compcomponentsCianiso标志表示如果各向异性是用Hill各向异性塑性Cyrdir直接成分Cyrshr剪切H订1各向杲性塑性元件CAMM的度量如果用曲线坐标使用曲线坐标时的公制标度因子CC***隐房X8(A-H0-Z)尺寸t(*)尺寸S(6),yrdir(3)、(3)、yrshrAMM(3)零二0.d0CC等效应力即为von米塞斯应力C如果(NDI。公式1)然后零二T(1)*T(1)*2.d0如果(
7、NDIo水利)然后零二2oOX(t(1)*T(1)+(2)XT(2)-(1)*T(2))如果(NDI。式)然后零二(1)*t(1)t(2)*t(2)t(3)*t(3)*T(1)*t(2)T(2)*T(3)T(3)*T(1)零二2oOX零endif如果(nshear。公式1)然后零二零+6oOXT(NDI+1)*T(NDI+1)如果(nshearo水利）然后零-零+6oOX(NDI+2))(T(NDI+1)*T(NDI+1)+T(NDI+2)*T如果(nshearo式）然后零-零+6oOX(T(NDI+1)*T(NDI+1)*T
8、(NDI+2)*T(NDI+2)+T(NDI+3)*T(NDI+3))endifC如果（零。It.OoDO）零二0.dO零二SQRT(0.5d0*零)返回结束
