统一格式的显式与隐式任意混合异步算法

动力学问题的有限元分析需要在每一时步求解系统信息,相对于静力学问题,其计算量要大得多.因而,提高计算效率,节省计算工作量是动力学求解方法研究的主要内容.该文针对大型复杂动力学系统的高效求解问题,提出了一种基于Newmark离散格式的显式、隐式任意混合异步算法,根据整体系统不同局部的物理力学特性和求解精度要求,在空间域及时间域内对动力学系统方程进行多尺度求解.该方法根据显式、隐式算法固有的信息传递机制,采取动态的可变边界处理方法,避免了异步边界上的误差积累;并通过对整体系统能量平衡的校验,动态地确定和修正仿真计算时步,可以有效地预防不稳定性的产生和发展.数值算例表明:该算法能在保持较高的计算精度的同时,极大地降低计算资源消耗,因而具有一定的实用价值.      

并行有限元网格生成方法及其应用

大型工程数值仿真中,在前处理阶段需要生成千万甚至亿量级的网格,传统的串行网格生成方法由于内存和时间的限制,难以处理如此规模的网格。针对此问题,本文提出了一种大规模网格并行生成方法。首先基于推进波前法对几何模型进行初始体网格划分,接着利用图论理论进行区域分解,并通过表面单元恢复保持其几何精度,然后通过分裂法进行网格的并行生成。将所述方法应用到实际大型工程数值仿真前处理阶段,结果表明所述方法可以获得较好的并行效率,同时所产生的网格质量可以满足后续计算需要。     

泥水平衡盾构机掘进并行数值模拟方法研究

以某大型越江隧道工程为背景,建立该工程场地土层的大规模三维非线性数值模型。同时基于共享内存式并行计算平台(SMP),利用刚度迁移法,完成了泥水平衡盾构机逐步掘进百米长距离,并穿越大堤这一复杂力学过程的大规模数值模拟。为模拟泥水平衡盾构机掘进各施工阶段,不同施工因素之间的复杂力学作用关系,提出了对各施工因素模拟的简化方法。具体包括:开挖面泥水压力对前方土体的作用、盾构机机头的超挖、盾构机机身对土体的作用、盾尾注浆材料的时空变化特性、盾构机后方百米长拖车对已铺成隧道的作用等施工因素。通过对计算结果与工程试验段实测数据对比验证的基础上,预测了地表沉降在盾构机掘进过程中的变化规律。通过对SMP并行计算平台下盾构法施工数值模型之并行计算效率的对比研究,给出了各并行算法的加速比、计算效率以及接触模型对并行计算效率的影响关系。     

列车与结构动态耦合分析的并行计算方法

在分析结构动态响应时考虑列车与结构的动态耦合作用,采用详细三维有限元模型会带来计算量太大的问题。本文采用并行计算方法,根据列车与结构动态耦合模型的计算特点,设计实现了列车结构耦合均衡的分区算法,并以两个工程应用为例,利用该方法对列车结构三维数值模型进行分区计算,结果表明该分区方法比递归坐标二分法有更好的并行效率。     

面向杆系结构的离散元异步长计算方法

为了提升离散元法处理连续介质问题的计算效率, 提出一种基于重叠颗粒的离散元分区异步长计算方法来处理连续介质问题。该方法采用颗粒分割将连续介质划分为若干子分区, 各分区采用velocity-Verlet积分格式求解运动方程。相邻分区通过重叠颗粒构成局部耦合区域, 边界数据传递过程中不涉及插值和截断过程。数值算例表明, 该方法在保证高精度的同时有效地降低了计算时间。在实际应用中可以有针对性地将连续介质划分为不同尺度颗粒的分区, 根据问题规模及分区颗粒尺度特性采用不同时间步长, 节省存储空间且大幅提升计算效率。     

GDR法在大型结构模态分析中的应用

有限元模态分析中，往往需要减缩矩阵阶数来减小求解规模，目前应用较多的是Guyan静力凝聚法和Kuhar动力凝聚法。本文从理论上分析一种不常用的方法一广义动力化简法(GDR法)，阐述该方法的基本原理、应用过程及其与前两种方法相比的优越性，并通过对城市客车的有限元模态计算加以验证。     

基于响应面模型的薄壁构件耐撞性优化设计

建立了薄壁构件的结构耐撞性优化模型，通过试验设计在设计空间中选择少量样本获得设计的响应特性，建立响应面模型，并应用Pareto遗传算法进行结构耐撞性的优化设计。     

结构静力有限元分层并行计算方法

根据分布式存储并行计算机的体系结构特点,提出了一种结构静力有限元分层并行计算方法. 该方法在两级分区两次缩聚策略的基础上不仅实现了大量数据的分布式存储,提高了数据的内存访问速率;而且实现了计算过程的三层并行,有效提高了通信效率;此外,它还进一步降低了界面方程的规模,大幅度减少了界面方程的求解时间. 因此,它能够充分利用分布式存储并行计算机的体系结构特点提升大规模并行计算效率. 最后通过典型数值算例验证了该方法的正确性和有效性.      

爆炸冲击波对玻璃幕墙破坏作用的多物质ALE有限元模拟

应用LS-DYNA显式程序,基于多物质的ALE有限元法,利用罚函数实现了爆炸冲击波与玻璃幕墙结构间的耦合作用,重现了玻璃幕墙防爆炸实验中爆炸荷载的产生及其与玻璃幕墙相互作用的三维动态过程。研究了爆炸冲击波作用下幕墙结构的动态响应行为,分析了幕墙内外层玻璃的破坏情况。研究表明,数值模拟结果与实验结果一致,这为玻璃幕墙的抗爆设计及改进提供了重要的参考依据。     

冲击载荷下柔性储液罐动态响应数值模拟及规律分析

针对柔性材料大变形、非线性等特点,对柔性储液罐的冲击响应规律进行研究,依据实验结果,采用ALE有限元方法对冲击情况下柔性罐状储液容器及其中流体的动态响应进行了数值模拟.所得结果与实验结果具有较好的一致性.在此基础上,深入分析了空投高度、装水量等对柔性储液罐动态响应的影响规律.重点关注了表面一些关键部位的变形和应力变化情...