一种大型结构特征值问题的并行解法

本文提出了一种求解大型结构固有频率与模态的并行解法。该方法在子空间迭代过程中,利用子结构刚度阵和质量阵并行进行凝聚,求得下一次的迭代基矢量,直到收敛。该算法在西安交通大学ELXS1-6400并行机上程序实现,计算结果表明能大幅度节省计算时间,同时也有效地节省了内存。     

结构模态分析的有限元并行模拟

以子结构模态综合分析为基础,提出一种求解大型结构特征值问题的并行解法．采用子结构模态综合算法,结构特征模态采用子空间迭代方式并行求解．这种子空间迭代法的子结构并行计算的实施是利用子结构的刚度阵和质量阵而不必完全组集系统刚度阵和质量阵求解综合系统的特征值问题．数值结果表明这种求解大型结构特征值问题的并行算法是可行有效的．     

大型结构特征值问题的混合粒度并行算法

本文提出一种求解大形结构特征值问题的粗细粒度混合并行算法：在子结构模态综合粗粒度并行算基础上,综合系统的特性值问题采用细粒度并行方式求解。细粒度并行包括子空间迭代法的子结构并行算法、雅可比分块并行计算的方法和一种Ｎｅｗｔｏｎ－Ｒａｐｈｏｎ迭代法在多处理器上任力均衡分配的有效策略。子空间迭代法的子结构并行计算的实施是利用子结构的刚度阵和质量阵而不必完全组集系统刚度阵和国求综合系统的特征值问题。利用雅     

动力子结构模态综合解法的并行算法

本文根据动力子结构模态综合法,提出了相应的并行算法。该算法有效地将整个结构分成独立的多个子结构,然后由多个CPU同时独立求各子结构的分支模态和进行各子结构的分支模态变换。再串行组集界面刚度和质量阵,并求解缩减后的整体方程。最后返回各子结构求结点位移对于(`ω~2`)的模态(φ),这一步也在各CPU内独立地同时进行。该方法在西安交通大学的ELXSI-6400并行机上程序实现,表明能有效地节省计算时间,为一种大型结构动力分析方法。     

大型结构有限元并行分析程序——PARNFAP

本文介绍了作者在西安交通大学ELXSI-6400并行机调试的大型结构有限元分析程序—PARNFAP。该程序利用子结构方法实现有限元计算的并行化,计算结果表明其具有良好的并行性能。     

基于多重多级动力子结构的瞬态分析方法

提出利用多重多级子结构技术与Newmark算法求解结构动力学方程的高精度算法.该算法利用静凝聚技术列式简单,在凝聚过程中并不引入任何近似的优点,采用子结构周游树技术,分别对每个子结构求解Newmark等效平衡方程,最后通过回代求解得到整体结构的响应.由于该算法考虑了子结构内部自由度对整体求解的贡献,算法实施不受子结构划分方式的限制,因此可以得到系统高阶模态对响应分析的影响.该算法计算精度与传统的全结构求解相当,计算效率高,消耗计算机资源少,且可构造为统一的多重多级子结构综合分析算法框架.数值算例验证了该算法的正确性和有效性.     

并行子结构法求解结构动力响应

本文利用直接积分法结合子结构技术和并行机特点提出一个求解复杂系统动响应的并行子结构动响应法。该方法利用子结构连接条件找到内力传递向量，从而确定各子结构动响应，算例表明该方法是有效的，并且具有较高的并行效率。     

追踪再正交Lanczos算法及其在SAP5程序中的实现

本文提出了适于求解大型复杂结构部分固有频率及其相应振型的追踪再正交Lanczos算法。该算法采用追踪再正交判定准则和截断Lanczos收敛性判定准则,从而减少了运算次数,提高了求解效率。数值结果表明,该算法与子空间迭代法相比收敛快、计算量小,是求解大型稀疏实对称矩阵特征值问题的有效算法。     

一种求解局部非线性结构稳态响应的方法

为了降低求解局部非线性结构稳态响应的计算量,基于子结构和阻抗缩聚提出了一种用于求解局部非线性结构稳态响应的计算方法.将局部非线性结构分解为线性子结构和非线性子结构,利用谐波平衡构造各个子结构的阻抗方程,对线性子结构进行缩聚,将局部非线性动力学方程转化为求解一组非线性代数方程组问题,通过迭代求解非线性代数方程组,求解系统的稳态响应.     

隐式重启动Arnoldi/Lanczos法的子区域并行算法

针对求解有限元分析的特征值问题,提出了一种隐式重启动Arnoldi/Lanczos方法的子区域并行算法。隐式重启动Arnoldi/Lanczos利用重启动技术以提高所需谱的收敛性,并能有效处理Krylov基形成问题、存储所需的内存问题、计算成本问题。并行算法中采取子区域接子区域方法、重叠和非重叠网格划分技术。采用压缩数据结构来储存系数矩阵。对Krylov的数值线性代数运算和隐式重启动法中的数值线性代数运算的并行化进行了研究。数值算例表明:该算法具有良好的适用性和效率,适合分布式储存体系的机群。     

An iterative parallel algorithm of finite element method

In this paper, a parallel algorithm with iterative form for solving finite element equation is presented. Based on the iterative solution of linear algebra equations, the parallel computational steps are introduced in this method. Also by using the weighted residual method and choosing the appropriate weighting functions, the finite element basic form of parallel algorithm is deduced. The program of this algorithm has been realized on the ELXSI-6400 parallel computer of Xi'an Jiaotong University. The computational results show the operational speed will be raised and the CPU time will be cut down effectively. So this method is one kind of effective parallel algorithm for solving the finite element equations of large-scale structures.     

A finite element parallel algorithm of the parametric variational principle for elastic contact problems

Basetl on the finite element solution of the parametric varialional principle of elastic con/del problem, a corresponding parallel algorithm has been created bv utilizing the specialities of parallel computer and the architecture of concurrent processing in this paper. In this algorithm. the parallelisms have heen realized in the processes of creation and assembly of stiffness matrix, of the static condensation, of the solution of stresses and in many other aspects. The programme of this algorithm has been realized on ELXSI-6400 parallel computer of Xi'an Jiaotong University. The results of computation show that the computational time can be saved efficiently and it is an effective parallel algorithm for the analyses of contact problems.     

Sensitivity analysis based on Lanczos algorithm in structural dynamics

IntroductionTheproblemforcomputingeigenpairderivativeshasoccupiedmanyresearchersinthepastseveraldecades.SurveypaperbyAdelmanandHaftka[1]isanexcellentintroductiontothetopicandalsoindicatevariousdisciplinesofapplications.Thereasonwhysomanypeopleareinterestedinthisproblemisthatsensitivitiesplayaveryimportantroleinstructuraldynamics,optimalanalysisandcontrolsystemdesign ,andmayalsobeanexpensiveandtime_consumingtask.Researchonsensitivityofmultipleeigenvalueshasbeenanareaofsignificantinterestinrece…     

大型非比例阻尼线性系统的地震反应复振型分析方法

对于大型非比例阻尼线性系统,当采用基于复振型的振型叠加方法进行动力反应分析时,按照通用算法(例如雅克比法)求解结构全部振型向量的计算工作量很大,甚至是不现实的.本文将经典阻尼系统中行之有效的Lanczos法和子空同迭代法加以推广和改进,给出了一组可对原来的方程进行自由度缩减的实向量基,然后与Foss变换相结合,得到一组实用的复向量基,使之适用于求解复杂非比例阻尼线性系统的任意低阶复振型和相应的复特征值,适用于任意扰力作用下的动力反应分析.理论推导和实例计算表明,本文所给出的复向量基概念清晰,计算效率高,能够适应对具有非比例阻尼特性的大型复杂结构进行动力分析的实际需要,其中包括地震作用下的时程分析和反应谱振型叠加分析.     

岩质圆形隧洞围岩应力场弹塑性新解

针对动态接触问题的有限元并行计算，提出了一种新的接触算法. 新算法引入局部拉氏 乘子技术来计算接触力. 由于同时考虑了无穿透的接触约束条件和相邻接触对的相互影响， 较之广泛使用的罚参数法，新算法使接触约束条件和系统平衡方程得到更充分的满足. 虽然 为提高接触计算精度而在局部采用了迭代技术，但算法仍然具有较高的效率，且与显式时间 积分方案完全相容. 此外，通过构造专门的区域分解方案，实现了将现有为串行程序开发的 搜索算法平滑移植到并行环境的目标. 数值算例表明，所提出的接触算法具有很好的并行性， 在保证了接触问题并行计算精度的同时，取得了满意的并行效率.     

基于SIMPLER算法的多重网格方法研究

以二维方腔顶盖驱动流为模型,将多重网格方法和SIMPLER算法进行耦合,对不同雷诺数下多重网格加速SIMPLER算法和SIMPLER算法的计算效率进行了对比,数值计算表明:多重网格加速SIMPLER算法不仅能够解决SIMPLER算法不能准确模拟较高雷诺数流场的问题,而且其计算效率远远高于SIMPLER算法.本文也对松弛因子的选取、多重网格实现形式以及网格层数对多重网格加速SIMPLER算法的影响进行了研究,从而为多重网格加速SIMPLER算法的实施提供了计算技术.     

Optimal control applied to water flow using second order adjoint method

This paper presents an optimal control applied to water flow using the first and second order adjoint equations. The gradient of the performance function with respect to control variables is analytically obtained by the first order adjoint equation. It is not necessary to compute the Hessian matrix directly using the second order adjoint equation. Two numerical studies have been performed to show the adaptability of the present method. The performance of the second order adjoint method is compared with that of the weighted gradient method, Broyden–Fletcher–Goldfarb–Shanno method and Lanczos method. The precise forms of the adjoint equations and the gradient to use for the minimisation algorithm are derived. The computation by the Lanczos method is shown as superior to those of the other methods discussed in this paper. The message passing interface library is used for the communication of parallel computing.     

基于并行计算和遗传算法的钢-UHPC华夫板组合梁优化设计

为实现钢-超高性能混凝土(UHPC)华夫板组合梁结构快速经济合理的设计，提出了基于并行计算与遗传算法的结构优化设计方法。通过Python建立了并行计算平台，使Abaqus和Python能够执行同步数值模拟和数据处理，以成本最小化为目标，采用遗传算法对钢-UHPC华夫板组合梁进行了优化，验证了所提方法的可行性。结果表明，遗传算法中密集的分析任务可以并行化并分配给不同的计算资源以提高计算效率；使用并行计算可以提高8.6倍的优化效率；并行计算和串行计算的CPU平均使用率分别为82%和18%。本文方法的成功应用可为其他类型结构的优化设计提供参考。     

基于Lanczos算法的模态重分析法及其在车身结构设计中的应用

模态重分析是指在结构修改之后不需要重新求解广义特征值方程,仅需要根据初始计算结果对修改后的问题进行求解,并能够在保证精度的前提下,提高计算速度。随着结构复杂度和修正量的增加,传统重分析方法的求解精度和稳定性随之下降。为此,利用初始结构模态分析结果,结合Lanczos算法和投影技术,采用缩减基方法求解修改结构的特征值和特征向量,使其同时具备了Lanczos向量快速收敛的优点和基于全局近似的缩减基向量的高精度。为了验证该方法的性能和准确性,对本文方法基于扩展基向量和瑞利-里兹分析的模态重分析法以及改进的单步摄动瑞利商逆迭代法进行了测试。测试结果表明,该方法具有最高的计算精度。同时,将该方法成功用于车架和车门的前期设计中,计算结果表明,该方法具备处理计算规模大、拓扑修改变化量大的结构分析问题的潜力。