弹性力学轴对称问题的有限元线法

给出了解弹性力学空间轴对称问题的有限元线法的基本理论。该法包括了２－４条结线的等参数单元，沿结线方向的两点边值问题采用插值矩阵法解之。算例表明，本法具有良好的收敛性和较高的计算精度。     

Sliding frictional contact between a rigid punch and a laterally graded elastic medium

Analytical and computational methods are developed for contact mechanics analysis of functionally graded materials (FGMs) that possess elastic gradation in the lateral direction. In the analytical formulation, the problem of a laterally graded half-plane in sliding frictional contact with a rigid punch of an arbitrary profile is considered. The governing partial differential equations and the boundary conditions of the problem are satisfied through the use of Fourier transformation. The problem is then reduced to a singular integral equation of the second kind which is solved numerically by using an expansion–collocation technique. Computational studies of the sliding contact problems of laterally graded materials are conducted by means of the finite element method. In the finite element analyses, the laterally graded half-plane is discretized by quadratic finite elements for which the material parameters are specified at the centroids. Flat and triangular punch profiles are considered in the parametric analyses. The comparisons of the results generated by the analytical technique to those computed by the finite element method demonstrate the high level of accuracy attained by both methods. The presented numerical results illustrate the influences of the lateral nonhomogeneity and the coefficient of friction on the contact stresses.     

Transient dynamic responses of a cracked solid subjected to in-plane loadings

Transient dynamic responses of an elastic cracked solid subjected to in-plane surface loadings are investigated in this study. Two vertical cracks, a surface-breaking crack and a sub-surface crack, are considered. The frequency responses of the plane strain problem are calculated by the computational mechanics combining the finite element method with the boundary integral equation. The finite element method is used for the near-field enclosing the crack, while the boundary integral equation is applied for the far-field to satisfy the Sommerfeld radiation condition. The transient responses are then obtained using fast Fourier transform. Surface displacements, crack opening displacements, and dynamic stress intensity factors are presented to show the significant effects of the cracks. The interaction between the elastic waves and the cracks as well as the mode conversion phenomena can be observed and understood through the numerical simulations.     

A review of boundary and finite element methods in fracture mechanics

Reviewed in this work are the methods of finite and boundary element as applied to solve fracture mechanics problems. The former requires the discretization of the interior of the domain while the latter involves computing an integral equation over the boundary of the domain. Applications of these methods are made to two-dimensional elastic crack problems. Efficiency and accuracy of different approaches are discussed and compared by examples. The boundary element procedure employing special Green's functions for the plane crack problem is shown to be superior. The correlation between the hybrid element formulations and boundary element regions embedded into a finite element model is also given.     

A novel hybrid finite element analysis of inplane singular elastic field around inclusion corners in elastic media

This paper deals with the inplane singular elastic field problems of inclusion corners in elastic media by an ad hoc hybrid-stress finite element method. A one-dimensional finite element method-based eigenanalysis is first applied to determine the order of singularity and the angular dependence of the stress and displacement field, which reflects elastic behavior around an inclusion corner. These numerical eigensolutions are subsequently used to develop a super element that simulates the elastic behavior around the inclusion corner. The super element is finally incorporated with standard four-node hybrid-stress elements to constitute an ad hoc hybrid-stress finite element method for the analysis of local singular stress fields arising from inclusion corners. The singular stress field is expressed by generalized stress intensity factors defined at the inclusion corner. The ad hoc finite element method is used to investigate the problem of a single rectangular or diamond inclusion in isotropic materials under longitudinal tension. Comparison with available numerical results shows the present method is an efficient mesh reducer and yields accurate stress distribution in the near-field region. As applications, the present ad hoc finite element method is extended to discuss the inplane singular elastic field problems of a single rectangular or diamond inclusion in anisotropic materials and of two interacting rectangular inclusions in isotropic materials. In the numerical analysis, the generalized stress intensity factors at the inclusion corner are systematically calculated for various material type, stiffness ratio, shape and spacing position of one or two inclusions in a plate subjected to tension and shear loadings.     

广义热弹模型的热力学基础与瞬态响应

微纳科技的快速发展与超短脉冲激光技术的广泛运用, 对描述微纳尺度超快热冲击的广义热传导及其热弹耦合理论提出迫切需求. 基于拓展热力学原理, 本文建立了考虑热传导双相滞后效应和高阶热流率的广义热弹耦合理论. 类比于力学领域黏弹性本构关系的串联、并联模型, 并受Green-Naghdi (GN)广义热传导模型启发, 本文提出了热学“弹性”单元和“黏性”单元模型, 并采用串联、并联方法实现了Cattaneo-Vernotte (CV)、GN、双相滞后(DPL)和Moore-Gibson-Thompson (MGT) 热传导模型的重构. 理论推导进一步表明, 本文新建模型对应于热学Burgers模型, 并得到了新模型中各相位滞后中松弛时间之间的比例关系. 运用拉普拉斯变换方法, 研究了一维结构受边界热冲击和移动热源作用下的瞬态响应, 计算结果表明: 新模型克服了热波速度无限大的悖论; 仅有边界热冲击载荷时, 新模型得到的响应结果均较大, 响应范围最小; 相比于无热源作用情形, 受移动热源作用时, 新模型会产生更大的峰值响应. 新模型与经典弹性理论耦合构建了广义热弹性理论, 运用该理论, 可以清晰观察到在热波和弹性波波前的应力突变. 理论方面, 本文推动了拓展热力学与连续介质力学的结合, 对于远离平衡态极端力学基础理论问题的研究具有启发意义; 应用方面, 本文研究结果可为激光等移动热源作用下材料的瞬态响应分析提供理论基础和数值方法.      

Finite element approximation of field dislocation mechanics

A tool for studying links between continuum plasticity and dislocation theory within a field framework is presented. A finite element implementation of the geometrically linear version of a recently proposed theory of field dislocation mechanics (J. Mech. Phys. Solids 49 (2001) 761; Proc. Roy. Soc. 459 (2003) 1343; J. Mech. Phys. Solids 52 (2004) 301) represents the main idea behind the tool. The constitutive ingredients of the theory under consideration are simply elasticity and a specification of dislocation velocity and nucleation. The set of equations to be approximated are non-standard in the context of solid mechanics applications. It comprises the standard second-order equilibrium equations, a first-order div-curl system for the elastic incompatibility, and a first-order, wave-propagative system for the evolution of dislocation density. The latter two sets of equations require special treatment as the standard Galerkin method is not adequate, and are solved utilizing a least-squares finite element strategy. The implementation is validated against analytical results of the classical elastic theory of dislocations and analytical results of the theory itself. Elastic stress fields of dislocation distributions in generally anisotropic media of finite extent, deviation from elastic response, yield-drop, and back-stress are shown to be natural consequences of the model. The development of inhomogeneity, from homogeneous initial conditions and boundary conditions corresponding to homogeneous deformation in conventional plasticity, is also demonstrated. To our knowledge, this work represents the first computational implementation of a theory of dislocation mechanics where no analytical results, singular solutions in particular, are required to formulate the implementation. In particular, a part of the work is the first finite element implementation of Kröner's linear elastic theory of continuously distributed dislocations in its full generality.     

广义节点有限元法

应用流形方法思想,通过引入广义节点的概念,对传统有限元方法进行改进,建立了可具有任意高阶多项式托值函数的广义节点有限元方法,计算结果表明,广义节点有限元方法较之传统有限元方法有较高的精度。     

A layerwise boundary integral equation model for layers and layered media

A hybrid method is presented for the analysis of layers, plates, and multilayered systems consisting of isotropic and linear elastic materials. The problem is formulated for the general case of a multilayered system using a total potential energy formulation. The layerwise laminate theory of Reddy is employed to develop a layerwise, two-dimensional, displacement-based, hybrid boundary element model that assumes piecewise continuous distribution of the displacement components through the system's thickness. A one-dimensional finite element model is used for the analysis of the multilayered system through its thickness, and integral Fourier transforms are used to obtain the exact solution for the in-plane problem. Explicit expressions are obtained for the fundamental solution of a typical infinite layer (element) assuming linear displacement distribution through its thickness. This fundamental solution is given in a closed form in the cartesian space, and it can be applied in the two-dimensional boundary integral equation model to analyze layered structures with finite dimensions. The proposed method provides a simple, efficient, and versatile model for a three-dimensional analysis of thick plates or multilayered systems.     

非连续变形计算力学模型的粘弹性分析方法Ⅰ:基本理论

基于粘弹性广义有限单元和接触力元，发展了适用于多体相互作用系统非连续变形分析的粘弹性数值分析方法，通过虚功原理，给出了其分区参变量最小势能原理，从而阐明了其理论基础。粘弹性广义有限单元的本构关系可由粘弹性退化为弹性或刚性，因此本文所提出的方法可对由刚体、弹性体和粘弹性体所构成的复杂多体系统在外荷载作用下的力学行为进行数值模拟，同时能够比本文精确地直接得到多体之间的接触应力。     

A new formula for the C-matrix in the Somigliana identity

By making use of a convenient decomposition of the fundamental tractions, a new formula for the C-matrix in the Somigliana identity for a three- or two-dimensional elastic isotropic body is derived. This kind of formula is more advantageous for analytical and computational C-matrix evaluations than the currently well-known formula. A general closed analytical formula of the C-matrix for the case of any finite number of tangent planes to the boundary of the body at a non-smooth boundary point, presented in the final section of this paper, demonstrates the usefulness of the new formula.     

结构工程科学中若干计算结构力学问题的研究展望

本文概述了计算结构力学的形成过程,评述了结构工程科学中以下７个需要重点研究的计算结构力学问题:数值方法基本理论问题（误差估计理论、网格自适应加密技术、多变量有限元理论和半解析数值方法）、工程结构优化设计、结构施工力学、计算机数值模拟和仿真技术、本构模型、计算机技术新发展的影响、计算机辅助设计。并论述了数值计算、理论和试验这三者之间相互依赖、相互促进、相互交叉又相互制约的辩证关系。      

热-机载荷下不规则夹杂的奇异性应力场分析

基于有限元特征分析法得到的夹杂角部场数值特征解开发了一种超级奇异单元模型,并将其与普通四节点单元紧密结合,用于热-机载荷下夹杂角端部的应力场分析。在数值计算中,考察了热-机载荷下不同弹性比和不同夹杂尺寸的应力强度因子,并将所得结果与文献解和传统有限元方法解比对。结果表明,本文方法对热-机耦合条件下的不规则夹杂角端部的热弹性应力分析极为有效,可避免局部网格的高度加密,并提高计算效率。模型在复合材料夹杂的局部强度问题分析方面具有很好的实用性。     

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

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

不同模量弹性问题理论及有限元法研究进展

随着科学技术的日益发展，对材料力学性质的研究提出了更高的要求，研制新型的材料以及挖掘材料自身特性的潜力，已成为新的研究动向．简述了不同模量弹性问题理论及其有限元法的研究与发展．利用等效的概念，对Ambartsumyan有限元计算模型、Jones有限元计算模型、张允真等有限元计算模型、叶志明等计算模型进行改进和探讨．通过不同模量弹性理论及其有限元方法在实际工程构件问题分析中的应用表明，若沿用相同模量弹性理论或有限元对有关问题进行计算，其结果将与采用不同模量模型材料所得的刚度和强度有较大的偏差。     

Numerical solution of axisymmetric cavity flows using the boundary element method

This paper presents a formulation of the boundary element method (BEM) for solution of axisymmetric cavity flow problems. The governing equation is written in terms of Stokes' stream function, requiring a new fundamental solution to be found. The iterative procedure for adjusting the free-surface position is similar to that used for planar cavity flows. Numerical results are compared with finite difference and finite element solutions, showing the robustness of the BEM model.     

Higher-order accurate implicit time integration schemes for transport problems

The present paper is concerned with the numerical solution of transient transport problems by means of spatial and temporal discretization methods. The generalized initial boundary value problem of various nonlinear transport phenomena like heat transfer or mass transport is discretized in space by p-finite elements. After finite element discretization, the resulting first-order semidiscrete balance has to be solved with respect to time. Next to the classical generalized-α integration method predicated on the Newmark approach and the evaluation at a generalized midpoint also implicit Runge–Kutta time integration schemes, are presented. Both families of finite difference-based integration schemes are derived for general first-order problems. In contrast to the above-mentioned algorithms, temporal discontinuous and continuous Galerkin methods evaluate the balance equation not at a selected time instant within the timestep, but in an integral sense over the whole time step interval. Therefore, the underlying semidiscrete balance and the continuity of the primary variables are weakly formulated within time steps and between time steps, respectively. Continuous Galerkin methods are obtained by the strong enforcement of the continuity condition as special cases. The introduction of a natural time coordinate allows for the application of standard higher-order temporal shape functions of the p-Lagrange type and the well-known Gau?–Legendre quadrature of associated time integrals. It is shown that arbitrary order accurate integration schemes can be developed within the framework of the proposed temporal p-Galerkin methods. Selected benchmark analyses of calcium diffusion demonstrate the properties of all three methods with respect to non-smooth initial or boundary conditions. Furthermore, the robustness of the present time integration schemes is also demonstrated for the highly nonlinear reaction–diffusion problem of calcium leaching, including the pronounced changes of the reaction term and non-smooth changes of Dirichlet boundary conditions of calcium dissolution.     

固体力学边界元分析程序包BESMAP——基本结构功能

要使边界元法象有限元法那样得到广泛应用，必须开发与有限元法相媲美的大型边界元法程序包。国外虽出现少数几个边界元法程序包，但由于结构、功能等方面的原因，还不能成为边界元法广泛应用的工具。针对这种情况，本文借鉴现有的先进的有限元法程序包，研制了大型固体力学边界元法程序包ＢＥＳＭＡＰ，ＢＥＳＭＡＰ在总体设计、计算能力、使用方便和功能齐全等方面作了研究和改进。     

算子自定义小波弹性板单元构造及自适应分析

提出一类用于分析弹性板问题的算子自定义小波弹性板单元构造方法。该方法的优点在于根据工程问题的求解需要灵活构造具有解耦特性的算子自定义小波基,使得系统多尺度刚阵具有沿对角线的强稀疏性,从而实现了该算法在每个尺度上独立、快速求解,系统方程的求解效率得到较大提高。建立多分辨Lagrange有限元空间和多尺度计算理论,提出基于稳定完备法的算子自定义小波弹性板单元构造方法及解耦条件。依据两尺度相对误差估计,提出自适应算子自定义小波有限元算法。数值算例证明,算子自定义小波弹性板单元具有求解精度与计算效率高等特点。     

修正的偶应力线弹性理论及广义线弹性体的有限元方法

以含偶应力的弹性理论为基础,考虑小变形情况下变形体的平动变形和旋转变形,提出关于偶应力与曲率张量的线性本构关系,建立一般弹性体的线性模型。为满足有限单元C1连续性要求,考虑转角为独立变量,利用罚方法引入约束条件,构造一般弹性体的约束变分形式。应用8节点48个自由度的实体等参元,建立一般弹性体力学响应分析的有限元方程。对悬臂梁的静力和动力分析表明,一般弹性体模型较之经典弹性力学更适合结构分析;较之Timoshenko梁模型,一般弹性体模型能够计及结构尺度对结构动力特性和动力响应造成的显著影响。