利索能及
我要发布
收藏
专利号: 2022108151729
申请人: 燕山大学
专利类型:发明专利
专利状态:已下证
更新日期:2026-06-16
缴费截止日期: 暂无
联系人

摘要:

权利要求书:

1.一种基于模式搜索法的双相钢Ohno‑Wang Ⅰ本构模型参数优化方法,其特征在于,其包括以下步骤:步骤1:建立双相钢材料的Ohno‑Wang Ⅰ本构模型;

所述双相钢材料的Ohno‑Wang Ⅰ本构模型包括材料流动法则、材料的屈服函数和材料的背应力方程组;

步骤2:根据双相钢材料的应力应变模型建立单轴拉伸曲线关系,获得本构模型参数初始值;

根据单轴拉伸曲线确定的求解双相钢材料本构模型参数初始值的方程组,如下所示:式中:γi表示第i个背应力分量对应的本构模型第一参数;ri表示第i个背应力分量对应的本构模型第二参数;σi和 分别为所选取的第i个单轴拉伸曲线上的点所对应的应力和塑性应变值;σ0表示单轴拉伸曲线上塑性应变值为0时的应力值;σmax表示单轴拉伸曲线应力峰值;M表示背应力分量个数;i表示背应力分量的编号;r1表示第一个背应力分量的模型第二参数;

步骤3:使用基于Isight软件搭建的双相钢材料本构模型优化平台,利用模式搜索法获取准确的本构模型参数;

步骤31:获取步骤2得到的本构模型参数初始值,搭建基于Isight软件的双相钢材料本构模型优化平台,求解本构模型参数的最优解;优化过程需要比较根据参数初始值产生的时间‑应力曲线和根据优化值产生的时间‑应力曲线,基于最小二乘法的目标函数计算公式如下:式中:L表示目标曲线和待优化曲线之间的应力差值平方之和;tj表示第j个时间‑应力曲线点中的时间,即曲线中横坐标; 表示第j个目标曲线中点的应力值,即曲线纵坐标;

表示第j个待优化曲线点的应力值;N表示时间‑应力曲线中的点个数;j表示目标曲线的段数编号;

步骤32:分别将本构模型参数的初始值赋值给x(1)作为优化过程的初始值;探测步长d,加速因子μ,μ>0,缩减因子β,n个探测方向e1,e2,e3,...,en,计算精度ε,令k=1;确定迭代参考值y:y=xk;

式中:y表示参数的迭代参考值;xk表示第k次迭代的过程参数;

步骤33:从本构模型参数的参考值开始计算,依次沿着方向e1,...,en,进行计算,并沿着正负两个方向进行;

沿正方向探测:若目标函数值满足L(y+de)<L(y),则迭代参考值更新为:y=y+de;

式中:d表示参数的探测步长;e表示参数的探测方向;

若L(y+de)≥L(y),进行负方向探测;

沿负方向探测:若目标函数值L(y‑de)<L(y),则迭代参考值更新如下所示:y=y‑de;

若L(y‑de)≥L(y),令y=y;

通过以上的探测,获得更新参考值y,第k+1次迭代的过程参数xk+1=y;

步骤34:沿着本构模型参数最优解的方向进行加速计算,过程如下:探测加速方向的计算公式如下所示:

pk=xk+1‑xk;

式中:pk表示计算过程的探测加速方向;

若目标函数值L(xk+1)<L(xk),迭代次数更新如下所示:式中:μ表示计算过程的探测加速因子;k表示计算过程的迭代次数;

否则,缩短步长d,即使d=dβ,判断是否满足迭代终止条件d≤ε,满足则终止迭代,不满足则转到步骤33继续进行探测;

步骤4:确定所有优化后的本构模型参数,建立准确的Ohno‑WangⅠ本构模型;

根据步骤3的优化过程,对M个本构模型第一参数和本构模型第二参数实现优化,得到全部的本构模型参数最优值,建立准确的双相钢材料Ohno‑WangⅠ本构模型实现对双相钢材料性质的准确分析。

2.根据权利要求1所述的基于模式搜索法的双相钢Ohno‑WangⅠ本构模型参数优化方法,其特征在于,所述步骤1中的材料流动法则,具体为:所述材料流动服从流动法则,如下所示:

式中:dεp表示材料的塑性应变增量;dλ表示材料的塑性乘子增量;f表示材料的屈服面函数;σ表示材料的应力值; 表示材料塑性应变演化的方向。

3.根据权利要求1所述的基于模式搜索法的双相钢Ohno‑WangⅠ本构模型参数优化方法,其特征在于,所述步骤1中材料的屈服函数,如下所示:所述屈服函数服从von Mises屈服准则,则屈服函数的获取方法如下所示:式中: 表示屈服函数的偏应力张量;α表示材料的背应力;I表示单位偏应力张量;tr(σ)表示偏应力张量的迹;σy表示材料的屈服应力,即屈服面的尺寸;

由vonMises屈服准则,塑性乘子增量与累积塑性应变增量相等,如下所示:式中:dp表示材料累积塑性应变增量。

4.根据权利要求1所述的基于模式搜索法的双相钢Ohno‑WangⅠ本构模型参数优化方法,其特征在于,所述步骤1中材料的背应力方程组,具体为:所述背应力方程组的公式如下所示:

式中:αi表示第i个背应力分量;dα表示材料的背应力增量;||αi||表示材料背应力分量的模,表示背应力分量值;H表示Heaviside阶跃函数;Ki表示第i个背应力分量的方向;mi表示第i个背应力分量的关键状态;<>表示Macauley括号,当x>0时=x,当x≤0时=0。

5.根据权利要求4所述的基于模式搜索法的双相钢Ohno‑WangⅠ本构模型参数优化方法,其特征在于,所述的Heavisied函数具体为:在Ohno‑Wang Ⅰ本构模型中Heavisied函数的计算结果与背应力的模及模型第二参数有关,获取方法如下所示:式中:||αi||为背应力分量的模;ri表示第i个背应力分量的模型第二参数。