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

摘要

本发明公开了一种平面应力状态下推进剂精细本构关系及数值化方法，包括以下步骤：建立平面应力状态下一般弹性材料的应力应变关系，基于应力应变关系，考虑粘弹性泊松比，建立球偏分解状态下的推进剂精细本构关系，建立推进剂精细本构关系中偏张量部分的增量关系，建立推进剂精细本构关系中球张量部分的增量关系。

说明书

技术领域

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

背景技术

推进剂属于典型的粘弹性材料，其泊松比、模量等力学参数与加载时间以及加载温度息息相关。一直以来，为了计算方便，将推进剂的泊松比处理成常数。仿真结果表明，泊松比的细微变化，将会严重影响固体火箭发动机结构完整性计算的结果。为此，推进剂的本构关系有必要考虑泊松比的粘弹性特性。

近年来，计算机仿真手段不断更新，已经有相关研究成果将推进剂的三维本构关系考虑粘弹性泊松比的效应，并且给出了详细的数值化方法。然而，对于二维平面应力问题，考虑粘弹性泊松比的推进剂本构关系尚未有具体的研究成果。

发明内容

为解决现有技术的不足，本发明的目的在于提供平面应力状态下推进剂精细本构关系及数值化方法，能对推进剂的粘弹性泊松比效应进行准确描述。

为了实现上述目标，本发明采用如下的技术方案：

一种平面应力状态下推进剂精细本构关系及数值化方法，包括以下步骤：

建立平面应力状态下一般弹性材料的应力应变关系；

基于应力应变关系，考虑粘弹性泊松比，建立球偏分解状态下的推进剂精细本构关系；

建立推进剂精细本构关系中偏张量部分的增量关系；

建立推进剂精细本构关系中球张量部分的增量关系。

进一步地，前述平面应力状态下弹性材料的应力应变关系表达式为：

式中，σ

进一步地，前述建立球偏分解状态下的推进剂精细本构关系，包括以下步骤：

建立平面应力条件下本构关系的球偏分解关系：

式中，S

在球偏分解关系中，考虑粘弹性泊松比，得到含粘弹性泊松比的推进剂本构关系：

式中，δ

进一步地，前述松弛模量的表达式为

式中，E

粘弹性泊松比表达式为：

式中，ν

进一步地，前述建立偏张量部分的增量关系的步骤包括：

对公式(1)，求解t

式中，ΔS

对

式中，

进一步地，前述建立球张量部分的增量关系的步骤包括：

对公式(2)，求解t

式中，Δσ

对

式中，

本发明所达到的有益效果：

本发明是针对现有的推进剂平面应力本构关系不能精确描述粘弹性泊松比效应的问题，在现有的弹性平面应力本构关系基础上，考虑粘弹性泊松比效应，建立平面应力条件下的推进剂精细本构关系。在现有的球偏分解状态下，弹性平面应力本构关系基础上构建的考虑粘弹性泊松比的推进剂平面应力本构关系，不仅可以真实地反映推进剂的粘弹性泊松比效应，还可以方便地开展本构关系的数值化计算。

具体实施方式

下面通过具体实施例对本发明技术方案做详细的说明，应当理解本申请实施例以及实施例中的具体特征是对本申请技术方案的详细的说明，而不是对本申请技术方案的限定，在不冲突的情况下，本申请实施例以及实施例中的技术特征可以相互组合。

本实施例公开了一种平面应力状态下推进剂精细本构关系及数值化方法，实施步骤包括：

1)建立平面应力状态下一般弹性材料的应力应变关系；

2)基于应力应变关系，考虑粘弹性泊松比，建立球偏分解状态下的推进剂精细本构关系；

3)建立推进剂精细本构关系中偏张量部分的增量关系；

4)建立推进剂精细本构关系中球张量部分的增量关系。

具体地，步骤1)中引入的描述平面应力状态下一般弹性材料的应力应变关系，如式(1) 所示；

其中，σ

将上式整理成矩阵的形式，

相应地，步骤2)，由以上式子整理可得弹性材料平面应力条件下本构关系的球偏分解关系

其中，S

进一步，步骤2)参考以上本构关系，可以得到平面应力下，考虑粘弹性泊松比的推进剂本构关系

其中，δ

松弛模量的表达式为

其中，E

粘弹性泊松比的表达式为

其中，ν

根据Stieltjes卷积定理，式(4)和式(5)可进一步表示成

步骤3)中本构方程偏张量部分表达式

基于式(8)，在t

基于式(8)，在t

上面两式相减可得推进剂精细本构关系中偏张量部分的增量关系

上述表达式在形式上比较复杂，为便于数值分析，将其改写成如下形式

其中，ΔS

其中，

整理可得

整理有

即

其中，

其中，

相应地，步骤4)中本构方程球张量部分表达式

基于式(9)，在t

基于式(9)，在t

上面两式相减可得推进剂精细本构关系中球张量部分的增量关系

上述表达式在形式上比较复杂，为便于数值分析，将其改写成如下形式

其中，Δσ

其中，

整理可得

其中，

其中，

本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。