自然形高空气球发放阶段动力学的研究.pdfVIP

2007年中国浮空器大会论文集 229 自然形高空气球发放阶段动力学研究 张向强1 姜鲁华2 王根华2 郑荣廷2 丁克尔2 高 衡1 (1．中国科学院光电研究院，北京100080；2．中国科学院高能物理研究所，北京100049) 摘要：探讨了如何获得充气后自然形球体的形状这一研究自然形高空气球发射阶段动力学的关键问题，通过选择适当的单元算 法、接触算法和设置流固耦合参数，初步给出了半充气自然形球体的形状，该问题的解决使得自然形高空气球发放阶段动力学研究 成为可能，为研究自然形球体的冲击响应和加强带和多层头部设计等问题奠定了基础。 关键词：高空气球；自然形气球；流固耦合 中圈分类号：V211 张量与变形张量的积，从而物体的任何运动都可以 0引言 分解为变形和转动，应力应变的度量采用满足应力 张量的客观率。 ’零压高空气球一般由极薄的聚乙烯膜(厚度 接触算法：接触中必然包含碰撞和接触面问的 滑移，导致接触面法向和切向的速度瞬时不连续， 20pan，弹性模量108量级)制成，它内充轻质气体 (氦气或氢气)，升限高(距离地面30—40kin)，载给离散方程的时间积分带来困难，罚函数法是目前 重大(1t到数吨)，具有重要的科学应用价值…。 常用的方法。很多软件提供多种接触算法，但如何 由于高空气球体积大(几万到几百万立方米)， 选择接触类型，尤其是众多接触参数的控制，以及 升限高，所以地面充气体积占气球容积的比例很小， 单元网格的大小需反复试验。与其他结构接触不同 以升限30km为例，地面充气体积只占气球体积的的是，高空气球的接触厚度极薄，弹性模量很低， 1％左右。 给接触算法参数的调节带来更大困难。 高空气球从地面充气到发放是高空气球动力学 单元与流体的耦合算法：结构变形(有限元中 行为非常复杂的阶段，涉及风载荷、冲击过载、球 膜和内含轻质气体的耦合运动。 ler方法)之间能量的传递是流固耦合的关键，目前 参考文献[2]给出了地面不同充气体积下球体的动力学分析领域多采用任意拉格朗日一欧拉方法 二维形状，本文采用有限元方法，以体积很小的零 (Arbitrary蛔ge 压高空气球为例，研究在无风状态下，三维高空气 球发放前的充气形状和冲击动力学问题。 元中，实现结构和流体之间的耦合¨J，ALE算法中 各参数的调节是获得充气后形状的关键。 研究方法和分析工具 本文选择的商业软件LS—Dyna基本满足上述条 件。需要特别指出的是，涉及接触碰撞或ALE算法， 从有限元的角度分析，高空气球充气是膜单元 对计算机硬件配置要求很高，受硬件限制，本文只 组成的封闭结构在重力、浮力外部的空气和内含的 能采用小体积气球进行方法研究。 轻质氦气共同作用下的耦合运动。动力学行为包括 膜单元大变形、大转动、球膜褶皱和自接触、冲击 2高空气球的形状 响应、球膜与内外异质流体之间的耦合行为。 研究该动力学问题，必须解决单元算法、接触 自然形气球是最普遍采用的优化方案，基于4 算法和结构单元与流体的耦合算法等关键问题。 个假设：(1)关于中心轴轴对称；(2)与中心轴正 单元算法：由于基于无限小变形假设的有限元 交的截面上，每点的圆周应力和子午方向的应力相 方法中，单元刚体转动导致变形率不为零，正确的 等，同时圆周向应力为零；(3)浮升气体和周围空 方法应该是连续介质力学中基于有限变形的算法， 气的密度相等；