第五章刚塑性有限元法基本理论与模拟方法.pptVIP

对于轴对称问题，当单元采用等参四边形单元时，式5-81中的 有下列展开式： 式中： (5-82) (5-83) (5-84) (5-85) 式中：k?2、k?1、k、k+1为下标i、j、k、m 的代号，0为下标时表示平均值。 B、C、D 分别由式2-28和式2-36中对应A 、B 、G 的给出。 5.5 罚函数法 5.5.1 求解方程的建立 刚塑性有限元法的一个基本假设是体积不变，罚函数法从这一点入手，引入一个很大的正数乘以体积应变速率的平方，即： (5-86) 将此项添加到泛函中去，于是可将式 5-11 改写如下： (5-87) 当 在每一点的变化率都接近零时，这个泛函就取得最小值，这时所对应的速度场就逼近真实速度场。可以按照处理拉格朗日乘子法同样的步骤导出罚函数法的计算公式。 罚函数法与拉格朗日乘子法的不同之处在最后一项，下面只对这一项进行分析。令这一项为 ，即：