CN101727541A - 伞衣织物透气性计算方法 - Google Patents

Abstract

一种伞衣织物透气性计算方法，该方法首次引用多孔介质边界条件来模拟降落伞伞衣织物的透气性，有效而真实地模拟了气流流过透气体的流场结构，解决了目前亟待解决的伞衣织物透气性的模拟问题。计算方法所采用计算步骤为：(1)使用伞衣材料的透气性实验数据拟合得到一个压力损失与透气速度之间的二次多项式关系；(2)设定伞衣厚度Δm，使用二次多项式的系数求出伞衣透气性参数：α和C₂；(3)使用参数Δm、α和C₂对多孔介质边界条件进行设置。经过以上设置，则可得到伞衣“压力损失-透气速度”的二次关系式。在对降落伞流场进行求解的过程中，即可根据来流速度v得出伞衣内外的压力差Δp。

Description

伞衣织物透气性计算方法

所用技术领域

本计算方法涉及降落伞伞衣的透气性研究，首次引用多孔介质边界条件来模拟降落伞伞衣织物的透气性，有效而真实地模拟了气流流过透气体的流场结构，解决了目前亟待解决的伞衣织物透气性的模拟问题。

背景技术

透气性是降落伞织物的重要指标之一，降落伞的气动特性很大程度上是由降落伞织物的透气性决定的。诸如伞衣充气时间、开伞动载、稳定性以及下降速度等降落伞的重要性能都与伞衣织物的透气性有关。但是，由于降落伞是个柔性透气体，它在工作过程中的变形、透气性都影响了伞衣周围压力场的变化，压力场的变化反过来对伞衣的形状有影响，由此形成了复杂的相互作用的过程。在降落伞开伞过程中，这一现象尤为复杂。因此建立何种合适的模型对其进行气动力研究一直是困扰降落伞理论向前发展的难点和盲点问题。最近有研究者采用商用软件对伞衣织物透气性进行了研究(贾贺，荣伟，陈国良，基于LS-DYNA的降落伞伞衣织物透气性参数仿真验证[J]，航天返回与遥感，2009，30(1)：15-19)但从文献中可看出，随着速度的减少计算结果与试验结果误差增大，当流速为20m/s时，二者误差达到了9％，进一步减小流速，误差会进一步增大，因此该方法在流速低时不可靠。

发明内容

为了解决目前亟待解决的伞衣织物透气性的模拟问题。首次引用多孔介质边界条件来模拟降落伞伞衣织物的透气性，有效而真实地模拟了气流流过透气体的流场结构，求得气流透过透气体时的压力损失。

本计算方法所采用技术方案是：

第一步：使用伞衣材料的透气性实验数据拟合得到一个压力损失与透气速度之间的二次多项式关系。

对多孔介质进行一维简化，“压力损失-透气速度”关系由以下公式给定：

$\mathrm{\Δp}=(\frac{\mathrm{\μ}}{\mathrm{\α}}v+\frac{1}{2}{C}_2{\mathrm{\ρv}}^2)\mathrm{\Δm}---\left(1\right)$

其中μ是流体粘性系数，ρ为流体密度，α是介质的透气系数，C₂是压力跃变系数，v是垂直于介质表面的速度大小(即透过伞衣表面的气流速度)，Δm是介质厚度(即伞衣厚度)。α、C₂和Δm是多孔介质边界条件中需要设置的参数。从公式(1)中可见，压力损失Δp表示为透气速度v的二次函数。

对伞衣材料进行透气性试验，获得不同气流速度v对应的压力损失Δp。对试验结果进行二次多项式插值拟合。得到v压力损失Δp之间的二次函数关系式：

Δp＝av²+bv           (2)

式中，a，b为已知的常数。比较(1)式和(2)式，可得：

$a=\frac{1}{2}{C}_2\mathrm{\ρ\Δm},$ $b=\frac{\mathrm{\μ}}{\mathrm{\α}}\mathrm{\Δm}---\left(3\right)$

第二步：设定伞衣厚度Δm，使用二次多项式的系数求出伞衣透气性参数：α和C₂。

透过伞衣的气体密度ρ和气体粘性系数μ可根据降落伞工作高度来定，因此ρ和μ是已知数。那么，根据伞衣材料，给定伞衣厚度Δm，则可根据(3)式求出伞衣透气性参数：α和C₂

第三步：使用参数Δm、α和C₂对多孔介质边界条件进行设置。对降落伞流场进行求解的过程中，即可根据来流速度v得出伞衣内外的压力差Δp。

边界条件进行设置即：把求得的α和C₂代入(1)中，给定伞衣厚度Δm，则可得到伞衣“压力损失-透气速度”的二次关系式。

本计算方法的有益效果是：模型更接近伞衣的物理模型，结果更可靠，计算过程简单。

附图说明

图1是二次多项式曲线拟合结果。

图2是透气性验证计算模型的边界条件和网格示意图。

图中：1.速度入口，2.多孔介质(伞衣)，3.压力出口

具体实施方式

下面以602锦丝66绸为例进行计算模型验证。

首先对多孔介质进行一维简化，“压力损失-透气速度”关系由以下公式给定：

$\mathrm{\Δp}=(\frac{\mathrm{\μ}}{\mathrm{\α}}v+C_2\frac{1}{2}{\mathrm{\ρv}}^2)\mathrm{\Δm}---\left(1\right)$

其中μ是流体粘性系数，ρ为流体密度，α是介质的透气系数，C₂是压力跃变系数，v是垂直于介质表面的速度大小，Δm是介质厚度。α、C₂和Δm是多孔介质边界条件中需要设置的参数。

从公式(1)中可见，压力损失Δp表示为透气速度v的二次函数。在实际计算中，流体粘性系数μ、流体密度ρ和介质厚度Δm为已知量，还有α和C₂两个参数需要确定。如果已知一组透气速度和压力损失的实验数据，则可以通过对实验数据进行二次多项式插值的方法确定α和C₂两个参数。表1中的第一二列为602锦丝66绸透气速度与压力损失的一组实验数据，采用过原点的二次多项式对这组数据进行拟合(如图1所示)可得：

Δp＝133.34v²+79.386v        (2)

比较式(2)和式(1)，由对应常数项相等可得：

$C_2\frac{1}{2}\mathrm{\ρ\Δm}=133.34---\left(3\right)$

$\frac{\mathrm{\μ}}{\mathrm{\α}}\mathrm{\Δm}=79.386---\left(4\right)$

标准大气条件下空气密度ρ为1.225kg/m³，粘性系数μ为1.726×10^-5kg/ms，设介质厚度Δm为0.0005m，则由以上两式计算可得：

C₂＝435395.9，α＝1.1270×10^-10

为了验证采用多孔介质边界条件模拟降落伞伞衣透气性的可靠性，制作了一个如图2所示的计算模型。该模型的计算域为二维长方形“通道”，侧边指定为周期边界，多孔介质边界设置在“通道”中间位置，相当于在通道中间加了一层“透气膜”，在速度入口指定空气流速，压力出口设置为标准大气压。使用所编写的程序对该模型进行计算直至收敛后，可以得到速度入口和压力出口之间的压力差，此压力差即为多孔介质边界模拟的“透气膜”在当前流速(即透气速度)下的压力损失。

根据之前对602锦丝66绸透气参数的计算，设定多孔介质边界的参数为：Δm＝0.0005m C₂＝435395.9，α＝1.1270×10^-10。将速度入口的空气流速分别设定为表1中的速度，即可计算得到相应的压力损失数值。压力损失计算值与实验值的对比如表1的第二三列所示，可以看出，在流速非常低的情况下，误差也不超过5％。表1中的压力损失计算值与使用之前得到的拟合公式Δp＝133.34v²+79.386v计算的结果吻合较好。因此，使用公式(1)可以较精确地将透气速度与压力损失之间的关系加入到流场数值模拟中。

以上验证表明，可以使用多孔介质边界条件模拟降落伞伞衣的透气特性，其具体设置步骤为：(1)使用伞衣材料的透气性实验数据拟合得到一个压力损失与透气速度之间的二次多项式关系；(2)设定伞衣厚度Δm，使用二次多项式的系数求出α和C₂；(3)使用参数Δm、α和C₂对多孔介质边界条件进行设置，则可得到伞衣“压力损失-透气速度”的二次关系式。在对降落伞流场进行求解的过程中，即可根据来流速度v得出伞衣内外的压力差Δp。

表1 602锦丝66绸透气速度与压力损失的计算值与实验值对比

空气流速m/s	压力损失Pa(实验值)	压力损失Pa(计算值)	误差
				0.972	195.64	203.14	3.6％
1.225	293.47	297.34	1.3％
				1.445	391.29	393.13	0.5％
1.61	489.11	483.44	1.1％
				1.835	588.94	594.65	0.9％

Claims (1)

1.一种伞衣织物透气性计算方法，其特征在于包括以下过程：

第一步：使用伞衣材料的透气性实验数据拟合得到一个压力损失与透气速度之间的二次多项式关系；

对多孔介质进行一维简化，“压力损失-透气速度”关系由以下公式给定：

$\mathrm{\Δp}=(\frac{\mathrm{\μ}}{\mathrm{\α}}v+\frac{1}{2}{C}_2{\mathrm{\ρv}}^2)\mathrm{\Δm}---\left(1\right)$

其中μ是流体粘性系数，ρ为流体密度，α是介质的透气系数，C₂是压力跃变系数，v是垂直于介质表面的速度大小(即透过伞衣表面的气流速度)，Δm是伞衣厚度；α、C₂和Δm是多孔介质边界条件中需要设置的参数；从公式(1)中可见，压力损失Δp表示为透气速度v的二次函数；

对伞衣材料进行透气性试验，获得不同气流速度v对应的压力损失Δp；对试验结果进行二次多项式插值拟合；得到v压力损失Δp之间的二次函数关系式：

Δp＝av²+bv        (2)

式中，a，b为已知的常数。比较(1)式和(2)式，可得：

$a=\frac{1}{2}{C}_2\mathrm{\ρ\Δm},$ $b=\frac{\mathrm{\μ}}{\mathrm{\α}}\mathrm{\Δm}---\left(3\right)$

第二步：设定伞衣厚度Δm，使用二次多项式的系数求出伞衣透气性参数：α和C₂；

透过伞衣的气体密度ρ和气体粘性系数μ可根据降落伞工作高度来定，因此ρ和μ是已知数；那么，根据伞衣材料，给定伞衣厚度Δm，则可根据(3)式求出伞衣透气性参数：α和C₂

第三步：使用参数Δm、α和C₂对多孔介质边界条件进行设置；对降落伞流场进行求解的过程中，即可根据来流速度v得出伞衣内外的压力差Δp；

边界条件进行设置即把求得的α和C₂代入(1)中，给定伞衣厚度Δm，则可得到伞衣“压力损失-透气速度”的二次关系式；

Images