壳厚度方向一点进入塑性，都有什么途径修正应力？

吃草的鸟

当一点材料进入塑性，如壳是平面应力问题，厚度方向一个积分点进入塑性，有几种选择修正应力啊？我编程的结果ＭＩＳＥＳ和软件对不上呢？具体流程如下：
ｄσ＝Ｃｄε
σ（ｎ＋１）＝σ（ｎ）＋ｄσ
σ（Ｍｉｓｅｓ）＝sqrt(σｔ｛１，－０．５，０；－０．５，１，０；０，０，３｝σ)
σｍ＝（σｘｘ＋σｙｙ＋σｚｚ）／３
Ｓ＝σ（ｎ＋１）－σｍＩ　！偏应力
ｉｆ（σ（Ｍｉｓｅｓ）＞σ（ｙｉｅｌｄ））
Ｓ＝Ｓ＊σ（ｙｉｅｌｄ）／σ（Ｍｉｓｅｓ）
ｅｎｄｉｆ
σ（ｎ＋１）＝σｍＩ＋Ｓ
各位大师，谁有相关的资料可以给点吗？还有其他的，利用迭代求解的，希望各位大师能多多指点，此问题困扰小弟好几天了！

tonnyw

作为入门可以先看看王瑁成的书是怎么谈论弹塑性的。对于编程，想看一下NONSAP
是如何处理平面应力/应变问题的弹塑性的。这两项都了解了，在来解决壳单元的问题。

tianyikillua

你所说的"修正应力"是什么意思？意思是说考虑材料的弹塑性响应？

吃草的鸟

tianyikillua 发表于 2013-4-28 19:38
你所说的"修正应力"是什么意思？意思是说考虑材料的弹塑性响应？

就是每一时间步开始进行弹性预测，如果材料超过屈服极限，此时，应力怎么处理？如下过程，我编的过程，但结果和软件对不上，弹性阶段没有问题的，但是材料进入塑性后，Mises曲线就开始对不上了~~~不解为什么？


ｄσ＝Ｃｄε
σ（ｎ＋１）＝σ（ｎ）＋ｄσ
σ（Ｍｉｓｅｓ）＝sqrt(σｔ｛１，－０．５，０；－０．５，１，０；０，０，３｝σ)
σｍ＝（σｘｘ＋σｙｙ＋σｚｚ）／３
Ｓ＝σ（ｎ＋１）－σｍＩ　！偏应力
ｉｆ（σ（Ｍｉｓｅｓ）＞σ（ｙｉｅｌｄ））
Ｓ＝Ｓ＊σ（ｙｉｅｌｄ）／σ（Ｍｉｓｅｓ）
ｅｎｄｉｆ
σ（ｎ＋１）＝σｍＩ＋Ｓ

tianyikillua

...啊我明白你要做什么了...

你用什么模型？Perfect plasticity？Isotropic 还是 Kinematic hardening？塑性的修正没那么容易（至少4~5行代码吧）...因为需要更新塑性形变，还要根据你硬化的类型（线性，非线性）解一个方程，如果非线性硬化还要用比如牛顿法解一个非线性方程...

吃草的鸟

tianyikillua 发表于 2013-4-29 00:55
...啊我明白你要做什么了...

你用什么模型？Perfect plasticity？Isotropic 还是 Kinematic hardening？塑 ...

我起初用的理想弹性-塑性线性强化，就是双线性，但是因为结果和软件对不上，就换成理想弹塑性了！可是结果依然不对！感觉这个已经很简单了！我给你参数啊，密度=7850，泊松比=0.3，屈服强度=235MPa,弹性模量=2.1e11!就是几个参数，理想弹塑性！测试模型为方形板边长0.5m，四周刚性固定，施加恒定的均布阶跃载荷1MPa！

吃草的鸟

tianyikillua 发表于 2013-4-29 00:55
...啊我明白你要做什么了...

你用什么模型？Perfect plasticity？Isotropic 还是 Kinematic hardening？塑 ...

方板中心点上表面的Mises应力曲线，感觉就是塑性应力修正有问题！

hillyuan

1.  Ｓ＝Ｓ＊σ（ｙｉｅｌｄ）／σ（Ｍｉｓｅｓ） wrong! We need do something concerning radial return mapping.
   And
2.  To make sure sigma(33)=0 for shell elements afterwards

吃草的鸟

hillyuan 发表于 2013-4-30 12:15
1.  Ｓ＝Ｓ＊σ（ｙｉｅｌｄ）／σ（Ｍｉｓｅｓ） wrong! We need do something concerning radial return  ...

我的程序，请帮忙！http://forum.simwe.com/forum.php ... p;extra=#pid3013979

吃草的鸟

tonnyw 发表于 2013-4-28 08:18
作为入门可以先看看王瑁成的书是怎么谈论弹塑性的。对于编程，想看一下NONSAP
是如何处理平面应力/应变问题 ...

谢谢版主你一直以来的帮助！我的问题的详细说明，请帮忙！http://forum.simwe.com/forum.php ... p;extra=#pid3013979