多波前并行处理的弹塑性子结构并行有限元

研究一种基于ＰＶＭ的弹塑性子结构并行有限元法。利用多个波前,对各子结构并行地进行静凝聚。再采用预条件共轭梯度法（ＰＣＧ）并行求解界面方程。算例表明该方法能获得较好的并行加速比,同时也能有癣地节省内存量。     

有限元方法在爆炸动力学中的应用

从动量和能量守恒方程出发，结合爆炸动力学的特点，建立适于模拟爆炸动力学过程的有限元方法，并对两个问题进行了模拟计算。计算结果与实验结果符合表明，通过选择合适的物质模型，有限元方法对爆炸动力学问题的模拟是可行的。     

球对称形式物质点法及其在水下爆炸中的应用

利用物质点法求解三维远场水下爆炸冲击响应问题时,爆炸冲击波到达结构物之前传播过程的计算会消耗大量的计算资源,影响求解效率。针对这一问题,并考虑到水下爆炸冲击波在自由场中的传播特点,提出了一维球对称形式的物质点法,对球形炸药水下爆炸过程进行了数值模拟,并与Cole公式以及DYNA计算结果进行了比较,结果吻合较好,验证了一维球对称形式物质点法的正确性。在此基础上,提出了基于物质点法的重映射算法,将爆炸冲击波在自由场中的传播过程在一维球对称模型中进行计算,并将计算结果映射到三维模型中,以便继续完成计算过程,有效避免了计算资源在冲击波传播过程计算上的消耗。     

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

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

有限元程序的一种并行化算法

传统的以串行计算机为物质基础的有限元算法分析未能满足科学研究和工程技术发展的需要．并行处理技术的发展为解决此类问题提供了新的契机．机群式并行计算机是当今较为流行的一种并行处理方式．根据有限元法的特点，提出了一种并行有限元算法，并在多处理机分布式内存并行计算机环境下完成了面向对象并行有限元程序的编制．算例表明该算法能大大提高其分析速度，而且算法简单，实现方便．     

水下目标在水下爆炸作用下冲击响应的时频特征

基于某水下目标的抗水下爆炸实验数据,利用小波变换良好的时频局部化性质对监测的水下目标内部装置冲击加速度信号进行时频特征分析,得到了不同频带上冲击信号振动分量的时间历程曲线和能量分布。结果表明,基于冲击信号小波时频特征分析的频带能量可以同时反映冲击振动的强度、频率和持续时间对目标损伤的影响,同时结合不同频带的时间历程曲线可以获得冲击波和二次压力波在不同频带上的分布和衰减的细节信息,由此可确定冲击波和二次压力波各自对目标毁伤的影响。     

装药量及水深对水下爆炸气泡动态特性的影响

以气泡体积加速度模型为基础研究水下爆炸气泡运动的初始条件,采用MSC.DYTRAN 非线性 有限元软件,结合开发的定义流场初始条件与边界条件的子程序,研究水下爆炸气泡运动特性,包括气泡的脉 动、坍塌以及射流等运动特性,并将气泡脉动体积计算结果与实验及边界积分方法计算结果进行对比,验证了 有限元模型的正确性与有效性。以此为基础,得到初始水深、装药量与气泡的脉动体积、最大半径、周期以及 射流速度之间的关系,计算结果与经验公式具有较好的一致性。得到一些有规律性的曲线,可为相关水下爆 炸气泡动态特性研究提供参考。     

铝纤维对黑索今水下爆炸性能的影响

将铝纤维添加到RDX中得到新型非理想炸药,并与RDX进行水下爆炸对比实验,得到2种炸药在不同位置的压力时程曲线,经过分析计算得到两者水下爆炸的冲击波压力峰值、冲量、冲击波能、气泡脉动周期以及气泡能。结果表明:距离药柱相同位置处,铝纤维炸药的压力峰值低于RDX,铝纤维炸药的冲击波冲量高于RDX,其差值受离药柱的距离影响不大。与RDX相比,铝纤维炸药的比冲击波能降低了2%~5.2%,比气泡能提高了9.4%~23.36%,总能量平均提高了3.5%。铝纤维炸药比气泡能与总能量的比值为55%~60%,高于RDX的50%~53%,其总能量与爆热比值为74%~84%,低于RDX的89%~95%。     

离心机水下爆炸缩比实验原理及数值研究

对离心机水下爆炸缩比实验方法进行探究,对离心机水下爆炸缩比实验的相似理论进行了推导,通过数值计算分析,探究了原模型、离心机缩比实验及常规缩比实验的冲击波载荷、气泡载荷以及气泡动力学行为。结果表明:常规缩比实验不能对气泡行为及垂直方向的近场载荷进行准确的预报,若要保证远场气泡脉冲峰值误差小于10%,则爆距需大于9.5倍气泡最大半径。而离心机缩比实验能够对原模型进行准确的预报,以小当量装药模拟大当量装药水下爆炸整个物理过程,且冲击波和气泡两个阶段完全相似。同时,水深也可以进行几何缩比,克服了常规缩比方法的缺陷。     

低能量导爆索水下爆炸冲击波特性实验研究

为了获得低能量导爆索水下爆炸冲击波特性,在水下爆炸试验容器中对复合铝纤维爆炸索进行水下爆炸实验。采用PCB压电型传感器测量水中爆炸冲击波压力脉冲,获得了复合铝纤维爆炸索冲击波压力峰值,计算了其冲击波比能和高压区冲击波比能。研究结果表明:复合铝纤维爆炸索冲击波能量的输出主要集中在特征时间内的高压区,占到冲击波比能的84%。通过复合铝纤维爆炸索水下爆炸冲击波特性的研究,为其下一步的推广应用奠定了基础。     

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.     

Prolonged simulation of near-free surface underwater explosion based on Eulerian finite element method

In the area of naval architecture and ocean engineering, the research about the underwater explosion problem is of great significance. To achieve prolonged simulation of near-free surface underwater explosion, the underwater explosion transient numerical model is established in this paper based on compressible Eulerian finite element method(EFEM). Compared with Geers–Hunter formula, EFEM is availably validated by simulating the free-field underwater explosion case. Then, the bubble pulsation and flow field dynamic characteristics of the cases with different underwater explosive depth are compared in this work. Lastly, the height of the water hump and the pressure of flow flied are analyzed quantitatively through the simulation results.     

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

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

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.     

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算法具有较高的并行加速比和并行效率,具有良好的性能,可有效地进行二维问题的随机有限元分析。     

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

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

近场水下爆炸气泡与双层破口结构的相互作用

针对双层结构在水中受到水下爆炸冲击这一问题,利用欧拉有限元数值模型对近场水下爆炸气泡与双层破口结构的相互作用机理进行了研究,分析了舱室涌流及流场演化等规律。首先,通过放电实验对数值模型进行了验证,结果表明,数值结果和实验结果吻合较好;然后,总结了不同破口尺寸、不同起爆位置和不同壳间水位条件下的耦合作用规律。在内部空气、流体惯性以及破口的联合作用下,气泡演化过程中会出现气泡分割现象。当内层破口尺寸系数小于0.5时,内舱室内会出现二次涌流现象,且涌流形态较细长;炸药起爆位置系数小于0.1时,自由液面处会出现破碎和重闭合现象;壳内水位对舱室涌流量的影响作用较为复杂,当水位满舱时,急速涌流会减少船艇的应急时间。

     

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.