求解固体力学问题的强?弱耦合形式单元微分法

单元微分法是一种新型强形式有限单元法. 与弱形式算法相比, 该算法直接对控制方程进行离散, 不需要用到数值积分. 因此该算法有较简单的形式, 并且其在计算系数矩阵时具有极高的效率. 但作为一种强形式算法, 单元微分法往往需要较多网格或者更高阶单元才能达到满意的计算精度. 与此同时, 对于一些包含奇异点的模型, 如在多材料界面、间断边界条件、裂纹尖端等处, 传统单元微分法往往得不到较精确的计算结果. 为了克服这些缺点, 本文提出了将伽辽金有限元法与单元微分法相结合的强?弱耦合算法, 即整体模型采用单元微分法的同时, 在奇异点附近或某些关键部件采用有限元法. 该策略在保留单元微分法高效率与简洁形式等优点的同时, 确保了求解奇异问题的精度. 在处理大规模问题时, 针对关键部件采用有限元法, 其他部件采用单元微分法, 可以在得到较精确结果的同时, 极大提高整体计算效率. 在本文中, 给出了两个典型算例, 一个是具有切口的二维问题, 一个是复杂的三维发动机问题. 针对这两个问题, 分析了该耦合算法在求二维奇异问题和三维大规模问题时的精度与效率.      

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.     

NEW ALGORITHM OF COUPLING ELEMENT-FREE GALERKIN WITH FINITE ELEMENT METHOD

IntroductionMeshless methods, as a special numerical method, originated from1970s. Since thediffuse element method was proposed by Nayroleset al.[1]in1992, the meshless methodshave received wide attentions in the mechanics area, and have shown obvious adv…     

自适应无网格法在生物涂层接触问题中的应用

自适应无网格法是针对有限元法无法求解或者不易求解的复杂问题,利用无网格法节点排布灵活、易于增删节点、便于自适应分析等优点发展起来的. 在对自适应无网格法理论基础和发展进行总结基础上,采用基于应变能梯度的自适应无网格—— 有限元耦合方法,对等离子喷涂制备的HA 生物涂层材料的无摩擦接触问题进行分析,对制备的两种不同厚度的生物涂层材料进行求解,分别给出了von Mises 应力分布云图. 结果表明,自适应无网格法能较好地应用于生物涂层接触问题中.     

A peridynamic implementation of crystal plasticity

This paper presents the first application of peridynamics theory for crystal plasticity simulations. A state-based theory of peridynamics is used (Silling et al., 2007) where the forces in the bonds between particles are computed from stress tensors obtained from crystal plasticity. The stress tensor at a particle, in turn, is computed from strains calculated by tracking the motion of surrounding particles. We have developed a quasi-static implementation of the peridynamics theory. The code employs an implicit iterative solution procedure similar to a non-linear finite element implementation. Peridynamics results are compared with crystal plasticity finite element (CPFE) analysis for the problem of plane strain compression of a planar polycrystal. The stress, strain field distribution and the texture formation predicted by CPFE and peridynamics were found to compare well. One particular feature of peridynamics is its ability to model fine shear bands that occur naturally in deforming polycrystalline aggregates. Peridynamics simulations are used to study the origin and evolution of these shear bands as a function of strain and slip geometry.     

弹性半空间二维出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性半空间自由波场时域计算的一维化有限元方法。首先利用Snell定律确定平面波沿水平方向的传播规律,在用有限元法对弹性半空间进行离散化时,竖向单元尺寸根据波动有限元模拟精度要求确定,而水平向有限元网格尺寸根据水平向波的传播规律和采用的离散时间步长确定,使得有限元离散模型中任意节点的运动可以用水平向相邻节点的运动表示,从而将二维有限元节点运动方程组化为一维的形式。求解此一维方程组,可得到弹性半空间中一列节点的运动,再根据行波的传播规律,可确定全空间自由波场。理论分析和数值算例表明,该方法具有较高的精度和良好的稳定性。     

近场动力学与有限元方法耦合求解热传导问题

求解含裂纹等不连续问题一直是计算力学的重点研究课题之一,以偏微分方程为基础的连续介质力学方法处理不连续问题时面临很大的困难. 近场动力学方法是一种基于积分方程的非局部理论,在处理不连续问题时有很大的优越性. 本文提出了求解含裂纹热传导问题的一种新的近场动力学与有限元法的耦合方法. 结合近场动力学方法处理不连续问题的优势以及有限元方法计算效率高的优势,将求解区域划分为两个区域,近场动力学区域和有限元区域. 包含裂纹的区域采用近场动力学方法建模,其他区域采用有限元方法建模. 本文提出的耦合方案实施简单方便,近场动力学区域与有限元区域之间不需要设置重叠区域. 耦合方法通过近场动力学粒子与其域内所有粒子(包括近场动力学粒子和有限元节点)以非局部方式连接,有限元节点与其周围的所有粒子以有限元方式相互作用. 将有限元热传导矩阵和近场动力学粒子相互作用矩阵写入同一整体热传导矩阵中,并采用Guyan缩聚法进一步减小计算量. 分别采用连续介质力学方法和近场动力学方法对一维以及二维温度场算例进行模拟,结果表明,本文的耦合方法具有较高的计算精度和计算效率. 该耦合方案可以进一步拓展到热力耦合条件下含裂纹材料和结构的裂纹扩展问题.      

A coupled IEM/FEM approach for solving elastic problems with multiple cracks

In this paper, the detailed two-dimensional infinite element method (IEM) formulation with infinite element (IE)–finite element (FE) coupling scheme for investigating mode I stress intensity factor in elastic problems with imbedded geometric singularities (e.g. cracks) is presented. The IE–FE coupling algorithm is also successfully extended to solve multiple crack problems. In this method, the domain of the primary problem is subdivided into two sub-domains modeled separately using the IEM for the multiple crack sub-domain, and the FEM for the uncracked sub-domain. In the IE sub-domain, the similarity partition concept together with the IEM formulation are employed to automatically generate a large number of infinitesimal elements, layer by layer, around the tip of each crack. All degrees of freedom related to the IE sub-domain, except for those associated with the coupling interface, are condensed and transformed to form a finite master IE for each crack with master nodes on sub-domain boundary only. All of the stiffness matrices constructed in the IE sub-domains are assembled into the system stiffness matrix for the FE sub-domain. The resultant FE solution with a symmetrical stiffness matrix, having the singularity effect of imbedded cracks in IEs, is required only for solving multiple crack problems.Using these efficient numerical techniques a very fine mesh pattern can be established around each crack tip without increasing the degree of freedom of the global FEM solution. One is easily allowed to conduct parametric analyses for various crack sizes without changing the FE mesh. Numerical examples are presented to show the performance of the proposed method and compared with the corresponding known results where available.     

An enhanced coupling of PD grids to FE meshes

The paper presents a new and simple way to couple FEM meshes to Peridynamic grids with different grid sizes. The excellent performance of the coupling technique is illustrated by means of 2D examples involving rigid body motions, uniform strain conditions and a case of crack branching. The method paves the way to an increased use of peridynamics within FEM software.     

近场动力学最小二乘和有限元耦合方法研究

准确高效地对损伤和断裂问题进行建模是计算力学中的关键研究课题之一。将近场动力学最小二乘在处理含裂纹等非连续问题上的优势和有限元计算效率高及便于施加边界条件的优势结合,提出了近场动力学最小二乘和有限元耦合方法。将裂纹及其可能扩展区域划分为近场动力学区域,边界及其他区域划分为有限元区域,并将其中的结点类型分为近场动力学结点和有限元结点。有限元结点仅与同单元中的其他结点产生作用,近场动力学结点则与其族内的所有结点产生作用。将以上的单元刚度矩阵和质量矩阵进行组装得到整体刚度矩阵和整体质量矩阵。本文的耦合方法数值实现简单有效,相对于键基和常规态基近场动力学,该耦合方法包含了应力和应变的概念,同时不受零能模式的影响。一维和二维静态和动态问题的研究,验证了本文的耦合方法的有效性和准确性。     

A 1D time-domain method for in-plane wave motions in a layered half-space

A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability. The project supported by the National Natural Science Foundation of China (50478014), the National 973 Program (2007CB714200) and the Beijing Natural Science Foundation (8061003). The English text was polished by Yunming Chen.     

Improved Coupling of Finite Shell Elements and 3D Boundary Elements

A strategy for the mixed-dimensional coupling of finite shell elements and 3D boundary elements is presented. The stiffness formulation for the boundary element domain is generated by the 3D symmetric Galerkin boundary element method and thus can be assembled to the global finite element formulation. Based on the equality of work at the coupling interface, coupling equations in an integral sense are derived for curved coupling interfaces and formulated as multipoint constraints in terms of kinematic quantities. Several examples show the highly accurate results compared to a strict kinematic coupling technique.     

功能梯度材料热传导问题的近场动力学模型

当材料中存在不连续性缺陷的时候,传统的基于连续介质理论的方法在研究热传导问题时存在诸多不便．我们基于近场动力学方法(Peridynamics)建立了功能梯度材料的热传导模型,并且研究了在温度荷载作用下功能梯度材料的温度场的变化表现．该文概述了PD方法应用于热传导问题时的详细理论基础,描述了其建模思路以及计算体系,给出了使用PD方法模拟结构承受温度荷载时的计算格式,讨论了不同梯度形式对功能梯度材料内热传导的影响．算例结果表明：通过PD模型所得到的温度场与解析解吻合较好,证明了PD方法在分析功能梯度材料热传导行为等问题时的可行性．     

基于绝对节点坐标法的平面梁有限变形下变形重构

现有的对有限变形条件下柔性结构变形重构的研究往往单纯基于曲率与应变间的几何关系, 同时忽略了被测体的纵向变形及其与弯曲变形的耦合效应. 为得到一种更加精确且能借助现有的力学工具进行应用方向扩展的变形重构方法, 以平面梁为对象, 借鉴变形重构逆有限元法的思想, 将平面梁的变形重构问题视作一类最优化问题. 首先, 通过引入绝对节点坐标法(absolute nodal coordinate formulation, ANCF)对柔性结构大变形下非线性的平面梁应变?位移关系进行精确描述, 构造了一种逆梯度缩减ANCF平面索梁单元. 然后, 对此逆ANCF单元进行改进, 在简化节点自由度的同时通过引入罚函数确保单元节点处的曲率连续性, 既保证了本问题的适定性, 也提升了最终解的精确性. 最后, 基于该单元利用Newton法构造了平面梁有限变形下变形重构问题的两种求解算法, 即逐单元算法和多单元整体算法, 以实现不同需求下的稳定求解. 数值仿真结果表明, 本方法在大变形条件下的变形重构误差小于1%, 而且在测点较少的情况下依然保持较高的精度, 同时验证了本方法的收敛性与计算效率.      

带源参数的二维热传导反问题的无网格方法

利用无网格有限点法求解带源参数的二维热传导反问题，推导了相应的离散方程. 与 其它基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不 需要划分网格，用配点法离散控制方程，可以直接施加边界条件，不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.     

近场动力学与有限元的混合建模方法

综述了近场动力学与有限元混合建模方法的研究进展,阐明了各种混合建模方法的基本原理与特点,并重点介绍本课题组在近场动力学与有限元方法混合建模方面的研究工作。现有近场动力学与有限元混合建模方法包括位移协调约束、力耦合、混合函数方法以及子模型方法等,除子模型方法外,都可归结为并行式多尺度分析方法,其基本思想是将计算结构划分为近场动力学子域、有限元子域以及两者的交界区域(或重叠区域、或界面单元、或过渡区域)。子模型方法可归结为显-显分析方法,先采用显式有限元进行整体分析,后采用近场动力学方法对重点区域进行分析。混合建模方法需要着重提高交界区域的计算精度,并且消除虚假力和虚假应力波问题。提出了通过力耦合的近场动力学与有限元混合建模的隐式分析方法,该方法不再设置重叠区,通过杆单元连接近场动力学子域与有限元子域,其中界面上的有限元结点不仅与其所在单元的其他结点发生作用,还通过杆单元与以其为圆心、一定半径的圆域内的其他物质点相互作用。研究表明,本文提出的混合模型和求解方法既能有效解决裂纹扩展等不连续问题,又可提高计算效率,为工程结构破坏问题的计算分析提供一种有效方法。     

Nodeless variable finite element method for heat transfer analysis by means of flux-based formulation and mesh adaptation

Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method. The solution accuracy is further improved by implementing an adaptive meshing technique to generate finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems are: (a) a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and (b) a transient heat conduction analysis in a long plate subjected to a moving heat source. The English text was polished by Yunming Chen.     

Coupling of finite volume and finite element subdomains using different time integrators

This paper proposes a coupling strategy that can be used for transient fluid–structure interaction. The objective of this paper is to propose a time integrator coupling strategy, which ensures good properties to couple typical solid and typical fluid time integrators in linear cases. It is evaluated on a 1‐D toy problem only dedicated to the study of the quality of time integrators coupling. The structure is discretized by the linear finite element method and solved in time by the Newmark scheme, whereas the finite difference and finite volume methods are used for the fluid subdomain. By projecting the fluid equations on the eigenvector, we obtain a compatibility relation that corresponds to the characteristic line that transfers the fluid information from the inside to the fluid–structure interface. With this appropriate compatibility relation and the solid equations, the interface status is predicted for the next time step, ensuring zero interface energy. Hence, the order of accuracy and the stability are preserved to the minimum level of the two parts, the structure and the fluid. Furthermore, the coupling strategy allows incompatible time steps. Some numerical results are obtained for the 1‐D linear problem, and a good agreement with the analytical solution has been found. Copyright © 2013 John Wiley & Sons, Ltd.     

近场波动模拟的一种应力人工边界

采用平面波和远场散射波混合透射，引入无限介质线弹性本构关系建立了一种应力人工 边界条件. 其优点在于边界结点反应与内部有限元结点反应采用相同的积分格式计算，有限 元积分方法稳定时不存在人工边界失稳问题. 数值算例表明：边界精度高于现有黏性边 界、黏弹性人工边界，以及一、二阶透射人工边界.     

《余能理论》引领中国变分原理的研究

正今年是钱令希院士诞辰100周年,带着崇敬的心情我们缅怀先生的一生。作为一名杰出的科学家,除了在很多研究工作中取得优秀的成果,钱令希先生的战略眼光更值得我们学习。1950年钱令希先生在中国科学杂志发表《余能理论》[1]。论文中钱令希先生引用了Westergaard 1941年关于余能原理的论文,特别引用了Westergaard的观点,认为余能方法没有受到与其价