流体与结构相互作用问题的FE-SPH耦合模拟

应用有限元(FE)-光滑粒子流体动力学(SPH)耦合法模拟了具有自由表面的不可压流体与结构的相互作用问题.流体和结构分别采用SPH法和有限元法同时求解,两者在交界面处的相互作用通过接触算法进行处理.为了避免隐式计算压力,通过引入人工压缩率,将不可压流体近似为人工可压缩流体.采用FE-SPH耦合法对弹性板在随时间变化的水压作用下的变形以及倒塌水柱冲击弹性结构两个问题进行了模拟.模拟结果与实验结果以及其他已有数值结果符合良好,说明FE-SPH耦合法用于流体与结构相互作用问题的模拟是可行和有效的.     

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

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

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

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

粘性流体中弹性板振动的有限元耦合问题

流体－结构耦合作用问题是工程中比较常见的问题，具有重要意义，由于流体计算的复杂性，迄今为止，大部分的流体－结构耦合分析都是建立在对流体充分简化的基础上，尤其是将流体视为无粘、无旋的理想流体。该文在近年来前人工作的基础上，发展了一种流体－弹性结构耦合计算模式。流体视为不可压、有粘的介质，流场没有自由表面。该文采用ＳＵ／ＰＧ方法形成流体的有限元方程，采用ＡＬＥ格式处理流体和结构之间的移动界面。采用预估     

非线性流体-刚体结构相互作用问题的一种数值模拟方法

给出了一种模拟非线性流体-刚体结构相互作用问题的数值方法.文中假定结构承受大的刚体运动,流体流动受非线性有粘或无粘的场方程支配并满足自由表面和两相耦合界面上的非线性边界条件,利用任意拉氏-欧氏(ALE)网格系统构造了数值模型.采用所探讨的多块数值格式,允许可动重造网格间有独立的相对运动,从而克服了流体网格与固体大运动匹配的困难.通过数值离散化,导出了描述非线性流固耦合动力学的数值方程并应用耦合迭代过程对其作了求解.通过算例,说明了所提出数值模型的应用.     

双层旋转壳流固耦合振动有限元分析

本文采用旋转壳体元和流体元分析了双层旋转壳体流固耦合振动特性,基于Novozhilov壳体理论、水弹性理论以及Hamilton变分原理推导出耦合系统运动方程,计算结果与实测值比较表明本分析方法有较好的精度。     

海洋地震工程流固耦合问题统一计算框架

海底地震动的模拟以及海洋工程结构的地震反应分析中,涉及到海水、饱和海床、弹性基岩、结构之间的相互耦合.传统的方法分别采用声波方程描述理想流体、Biot方程描述饱和海床、弹性波方程描述基岩和结构,分别考虑相互之间的耦合,十分不便.本文基于理想流体、固体分别为饱和多孔介质的特殊情形(孔隙率分别为1和0),由饱和多孔介质的Biot方程可退化得到理想流体的声波方程和固体的弹性波方程.然后,以饱和多孔介质方程为基础,经集中质量有限元离散,考虑不同孔隙率的饱和多孔介质之间耦合的一般情形,建立了该耦合情形的求解方法.进一步论证了该一般情形的耦合计算方法可分别退化到流体与固体、流体与饱和多孔介质、固体与饱和多孔介质之间的耦合计算,从而将流体、固体、饱和多孔介质间的耦合问题纳入到统一计算框架,并编制了相应的三维并行分析程序.以P-SV波垂直入射时,半无限层状海水-饱和海床、海水-弹性基岩、海水-饱和海床-弹性基岩三种情形的动力分析为例,采用统一计算框架结合透射边界条件进行求解,并与传递矩阵方法得到的解进行对比,验证了该统一计算框架的有效性以及并行计算的可行性.      

密闭腔体声-结构耦合系统的动力灵敏度分析

以密闭空腔为对象,开展了声-结构耦合系统的动力分析和灵敏度计算,为系统性态优化设计提供理论和算法基础。分别把结构和声场进行离散化,推导了声-结构耦合系统的有限元方程,求解了耦合系统的频率和声压级响应。在此基础上,以结构尺寸为设计变量,计算了耦合系统的固有频率和声压级响应的灵敏度,解决了声-结构耦合系统动力灵敏度的数值算法问题。     

声场-结构耦合系统声压约束下板重量优化设计研究

为了减小振动板重量并使其振动噪声满足设计要求,基于声场-结构耦合有限元模型,提出了使结构重量最小化、满足声压约束的优化设计模型。用有限元直接法计算耦合系统的声响应及域点声压对结构设计变量的灵敏度,用可行方向法对其进行了优化设计研究。以一矩形平板和立方体声场耦合系统为实例,以壳厚度为设计变量,给出了域点声压对结构设计变量的灵敏度、优化前和优化后声压级对比图及最优的设计变量值。经过优化设计,在1Hz~200Hz频段内域点声压最大值降低6分贝,结构重量减少1.88千克。结果表明,声压约束下结构重量最小化设计,通过对结构重量的重新分布,能同时满足结构轻量化及噪声指标的要求。     

一种新的并行自适应网格有限元算法

给出了一种新的适用于流体力学问题的并行自适应有限元算法。首先，基于初始稀网格上获得的事后误差估算值，应用反复谱对剖分方法对初网格进行划分，使各子域上总体误差近似相等，从而解决并行自适应计算中的负载平衡问题。然后在各处理器上独立地求解整体问题，并进行指定子域上的网格自适应处理。最后将各子域上的自适应网格组合成一个整体网格，应用基于粘接元技术的区域分裂法在该网格上获得最终解。文末给出了数值实验结果。     

A parallel two-level finite element method for the Navier-Stokes equations

Based on domain decomposition, a parallel two-level finite element method for the stationary Navier-Stokes equations is proposed and analyzed. The basic idea of the method is first to solve the Navier-Stokes equations on a coarse grid, then to solve the resulted residual equations in parallel on a fine grid. This method has low communication complexity. It can be implemented easily. By local a priori error estimate for finite element discretizations, error bounds of the approximate solution are derived. Numerical results are also given to illustrate the high efficiency of the method.     

非线性动力有限元重叠区域分裂的隐式并行算法

针对大规模结构非线性瞬态动力分析非常耗时,提出了相应的并行算法。该算法采用无条件稳定的Ne-wmark-β方法(平均加速技术)进行时间积分,并结合区域分裂技术进行分析。它不同于已有的采用非重叠区域的并行算法,而是采用重叠区域的并行算法。对给定结构有限元分析的质量、阻尼、刚度矩阵进行分裂可推出重叠区域分裂算法的计算公式。为改善每一步的求解,采用预估和校正子方案。编写了该算法的程序,在工作站机群上实现了数值算例,验证了算法的性能。计算结果表明该算法优于非重叠区域分裂算法。     

A parallel solver based on the dual Schur decomposition of general finite element matrices

A parallel solver based on domain decomposition is presented for the solution of large algebraic systems arising in the finite element discretization of mechanical problems. It is hybrid in the sense that it combines a direct factorization of the local subdomain problems with an iterative treatment of the interface system by a parallel GMRES algorithm. An important feature of the proposed solver is the use of a set of Lagrange multipliers to enforce continuity of the finite element unknowns at the interface. A projection step and a preconditioner are proposed to control the conditioning of the interface matrix. The decomposition of the finite element mesh is formulated as a graph partitioning problem. A two-step approach is used where an initial decomposition is optimized by non-deterministic heuristics to increase the quality of the decomposition. Parallel simulations of a Navier–Stokes flow problem carried out on a Convex Exemplar SPP system with 16 processors show that the use of optimized decompositions and the preconditioning step are keys to obtaining high parallel efficiencies. Typical parallel efficiencies range above 80%. © 1998 John Wiley & Sons, Ltd.     

随机有限元分析的二级区域分解并行求解算法

采用蒙特卡罗模拟(MCS)和加权积分法对二维问题进行随机有限元分析。尽管MCS方法对任何有确定解的问题都具有求解精度高的优点,但由于求解所需的计算量巨大使其应用受到限制。利用并行求解技术可有效地处理这种密集型计算问题。基于有限元分裂对接法(FETI)的并行特性并利用预处理共轭梯度法(PCG)的求解高效性,结合整体子区域实现(GSI-PCG)和FETI法,提出二级求解算法,并在工作站机群上实现了数值算例。算例计算结果表明本文GSI(PCG)-FETI算法具有较高的并行加速比和并行效率,具有良好的性能,可有效地进行二维问题的随机有限元分析。     

A parallel finite element sliding mesh technique for the simulation of viscous flows in agitated tanks

A parallel sliding mesh algorithm for the finite element simulation of viscous fluid flows in agitated tanks is presented. Lagrange multipliers are used at the sliding interfaces to enforce the continuity between the fixed and moving subdomains. The novelty of the method consists of the coupled solution of the resulting velocity–pressure‐Lagrange multipliers system of equations by an ILU(0)‐QMR solver. A penalty parameter is introduced for both the interface and the incompressibility constraints to avoid pivoting problems in the ILU(0) algorithm. To handle the convective term, both the Newton–Raphson scheme and the semi‐implicit linearization are tested. A penalty parameter is introduced for both the interface and the incompressibility constraints to avoid the failure of the ILU(0) algorithm due to the lack of pivoting. Furthermore, this approach is versatile enough so that it allows partitioning of sliding and fixed subdomains if parallelization is required. Although the sliding mesh technique is fairly common in CFD, the main advantage of the proposed approach is its low computational cost due to the inexpensive and parallelizable calculations that involve preconditioned sparse iterative solvers. The method is validated for Couette and coaxial stirred tanks. Copyright © 2011 John Wiley & Sons, Ltd.     

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 DOMAIN DECOMPOSITION ALGORITHM WITH FINITE ELEMENT-BOUNDARY ELEMENT COUPLING

A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.     

基于PVM的网络并行有限元面向对象的设计

子结构是有限元并行计算常用的一种方法,本文采用面向对象的方法,首先对子结构进行了面向对象的设计,得到了其类层次结构图;然后针对工作站网络有限元并行计算环境。提出了基于PVM消息传递平台上的Shadow—Mirror数据传输模型,该模型在有限元并行计算数据传输时,充分发挥数据面向对象的特性,采用设置数据缓冲区、短消息合并等方法以缩短数据通信时间,并据此编制了相应的程序。计算结果表明,使用文中提出的面向对象的Shadow—Mirror数据传输模型可以得到较为理想的并行加速比,而且随着问题规模增大,并行加速比增高。本文研究内容为进一步开展基于工作站网络的并行有限元研究提供了一个可参考的基础。     

Parallel finite element algorithm based on full domain partition for stationary Stokes equations

Based on the full domain partition, a parallel finite element algorithm for the stationary Stokes equations is proposed and analyzed. In this algorithm, each subproblem is defined in the entire domain. Majority of the degrees of freedom are associated with the relevant subdomain. Therefore, it can be solved in parallel with other subproblems using an existing sequential solver without extensive recoding. This allows the algorithm to be implemented easily with low communication costs. Numerical results are given showing the high efficiency of the parallel algorithm.