储液容器跌落碰撞的并行仿真计算研究
|
摘要
基于多物质的ALE有限元法,应用通用动力学软件LS-DYNA完成了空投柔性储液容器跌落碰撞过程的数值仿真计算.容器采用Lagrangian描述,溶液采用ALE描述,容器与溶液间的耦合作用及容器与地面的碰撞接触均采用罚函数法得以实现.在分析空投容器力学行为的同时,对该问题计算的并行性进行了讨论.针对储液容器跌落碰撞过程中,即存在十分耗时的非线性流体与结构耦合计算又具有碰撞接触计算的特点,在递归坐标二分分区方法的基础上,提出了负载均衡的分区策略,通过程序实现该分区过程.结果表明,相对于坐标递归二分分区方法,负载均衡的分区方法使得各分区的负载更加均衡,具有更好的加速比与并行效率.
The simulation of the three-dimensional flexible fluid-filled container in the crash was achieved through multi-material Arbitrary Lagrangian-Eulerian(ALE) finite element method because of its ability to control mesh geometry independently from geometry.The ALE description was adopted for the fluid domain,while for the structural domain the Lagrangian formulation was adopted.The computation of the fluid-structure interaction and the crash contact between the container and the ground was realized by the penalty-based coupling method.Then the dynamic behavior of the container in the crash was analyzed and the parallelism was discussed because the computation of the fluid-structure interaction and the crash contact was quite time-consuming.Based on the recursive coordinate bisection(RCB),the Load balanced bisection(LBB) was designed and implemented according to the time-consuming characteristics of the fluid-filled container in the crash.The result indicated,in comparison with RCB method,the LBB method had improved the speedup and the parallel efficiency.
The simulation of the three-dimensional flexible fluid-filled container in the crash was achieved through multi-material Arbitrary Lagrangian-Eulerian(ALE) finite element method because of its ability to control mesh geometry independently from geometry.The ALE description was adopted for the fluid domain,while for the structural domain the Lagrangian formulation was adopted.The computation of the fluid-structure interaction and the crash contact between the container and the ground was realized by the penalty-based coupling method.Then the dynamic behavior of the container in the crash was analyzed and the parallelism was discussed because the computation of the fluid-structure interaction and the crash contact was quite time-consuming.Based on the recursive coordinate bisection(RCB),the Load balanced bisection(LBB) was designed and implemented according to the time-consuming characteristics of the fluid-filled container in the crash.The result indicated,in comparison with RCB method,the LBB method had improved the speedup and the parallel efficiency.
引文
[1]Belytschko T,Kennedy J M,Schoeberle D F.Quasi-eulerian finite element formulation for fluid-structure interaction[J].Journal of Pressure Vessel Technology,1980,102(4):62-69
[2]Nitikitpaiboon C,Bathe K J.An arbitrary Lagrangian-eulerian velocity potential formulation for fluid-structure[J].Computers and Structures,1993,47(4):871-891
[3]Marco A,Luigi ML,Maurizio T.Fluid-structure interaction of water filled tanks during the impact with the ground[J].International Journal of Impact Engineering,2005,31:235-254
[4]Johan Z,Mathur K K,Johnson S L,et al.An efficient communications strategy for finite element methods on theConnection Machine CM-5 system[J].Computer Methods in Applied Mechanics and Engineering,1994,113:363-387
[5]Ray S E,Wren G P,Tezduyar T E.Parallel implementations of a finite element formulation for fluid-structureinteractions in interior flows[J].Parallel Computing,1997,23:1279-1292
[6]Lama H.Review of domain decomposition for the implementation of FEM on massively parallel computers[J].IEEEAntennas and Propagation Magazine,1995,37(1):93-98
[7]岳宝增,刘延柱,王照琳.求解液体大幅晃动问题的数值方法评述[J].上海交通大学学报,1999,33(6):760-763Yue Baozeng,Liu Yanzhu,Wang Zhaolin.Reviewof the numerical methods for simulating large-scale amplitude sloshingof liquid[J].Chinese Journal of Shanghai Jiao tong University,1999,33(6):760-763
[8]Wang X W,Guo L,Ge W,et al.Parallel implementation of macro-scale pseudo-particle simulation for particle-filledsystems[J].Computers and Chemical Engineering,2005,29:1543-1553
[9]郭毅之,金先龙,定峻宏,等.并行数值仿真技术在盾构隧道地震相应分析中的应用[J].应用基础与工程科学学报,2005,13(1):43-50Guo Yizhi,Jin Xianlong,Ding Junhong,et al.Parallel numerical simulation techniques in seismic response analysis ofshield tunnel[J].Journal of Basic Science and Engineering,2005,13(1):43-50
[2]Nitikitpaiboon C,Bathe K J.An arbitrary Lagrangian-eulerian velocity potential formulation for fluid-structure[J].Computers and Structures,1993,47(4):871-891
[3]Marco A,Luigi ML,Maurizio T.Fluid-structure interaction of water filled tanks during the impact with the ground[J].International Journal of Impact Engineering,2005,31:235-254
[4]Johan Z,Mathur K K,Johnson S L,et al.An efficient communications strategy for finite element methods on theConnection Machine CM-5 system[J].Computer Methods in Applied Mechanics and Engineering,1994,113:363-387
[5]Ray S E,Wren G P,Tezduyar T E.Parallel implementations of a finite element formulation for fluid-structureinteractions in interior flows[J].Parallel Computing,1997,23:1279-1292
[6]Lama H.Review of domain decomposition for the implementation of FEM on massively parallel computers[J].IEEEAntennas and Propagation Magazine,1995,37(1):93-98
[7]岳宝增,刘延柱,王照琳.求解液体大幅晃动问题的数值方法评述[J].上海交通大学学报,1999,33(6):760-763Yue Baozeng,Liu Yanzhu,Wang Zhaolin.Reviewof the numerical methods for simulating large-scale amplitude sloshingof liquid[J].Chinese Journal of Shanghai Jiao tong University,1999,33(6):760-763
[8]Wang X W,Guo L,Ge W,et al.Parallel implementation of macro-scale pseudo-particle simulation for particle-filledsystems[J].Computers and Chemical Engineering,2005,29:1543-1553
[9]郭毅之,金先龙,定峻宏,等.并行数值仿真技术在盾构隧道地震相应分析中的应用[J].应用基础与工程科学学报,2005,13(1):43-50Guo Yizhi,Jin Xianlong,Ding Junhong,et al.Parallel numerical simulation techniques in seismic response analysis ofshield tunnel[J].Journal of Basic Science and Engineering,2005,13(1):43-50
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心