MSC.MARC 等强硬化弹塑性本构源程序(国外英语资料).doc

MSC.MARC 等强硬化弹塑性本构源程序 C ===================================================== C清华大学研究生精品课程《钢筋混凝土有限元》教学程序 C等强双线性硬化本构模型 C各向同性哈丁塑料子程序模拟 C陆新征江见鲸 C清华大学土木工程系，北京，100084 C路新郑江见鲸，土木工程系。清华大学 c最后修订：2009年9月2003。 C主要参考文献： C 1。马克体积D用户子程序和特殊程序” C 2。江见鲸”钢筋混凝土结构非线性有限元分析” C C ===================================================== C子程序说明： C：得到等效应力零 C：得到当前屈服应力产量 C wkslp：得到强化模量 C公司：得到流动方向 子公司（弹起，辛格，SC，T，ngens，NDI，nshear，N，NNN，层） C *** C c提供通用可塑性选项的流规则。 C流动法则 C C随机应力阵列！应力矩阵 C从流动方向阵列！流动向量 静水压力！静水压力 未使用 C ngens应力分量的数目 C NDI（3）直接应力分量数 C nshear（3）数量的剪应力分量 n元素数 C NNN积分点数 层层数 C C *** 隐房×8（A-H O-Z）DP 维斯托特（*），辛格（*） 26我= 1，NDI 26从（我）= 0.5 *（3 *斯托特（I）- SC） 如果（nshear。eq.0）去30 F3＝3。 28我= 1，nshear J = NDI +我 28正弦（j）= F3 *斯托特（J） 30继续 返回 结束 子程序wkslp（M、NN、KC、垫、坡度、ebarp，eqrate，stryt，DT，第一次） 隐房×8（A-H O-Z）DP C *** C c用户子程序定义加工硬化斜率。 C C硬化子程序 c是当前的用户元素号。 神经网络是积分点数 C是层数（非分层元素＝1）。 C垫子是当前的材质id c边坡工作硬化斜坡应在此例程中定义 C ebarp是等效塑性应变 C eqrate是等效塑性应变率 C stryt电流屈服应力，也可以在这个程序的定义 是当前的总温度。 C第一标志区分第十个周期特性 C ORNL选项 C 可以获得内部元素号薄荷。 C薄荷= ielint（M） C C *** 20e3斜率！强化模量为ES = 2gpa stryt = 210 + ebarp * 20e3！强化后的屈服应力为210mpa + eps_eq_pl * ES 返回 结束 真正的8功能分子（M，NN，KC，你，第一次，DT，eplas，同时，该， *用来存储） C *** C得到屈服应力 C 求屈服应力的电流值 C 碳元素指数 C你初始屈服应力 C第一旗第十周期产量 C DT温度 C eplas等效塑性应变 C，等效塑性应变率 C用来存储表携产量 C C *** 隐房×8（A-H O-Z） C *** C C 效率=你 C C C用户子程序wkslp C调用做功强化子程序得到强化模量和屈服应力 C stryt = 0.d0 电话wkslp（M、NN、KC、垫、坡度、eplas，，stryt，DT，第一次） 效率= stryt 返回 结束 真正的8函数零点（NDI，nshear，T，IORT，ianiso，yrdir，yrshr，AMM，A0） C *** C C得到等效应力 求应力张量的等效应力。 C 应力的直接成分C NDI数 C nshear数量的剪应力分量 c（I）应力分量I C中说明如果使用曲线坐标 C = 0放疗没有使用曲线坐标 C = 1 T获得IORT混合组件 C = 2假定IORT已经混合compcomponents C ianiso标志表示如果各向异性是用 Hill各向异性塑性C yrdir直接成分 C yrshr剪切Hill各向异性塑性元件 C AMM的度量如果用曲线坐标 使用曲线坐标时的公制标度因子 C C *** 隐房×8（A-H O-Z） 尺寸t（*） 尺寸S（6），yrdir（3）、（3）、yrshr AMM（3） 零= 0.d0 C C等效应力即为von米塞斯应力 C 如果（NDI。公式1）然后 零= T（1）* T（1）* 2.d0 如果（NDI。水利）然后 零= 2。0×（t（1）* T（1）+（2）×T（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） 零= 2。0×零 endif 如果（nshear。公式1）然后 零=零+ 6。0×T（NDI + 1）* T（NDI + 1） 如果