船舶力学

角度/°	屈服强度$ {\sigma }_{\mathrm{s}}$ /MPa	抗拉强度$ {\sigma }_{\mathrm{b}} $ /MPa	伸长率A /%	厚向异性系数R
0	348.9	636.0	31.5	0.887
15	348.4	628.8	30.6	0.868
30	354.4	641.3	31.4	0.897
45	361.3	634.4	31.3	0.945
60	357.1	652.5	31.8	0.847
75	357.0	639.0	30.8	0.808
90	359.7	639.3	32.9	0.859

角度/°	屈服强度$ {\sigma }_{\mathrm{s}}$ /MPa	抗拉强度$ {\sigma }_{\mathrm{b}} $ /MPa	伸长率A /%	厚向异性系数R
0	348.9	636.0	31.5	0.887
15	348.4	628.8	30.6	0.868
30	354.4	641.3	31.4	0.897
45	361.3	634.4	31.3	0.945
60	357.1	652.5	31.8	0.847
75	357.0	639.0	30.8	0.808
90	359.7	639.3	32.9	0.859

参数	$ {\sigma }_{0} $	K	Q	b	a
Hockett-Sherby+Linear	983.835	131.570	7.387	0.663	−150.245

参数	$ {\sigma }_{0} $	K	Q	b	a
Hockett-Sherby+Linear	983.835	131.570	7.387	0.663	−150.245

	F	G	H	L	M	N
Hill48-S	0.4889	0.5479	0.4521	1.5	1.5	1.3472
Hill48-R	0.5472	0.5299	0.4701	1.5	1.5	1.5565

	F	G	H	L	M	N
Hill48-S	0.4889	0.5479	0.4521	1.5	1.5	1.3472
Hill48-R	0.5472	0.5299	0.4701	1.5	1.5	1.5565

$ \alpha 1 $	$ \alpha 2 $	$ \alpha 3 $	$ \alpha 4 $	$ \alpha 5 $	$ \alpha 6 $	$ \alpha 7 $	$ \alpha 8 $
1.069	0.918	1.186	1.025	1.073	1.355	0.981	0.773

$ \alpha 1 $	$ \alpha 2 $	$ \alpha 3 $	$ \alpha 4 $	$ \alpha 5 $	$ \alpha 6 $	$ \alpha 7 $	$ \alpha 8 $
1.069	0.918	1.186	1.025	1.073	1.355	0.981	0.773

参数	$ c_{12}^{'} $	$ c_{13}^{'} $	$ c_{21}^{'} $	$ c_{23}^{'} $	$ c_{31}^{'} $	$ c_{32}^{'} $	$ c_{44}^{'} $	$ c_{55}^{'} $	$ c_{66}^{'} $
实验值	0.8689	0.7209	0.5009	0.9520	0.6541	1.2411	0.9357	1	1

参数	$ c_{12}^{''} $	$ c_{13}^{''} $	$ c_{21}^{''} $	$ c_{23}^{''} $	$ c_{31}^{''} $	$ c_{32}^{''} $	$ c_{44}^{''} $	$ c_{55}^{''} $	$ c_{66}^{''} $
实验值	0.6121	1.1445	0.9676	0.7189	1.2240	0.4174	0.9854	1	1

参数	$ c_{12}^{'} $	$ c_{13}^{'} $	$ c_{21}^{'} $	$ c_{23}^{'} $	$ c_{31}^{'} $	$ c_{32}^{'} $	$ c_{44}^{'} $	$ c_{55}^{'} $	$ c_{66}^{'} $
实验值	0.8689	0.7209	0.5009	0.9520	0.6541	1.2411	0.9357	1	1

参数	$ c_{12}^{''} $	$ c_{13}^{''} $	$ c_{21}^{''} $	$ c_{23}^{''} $	$ c_{31}^{''} $	$ c_{32}^{''} $	$ c_{44}^{''} $	$ c_{55}^{''} $	$ c_{66}^{''} $
实验值	0.6121	1.1445	0.9676	0.7189	1.2240	0.4174	0.9854	1	1

方向/$ {^{\circ}} $	实验值	Hill48-S	Hill48-R	Yld2000-2d	Yld2004-18p
0	1	0.9903	1.0156	1.0004	0.9990
15	0.9987	1.0030	1.0101	1.0068	0.9997
30	1.0158	1.0184	0.9990	1.0217	1.0170
45	1.0354	1.0316	0.9923	1.0351	1.0331
60	1.0236	1.0345	0.9949	1.0388	1.0248
75	1.0233	1.0301	1.0027	1.0345	1.0244
90	1.0309	1.0271	1.0070	1.0314	1.0298

方向/$ {^{\circ}} $	实验值	Hill48-S	Hill48-R	Yld2000-2d	Yld2004-18p
0	1	0.9903	1.0156	1.0004	0.9990
15	0.9987	1.0030	1.0101	1.0068	0.9997
30	1.0158	1.0184	0.9990	1.0217	1.0170
45	1.0354	1.0316	0.9923	1.0351	1.0331
60	1.0236	1.0345	0.9949	1.0388	1.0248
75	1.0233	1.0301	1.0027	1.0345	1.0244
90	1.0309	1.0271	1.0070	1.0314	1.0298

方向/$ {^{\circ}} $	实验值	Hill48-S	Hill48-R	Yld2000-2d	Yld2004-18p
0	0.887	0.825	0.887	0.888	0.887
15	0.868	0.814	0.903	0.902	0.862
30	0.897	0.795	0.934	0.931	0.893
45	0.945	0.799	0.945	0.948	0.957
60	0.847	0.842	0.920	0.923	0.852
75	0.808	0.898	0.879	0.880	0.814
90	0.859	0.925	0.859	0.861	0.858

方向/$ {^{\circ}} $	实验值	Hill48-S	Hill48-R	Yld2000-2d	Yld2004-18p
0	0.887	0.825	0.887	0.888	0.887
15	0.868	0.814	0.903	0.902	0.862
30	0.897	0.795	0.934	0.931	0.893
45	0.945	0.799	0.945	0.948	0.957
60	0.847	0.842	0.920	0.923	0.852
75	0.808	0.898	0.879	0.880	0.814
90	0.859	0.925	0.859	0.861	0.858

节点数	单元数	内存	软件	操作系统	处理器
18227	15360	8 GB	ABAQUS/Explicit	64-bit Windows 10	Intel（R） Core（TM） i5-12600KF 3.70 GHz

节点数	单元数	内存	软件	操作系统	处理器
18227	15360	8 GB	ABAQUS/Explicit	64-bit Windows 10	Intel（R） Core（TM） i5-12600KF 3.70 GHz

弯曲半径	屈服模型
$ r $ = 500 mm	4 min 20 s	7 min 28 s	10 min 17 s
$ r $ = 200 mm	4 min 47 s	8 min 03 s	11 min 16 s

弯曲半径	屈服模型
$ r $ = 500 mm	4 min 20 s	7 min 28 s	10 min 17 s
$ r $ = 200 mm	4 min 47 s	8 min 03 s	11 min 16 s

E40钢各向异性行为的实验和模拟研究

PDF下载

蔡一杰 ¹^,² , 张文杰 ¹^,² , 张文谦 ¹^,² , 施凤燕 ³ , 胡勇 ³ , 杨青云 ⁴

船舶力学 | 结构力学 2026,30(5): 804-815

收起

船舶力学 |结构力学 2026 , 30 (5) : 804 -815

E40钢各向异性行为的实验和模拟研究

全屏

蔡一杰¹^,², 张文杰¹^,², 张文谦¹^,², 施凤燕³, 胡勇³, 杨青云⁴

作者信息

^1.湖北工业大学机械工程学院，武汉 430068

^2.现代制造质量工程湖北省重点实验室，武汉 430068

^3.武汉理工大学船海与能源动力工程学院，武汉 430063

^4.中国船舶集团有限公司第七二五研究所，河南洛阳 471023

通讯作者:

张文谦（1990–），男，博士，副教授，通讯作者，E-mail： wenqian_zh@hbut.edu.cn

作者简介:

蔡一杰（1987–），男，博士，副教授，E-mail：yijie.cai@hbut.edu.cn

Experimental and simulation studies on the anisotropic behavior of E40 steel

Yi-jie CAI¹^,², Wen-jie ZHANG¹^,², Wen-qian ZHANG¹^,², Feng-yan SHI³, Yong HU³, Qing-yun YANG⁴

Affiliations

^1.School of Mechanical Engineering, Hubei University of Technology, Wuhan 430068, China

^2.Key Laboratory of Modern Manufacturing Quality Engineering in Hubei Province, Wuhan 430068, China

^3.Ocean and Energy Power Engineering, Wuhan University of Technology, Wuhan 430063, China

^4.Luoyang Ship Material Research Institute, Luoyang 471023, China

出版时间: 2026-05-15 doi: 10.3969/j.issn.1007-7294.2026.05.012

文章导航

摘要

收起

本文通过单向拉伸实验研究了E40钢在拉伸状态下的塑性变形行为，结果表明其具有显著的各向异性特征。同时基于实验结果，对Hill48、Yld2000-2d与Yld2004-18p三种屈服模型进行了参数标定，并采用修正型Hockett-Sherby模型描述其应变硬化行为。最后基于ABAQUS软件开发了VUMAT子程序，结合板料冲压实验系统评估不同屈服模型的工程适用性。数值模拟与实验结果对比显示，在弯曲半径$ r $ = 200 mm和$ r $ = 500 mm的工况下，Yld2004-18p模型的预测误差小于4%，其预测精度较Hill48模型提高约3倍。研究表明，Yld2004-18p模型能更准确地刻画E40钢的各向异性响应，为船舶构件精密成形仿真提供了有效的本构建模策略。

关键词

屈服模型 / 各向异性 / E40钢 / 冲压成形

Abstract

收起

The plastic deformation behavior of E40 steel in the tensile state was investigated by unidirectional tensile experiments, and the results showed that it has significant anisotropic characteristics. Based on the experimental results, three yield models, Hill48, Yld2000-2d, and Yld2004-18p, were parametrically calibrated, and their strain hardening behaviors were described using a modified Hockett-Sherby model. A VUMAT subroutine was developed and implemented using ABAQUS software and combined with a sheet metal stamping experimental system to evaluate the engineering applicability of different yield models. Comparison between numerical simulations and experimental results shows that under bending radii of $ r $ = 200 mm and $ r $ = 500 mm, the prediction error of the Yld2004-18p model is less than 4%, with its prediction accuracy approximately three times higher than that of the Hill48 model. The study shows that the Yld2004-18p model can more accurately portray the anisotropic response of E40 steel, which provides an effective constitutive modeling strategy for the simulation of precision forming of ship components.

Key words

yield model / anisotropy / E40 steel / stamping forming

引用本文

蔡一杰, 张文杰, 张文谦, 施凤燕, 胡勇, 杨青云. E40钢各向异性行为的实验和模拟研究. 船舶力学, 2026 , 30 (5) : 804 -815 . DOI: 10.3969/j.issn.1007-7294.2026.05.012

Yi-jie CAI, Wen-jie ZHANG, Wen-qian ZHANG, Feng-yan SHI, Yong HU, Qing-yun YANG. Experimental and simulation studies on the anisotropic behavior of E40 steel[J]. Journal of Ship Mechanics, 2026 , 30 (5) : 804 -815 . DOI: 10.3969/j.issn.1007-7294.2026.05.012

正文

收起

0　引　言

收起

随着海洋工程与船舶制造对高性能结构材料需求的不断提升，高强度船用钢在大型船体结构中的应用日益广泛。其中，E40钢作为一种在低温环境下仍具备优良强度与韧性的高强钢，广泛应用于船体外板、甲板结构以及高应力关键承载区域，在承受复杂海洋载荷与波浪冲击方面发挥着重要作用^[1−2]。在此背景下，船体板材的成形精度直接关系到结构安全性与装配质量。冷弯成形作为一种兼具高精度与低能耗的加工技术^[3]，相较于传统火工成形无需经历加热与冷却过程，从而有效避免了材料性能劣化与组织变化问题，已成为船舶构件制造的关键工艺。然而，船用钢在成形过程中的回弹效应较为明显，会导致零件几何尺寸偏差及局部应力集中，从而影响成形精度和后续加工质量^[4−5]。因此，作为数值仿真理论的核心，本构模型对真实变形过程的精确描述不仅是保证材料成形有限元模拟准确性的前提，也直接决定了工艺优化的成效。本文通过解析E40钢各向异性特征及其与回弹行为的关联机制，为有限元仿真优化及成形工艺改进提供理论支撑。

目前，国内外学者围绕金属材料的各向异性已开展了广泛研究，多种描述材料各向异性的模型被提出，包括基于介观尺度的Bishop-Hill模型^[6]、粘塑性自洽模型VPSC（Visco-Plastic Self-Consistent Model）^[7]、晶体塑性有限元法CPFEM（Crystal Plasticity Finite Element Method）^[8]等，但其高昂的计算成本制约了工程应用。工程实践中，基于各向异性屈服模型与硬化定律耦合的唯象模型凭借计算效率优势成为主流方案。经典的Hill48屈服模型奠定了正交各向异性理论框架^[9]，Chen等^[10]基于该模型构建的板材本构模型经反向深冲实验验证,显示良好精度。但Mu等^[11]通过系统分析指出，Hill48模型的预测有效性取决于材料各向异性演化的单调性特征，且与

$ R $

值参数无显著相关性。为突破该限制，Barlat等^[12−15]通过对应力张量进行线性变换,相继提出了Yld89、Yld91、Yld2000-2d、Yld2004-18p等各向异性屈服模型，其多参数架构显著提升了复杂应力状态下的各向异性表征能力。张玮等^[16]对船用E36钢的力学特性研究表明，其非线性应变硬化效应会显著改变材料的塑性流动规律，这要求本构建模需实现各向异性与硬化行为的精确耦合。Chen等^[17]通过整合传统实验与有限元逆向识别（FEMU）方法，构建了AA5086铝合金的Yld2000-2d参数体系，并通过深冲成形实验验证了模型的可靠性。Wang等^[18]进一步构建了基于Yld2004-18p的三维非关联本构模型，结合圆杯拉伸与扩孔实验，验证了模型在三维应力下对板材各向异性变形行为的高保真预测能力。

为了提升船用E40钢板冷弯成形回弹预测精度，本研究基于单轴拉伸实验分别实现对Hill48、Yld2000-2d、Yld2004-18p三种屈服模型的参数标定，并结合修正型Hockett-Sherby硬化定律构建了E40钢各向异性本构模型。为验证模型的适用性，进一步开展了冲压成形实验，并通过用户子程序VUMAT将本构模型嵌入ABAQUS平台进行数值模拟，系统评估了三种各向异性塑性模型的预测精度，为船用E40钢精密成形工艺参数优化提供理论支撑。

1　E40钢力学性能测试

收起

1.1　单向拉伸实验

为了获取E40钢的力学性能参数，依据GB/T 228.1-2021《金属材料拉伸试验第1部分：室温试验方法》标准规范制备测试样本^[19]。其中，轧制方向（Rolling Direction, RD）定义为板材在轧制过程中的主变形方向，横向（Transverse Direction, TD）为与RD垂直的板面方向。本文中“角度”定义为试样拉伸方向与轧制方向（RD）之间的夹角。以RD为0°基准方向，沿RD至TD方向按15°间隔选取不同取向（0°~90°）制备7组拉伸试样。文中所有涉及的

$ R $

值均表示厚向异性系数，用于表征材料在塑性变形过程中宽度方向塑性应变与厚度方向塑性应变之比。不同下标（如

$ {R}_{\text{RD}} $

、

$ {R}_{\text{TD}} $

、

$ {R}_{\theta } $

、

$ {R}_{\mathrm{b}} $

等）仅表示取向或加载状态的差异，其物理含义保持一致。试样长度为180 mm，几何形状为扁平狗骨形。单向拉伸实验在YYF-50型拉伸试验机上执行，设置位移速率为2 mm/min，同步采用位移传感器采集载荷位移数据。图1展示了拉伸试样尺寸参数与实验装置。

1.2　实验结果分析

通过单向拉伸实验获取E40钢的力学性能参数，结果如图2和表1所示。实验结果表明，E40钢在不同取向上的力学性能存在一定差异，表现出明显的各向异性特征。该各向异性特征对板材成形极限与回弹控制具有双重影响，在工程实践中需要建立包含各向异性参数的成形工艺优化模型。

依据GB/T 5027-2007^[20]，当应变达10%时，采集横向应变与厚向应变并通过最小二乘法建立回归模型（R²

$ \geq $

0.98），其斜率绝对值为厚向异性系数

$ R $

值，如图3所示。

基于RD方向的应力–应变曲线来确定材料的硬化准则。E40钢的加工硬化过程通常包括一个初始快速硬化阶段和后续的线性硬化阶段。因此，为了更准确地描述材料的塑性流动特征，本研究对Hockett-Sherby函数进行了修正，通过引入线性项以同时表征初始快速加工硬化阶段和随后的线性硬化行为，如公式（1）所示

$ \sigma ={\sigma }_{0}+K\varepsilon +Q(1-{{\mathrm{e}}}^{-b{{\varepsilon }^{a}}}) $

(1)

公式（1）中的

$ {\sigma }_{0} $

表示初始屈服应力，

$ \varepsilon $

为等效塑性应变，

$ {\sigma }_{0} $

、

$ K $

、

$ Q $

、

$ b $

、

$ a $

表示材料参数，其具体数值通过单轴拉伸试验逆向识别获得（见表2）。由图4可知，实验数据与拟合曲线吻合较好，说明所采用的硬化准则能够较为准确地预测E40钢板材的塑性变形情况。

2　各向异性塑性模型标定

收起

2.1　Hill48屈服模型

在平面应力条件下，Hill48屈服模型表达式为

$ f=\left(G+H\right)\sigma _{xx}^{2}-2H{\sigma }_{xx}{\sigma }_{yy}+\left(F+H\right)\sigma _{yy}^{2}+2N\sigma _{xy}^{2}=\overline{{\sigma }^{2}} $

(2)

式中：

$ f $

为屈服函数，用于描述材料在各向异性应力状态下的屈服行为；

$ \overline{{\sigma }}$

表示等效应力，

$ \overline{{\sigma }^{2}} $

由应力分量通过Hill屈服准则计算得到。F、G、H和N是表征材料各向异性的参数。有两种求解方法，一种是采用屈服应力进行求解，表示为Hill48-S，如下式所示

$ \begin{array}{l}2F={\left(\dfrac{{\sigma }_{\rm{RD}}}{{\sigma }_{\rm{TD}}}\right)}^{2}+{\left(\dfrac{{\sigma }_{\rm{RD}}}{{\sigma }_{\rm{b}}}\right)}^{2}-1\\2G={\left(\dfrac{{\sigma }_{\rm{RD}}}{{\sigma }_{\rm{b}}}\right)}^{2}+1-{\left(\dfrac{{\sigma }_{\rm{RD}}}{{\sigma }_{\rm{TD}}}\right)}^{2}\\2H=1+{\left(\dfrac{{\sigma }_{\rm{RD}}}{{\sigma }_{\rm{TD}}}\right)}^{2}-{\left(\dfrac{{\sigma }_{\rm{RD}}}{{\sigma }_{\rm{b}}}\right)}^{2}\\2N={\left(\dfrac{2{\sigma }_{\rm{RD}}}{{\sigma }_{\rm{DD}}}\right)}^{2}-{\left(\dfrac{{\sigma }_{\rm{RD}}}{{\sigma }_{\rm{b}}}\right)}^{2}\\{\sigma }_{\rm{b}}=\dfrac{{\sigma }_{\rm{RD}}+2{\sigma }_{\rm{DD}}+{\sigma }_{\rm{TD}}}{4}\end{array} $

(3)

式中：

$ {\sigma }_{\mathrm{b}} $

表示等双拉应力状态下的屈服应力，可由0°、45°和90°方向的单轴屈服应力计算得到^[21]。

第二种方法为采用

$ R $

值进行求解，表示为Hill48-

$ R $

，如下式所示

$ \begin{array}{l}F=\dfrac{{R}_{\rm{RD}}}{\left(1+{R}_{\rm{RD}}\right){R}_{\rm{TD}}}\\G=\dfrac{1}{1+{R}_{\rm{RD}}}\\H=\dfrac{{R}_{\rm{RD}}}{1+{R}_{\rm{RD}}}\\N=\dfrac{\left(1+2{R}_{\rm{DD}}\right)\left({R}_{\rm{RD}}+{R}_{\rm{TD}}\right)}{2\left(1+{R}_{\rm{RD}}\right){R}_{\rm{DD}}}\end{array} $

(4)

式中：

$ {R}_{\text{RD}} $

、

$ {R}_{\text{DD}} $

和

$ {R}_{\text{TD}} $

分别表示沿轧制方向（0°）、45°方向及横向（90°）的厚向异性系数，其数值由表1中的单轴拉伸试验结果获得。

根据表1中的各向异性参数求得Hill48屈服模型的参数，结果见表3。

2.2　Yld2000-2d屈服模型

Yld2000-2d是由Barlat等^[14]所提出的一种非二次的各向异性屈服模型，对应的等效应力形式如下

$ {\sigma }_{{\mathrm{eq}}}={\left\{\dfrac{{\left| X_{1}^{'}-X_{2}^{'}\right| }^{{\rm{M}}}+{\left| 2X_{2}^{''}+X_{1}^{''}\right| }^{{\rm{M}}}+{\left| 2X_{1}^{''}+X_{2}^{''}\right| }^{{\rm{M}}}}{2}\right\}}^{1/{\rm{M}}} $

(5)

式中：M为与晶体结构有关的常数，对于E40钢，M=6，

$ X_{i}^{'} $

和

$ X_{i}^{''} $

（

$ i $

=1,2）为二阶张量

$ {{\boldsymbol{X}}}^{'} $

和

$ {{\boldsymbol{X}}}^{''} $

的主值，其中

$ i=1 $

对应“+”号，

$ i=2 $

对应“−”号。形式如下

$ \begin{cases} X_{i}^{'}=\dfrac{1}{2}\left(X_{11}^{'}+X_{22}^{'}\pm \sqrt{\left(X_{11}^{'}-X_{22}^{'}\right)+4X_{12}^{'2}}\right)\\X_{i}^{''}=\dfrac{1}{2}\left(X_{11}^{''}+X_{22}^{''}\pm \sqrt{\left(X_{11}^{''}-X_{22}^{''}\right)+4X_{12}^{''2}}\right)\end{cases} $

(6)

$ {\boldsymbol{X}}^{{'}}={\boldsymbol{L}}^{{'}}\boldsymbol{\sigma };\;\;{\boldsymbol{X}}^{{''}}={\boldsymbol{L}}^{{''}}\boldsymbol{\sigma } $

(7)

式中：

$ {\boldsymbol{\sigma}} $

为柯西应力张量,

$ {X}^{'} $

和

$ {X}^{''} $

为线性转换向量，其形式如下

$ \left[\begin{array}{c}L_{11}^{'}\\L_{12}^{'}\\L_{21}^{'}\\L_{22}^{'}\\L_{66}^{'}\end{array}\right]=\dfrac{1}{3}\left[\begin{array}{ccc}2 & 0 & 0\\-1 & 0 & 0\\0 & -1 & 0\\0 & 2 & 0\\0 & 0 & 3\end{array}\right]\left[\begin{array}{c}{\alpha }_{1}\\{\alpha }_{2}\\{\alpha }_{2}\end{array}\right] , \left[\begin{array}{c}L_{11}^{''}\\L_{12}^{''}\\L_{21}^{''}\\L_{22}^{''}\\L_{66}^{''}\end{array}\right]=\dfrac{1}{9}\left[\begin{array}{ccccc}-2 & 2 & 8 & -2 & 0\\-1 & -4 & -4 & 4 & 0\\0 & -4 & -4 & 1 & 0\\-2 & 8 & 2 & -2 & 0\\0 & 0 & 0 & 0 & 9\end{array}\right]\left[\begin{array}{c}{\alpha }_{3}\\{\alpha }_{4}\\{\alpha }_{5}\\{\alpha }_{6}\\{\alpha }_{8}\end{array}\right] $

(8)

式中：

$ \alpha 1\sim \alpha 8 $

为材料的八个各向异性参数。

相较于仅含四个各向异性参数的Hill48屈服模型，Yld2000-2d包含八个各向异性参数，其扩展的参数空间显著增强了非对称屈服行为的表征能力。基于表1中不同取向单轴拉伸试验获得的屈服应力及厚向异性系数，采用粒子群优化算法进行参数优化，以最小化方程（9）中定义的误差函数，得出Yld2000-2d的各个参数见表4。

$ E({\alpha }_{1},\cdots,{\alpha }_{8})=\sum\limits_{i}\,{\left(\dfrac{{\sigma }_{i}}{\sigma _{i}^{\exp }}-1\right)}^{2}+\sum\limits_{i}\,{\left(\dfrac{{R}_{i}}{R_{i}^{\exp }}-1\right)}^{2} $

(9)

式中：

$ \sigma _{i}^{\exp } $

和

$ R_{i}^{\exp } $

分别代表第

$ i $

个方向下测得的屈服应力和厚向异性系数的实验数据，

$ {\sigma }_{i} $

和

$ {R}_{i} $

则表示Yld2000-2d屈服模型在第

$ i $

个方向上的预测值。

2.3　Yld2004-18p屈服模型

针对Yld2000-2d屈服模型在复杂应力状态下的表征局限，Barlat等^[15]基于高阶张量变换理论提出了具有18个各向异性参数的Yld2004-18p屈服模型。该模型通过引入双线性变换张量对应力偏张量

$ {\boldsymbol{S }}$

进行线性变换。该屈服模型如下

$ \varphi (\sigma )=\sum\limits_{i=1}^{3}\,\sum\limits_{j=1}^{3}\,|\tilde{S}_{i}^{'}-\tilde{S}_{j}^{''}{|}^{a}=4{\overline{\sigma }}^{a} $

(10)

对于E40钢，

$ a $

=6。

$ \sigma $

为流动应力（等效应力），

$ \tilde{S}_{1}^{'} $

和

$ \tilde{S}_{1}^{''} $

为张量

$ {\tilde{{\boldsymbol{S}}}}^{\mathbf{'}}、{\tilde{{\boldsymbol{S}}}}^{\mathbf{''}} $

的主应力值，它们由应力偏张量

$ {\boldsymbol{S }} $

的两次线性变换定义。如下式

$ {\tilde{\boldsymbol{S}}}^{\boldsymbol{'}}={\boldsymbol{C}}^{\boldsymbol{'}}\boldsymbol{S}={\boldsymbol{C}}^{\boldsymbol{'}}\boldsymbol{T}{\sigma }={\boldsymbol{L}}^{\boldsymbol{'}}{\sigma } $

(11)

$ {\tilde{\boldsymbol{S}}}^{{'}{'}}={\boldsymbol{C}}^{{'}{'}}\boldsymbol{S}={\boldsymbol{C}}^{{'}{'}}\boldsymbol{T}{\sigma }={\boldsymbol{L}}^{{'}{'}}{\sigma } $

(12)

公式（11）、（12）中，

$ {{\boldsymbol{C}}}^{'} $

和

$ {{\boldsymbol{C}}}^{''} $

是应力张量线性变换矩阵，如下式

$ {\boldsymbol{C}}^{\boldsymbol{'}}=\left[\begin{array}{cccccc}0 & -c_{12}^{'} & -c_{13}^{'} & 0 & 0 & 0\\-c_{21}^{'} & 0 & -c_{23}^{'} & 0 & 0 & 0\\-c_{31}^{'} & -c_{32}^{'} & 0 & 0 & 0 & 0\\0 & 0 & 0 & c_{44}^{'} & 0 & 0\\0 & 0 & 0 & 0 & c_{55}^{'} & 0\\0 & 0 & 0 & 0 & 0 & c_{66}^{'}\end{array}\right] $

(13)

$ {\boldsymbol{C}}^{{''}}=\left[\begin{array}{cccccc}0 & -c_{12}^{''} & -c_{13}^{''} & 0 & 0 & 0\\-c_{21}^{''} & 0 & -c_{23}^{''} & 0 & 0 & 0\\-c_{31}^{''} & -c_{32}^{''} & 0 & 0 & 0 & 0\\0 & 0 & 0 & c_{44}^{''} & 0 & 0\\0 & 0 & 0 & 0 & c_{55}^{''} & 0\\0 & 0 & 0 & 0 & 0 & c_{66}^{''}\end{array}\right] $

(14)

$ {\boldsymbol{T}} $

也是线性变换，如下式所示

$ \boldsymbol{T}=\dfrac{1}{3}\left[\begin{array}{cccccc}2 & -1 & -1 & 0 & 0 & 0\\-1 & 2 & -1 & 0 & 0 & 0\\-1 & -1 & 2 & 0 & 0 & 0\\0 & 0 & 0 & 3 & 0 & 0\\0 & 0 & 0 & 0 & 3 & 0\\0 & 0 & 0 & 0 & 0 & 3\end{array}\right] $

(15)

为了确定Yld2004-18p的各向异性系数，基于表1中不同取向单轴拉伸试验获得的力学性能数据，采用梯度下降算法对屈服模型参数进行反演标定，以最小化误差函数。最终得到的模型参数如表5所示。

在本工作中，误差函数定义为

$ E\left(c_{ij,}^{'}c_{ij}^{''}\right)=\sum\limits_{p}\,{w}_{p}{\left(\dfrac{\sigma _{p}^{{\mathrm{pr}}}}{\sigma _{p}^{{\mathrm{ex}}}}-1\right)}^{2}+\sum\limits_{q}\,{w}_{q}{\left(\dfrac{r_{q}^{{\mathrm{pr}}}}{r_{q}^{{\mathrm{ex}}}}-1\right)}^{2} $

(16)

式中：

$ p $

为不同方向单轴拉伸或单轴剪切时可用的实验流动应力个数，

$ q $

为实验

$ R $

值的个数。上式中，上标ex或pr表示对应值是实验值或预测值。误差函数中的每一项都用

$ w $

加权。该权重可以用来区分流动应力或

$ r $

值。将平面内流动应力的权重、单轴拉伸时的

$ R $

值，以及其他流动应力的权重分别设为1.00、0.10和0.01。Yld2004-18p屈服模型中包含14个代表板料面内性能的参数。由于板料在厚度方向上不存在弯曲，

$ c_{55}^{'} $

、

$ c_{55}^{''} $

和

$ c_{66}^{'} $

、

$ c_{66}^{''} $

表征板料面外的剪切属性设置为1。基于表1的试验结果，得出

$ c_{ij}^{'} $

和

$ c_{ij}^{''} $

的具体参数见表5。

2.4　塑性各向异性参数

根据Hill正交各向异性塑性理论^[9]，用于表征板材的面内各向异性特征。该参数与表1中所列

$ R $

值具有相同物理含义，均表示宽度方向塑性应变与厚度方向塑性应变之比。其表达式为

$ R=\dfrac{{\mathrm{d}}{\varepsilon }_{b}}{{\mathrm{d}}{\varepsilon }_{t}} $

(17)

对于沿x轴（RD）的单轴拉力，

$ R $

值可进一步表示为

$ {R}_{\theta }=-\dfrac{\dfrac{\partial \varphi }{\partial {\sigma }_{xx}}{\sin}^{2} \theta +\dfrac{\partial \varphi }{\partial {\sigma }_{yy}}{\cos}^{2} \theta -\dfrac{\partial \varphi }{\partial {\sigma }_{xy}}\sin 2\theta }{\left(\dfrac{\partial \varphi }{\partial {\sigma }_{xx}}+\dfrac{\partial \varphi }{\partial {\sigma }_{yy}}\right)} $

(18)

式中：

$ \varphi $

表示屈服函数，其具体形式由所选各向异性屈服模型（Hill48、Yld2000-2d或Yld2004-18p）确定。需说明的是，前文中的

$ f $

与此处的

$ \varphi $

均表示屈服函数，仅符号表示不同，其物理含义一致。式中各偏导数项基于关联流动法则由屈服函数对应力分量求导得到，用于描述塑性应变增量方向。

$ \theta $

为单轴拉伸方向相对于轧制方向（RD）的取向角。

关于正交各向异性坐标的单轴应力张量与柱坐标的关系为

$ {\boldsymbol{\sigma }}_{\boldsymbol{\theta }}({\boldsymbol{Y}}_{\boldsymbol{\theta }},\boldsymbol{\theta })=\left[\begin{matrix}{\sigma }_{11} & {\sigma }_{12} & 0\\{\sigma }_{21} & {\sigma }_{22} & 0\\0 & 0 & 0\end{matrix}\right]=\left[\begin{array}{ccc}{Y}_{\theta }{\cos}^{2} \theta & {Y}_{\theta }\sin \theta \cos \theta & 0\\{Y}_{\theta }\sin \theta \cos \theta & {Y}_{\theta }{\sin}^{2} \theta & 0\\0 & 0 & 0\end{array}\right] $

(19)

式中：

$ {Y}_{\theta } $

表示沿取向角

$ \theta $

方向的单轴拉伸屈服应力。

3　不同屈服模型对E40钢各向异性的预测能力

收起

分别采用Hill48、Yld2000-2d和Yld2004-18p屈服模型计算得到不同加载方向时归一化屈服应力的理论值，并将其与实验结果进行比较，结果如表6和表7所示。

根据公式（20）计算不同屈服模型预测屈服应力和

$ R $

值的误差

$ E $

，其结果如图5所示。数据分析表明，Yld2004-18p屈服模型预测精度最高，屈服应力与

$ R $

值分别实现0.12%和0.55%的最低平均偏差；Hill48屈服模型预测精度最低，对应误差分别为2.24%和8.5%；Yld2000-2d屈服模型预测精度中等，对应偏差分别为0.58%和3.75%。

$ {E}_{1}=\left| \dfrac{{\sigma }_{{\mathrm{pred}}}-{\sigma }_{{\mathrm{exp}}}}{{\sigma }_{{\mathrm{exp}}}}\right| \times 100\% $

(20)

式中：

$ {\sigma }_{\exp } $

为实验值，

$ {\sigma }_{\text{pred}} $

为理论预测值。

图6给出了不同屈服模型预测结果与实验数据的对比情况，其中图6（a）为屈服应力随取向角变化的分布，图6（b）为厚向异性系数

$ R $

随取向角的变化规律。从图6中可以看出，Yld2004-18p屈服模型在非45°轧制取向区间展现出较好的吻合度，说明其复杂参数体系能较好捕捉多角度各向异性特征。实验数据点与Yld2000-2d屈服模型的理论预测曲线有一定误差，与Hill48的理论预测曲线误差最大。该现象源于Hill48模型受限于其各向异性参数采集维度不足，对非对称各向异性行为的适应性较差。

图7为不同屈服模型下的屈服轨迹与实验值的对比，Hill48屈服模型因参数简化及数学形式限制，屈服轨迹与实验值偏差显著；而Yld2000-2d屈服模型轨迹形状较接近实验曲线，尤其在剪切主导区，吻合度提升明显，但在纯剪切区偏差较大，可能是由于二维平面假设导致对纯剪切和双轴拉伸的耦合效应捕捉不足。综上所述，Yld2004-18p屈服模型对E40钢的各向异性行为具有更优越的表征能力。

4　模型验证分析

收起

为进一步验证模型预测效果，加工了两块尺寸为100 mm×30 mm×10 mm的长方体坯料，分别采用半径为200 mm和500 mm的模具进行冲压实验。

基于ABAQUS平台构建了E40钢冷弯成形显–隐式顺序耦合仿真模型。如图8所示，模具简化为刚体，板材采用C3D8R单元，弹性模量E = 210 GPa，泊松比

$ \nu $

= 0.3，摩擦系数设为0.1。成形阶段使用ABAQUS/Explicit模块进行显式动态分析：设置位移载荷驱动上模下压，采用自主开发的VUMAT子程序分别将Hill48、Yld2000-2d和Yld2004-18p各向异性屈服模型及硬化准则嵌入材料属性中，模拟r = 200 mm/500 mm两种弯曲半径下的非线性变形；回弹阶段转入ABAQUS/Standard模块执行隐式静态分析，将成形终态应力场作为初始条件，通过完全释放边界约束，精确模拟卸载后的弹性回复行为。该方法通过显式算法保障大变形计算的收敛性，采用隐式求解确保回弹量的计算精度。数值分析的基本条件如表8所示。

表9比较了Hill48、Yld2000-2d和Yld2004-18p三种屈服模型在冲压模拟中的计算时间。可以看出，Hill48因形式简单、参数少而计算耗时最短，Yld2000-2d引入更多面内各向异性参数计算耗时显著增加，Yld2004-18p由于阶数高且需求解主应力，计算时间最长。此外，小弯曲半径（

$ r $

= 200 mm）下局部曲率大，显式积分求解时间步减小，导致计算时间反而略长于大半径工况。

图9展示了在不同弯曲半径（

$ r $

= 200 mm与

$ r $

= 500 mm）条件下，沿板材上表面长度方向（X方向）的回弹变形法向位移误差分布情况，其中误差定义为模型预测结果与实验测量结果在板材表面法向方向上的位移差。所有模型的回弹误差在X = 0附近均表现出对称分布。在变形程度较大时（

$ r $

= 200 mm），Hill48与Yld2000-2d在边缘区域误差较大，峰值超过1 mm，预测精度不理想；而Yld2004-18p模型在整个区域内误差均控制在5%以内，表现出更高的预测一致性和空间均匀性。随着变形程度减小（

$ r $

= 500 mm），所有模型的误差明显减小，其中Yld2004-18p误差小于3%，为最优表现。Hill48模型精度中等，边缘区域误差略有波动，而Yld2000-2d虽在整体应变较缓区域有所改善，但边缘误差仍较明显。

上述预测差异主要源于各屈服模型对材料各向异性和非线性塑性行为的描述能力差异。Hill48模型由于采用二次屈服函数，难以准确捕捉材料在复杂应力状态下的真实变形特征，尤其在小变形条件下，其相对误差可达15%。Yld2000-2d虽然通过引入参数扩展改善了各向异性描述能力，但在应力梯度较高区域依然存在预测偏差。相比之下，Yld2004-18p模型通过更高阶的屈服函数形式和优化参数设置，能更准确地描述E40钢板材的塑性应变演化和厚向回弹特性，在不同工况下均表现出优异的成形精度与工程适用性。

5　结　论

收起

本文围绕E40钢板材在拉伸状态下的各向异性行为及其冲压回弹预测，结合单向拉伸实验、屈服模型参数标定及数值模拟，对比分析了Hill48、Yld2000-2d和Yld2004-18p三种各向异性屈服模型的适用性与预测精度，主要结论如下：

（1） E40钢在拉伸状态下表现出显著的各向异性特征，实验结果显示，不同方向屈服应力最大相差13  MPa，伸长率相差2.3%，厚向异性系数均小于1。

（2）基于实验数据对Hill48、Yld2000-2d和Yld2004-18p三种屈服模型进行了参数标定，并采用修正型Hockett-Sherby模型描述材料的应变硬化行为。结果表明，Yld2004-18p模型在主应力空间下的屈服面轨迹与实验吻合度最高。

（3）冲压回弹实验结果验证了不同屈服模型的成形预测精度。Hill48和Yld2000-2d在厚度减薄率及回弹模拟中的最大相对误差分别约为10%和15%，Yld2004-18p在小半径（

$ r $

= 200 mm）和大半径（

$ r $

= 500 mm）工况下预测误差均小于4%，较Hill48模型精度提高约3倍，能更准确地反映E40钢的塑性变形与回弹特性，具有良好的工程适用性。

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

Chen

, Guo

, Zhang

, et al. Effect of water flow on the microstructure, mechanical performance, and cracking susceptibility of underwater wet welded Q235 and E40 steel[J]. Journal of Materials Processing Technology, 2020, 277: 116435.

Wang

, Ma

, Liu

, et al. Influence of heat input on microstructure and corrosion resistance of underwater wet-welded E40 steel joints[J]. Journal of Materials Engineering and Performance, 2020, 29(11): 1‒10.

Mutafi

, Irwan

J M

, Yidris

, et al. An in-depth comparative FEA on through-thickness residual stresses in press-braked cold-formed steel section[J]. Results in Engineering, 2024, 22: 102124.

Wang

, Sun

, Liu

, et al. Research on springback mechanism of TA2/Q235B composite plate in roll forming[J]. Journal of Manufacturing Processes, 2024, 119: 715‒728.

梁棋钰, 张龙, 朱凌. 考虑部分卸载效应的多压头凹形成形及回弹预测(英文)[J]. 船舶力学, 2024, 28(12): 1953‒1969.

Liang

Q Y

, Zhang

, Zhu

. Forming and springback prediction of strips under multi-square punch concave forming process considering partial-unloading effects[J]. Journal of Ship Mechanics, 2024, 28(12): 1953‒1969. (in Chinese)

Juul

N Y

, Oddershede

, Beaudoin

, et al. Measured resolved shear stresses and Bishop-Hill stress states in individual grains of austenitic stainless steel[J]. Acta Materialia, 2017, 141: 388‒404.

Zhou

, Wang

, et al. Dislocation behavior in a polycrystalline Mg-Y alloy using multi-scale characterization and VPSC simulation[J]. Journal of Materials Science & Technology, 2022, 98: 87‒98.

Zhou

, Ben

Bettaieb M

, Abed-meraim

. Investigation of the effect of morphological and crystallographic textures on the ductility limits of thin metal sheets using a CPFEM-based approach[J]. European Journal of Mechanics-A/Solids, 2024, 106: 105293.

Hill

. A theory of the yielding and plastic flow of anisotropic metals[J]. Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences, 1948, 193(1033): 281‒297.

Chen

, Zhao

, Shi

, et al. Nonlinear kinematic hardening constitutive model based on Hill48 yield criterion and its application in reverse deep drawing[J]. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2022, 44(10).

, Zhao

, Meng

, et al. Applicability of Hill48 yield model and effect of Anisotropic Parameter Determination Methods on anisotropic prediction[J]. Journal of Materials Engineering and Performance, 2022, 31(3): 2023‒2042.

Barlat

, Lian

. Plastic behavior and stretchability of sheet metals. Part I: A yield function for orthotropic sheets under plane stress conditions[J]. International Journal of Plasticity, 1989, 5(1): 51‒66.

Barlat

, Lege

D J

, Brem

J C

. A six-component yield function for anisotropic materials[J]. International Journal of Plasticity, 1991, 7(7): 693‒712.

Barlat

, Brem

J C

, Yoon

J W

, et al. Plane stress yield function for aluminum alloy sheets—part 1: Theory[J]. International Journal of Plasticity, 2003, 19(9): 1297‒1319.

Barlat

, Aretz

, Yoon

J W

, et al. Linear transfomation-based anisotropic yield functions[J]. International Journal of Plasticity, 2005, 21(5): 1009‒1039.

张玮, 毛柳伟, 杨坤, 等. 船用E36钢的动态力学行为及本构关系(英文)[J]. 船舶力学, 2022, 26(6): 912‒918.

Zhang

, Mao

L W

, Yang

, et al. Dynamic mechanical behavior and constitutive relation of shipbuilding E36 steel[J]. Journal of Ship Mechanics, 2022, 26(6): 912‒918. (in Chinese)

Chen

, Wang

, Chu

, et al. Calibration of Yld2000-2D anisotropy yield criterion with traditional testing and inverse identification strategies[J]. Materials, 2023, 16(21): 6904.

Wang

, Fan

, Yoshimura

, et al. A non-associated constitutive model based on Yld2004-18p yield criterion and its applications on sheet metal forming analysis[J]. Scientific Reports, 2025, 15(1).

GB/T 228.1-2021. 金属材料拉伸试验第1部分: 室温试验方法[S].

GB/T 228.1-2021. Metallic materials-tensile testing-Part I: Method of test at room temperature[S]. (in Chinese)

GB/T 5027-2007. 金属薄板和薄带塑性应变比(r值)试验方法[S].

GB/T 5027-2007. Metallic materials sheet and strip determination of plastic strain ratio(r-values) [S]. (in Chinese)

Fuhui

, Sebastian

, Junhe

. Forming limit prediction by the Marciniak-Kuczynski model coupled with the evolving non-associated Hill48 plasticity model[J]. Journal of Materials Processing Tech., 2019, 287.

2026年第30卷第5期

PDF下载

163

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.1007-7294.2026.05.012

接收时间：2025-07-17
首发时间：2026-07-07
出版时间：2026-05-15

补充材料

相关文章

文章信息

作者

出版历史

收稿日期：2025-07-17

基金

作者信息

^1.湖北工业大学机械工程学院，武汉 430068

^2.现代制造质量工程湖北省重点实验室，武汉 430068

^3.武汉理工大学船海与能源动力工程学院，武汉 430063

^4.中国船舶集团有限公司第七二五研究所，河南洛阳 471023

通讯作者:

张文谦（1990–），男，博士，副教授，通讯作者，E-mail： wenqian_zh@hbut.edu.cn

参考文献

分享链接

https://castjournals.cast.org.cn/joweb/cblx/CN/10.3969/j.issn.1007-7294.2026.05.012

分享至

全文二维码

扫描看全文

引用本文

BibTeX

本文的引用情况

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT

角度/°	屈服强度$ {\sigma }_{\mathrm{s}}$ /MPa	抗拉强度$ {\sigma }_{\mathrm{b}} $ /MPa	伸长率A /%	厚向异性系数R
0	348.9	636.0	31.5	0.887
15	348.4	628.8	30.6	0.868
30	354.4	641.3	31.4	0.897
45	361.3	634.4	31.3	0.945
60	357.1	652.5	31.8	0.847
75	357.0	639.0	30.8	0.808
90	359.7	639.3	32.9	0.859

角度/°

屈服强度

$ {\sigma }_{\mathrm{s}}$

/MPa

抗拉强度

$ {\sigma }_{\mathrm{b}} $

/MPa

伸长率A /%

厚向异性系数R

348.9

636.0

31.5

0.887

348.4

628.8

30.6

0.868

354.4

641.3

31.4

0.897

361.3

634.4

31.3

0.945

357.1

652.5

31.8

0.847

357.0

639.0

30.8

0.808

359.7

639.3

32.9

0.859

参数	$ {\sigma }_{0} $	K	Q	b	a
Hockett-Sherby+Linear	983.835	131.570	7.387	0.663	−150.245

参数

$ {\sigma }_{0} $

Hockett-Sherby+Linear

983.835

131.570

7.387

0.663

−150.245

	F	G	H	L	M	N
Hill48-S	0.4889	0.5479	0.4521	1.5	1.5	1.3472
Hill48-R	0.5472	0.5299	0.4701	1.5	1.5	1.5565

Hill48-S

0.4889

0.5479

0.4521

1.5

1.3472

Hill48-R

0.5472

0.5299

0.4701

1.5

1.5565

$ \alpha 1 $	$ \alpha 2 $	$ \alpha 3 $	$ \alpha 4 $	$ \alpha 5 $	$ \alpha 6 $	$ \alpha 7 $	$ \alpha 8 $
1.069	0.918	1.186	1.025	1.073	1.355	0.981	0.773

$ \alpha 1 $

$ \alpha 2 $

$ \alpha 3 $

$ \alpha 4 $

$ \alpha 5 $

$ \alpha 6 $

$ \alpha 7 $

$ \alpha 8 $

1.069

0.918

1.186

1.025

1.073

1.355

0.981

0.773

参数	$ c_{12}^{'} $	$ c_{13}^{'} $	$ c_{21}^{'} $	$ c_{23}^{'} $	$ c_{31}^{'} $	$ c_{32}^{'} $	$ c_{44}^{'} $	$ c_{55}^{'} $	$ c_{66}^{'} $
实验值	0.8689	0.7209	0.5009	0.9520	0.6541	1.2411	0.9357	1	1

参数	$ c_{12}^{''} $	$ c_{13}^{''} $	$ c_{21}^{''} $	$ c_{23}^{''} $	$ c_{31}^{''} $	$ c_{32}^{''} $	$ c_{44}^{''} $	$ c_{55}^{''} $	$ c_{66}^{''} $
实验值	0.6121	1.1445	0.9676	0.7189	1.2240	0.4174	0.9854	1	1

参数

$ c_{12}^{'} $

$ c_{13}^{'} $

$ c_{21}^{'} $

$ c_{23}^{'} $

$ c_{31}^{'} $

$ c_{32}^{'} $

$ c_{44}^{'} $

$ c_{55}^{'} $

$ c_{66}^{'} $

实验值

0.8689

0.7209

0.5009

0.9520

0.6541

1.2411

0.9357

参数

$ c_{12}^{''} $

$ c_{13}^{''} $

$ c_{21}^{''} $

$ c_{23}^{''} $

$ c_{31}^{''} $

$ c_{32}^{''} $

$ c_{44}^{''} $

$ c_{55}^{''} $

$ c_{66}^{''} $

实验值

0.6121

1.1445

0.9676

0.7189

1.2240

0.4174

0.9854

方向/$ {^{\circ}} $	实验值	Hill48-S	Hill48-R	Yld2000-2d	Yld2004-18p
0	1	0.9903	1.0156	1.0004	0.9990
15	0.9987	1.0030	1.0101	1.0068	0.9997
30	1.0158	1.0184	0.9990	1.0217	1.0170
45	1.0354	1.0316	0.9923	1.0351	1.0331
60	1.0236	1.0345	0.9949	1.0388	1.0248
75	1.0233	1.0301	1.0027	1.0345	1.0244
90	1.0309	1.0271	1.0070	1.0314	1.0298

方向/

$ {^{\circ}} $

实验值

Hill48-S

Hill48-R

Yld2000-2d

Yld2004-18p

0.9903

1.0156

1.0004

0.9990

0.9987

1.0030

1.0101

1.0068

0.9997

1.0158

1.0184

0.9990

1.0217

1.0170

1.0354

1.0316

0.9923

1.0351

1.0331

1.0236

1.0345

0.9949

1.0388

1.0248

1.0233

1.0301

1.0027

1.0345

1.0244

1.0309

1.0271

1.0070

1.0314

1.0298

方向/$ {^{\circ}} $	实验值	Hill48-S	Hill48-R	Yld2000-2d	Yld2004-18p
0	0.887	0.825	0.887	0.888	0.887
15	0.868	0.814	0.903	0.902	0.862
30	0.897	0.795	0.934	0.931	0.893
45	0.945	0.799	0.945	0.948	0.957
60	0.847	0.842	0.920	0.923	0.852
75	0.808	0.898	0.879	0.880	0.814
90	0.859	0.925	0.859	0.861	0.858

方向/

$ {^{\circ}} $

实验值

Hill48-S

Hill48-R

Yld2000-2d

Yld2004-18p

0.887

0.825

0.887

0.888

0.887

0.868

0.814

0.903

0.902

0.862

0.897

0.795

0.934

0.931

0.893

0.945

0.799

0.945

0.948

0.957

0.847

0.842

0.920

0.923

0.852

0.808

0.898

0.879

0.880

0.814

0.859

0.925

0.859

0.861

0.858

节点数	单元数	内存	软件	操作系统	处理器
18227	15360	8 GB	ABAQUS/Explicit	64-bit Windows 10	Intel（R） Core（TM） i5-12600KF 3.70 GHz

节点数

单元数

内存

软件

操作系统

处理器

18227

15360

8 GB

ABAQUS/Explicit

64-bit Windows 10

Intel（R） Core（TM） i5-12600KF 3.70 GHz

弯曲半径	屈服模型
$ r $ = 500 mm	4 min 20 s	7 min 28 s	10 min 17 s
$ r $ = 200 mm	4 min 47 s	8 min 03 s	11 min 16 s

弯曲半径

屈服模型

Hill48

Yld2000-2d

Yld2004-18p

$ r $

= 500 mm

4 min 20 s

7 min 28 s

10 min 17 s

$ r $

= 200 mm

4 min 47 s

8 min 03 s

11 min 16 s