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

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

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

本计算方法所采用技术方案是：
第一步：使用伞衣材料的透气性实验数据拟合得到一个压力损失与透气速度之间 的二次多项式关系。

对多孔介质进行一维简化，"压力损失-透气速度"关系由以下公式给定：
其中ii是流体粘性系数，P为流体密度，a是介质的透气系数，C2是压力跃变系 数，v是垂直于介质表面的速度大小（即透过伞衣表面的气流速度），Am是介质厚度（即 伞衣厚度）。a、G和Am是多孔介质边界条件中需要设置的参数。从公式（1)中可见，压 力损失A p表示为透气速度v的二次函数。

对伞衣材料进行透气性试验，获得不同气流速度v对应的压力损失A p。对试验结 果进行二次多项式插值拟合。得到v压力损失Ap之间的二次函数关系式： [0011] Ap = av2+bv (2)

式中，a，b为已知的常数。比较（1)式和（2)式，可得：[0013] <formula>formula see original document page 4</formula>[0014] 第二步：设定伞衣厚度Am，使用二次多项式的系数求出伞衣透气性参数：a和 C2。

透过伞衣的气体密度P和气体粘性系数可根据降落伞工作高度来定，因此P 和P是已知数。那么，根据伞衣材料，给定伞衣厚度Am，则可根据（3)式求出伞衣透气性 参数：a和C2

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

边界条件进行设置即：把求得的a和C2代入（1)中，给定伞衣厚度Am，则可得到 伞衣"压力损失_透气速度"的二次关系式。

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

具体实施方式
下面以602锦丝66绸为例进行计算模型验证。
首先对多孔介质进行一维简化，"压力损失_透气速度"关系由以下公式给定： ," 1 '
其中是流体粘性系数，P为流体密度，a是介质的透气系数，C2是压力跃变系 数，v是垂直于介质表面的速度大小，Am是介质厚度。a 、 （:2和Am是多孔介质边界条件 中需要设置的参数。
从公式（1)中可见，压力损失Ap表示为透气速度v的二次函数。在实际计算中， 流体粘性系数P 、流体密度P和介质厚度Am为已知量，还有a和G两个参数需要确定。 如果已知一组透气速度和压力损失的实验数据，则可以通过对实验数据进行二次多项式插 值的方法确定a和C2两个参数。表1中的第一二列为602锦丝66绸透气速度与压力损 失的一组实验数据，采用过原点的二次多项式对这组数据进行拟合（如图1所示）可得： [0027] Ap = 133. 34v2+79. 386v (2)

比较式（2)和式（1)，由对应常数项相等可得： 1
标准大气条件下空气密度P为1. /m?粘性系数y为1. 726 X 10—5kg/ms，设 介质厚度A m为0. 0005m，则由以上两式计算可得：[0032] C2 = 435395. 9， a = 1. 1270X 10—10
为了验证采用多孔介质边界条件模拟降落伞伞衣透气性的可靠性，制作了 一个如 图2所示的计算模型。该模型的计算域为二维长方形"通道"，侧边指定为周期边界，多孔介 质边界设置在"通道"中间位置，相当于在通道中间加了一层"透气膜"，在速度入口指定空 气流速，压力出口设置为标准大气压。使用所编写的程序对该模型进行计算直至收敛后，可 以得到速度入口和压力出口之间的压力差，此压力差即为多孔介质边界模拟的"透气膜"在 当前流速（即透气速度）下的压力损失。

根据之前对602锦丝66绸透气参数的计算，设定多孔介质边界的参数为：Am = 0.0005m C2 = 435395.9， a = 1. 1270 X 10—1CI。将速度入口的空气流速分别设定为表1中 的速度，即可计算得到相应的压力损失数值。压力损失计算值与实验值的对比如表l的第 二三列所示，可以看出，在流速非常低的情况下，误差也不超过5% 。表1中的压力损失计算 值与使用之前得到的拟合公式Ap= 133.34v2+79.386v计算的结果吻合较好。因此，使用 公式（1)可以较精确地将透气速度与压力损失之间的关系加入到流场数值模拟中。 [0035] 以上验证表明，可以使用多孔介质边界条件模拟降落伞伞衣的透气特性，其具体
设置步骤为：
（l)使用伞衣材料的透气性实验数据拟合得到一个压力损失与透气速度之间 的二次多项式关系；
（2)设定伞衣厚度Am，使用二次多项式的系数求出a和C2 ;
(3)使用 参数Am、 a和（:2对多孔介质边界条件进行设置，则可得到伞衣"压力损失-透气速度"的 二次关系式。在对降落伞流场进行求解的过程中，即可根据来流速度v得出伞衣内外的压 力差Ap。

更多关于:透气性测试仪