Abaqus 采用YLD2000-2D屈服准则的UMAT子程序
2020年5月23日 浏览:2192 评论:12 收藏:3
Barlat在2003年提出了专门针对平面应力问题的各向异性屈服准则,该屈服准则对于各向异性材料具有很高的精度,得到了广泛的应用。
YLD2000-2D屈服面示意图
Yld2000-2d屈服准则由下式给出
![1.png 1.png](https://img.jishulink.com/upload/202005/d09ebf7e85054d3cbb3f855eabc08a51.png)
其中
![2.png 2.png](https://img.jishulink.com/upload/202005/046b97da7b944074922a37dfcc0e109a.png)
矩阵X′和X″的元素分别由柯西应力的下列线性变换获得
![3.png 3.png](https://img.jishulink.com/upload/202005/783f8b5181b842e2b8c01d7a7ec660de.png)
L′和L″的分量由下式求得
![45.png 45.png](https://img.jishulink.com/upload/202005/a0f9625a82f1406c97519721e721643c.png)
积分算法采用径向返回算法,该方法是稳健和精确的。
![6.png 6.png](https://img.jishulink.com/upload/202005/9203f5614ddf4e7da75797472ba7aea4.png)
![7.png 7.png](https://img.jishulink.com/upload/202005/111ccbc02bbf418089087f04aa9e6bcd.png)
当弹性试算超出屈服面时,则需要进行塑性修正
![8.png 8.png](https://img.jishulink.com/upload/202005/206a9ceb7f6346c28690103040afc312.png)
使其满足
![9.png 9.png](https://img.jishulink.com/upload/202005/4bf1bafab5174cc6a2840e41b7f327d3.png)
公式9可以通过牛顿法进行迭代求解。
计算的应力应变曲线如下图所示
![YLD2000.gif YLD2000.gif](https://img.jishulink.com/upload/202005/32713889765d4c0dbda014fe0feb5e2c.gif)
B, F. Barlat A , et al. "Plane stress yield function for aluminum alloy sheets—part 1: theory." International Journal of Plasticity 19. 9(2003):1297-1319.
王海波, 万敏, 阎昱,等. 屈服准则在有限元软件中实现的正确性验证[J]. 固体力学学报, 2010, 031(002):173-180.
![无标题.png 无标题.png](https://img.jishulink.com/upload/202005/53a7c559356040829ce826c7a2c554f4.png)
点赞 19 评论 12 收藏 3
查看更多评论 >