一种复杂工况下通用弹粘塑性本构关系构建方法及介质技术

本发明专利技术公开了一种复杂工况下通用弹粘塑性本构关系构建方法及介质，方法包括：S1：构建非耦合式的本构关系的基本模型，基本模型中至少引入四个影响因子；S2：根据不同的影响因子，确定对应的独立函数关系；S3：根据影响因子和对应的独立函数关系，获得最终的非耦合式的本构关系。本发明专利技术引入了多个影响因子建立了非耦合式的本构关系模型，能够精确描述多种应力状态下力学行为，同时保证了仿真计算效率。同时保证了仿真计算效率。同时保证了仿真计算效率。

【技术实现步骤摘要】
一种复杂工况下通用弹粘塑性本构关系构建方法及介质


[0001]本专利技术涉及弹塑性本构关系建模
，特别涉及一种复杂工况下通用弹粘塑性本构关系构建方法及介质。

技术介绍

[0002]CAE(Computer Aided Engineering)是用计算机辅助求解复杂工程和产品结构强度、刚度、屈曲稳定性、动力响应、热传导、三维多体接触、弹塑性等力学性能的分析计算以及结构性能的优化设计等问题的一种近似数值分析方法。
[0003]当前CAE仿真采用的传统塑性标量本构中，冯
·
米塞斯等效应力一般为等效塑性应变的单因子函数，而单组函数参数无法精确描述多种应力状态下力学行为。
[0004]动载荷是指随时间作明显变化的载荷，即具有较大加载速率的载荷，在工程实际中，高速运行的构件，由加速度引起的载荷一般称为动载荷。
[0005]目前动载荷下需要多条单轴拉伸测试获得的应力
‑
应变曲线用于插值，并且认为该插值曲线全段应变率恒定，然而这与实际不符，因为颈缩后应变率会极剧增加。
[0006]基于现有的弹粘塑性本构关系，在当考虑多因子本构时，又会因计算层级增加，会严重影响计算效率，不利于工业应用。

技术实现思路

[0007]为解决上述问题，本专利技术提供了一种复杂工况下通用弹粘塑性本构关系构建方法及介质，建立了非耦合式的本构关系模型，解决考虑多因子本构时，计算效率低的问题。
[0008]本专利技术提供了一种复杂工况下通用弹粘塑性本构关系构建方法，具体技术方案如下：
[0009]S1：构建非耦合式的本构关系的基本模型，所述基本模型中包括由若干不同影响因子构建的独立函数关系；所述基本关系模型中至少引入四个影响因子；
[0010]S2：根据不同的所述影响因子，确定对应的独立函数关系；
[0011]S3：根据所述影响因子和对应的独立函数关系，获得最终的非耦合式的本构关系。
[0012]进一步的，所述非耦合式的本构关系的基本模型，表示如下：
[0013][0014]其中，表示von Mises等效应力(冯
·
米塞斯等效应力)，a
n
表示不同的影响因子，f
n
()表示对应的独立函数关系。
[0015]进一步的，所述影响因子包括等效塑性应变、应力三轴度、归一化洛德角参数和等效塑性应变率，包含等效塑性应变、应力三轴度、归一化洛德角参数和等效塑性应变率，四个影响因子的所述非耦合式的本构关系，表示如下：
[0016][0017]其中，表示von Mises等效应力(冯
·
米塞斯等效应力)，表示等效塑性应变，η
表示应力三轴度，表示归一化洛德角参数，表示等效塑性应变率，表示等效塑性应变对应的独立函数关系，f2(η)表示应力三轴度影响因子对应的独立函数关系，表示归一化洛德角参数影响因子对应的独立函数关系，表示等效塑性应变率影响因子对应的独立函数关系。
[0018]进一步的，所述影响因子还包括温度，包含等效塑性应变、应力三轴度、归一化洛德角参数、等效塑性应变率和温度，五个影响因子的所述非耦合式的本构关系，表示如下：
[0019][0020]其中，T表示温度，f5(T)表示温度影响因子对应的独立函数关系。
[0021]进一步的，所述独立函数关系为单值函数或由数对序列构成的曲线关系。
[0022]进一步的，所述应力三轴度包括第一应力三轴度和第二应力三轴度，所述第一应力三轴度，表示如下：
[0023][0024]所述第二应力三轴度，表示如下：
[0025][0026]其中，p表示静水压强，表示von Mises等效应力(冯
·
米塞斯等效应力)。
[0027]进一步的，所述等效塑性应变率影响因子对应的独立函数关系，包括：
[0028][0029]其中，V1为本构参数。
[0030]进一步的，所述温度影响因子对应的独立函数关系，包括：
[0031][0032]或
[0033][0034]其中，T
ref
为软化参考温度，T
melt
为熔化温度，m为软化凸凹性因子，T
0.5
为材料强度为参考强度一半时对应的温度，g3(T)＝exp(
‑
x)/(1+exp(
‑
x))，x＝(T
‑
T
0.5
)/r
scale
，r
scale
为缩放因子。
[0035]本专利技术还提供了一种存储介质，所述存储介质上存储有复杂工况下通用弹粘塑性本构关系的构建程序以及针对本构关系模型设定的仿真计算方法，所述复杂工况下通用弹粘塑性本构关系的构建程序被处理器执行时实现上述所述的复杂工况下通用弹粘塑性本构关系构建方法的步骤。
[0036]本专利技术的有益效果如下：
[0037]本专利技术提供了非耦合式的本构关系的基本模型，引入了多个影响因子，建立了非耦合式的本构关系模型，通过引入多组独立函数关系实现精确描述多种应力状态下力学行为，同时基于该非耦合式的本构关系模型，不用设计过多的计算方式方法，能够在考虑多因子本构时，解决仿真计算效率低的问题。
附图说明
[0038]图1是本专利技术方法流程示意图。
具体实施方式
[0039]在下面的描述中对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0040]在本专利技术实施例的描述中，需要说明的是，指示方位或位置关系为基于附图所示的方位或位置关系，或者是该专利技术产品使用时惯常摆放的方位或位置关系，或者是本领域技术人员惯常理解的方位或位置关系，或者是该专利技术产品使用时惯常摆放的方位或位置关系，仅是为了便于描述本专利技术和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本专利技术的限制。此外，术语“第一”、“第二”仅用于区分描述，而不能理解为指示或暗示相对重要性。
[0041]在本专利技术实施例的描述中，还需要说明的是，除非另有明确的规定和限定，术语“设置”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是直接连接，也可以通过中间媒介间接连接。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本专利技术中的具体含义。
[0042]实施例1
[0043]本专利技术的实施例1公开了一种复杂工况下通用弹粘塑性本构关系构建方法，如图1所示，具体步骤流程如下：
[0044]S1：构建非耦合式的本构关系的基本模型，所述基本模型中包括由若干不同影响因子构建的独立函数关系，表示如下：
[0045][0046]其中，表本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种复杂工况下通用弹粘塑性本构关系构建方法，其特征在于，包括：S1：构建非耦合式的本构关系的基本模型，所述基本模型中包括由若干不同影响因子构建的独立函数关系；S2：根据不同的所述影响因子，确定对应的独立函数关系；S3：根据所述影响因子和对应的独立函数关系，获得最终的非耦合式的本构关系。2.根据权利要求1所述的复杂工况下通用弹粘塑性本构关系构建方法，其特征在于，所述非耦合式的本构关系的基本模型，表示如下：其中，表示等效应力，a
n
表示不同的影响因子，f
n
()表示对应的独立函数关系。3.根据权利要求1所述的复杂工况下通用弹粘塑性本构关系构建方法，其特征在于，所述影响因子包括等效塑性应变、应力三轴度、归一化洛德角参数和等效塑性应变率，包含等效塑性应变、应力三轴度、归一化洛德角参数和等效塑性应变率，四个影响因子的所述非耦合式的本构关系，表示如下：其中，表示von Mises等效应力，表示等效塑性应变，η表示应力三轴度，表示归一化洛德角参数，表示等效塑性应变率，表示等效塑性应变对应的独立函数关系，f2(η)表示应力三轴度影响因子对应的独立函数关系，表示归一化洛德角参数影响因子对应的独立函数关系，表示等效塑性应变率影响因子对应的独立函数关系。4.根据权利要求3所述的复杂工况下通用弹粘塑性本构关系构建方法，其特征在于，所述影响因子还包括温度，包含等效塑性应变、应力三轴度、归一化洛德角参数、等效塑性应变率和温度，五个影响因子的所述非耦合式的本构关系，表示如下：其中，T表示温度，f5(T)表示温度影响因子对应的独立函数关系。5.根据权利要求1...

【专利技术属性】
技术研发人员：高峰，范吉富，梁宾，刘应波，姜子涵，王腾腾，黄荣亚，陈悟果，赵岩，王扬卫，
申请(专利权)人：北京理工大学重庆创新中心，
类型：发明
国别省市：