【摘 要】
：
开展有限元并行计算可以扩大计算规模、加快计算速度，但是由于开展有限元并行计算需要对并行计算技术有比较深入的了解，从而使得有限元并行计算的应用与推广受到限制。Aztec是一个并行线性解法器，提供多种方程组迭代求解方法和预条件子，可移植可扩展科学计算工具箱PETSc 提供了高性能求解偏微分方程组的大量对象和解法库。基于这两个软件开展有限元并行计算，无需复杂的数据分布和通信，降低了开发难度和成本，缩短了
【机 构】
：
河海大学土木工程学院,南京 210098 中国水利水电科学研究院,北京 100044
论文部分内容阅读
开展有限元并行计算可以扩大计算规模、加快计算速度，但是由于开展有限元并行计算需要对并行计算技术有比较深入的了解，从而使得有限元并行计算的应用与推广受到限制。Aztec是一个并行线性解法器，提供多种方程组迭代求解方法和预条件子，可移植可扩展科学计算工具箱PETSc 提供了高性能求解偏微分方程组的大量对象和解法库。基于这两个软件开展有限元并行计算，无需复杂的数据分布和通信，降低了开发难度和成本，缩短了开发周期。本文首先结合Aztec，开发了支持自由格式和一维变带宽存储格式有限元方程组求解并行接口，能方便的为用户调用；并采用不同算例对各种求解方法进行了测试和分析。然后，给出了基于PETS的结构有限元并行计算实现方法，包括有限元方程组的并行形成和并行求解的实现，并根据PETSc 的特点，提出了提高计算性能的优化措施，即数据局部化和存储预分配，数值试验表明实现方法可行，优化措施效果明显。
其他文献
文章以高孔隙率的烧结纤维多孔金属为研究对象，提出了材料在高声压级环境下吸声性能的优化方法。首先从纤维多孔金属的微观结构出发，分析了在线性声压级条件下材料各基本参数（孔隙率、纤维直径，厚度等）的变化对法向声入射表面阻抗的影响，并给出了如何通过优化各参数来提高烧结纤维多孔金属在线性声压级条件下吸声性能的方法；然后分析在高声压级条件下声压级变化对烧结纤维多孔金属的表面阻抗的影响，提出了进一步提高烧结纤维
本文研究了直杆点阵型夹芯结构和斜杆点阵型夹芯结构的弯曲性能，并对比了这两种结构的力学性能。首先采用一阶剪切变形理论，对夹芯结构的弯曲性能从理论上加以预测。然后实验和数值模拟的方法对理论结果加以验证，对比了三者之间的不同。结果表明，斜杆点阵型夹芯结构的弯曲刚度和强度都明显大于直杆点阵型夹芯结构，通过提高芯子与面板之间结点的强度，可以较大范围内提高夹芯结构的强度。
本文实验测试了孔隙率为0.92,孔密度不同的三种通孔铝泡沫热沉在轴流风扇冲击射流下的整体热阻、表面温差等性能表征参数，并与传统翅片式热沉的相应结果进行了比较分析。结果表明：铝泡沫热沉能使加热模块表面温度的分布明显变得均匀，更有利于延长电子元器件的使用寿命。且铝泡沫热沉可在保持整体换热性能与翅片式热沉相当的条件下，使散热装置的体积和重量均减少50％。另外，实验也显示出铝泡沫热沉在电子元器件冷却应用上
本文使用Daubechies（DB）小波求解二维平面弹性问题。小波尺度函数直接用于函数近似，这种方法不需要网格，利用无网格法的思想，建立了求解二维弹性问题的公式。为处理一般的边界并提高求解效率和精度，提出了一种积分估计方案。二维弹性问题的数值算例显示该方法的稳定和有效性，并有潜力求解更复杂问题。
本文基于薄板的稳定性理论，对挠度函数采用Hermite径向基函数（HRBF ）进行插值，研究无网格Galerkin法在板屈曲问题中的应用HRBF方法利用挠度及其导数进行插值，具有C1连续性，而普通的形函数只具有C0连续性，故HRBF法非常适合对光滑性要求较高的板弯曲问题。而基于Galerkin法的无网格法实现过程与有限元法类似，具有精度高、稳定性好等优点。根据本文构造的无网格法数值计算了矩形薄板在
介绍了一种不需要内部网格计算导热系数随温度变化的热传导问题的边界元算法.该算法是建立在一种能将任何域积分转换成边界积分的径向积分(RIM)转换法基础上,首先用对应各向同性问题的基本解来建立以正规化位移表示的非线形问题的标准积分方程,然后用径向积分转换法将出现在积分方程中的域积分转换成边界积分,从而形成不需要使用内部网格来计算域积分的纯边界元算法.在利用径向积分法进行积分转换的过程中，由于存在与温度
采用径向积分边界单元法对多功能复合涂层结构进行热应力分析。首先使用加权余量法建立非线性热传导和热应力问题的边界域积分方程，然后用径向积分法将出现在积分方程中的域积分转换成等效的边界积分。热学与力学物性参数可以是温度的函数，也可以是空间的函数。采用不等间隔单元子分技术，消除结构多尺度问题中的积分近奇异性；采用三步变量凝聚多区域边界元技术，求解由多种材料组成的复合结构问题。论文给出了数值算例，验证所述
对Navier-Stokes 方程非结构化网格有限体积法离散方法进行了研究，给出了流动边界条件处理方法。将聚合代数多重网格方法成功应用到基于压力的流场计算过程，通过方程聚合方法形成粗化子集，设计了一种高效的聚合算法以及各算子计算方法。从线性方程组解法层面上，利用代数多重网格方法对压力修正方程求解，而对速度以及其它被动标量则采用CGS（Conjugate Gradient Square）方法，两类线
基于建筑玻璃幕墙防爆炸实验，建立了包括炸药、空气、幕墙防爆炸测试系统、地面等多个对象在内的精细三维分析模型。在高性能计算平台上应用LS-DYNA显式程序，采用多物质的ALE 有限元法实现爆炸冲击波与玻璃幕墙结构间的耦合作用，重现了玻璃幕墙防爆炸实验中爆炸荷载的产生及与玻璃幕墙相互作用的三维动态过程。研究了爆炸冲击波作用下建筑幕墙结构的动态响应行为，分析了幕墙玻璃的应力分布情况，得到了与实验一致的结
转向架是轨道车辆非常重要的部件，它承受车辆的重量及各方向的载荷，因此有必要对其进行有限元分析及计算。本文应用Altail公司的Hypermesh8.0中的接触模块来定义轨道车辆转向架轴箱与轴的接触状态，在车轮和车轴的接触处施加全约束，并利用大型有限元分析软件Ansys进行强度校核。结果表明，该方法比以往直接在轴箱套筒内壁施加全约束的计算方法更符合轴箱的实际受力状态，亦可广泛应用于工程中的其他接触问