Hybrid graded element model for transient heat conduction in functionally graded materials

This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest’s algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.     

Fundamental-solution-based finite element model for plane orthotropic elastic bodies

A new hybrid finite element formulation is presented for solving two-dimensional orthotropic elasticity problems. A linear combination of fundamental solutions is used to approximate the intra-element displacement fields and conventional shape functions are employed to construct elementary boundary fields, which are independent of the intra-element fields. To establish a linkage between the two independent fields and produce the final displacement-force equations, a hybrid variational functional containing integrals along the elemental boundary only is developed. Results are presented for four numerical examples including a cantilever plate, a square plate under uniform tension, a plate with a circular hole, and a plate with a central crack, respectively, and are assessed by comparing them with solutions from ABAQUS and other available results.     

Dual reciprocity hybrid radial boundary node method for Winkler and Pasternak foundation thin plate

An efficient dual reciprocity hybrid radial boundary node method is developed for the analysis of Winkler and Pasternak foundation thin plate, in which a hybrid displacement variational principle, radial point interpolation method (RPIM) and dual reciprocity method (DRM) are combined. Firstly, the hybrid displacement variational principle is developed, in which the domain variables are interpolated by two groups of symmetric fundamental solutions, while the boundary variables are interpolated by RPIM instead of the traditional moving least square, and the shape function obtained by RPIM satisfies the delta function property, so boundary conditions can be applied directly. Besides, DRM is exploited to evaluate the particular solutions of inhomogeneous terms, which can be used to transform the domain integrals arising from the inhomogeneous term into equivalent boundary integrals. Finally, some additional equations based on the DRM theory are proposed to overcome the problem that the boundary integral equations are not enough to solve all variables. This method has the advantages of both no element mesh of meshless method and dimensionality reduction of boundary element method. Numerical examples of Winkler and Pasternak foundation plates are given to illustrate that the present method is effective, accurate and it can be further expanded into practical engineering.     

功能梯度材料结构分析的半解析数值方法研究

一种典型的半解析数值方法——线法被引入功能梯度材料的结构分析。首先推导了功能梯度材料位移形式的平衡方程和边界条件,然后阐述了线法功能梯度材料结构分析的基本步骤和数值原理。该方法的基本思想是通过有限差分将问题的控制方程半离散为定义在沿梯度方向离散节线上的常微分方程组,然后应用B样条函数Gauss配点法求解该常微分方程组得到问题的解答。为演示线法在功能梯度材料结构分析中的应用,给出了线性梯度和指数梯度功能梯度材料板分别受恒定位移、均匀拉伸载荷和弯曲载荷作用的数值算例。与相应问题解析解和其他数值方法的比较表明,线法的计算结果具有很高的精度,而且不需要任何特殊的考虑就能够有效模拟材料内部物性参数的连续变化,也无需事先选取满足特定条件的待定场函数,是一种非常适合功能梯度材料结构形式和材料特点的半解析数值方法。     

周期分布多边形夹杂奇异性应力干涉问题的研究

采用一种新型的杂交元模型和一种单胞模型来解决周期分布多边形夹杂角部的奇异性应力相互干涉的问题。新型杂交元模型是基于广义Hellinger-Reissner变分原理建立的,其中奇异性应力场分量和位移场分量是采用有限元特征分析法的数值特征解得到的。使用当前的新型杂交元模型,只需要在夹杂角部邻域的周界上划分一维单元,避免了像传统有限元模型那样需要划分高密度二维单元。文中给出了代表奇异性应力场强度的夹杂角部广义应力强度因子数值解,并考虑材料属性、夹杂尺寸和夹杂位置关系的影响。算例中,考虑了夹杂和基体完全接合的情况,并给出了考核例。结果表明:当前模型能得到高精度数值解,且收敛性好;与传统有限元法和积分方程方法相比,该模型更具有通用性,为非均质材料的细观力学分析打下了基础。     

HYBRID FEM WITH FUNDAMENTAL SOLUTIONS AS TRIAL FUNCTIONS FOR HEAT CONDUCTION SIMULATION

A new type of hybrid finite element formulation with fundamental solutions as internal interpolation functions, named as HFS-FEM, is presented in this paper and used for solving two dimensional heat conduction problems in single and multi-layer materials. In the proposed approach, a new variational functional is firstly constructed for the proposed HFS-FE model and the related existence of extremum is presented. Then, the assumed internal potential field constructed by the linear combination of fundamental solutions at points outside the elemental domain under consideration is used as the internal interpolation function, which analytically satisfies the governing equation within each element. As a result, the domain integrals in the variational functional formulation can be converted into the boundary integrals which can significantly simplify the calculation of the element stiffness matrix. The independent frame field is also introduced to guarantee the inter-element continuity and the stationary condition of the new variational functional is used to obtain the final stiffness equations. The proposed method inherits the advantages of the hybrid Trefftz finite element method （HT-FEM） over the conventional finite element method （FEM） and boundary element method （BEM）, and avoids the difficulty in selecting appropriate terms of T-complete functions used in HT-FEM, as the fundamental solutions contain usually one term only, rather than a series containing infinitely many terms. Further, the fundamental solutions of a problem are, in general, easier to derive than the T-complete functions of that problem. Finally, several examples are presented to assess the performance of the proposed method, and the obtained numerical results show good numerical accuracy and remarkable insensitivity to mesh distortion.     

A preparation method of functionally graded materials with phase transition under shock loading

The propagation of phase boundary in a material undergoing shock induced irreversible phase transition is studied theoretically using a model based on simple-mixture rule. It is found that along with the decay of the phase boundary, a functionally graded material (FGM) forms in the mixed-phase region. Such FGMs are composed of parent phase and product phase, and the composition and physical properties are changing continuously without apparent macro-interfaces. The effect of stress boundary conditions on formation of the FGM is investigated in detail with a numerical method. The possibility of producing FGMs with impact method is proposed and the limit of this method is discussed.     

动态载荷下功能梯度复合材料的圆币形裂纹问题

研究了动态载荷下功能梯度材料中的圆币形裂纹问题．假设材料为横观各向同性,并且含有多个垂直于厚度方向的裂纹,材料参数沿轴向（与裂纹面垂直的方向）为变化的,沿该方向将材料划分为许多单层,各单层材料参数为常数,利用Hankel变换祛,在Laplace域内推导出了控制问题的对偶积分方程组．利用Laplace数值反演,得出了裂纹尖端的动态应力强度因子和能量释放率．研究了含两个裂纹的功能梯度接头结构,分析了材料非均匀性参数对应力强度因子和能量释放率的影响．     

VARIATIONAL PRINCIPLES FOR HYDRODYNAMIC IMPACT PROBLEMS

We first establish the rigorous field equations of the two continuous stages before and after entering water. Then correspondently, we obtain the specific variational principles, bounded theorems, and boundary integral equations of the second stage problems. The existence of solutions are proved and the scheme of solving the solutions are provided. Finally, as a numerical example, the ship's wave resistence problem is used to demonstrate the specific application of the second stage problems and its accuracy. Then we provide a rigorous and sound theoretical basis of variational finite element method and boundary element method for calculating the accurately fundamental equations.     

A boundary integral method for computing forces on particles in unsteady Stokes and linear viscoelastic fluids

Accurately characterizing the forces acting on particles in fluids is of fundamental importance for understanding particle dynamics and binding kinetics. Conventional asymptotic solutions may lead to poor accuracy for neighboring particles. In this paper, we develop an accurate boundary integral method to calculate forces exerted on particles for a given velocity field. We focus our study on the fundamental two‐bead oscillating problem in an axisymmetric frame. The idea is to exploit a correspondence principle between the unsteady Stokes and linear viscoelasticity in the Fourier domain such that a unifying boundary integral formulation can be established for the resulting Brinkman equation. In addition to the dimension reduction vested in a boundary integral method, our formulation only requires the evaluation of single‐layer integrals, which can be carried out efficiently and accurately by a hybrid numerical integration scheme based on kernel decompositions. Comparison with known analytic solutions and existing asymptotic solutions confirms the uniform third‐order accuracy in space of our numerical scheme. Copyright © 2016 John Wiley & Sons, Ltd.     

A meshfree method for transverse vibrations of anisotropic plates

A meshfree approach, called displacement boundary method, for anisotropic Kirchhoff plate dynamic analysis is presented. This method is deduced from a variational principle, which uses a modified hybrid functional involving the generalized displacements and generalized tractions on the boundary and the lateral deflection in the domain as independent variables. The discretization process is based on the employment of the fundamental solutions of the static problem operator for the expression of the variables involved in the functional. The stiffness and mass matrices obtained for the dynamic model are frequency-independent, symmetric and positive definite and their computation involves boundary integrals of regular kernels only. Due to its features, the final resolving system can be solved with the classical approaches by using standard numerical procedures. To assess the formulation, the free vibrations of some anisotropic plates were calculated and the results compared with those obtained using other solution techniques. The present results are in good agreement with those found in the literature showing the accuracy and effectiveness of the proposed approach.     

黏弹性功能梯度材料裂纹问题的有限元方法

针对组分材料体积含量任意分布的黏弹性功能梯度材料裂纹问题建立有限元分析途径. 通过Laplace变换,将黏弹性问题转化到象空间中求解,基于反映材料非均匀的梯度单元和裂纹尖端奇异特性的奇异单元计算象空间中的位移、应力和应变场,应用虚拟裂纹闭合方法得到应变能释放率,分别由应力和应变能释放率确定应力强度因子. 给出这些断裂参量在物理空间和象空间之间的对应关系,由数值逆变换求出其在物理空间的相应值. 文中分析两端均匀受拉的黏弹性边裂纹板条,首先针对松弛模量表示为空间函数和时间函数乘积的特殊梯度材料进行计算,结合对应原理验证方法的有效性. 然后分析组分材料体积含量具有任意梯度分布的情形,由Mori-Tanaka方法预测象空间中的等效松弛模量. 计算结果表明,蠕变加载条件下,应变能释放率随时间增加,其增大程度与黏弹性组分材料体积含量相关. 由于梯度材料的非均匀黏弹性性质,产生应力重新分布,导致应力强度因子随时间变化,其变化范围与组分材料的体积含量分布方式有关.     

敷设多孔吸声材料声腔特征值分析的径向积分边界元法

由于Helmholtz方程的基本解是频率的函数,因此传统边界元法在处理声场特征值问题时具有天生的缺陷。本文采用Laplace方程基本解生成积分方程,通过径向积分法将在此过程中产生的域积分项转化为边界积分。此方法克服了传统边界元法系数矩阵对频率的依赖,同时克服了特解积分法对特解的依赖,并通过对表面声导纳的多项式逼近,将敷设多孔吸声材料声腔特征值问题转化为矩阵多项式,从而避免了复杂的非线性求解。通过数值算例验证了算法的有效性。     

Dual variational formulation for Trefftz finite element method of elastic materials

A dual variational principle is presented for Trefftz finite element analysis. The proof of the stationary conditions of the variational functional and the theorem on the existence of extremum are provided in this paper. They are boundary displacement condition, surface traction condition and interelement continuity condition. Based on the assumed intraelement and frame fields, element stiffness matrix equation is obtained which can easily be implemented into computer programs for numerical analysis with Trefftz finite element method. Two numerical examples are considered to illustrate the effectiveness and applicability of the proposed element model.     

三维势流场的比例边界有限元求解方法

比例边界有限元法(SBFEM)是线性偏微分方程的一种新的数值求解方法。该方法只对计算域边界利用Galerkin方法进行数值离散,相对于有限元方法(FEM)减少了一个空间坐标的维数,而在减少的空间坐标方向利用解析方法进行求解;相对于边界元法(BEM),比例边界有限元方法不需要基本解,避免了奇异积分的计算,所以它结合了有限元和边界元方法的优点。本文建立了利用比例边界有限元法求解三维Laplace方程的数值模型并用于计算三维物体周围的水流场,将计算结果与解析解和边界元方法进行了对比,结果表明此方法可以很好地模拟水流场,且具有较高的计算精度。     

压电材料平面裂纹尖端场的杂交应力有限元分析

基于复势理论和杂交变分原理建立了一种适用于力电耦合分析的杂交应力有限元模 型. 给出了建立刚度矩阵的主要公式和推导过程，单元内的位移场和应力场采用满足平 衡方程的复变函数级数解，假设的复变函数级数解事先精确满足裂纹的无应力和电位移法向 分量为零的条件，单元外边界的位移场假设按抛物线变化, 单元的刚度矩阵采用Gauss积分的方法得出. 通过对力电耦合裂尖场的数值计算验证了程序 的正确性和单元的有效性，同时也用所得结果校验了理论解.     

A high order control volume finite element procedure for transient heat conduction analysis of functionally graded materials

A robust unstructured control-volume finite element method is presented for the solution of two-dimensional transient heat conduction in functionally graded materials (FGMs) with isotropic properties. The material properties at a point in the domain vary exponentially to spatial coordinates. A triangular mesh is chosen for spatial discretization and a fully implicit scheme is adapted for time discretization. Several problems are investigated and the results are successfully validated by using analytical and other numerical solutions available in the literature.     

A HYBRID METHOD FOR ANALYZING THE DYNAMIC RESPONSES OF CAVITIES OR SHELLS BURIED IN AN ELASTIC HALF-PLANE

In this paper, based on a variational formalism which originally proposed by Mei [1] for infinite elastic medium and extended by Yeh, et al. [2,3] for elastic half-plane, a hybrid method which combines the finite element and series expansion method is implemented to solve the diffraction of plane waves by a cavity buried in an elastic half-plane. The finite domain which encloses all inhomogeneities including the cavity can be easily formulated by finite element methods. The unknown boundary data obtained by subtracting the known free fields from the total fields which include the boundary nodal displacements and tractions at the interface between the finite domain and the surrounding elastic half-plane are not independent of each other and can be correlated through aseries repre sentation. Due to the continuity condition at the interface, the same series representation is still valid for the exterior elastic half-plane to represents the scattered wave. The unknown coefficients of this series are treated as generalized coordinates and can be easily formulated by the same variational principle. The expansion function of the series is composed of basis function. Each basis function is constructed from the basis function for an infinite plane by superimposing an additional homogeneous reflective term to satisfy both traction free conditions at ground surface and radiation conditions at infinity. The numerical results are made against those obtained by boundary element methods, and good agreements are found.     

Dynamic response for functionally graded materials with penny-shaped cracks

This paper provides a method for studying the penny-shaped cracks configuration in functionally graded material(FGM) structures subjected to dynamic or steady loading. It is assumed that the FGMs are transversely isotropic and all the material properties only depend on the axial coordinatez. In the analysis, the elastic region is treated as a number of layers. The material properties are taken to be constants for each layer. By utilizing the Laplace transform and Hankel transform technique, the general solutions for the layers are derived. The dual integral equations are then obtained by introducing the mechanical boundary and layer interface conditions via the flexibility/stiffness matrix approach. The stress intensity factors are computed by solving dual integral equations numerically in Laplace transform domain. The solution in time domain is obtained by utilizing numerical Laplace inverse. The main advantage of the present model is its ability for treating multiple crack configurations in FGMs with arbitrarily distributed and continuously varied material properties by dividing the FGMs into a number of layers with the properties of each layer slightly different from one another. This work was supported by Failure Mechanics Laboratory of State Education Commission and the Post-doctor Research Fund of China.     

Two-dimensional modeling of heterogeneous structures using graded finite element and boundary element methods

In the present work, graded finite element and boundary element methods capable of modeling behaviors of structures made of nonhomogeneous functionally graded materials (FGMs) composed of two constituent phases are presented. A numerical implementation of Somigliana’s identity in two-dimensional displacement fields of the isotropic nonhomogeneous problems is presented using the graded elements. Based on the constitutive and governing equations and the weighted residual technique, effective boundary element formulations are implemented for elastic nonhomogeneous isotropic solid models. Results of the finite element method are derived based on a Rayleigh–Ritz energy formulation. The heterogeneous structures are made of combined ceramic–metal materials, in which the material properties vary continuously along the in-plane or thickness directions according to a power law. To verify the present work, three numerical examples are provided in the paper.