平面应力状态下推进剂精细本构关系及数值化方法

1.本发明涉及平面应力状态下推进剂精细本构关系及数值化方法，属于推进剂粘弹性本构关系技术领域。


背景技术：

2.推进剂属于典型的粘弹性材料，其泊松比、模量等力学参数与加载时间以及加载温度息息相关。一直以来，为了计算方便，将推进剂的泊松比处理成常数。仿真结果表明，泊松比的细微变化，将会严重影响固体火箭发动机结构完整性计算的结果。为此，推进剂的本构关系有必要考虑泊松比的粘弹性特性。
3.近年来，计算机仿真手段不断更新，已经有相关研究成果将推进剂的三维本构关系考虑粘弹性泊松比的效应，并且给出了详细的数值化方法。然而，对于二维平面应力问题，考虑粘弹性泊松比的推进剂本构关系尚未有具体的研究成果。


技术实现要素：

4.为解决现有技术的不足，本发明的目的在于提供平面应力状态下推进剂精细本构关系及数值化方法，能对推进剂的粘弹性泊松比效应进行准确描述。
5.为了实现上述目标，本发明采用如下的技术方案：
6.一种平面应力状态下推进剂精细本构关系及数值化方法，包括以下步骤：
7.建立平面应力状态下一般弹性材料的应力应变关系；
8.基于应力应变关系，考虑粘弹性泊松比，建立球偏分解状态下的推进剂精细本构关系；
9.建立推进剂精细本构关系中偏张量部分的增量关系；
10.建立推进剂精细本构关系中球张量部分的增量关系。
11.进一步地，前述平面应力状态下弹性材料的应力应变关系表达式为：
[0012][0013]
式中，σ
11
、σ
22
以及σ
12
分别代表x、y以及xy方向的应力，ε
11
、ε
22
、ε
12
分别代表x，y， xy方向的应变，e、ν和g分别代表弹性材料的模量、泊松比和剪切模量。
[0014]
进一步地，前述建立球偏分解状态下的推进剂精细本构关系，包括以下步骤：
[0015]
建立平面应力条件下本构关系的球偏分解关系：
[0016][0017][0018]
式中，s
ij
和e
ij
分别表示弹性材料的偏应力以及偏应变张量，σ
kk
和ε
kk
分别表示弹性
材料的球应力以及球应变张量；
[0019]
在球偏分解关系中，考虑粘弹性泊松比，得到含粘弹性泊松比的推进剂本构关系：
[0020][0021][0022][0023]
式中，δ
ij
表示本构关系，σ
ij
()、s
ij
()和σ
kk
()分别表示推进剂的应力张量、偏应力张量以及球应力张量，e
ij
()和ε
kk
()分别表示推进剂的偏应变张量以及球应变张量，e()、ν()分别表示推进剂的松弛模量、粘弹性泊松比，t、θ以及ξ分别代表加载时间、考虑泊松比时温效应的缩减时间以及考虑松弛模量时温效应的缩减时间，τ、θ
′
以及ξ
′
分别代表t、θ以及ξ的被积数。
[0024]
进一步地，前述松弛模量的表达式为
[0025][0026]
式中，en和分别表示第n项松弛模量的两个参数，ne表示松弛模量prony级数的项数， e0为初始松弛模量；
[0027]
粘弹性泊松比表达式为：
[0028][0029]
式中，ν
∞
表示平衡泊松比，νn和分别表示第n项粘弹性泊松比的两个参数，
nν
表示粘弹性泊松比prony级数的项数。
[0030]
进一步地，前述建立偏张量部分的增量关系的步骤包括：
[0031]
对公式(1)，求解t
m+1
时刻和tm时刻的增长量得
[0032][0033]
式中，δs
ij
()、以及分别表示偏应力增量、偏应力增量第一分量、偏应力增量第二分量、偏应力增量第三分量、偏应力增量第四分量以及偏应力增量第五分量，其对应表达式分别为：
[0034][0035][0036][0037][0038]
[0039]
对以及求解，带入公式(3)得
[0040][0041][0042][0043][0044][0045][0046][0047][0048][0049]
式中，γ
ν
()、γe()、以及分别表示偏应力张量辅助第一变量、偏应力张量辅助第二变量、偏应力张量辅助第三变量、偏应力张量辅助第四变量、偏应力张量辅助第五变量、偏应力张量辅助第六变量、偏应力张量辅助第七变量以及偏应力张量辅助第八变量。
[0050]
进一步地，前述建立球张量部分的增量关系的步骤包括：
[0051]
对公式(2)，求解t
m+1
时刻和tm时刻的增长量得
[0052][0053]
式中，δσ
kk
()、以及分别表示球应力增量、球应力增量第一分量、球应力增量第二分量、球应力增量第三分量、球应力增量第四分量以及球应力增量第五分量，其对应表达式分别为：
[0054][0055][0056][0057]
[0058][0059]
对以及求解，带入公式(4)得
[0060][0061]
式中，分别表示球应力张量辅助第一变量、球应力张量辅助第二变量，其对应表达式分别为：
[0062][0063]
本发明所达到的有益效果：
[0064]
本发明是针对现有的推进剂平面应力本构关系不能精确描述粘弹性泊松比效应的问题，在现有的弹性平面应力本构关系基础上，考虑粘弹性泊松比效应，建立平面应力条件下的推进剂精细本构关系。在现有的球偏分解状态下，弹性平面应力本构关系基础上构建的考虑粘弹性泊松比的推进剂平面应力本构关系，不仅可以真实地反映推进剂的粘弹性泊松比效应，还可以方便地开展本构关系的数值化计算。
具体实施方式
[0065]
下面通过具体实施例对本发明技术方案做详细的说明，应当理解本技术实施例以及实施例中的具体特征是对本技术技术方案的详细的说明，而不是对本技术技术方案的限定，在不冲突的情况下，本技术实施例以及实施例中的技术特征可以相互组合。
[0066]
本实施例公开了一种平面应力状态下推进剂精细本构关系及数值化方法，实施步骤包括：
[0067]
1)建立平面应力状态下一般弹性材料的应力应变关系；
[0068]
2)基于应力应变关系，考虑粘弹性泊松比，建立球偏分解状态下的推进剂精细本构关系；
[0069]
3)建立推进剂精细本构关系中偏张量部分的增量关系；
[0070]
4)建立推进剂精细本构关系中球张量部分的增量关系。
[0071]
具体地，步骤1)中引入的描述平面应力状态下一般弹性材料的应力应变关系，如式(1) 所示；
[0072][0073]
其中，σ
11
、σ
22
以及σ
12
分别代表x、y以及xy方向的应力，ε
11
、ε
22
、ε
33
、ε
12
、ε
23
以及ε
13
分
别代表x，y，z、xy、yz以及xz方向的应变。e、ν和g分别代表弹性材料的模量、泊松比和剪切模量。
[0074]
将上式整理成矩阵的形式，
[0075][0076]
相应地，步骤2)，由以上式子整理可得弹性材料平面应力条件下本构关系的球偏分解关系
[0077][0078]
其中，s
ij
和e
ij
分别表示弹性材料的偏应力以及偏应变张量，σ
kk
和ε
kk
分别表示弹性材料的球应力以及球应变张量。
[0079]
进一步，步骤2)参考以上本构关系，可以得到平面应力下，考虑粘弹性泊松比的推进剂本构关系
[0080][0081][0082][0083]
其中，δ
ij
表示本构关系，σ
ij
()、s
ij
()和σ
kk
()分别表示推进剂的应力张量、偏应力张量以及球应力张量，e
ij
()和ε
kk
()分别表示推进剂的偏应变张量以及球应变张量，e()、ν()分别表示推进剂的松弛模量、粘弹性泊松比。t、θ以及ξ分别代表加载时间，考虑泊松比时温效应的缩减时间以及考虑松弛模量时温效应的缩减时间，τ、θ
′
以及ξ
′
分别代表t、θ以及ξ的被积数。
[0084]
松弛模量的表达式为
[0085][0086]
其中，en和分别表示第n项松弛模量的两个参数，ne表示松弛模量prony级数的项数。 e0为初始松弛模量。
[0087]
粘弹性泊松比的表达式为
[0088][0089]
其中，ν
∞
表示平衡泊松比，νn和分别表示第n项粘弹性泊松比的两个参数，n
ν
表示粘弹性泊松比prony级数的项数，代表初始泊松比。
[0090]
根据stieltjes卷积定理，式(4)和式(5)可进一步表示成
[0091][0092][0093]
步骤3)中本构方程偏张量部分表达式
[0094]
基于式(8)，在tm时刻有如下形式
[0095][0096]
基于式(8)，在t
m+1
时刻有
[0097][0098]
上面两式相减可得推进剂精细本构关系中偏张量部分的增量关系
[0099][0100]
上述表达式在形式上比较复杂，为便于数值分析，将其改写成如下形式
[0101][0102]
其中，δs
ij
()、以及分别表示偏应力增量、偏应力增量第一分量、偏应力增量第二分量、偏应力增量第三分量、偏应力增量第四分量以及偏应力增量第五分量。
[0103]
其中，
[0104][0105][0106][0107][0108][0109]
整理可得
[0110][0111][0112][0113]
[0114]
整理有
[0115][0116]
即
[0117][0118][0119]
其中，γ
ν
()、以及分别表示偏应力张量辅助第一变量、偏应力张量辅助第二变量、偏应力张量辅助第三变量以及偏应力张量辅助第四变量。是与粘弹性泊松比相关的张量，n取值范围为1至n
ν
。γ
ν
()是与粘弹性泊松比相关的标量。是与粘弹性泊松比相关的标量，n取值范围为1至n
ν
。
[0120][0121]
其中，γe()、以及分别表示偏应力张量辅助第五变量、偏应力张量辅助第六变量、偏应力张量辅助第七变量以及偏应力张量辅助第八变量。是与松弛模量相关的张量， n取值范围为1至ne。γe()是与松弛模量相关的标量。是与松弛模量相关的标量，n取值范围为1至ne。
[0122]
相应地，步骤4)中本构方程球张量部分表达式
[0123]
基于式(9)，在tm时刻有如下形式
[0124][0125]
基于式(9)，在t
m+1
时刻有
[0126][0127]
上面两式相减可得推进剂精细本构关系中球张量部分的增量关系
[0128][0129]
上述表达式在形式上比较复杂，为便于数值分析，将其改写成如下形式
[0130][0131]
其中，δσ
kk
()、以及分别表示球应力增量、球应力增量第一分量、球应力增量第二分量、球应力增量第三分量、球应力增量第四分量以及球应力增量第五分量。
[0132]
其中，
[0133][0134][0135][0136][0137][0138]
整理可得
[0139][0140][0141][0142][0143][0144]
其中，
[0145][0146]
其中，分别表示球应力张量辅助第一变量、球应力张量辅助第二变量。
[0147]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序
产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0148]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0149]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0150]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0151]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。