SPH 统一算法对自由流体冲击弹性结构体问题模拟.pdfVIP

振　动　与　冲　击 第３４卷第１６期 ＪＯＵＲＮＡＬ ＯＦ ＶＩＢＲＡＴＩＯＮ ＡＮＤ ＳＨＯＣＫ Ｖｏｌ．３４ Ｎｏ．１６ ２０１５　  ＳＰＨ统一算法对自由流体冲击弹性结构体问题模拟 沈雁鸣，何　琨，陈坚强，袁先旭 （中国空气动力研究与发展中心，四川　绵阳　６２１０００） 　　摘　要：对自由流体冲击弹性结构体涉及的自由界面运动、弹塑性变形、流固耦合作用等问题进行模拟时要求能 有效计算介质运动变形与流固间耦合作用。利用光滑粒子流体动力学方法将流体、弹性体视为连续介质，用统一的Ｎ－Ｓ 方程与虚拟压缩状态方程求解。基于初始粒子类型划分提出新的介质接触算法，通过在结构体表面增加结构体表面粒 子，使流体与弹性体介质间相互作用实现稳定传递。计算结果表明，所用ＳＰＨ统一算法与接触算法能准确有效模拟自由 流体冲击弹性结构体变形回弹过程。 关键词：流固耦合；变形回弹；ＳＰＨ方法；接触算法 中图分类号：Ｏ３５３．４　　　文献标志码：Ａ ＤＯＩ：１０．１３４６５／ｊ．ｃｎｋｉ．ｊｖｓ．２０１５．１６．０１１ Ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ ｏｆｆｒｅｅ ｓｕｒｆａｃｅ ｆｌｏｗ ｉｍｐａｃｔｉｎｇ ｅｌａｓｔｉｃ ｓｔｒｕｃｔｕｒｅ ｗｉｔｈ ＳＰＨ ｕｎｉｆｏｒｍ ｍｅｔｈｏｄ ＳＨＥＮ Ｙａｎｍｉｎｇ，ＨＥＫｕｎ，ＣＨＥＮＪｉａｎｑｉａｎｇ，ＹＵＡＮＸｉａｎｘｕ （Ｃｈｉｎａ Ａｅｒｏｄｙｎａｍｉｃｓ Ｒｅｓｅａｒｃｈ ａｎｄ Ｄｅｖｅｌｏｐｍｅｎｔ Ｃｅｎｔｅｒ，Ｍｉａｎｙａｎｇ６２１０００，Ｃｈｉｎａ） 　　Ａｂｓｔｒａｃｔ：　Ｆｒｅｅ ｓｕｒｆａｃｅ ｆｌｏｗ ｉｍｐａｃｔｉｎｇ ｅｌａｓｔｉｃ ｓｔｒｕｃｔｕｒｅ ｉｓ ａ ｃｏｍｍｏｎ ｐｒｏｂｌｅｍ，ｕｓｕａｌｌｙ ｒｅｆｅｒｒｉｎｇ ｔｏ ｉｎｔｅｒｆａｃｅ ｔｒａｃｋ， ｅｌａｓｔｉｃ ｏｒ ｐｌａｓｔｉｃ ｓｔｒｕｃｔｕｒａｌ ｄｅｆｏｒｍａｔｉｏｎ ａｎｄ ｆｌｕｉｄｓｔｒｕｃｔｕｒｅ ｉｎｔｅｒａｃｔｉｏｎ．Ｔｏ ｓｉｍｕｌａｔｅ ｔｈｅｓｅ ｐｒｏｂｌｅｍｓ，ｉｔ ｈａｓ ｔｏ ｃａｐｔｕｒｅ ｔｈｅ ｄｅｆｏｒｍａｔｉｏｎ ｏｆ ｉｎｔｅｒｆａｃｉａｌ ｆｌｏｗ ａｎｄ ｃｏｍｐｕｔｅ ｔｈｅ ｆｌｕｉｄ ｓｔｒｕｃｔｕｒｅ ｃｏｕｐｌｉｎｇ ｉｎｔｅｒａｃｔｉｏｎ ａｃｃｕｒａｔｅｌｙ．Ａ ｓｍｏｏｔｈｉｎｇ ｐａｒｔｉｃｌｅ ｈｙｄｒｏｄｙｎａｍｉｃｓ （ＳＰＨ）ｕｎｉｆｏｒｍ ｍｅｔｈｏｄ ｗｈｅｒｅ ｂｏｔｈ ｆｌｕｉｄ ａｎｄ ｓｏｌｉｄ ｐｈａｓｅｓ ａｒｅ ｄｅｓｃｒｉｂｅｄ ｂｙ ＳＰＨ ｗａｓ ｐｒｅｓｅｎｔｅｄ ａｎｄ ｔｈｅ ｕｎｉｆｏｒｍ ｃｏｎｔｉｎｕｏｕｓ ｍｅｄｉｕｍ ｇｏｖｅｒｎｉｎｇ ｅｑｕａｔｉｏｎｓ ａｎｄ ｗｅｅｋ ｐｒｅｓｓｕｒｅ ｓｔａｔｅ ｅｑｕａｔｉｏｎ ｗｅｒｅ ｓｏｌｖｅｄ．Ａ ｎｅｗ ｃｏｎｔａｃｔ ａｌｇｏｒｉｔｈｍ ｗｉｔｈ ｔｈｅ ｓｔｒｕｃｔｕｒｅ ｓｕｒｆａｃｅ ｐａｒｔｉｃｌｅｓ ｔｒｅａｔｅｄ ａｓ ｖｉｒｔｕａｌ ｐａｒｔｉｃｌｅｓ ｗａｓ ｄｅｖｅｌｏｐｅｄ．Ｔｈｅ ｖｉｒｔｕａｌ ｐａｒｔｉｃｌｅｓ ｗｅｒｅ ｄｅｆｉｎｅｄ ｉｎ ｉｎｉｔｉａｌ ｐａｒｔｉｃｌｅ ｃｏｎｆｉｇｕｒａｔｉｏｎ，ｃｏｎｔａｉｎｉｎｇ ｂｏｔｈ ｆｌｕｉｄ ａｎｄ ｓｔｒｕｃｔｕｒｅ ｐａｒｔｉｃｌｅ ｃｈａｒａｃｔｅｒｓ ａｎｄ ｐｒｏｄｕｃｉｎｇ ｒｅｐｕｌｓｉｖｅ ｆｏｒｃｅ ｔｏ ｐｒｅｖｅｎｔ ｐｅｎｅｔｒａｔｉｏｎ．Ｃｏｍｐａｒｉｎｇｔｈｅ ｎｕｍｅｒｉｃａｌｒｅｓｕｌｔｓｗｉｔｈ ｅｘｐｅｒｉｍｅｎｔａｌｏｎｅｓ，ｉｔｉｓｄｅｍｏｎｓｔｒａｔｅｄｔｈａｔｔｈｅ ＳＰＨｕｎｉｆｏｒｍｍｅｔｈｏｄｗｉｔｈ ｔｈｅ ｎｅｗ ｃｏｎｔａｃｔ ａｌｇｏｒｉｔｈｍ ｃ