沿板平面材料梯度变化功能材料板件分析

本文发展一种功能梯度构件分析的细观元法.细观元法在构件的常规有限元内部设置密集观细单元以反映材料特性变化,又通过协调条件将各细观元结点自由度转换为同一常规有限元自由度,再上机计算.这种细观元法既能充分反映材料功能梯度沿各方向任意变化特性;而其计算单元又和常规有限元一样,是一种针对功能梯度结构分析的有效数值方法.现有功能梯度板件分析中无论对不同形状还是不同边界的功能梯度构件,其材料特性均沿板厚度方向梯度变化,本文用细观元法进行计算与分析,给出了目前尚未得到的沿板平面方向功能梯度变化构件的力学量三维分布形态.     

具有不同功能梯度分布函数的板件的三维分析

发展了一种新颖的功能梯度结构分析的细观元法。细观元法对结构的常规有限单元内部设置密集细观单元以反映材料特性梯度变化,又通过协调条件将各细观元结点自由度转换为同一常规有限元自由度,再上机计算。这种细观元法既能充分反映材料功能梯度变化特性;而其计算单元又和常规有限元一样,是一种针对功能梯度结构分析的有效数值方法。文中通过细观元技术进行计算、分析,给出了具有不同功能梯度分布函数板件的力学量三维分布形态。     

功能梯度材料构件三维分析的细观元模型

提出一种新颖的功能梯度构件分析的细观元法,给出了方法模型、基本算式及特点与功能。细观元法对构件的常规有限单元内部设置密集细观单元以反映材料特性梯度变化,又通过协调条件将各细观元结点自由度转换为同一常规有限元自由度,再上机计算。这种细观元法既能充分反映材料功能梯度及组分变化特性,而其计算单元与自由度又与常规有限元一样,是一种针对功能梯度构件分析的有效数值方法。算例表明了细观元法对不同情况下功能梯度构件分析的适应性与精度。     

空隙分布对功能梯度板件宏观响应的影响

提出了一种新的功能梯度结构分析的细观元法. 细观元法对结构的常 规有限单元内部设置密集细观单元以反映材料细观构造,又通过协调条件将各细观元结点自 由度转换为同一常规有限元自由度. 此方法实现了材料细观结构到构件宏观响 应的直接过渡分析. 用细观元法来分析细观构造上空隙存在对宏观响应影响,给出了含 有不同分布与形状、位置空隙的功能梯度板件的力学量三维分布形态,实现功能梯度材料结 构宏、细观间跨尺度分析.     

AXISYMMETRIC BENDING OF TWO-DIRECTIONAL FUNCTIONALLY GRADED CIRCULAR AND ANNULAR PLATES

Assuming the material properties varying with an exponential law both in the thick- ness and radial directions,axisymmetric bending of two-directional functionally graded circular and annular plates is studied using the semi-analytical numerical method in this paper.The de- flections and stresses of the plates are presented.Numerical results show the well accuracy and convergence of the method.Compared with the finite element method,the semi-analytical nu- merical method is with great advantage in the computational efficiency.Moreover,study on ax- isymmetric bending of two-directional functionally graded annular plate shows that such plates have better performance than those made of isotropic homogeneous materials or one-directional functionally graded materials.Two-directional functionally graded material is a potential alter- native to the one-directional functionally graded material.And the integrated design of materials and structures can really be achieved in two-directional functionally graded materials.     

功能梯度材料板件三维分析的半解析梯度有限元法

将半解析有限元与梯度有限元相结合,形成一种半解析梯度有限元来求解功能梯度材料板件问题。该方法兼有有限元法的适应性强、程序统一,半解析有限元法的节省单元与计算工作量,梯度有限元法的适应构件内部材料性能任意梯度分布等特点,并实现用一维数值计算给出构件三维分析结果。算例分析表明了方法的精度、功能与上述特点,充分揭示了功能梯度材料板件力学响应的三维形态。半解析梯度有限元法可推广应用到其他功能梯度材料面结构的各类分析中。     

基于金相图片识别功能梯度材料参数的梯度分布函数

功能梯度材料具有复杂的细观结构.其内部构造远比匀质材料复杂.以目前的实验条件,测量功能梯度材料参数分布函数一般是十分困难的.该文提出一种细观元方法,直接利用功能梯度材料金相图片信息,来进行功能梯度构件宏观响应的跨尺度力学分析.另在功能梯度材料结构静、动力正演分析基础上,根据正演所求得的挠度或固有频率值,采用混和罚函数法来反演识别功能梯度材料参数的宏观梯度分布函数.     

基于高阶剪切变形理论的四边形求积元板单元及其应用

近年来由各类新型复合材料或功能梯度材料构成的板结构在工程领域得到了广泛应用,其显著特点是材料性能沿板厚变化.为合理考虑横向剪切应变,许多学者基于Reddy高阶剪切变形理论,构建了不同的有限元单元对该类板结构进行分析,但其中满足$C^{1}$连续条件的单元相对较少.本文基于Reddy高阶剪切变形理论,采用求积元方法,建立了$C^{1}$连续的四边形板单元.利用该单元对均质材料、复合材料、功能梯度材料构成的等厚度矩形板、变厚度矩形板及等厚度斜板的线弹性弯曲和自由振动问题进行了计算分析,并与现有文献中的相应计算结果进行了对比.研究表明:基于高阶剪切变形理论的四边形求积元板单元具有较高的计算效率和良好的适应性,文中各类材料构成的等变厚度矩形板及等厚度斜板均只需1个单元即可得到理想的计算结果.对于等/变厚度矩形板,可仅使用9$\times$9个积分点,而对于等厚度斜板,随着斜角的增大,所需积分点的数目逐渐增多至15$\times $15.该四边形求积元板单元可进一步用于新型复合材料板的非线性分析.      

Three-dimensional stress and free vibration analyses of functionally graded plates with circular holes by the use of the graded finite element method

Static and free vibration analyses of plates with circular holes are performed based on the three-dimensional theory of elasticity. The plates are made of a functionally graded material (FGM), and the volume fractions of the constituent materials vary continuously across the plate. The effective properties of the FGM plate are estimated by using the Mori–Tanaka homogenization method. A graded finite element method based on the Rayleigh–Ritz energy formulation is used to solve the problem. Effects of different volume fractions of the materials and hole sizes on the behavior of FGM plates under uniaxial tension are investigated. Natural frequencies of a fully clamped FGM plate with a circular cutout are derived. The results obtained are compared with available experimental data.     

Static analysis of simply supported functionally graded and layered magneto-electro-elastic plates

In this article, static analysis of functionally graded, anisotropic and linear magneto-electro-elastic plates have been carried out by semi-analytical finite element method. A series solution is assumed in the plane of the plate and finite element procedure is adopted across the thickness of the plate such a way that the three-dimensional character of the solution is preserved. The finite element model is derived based on constitutive equation of piezomagnetic material accounting for coupling between elasticity, electric and magnetic effect. The present finite element is modeled with displacement components, electric potential and magnetic potential as nodal degree of freedom. The other fields are calculated by post-computation through constitutive equation. The functionally graded material is assumed to be exponential in the thickness direction. The numerical results obtained by the present model are in good agreement with available functionally graded three-dimensional exact benchmark solutions given by Pan and Han [Pan, E., Han, F., in press. Green’s function for transversely isotropic piezoelectric functionally graded multilayered half spaces. Int. J. Solids Struct.]. Numerical study includes the influence of the different exponential factor, magneto-electro-elastic properties and effect of mechanical and electric type of loading on induced magneto-electro-elastic fields. In addition further study has been carried out on non-homogeneous transversely isotropic FGM magneto-electro-elastic plate available in the literature [Chen, W.Q., Lee, K.Y., Ding, H.J., 2005. On free vibration of non-homogeneous transversely isotropic magneto-electro-elastic plates].     

Mechanical buckling and free vibration of thick functionally graded plates resting on elastic foundation using the higher order B-spline finite strip method

In this study, the mechanical buckling and free vibration of thick rectangular plates made of functionally graded materials (FGMs) resting on elastic foundation subjected to in-plane loading is considered. The third order shear deformation theory (TSDT) is employed to derive the governing equations. It is assumed that the material properties of FGM plates vary smoothly by distribution of power law across the plate thickness. The elastic foundation is modeled by the Winkler and two-parameter Pasternak type of elastic foundation. Based on the spline finite strip method, the fundamental equations for functionally graded plates are obtained by discretizing the plate into some finite strips. The results are achieved by the minimization of the total potential energy and solving the corresponding eigenvalue problem. The governing equations are solved for FGM plates buckling analysis and free vibration, separately. In addition, numerical results for FGM plates with different boundary conditions have been verified by comparing to the analytical solutions in the literature. Furthermore, the effects of different values of the foundation stiffness parameters on the response of the FGM plates are determined and discussed.     

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.     

功能梯度板动力分析的光滑有限元法研究

为提高功能梯度板动响应问题求解精度,基于一阶剪切变形板理论,提出了求解功能梯度板自由振动问题的Cell-Based 光滑有限元格式．功能梯度板材料属性沿厚度方向呈梯度连续变化,计算系统刚度矩阵时在光滑域内进行光滑梯度操作,可提高求解精度．采用Cell-Based 光滑有限元法,讨论了长厚比、形状因子和边界条件对两种典型功能梯度方板无量纲频率参数的影响,并与FEM法的解和文献中的解做了对比．数值算例的结果表明,光滑梯度操作可改善有限元系统的刚度,Cell-Based 光滑有限元法的计算结果更加逼近真实解,从而为功能梯度材料的进一步应用打下基础．     

基于等几何有限元法的功能梯度微板热力耦合屈曲预测

基于修正偶应力理论和Kirchhoff板理论，建立了功能梯度微板热力耦合屈曲等几何有限元模型。该模型仅包含一个材料尺度参数，能够描述尺度效应现象，且满足修正偶应力理论的高阶连续性要求。基于虚功原理推导了功能梯度微板热力耦合屈曲等几何有限元方程。通过对板的典型算例分析，讨论了材料尺度参数、边长比及梯度指数对板稳定性的影响。结果表明，本文模型预测的屈曲载荷总是大于宏观理论的结果，即捕捉到了尺度效应现象；随着临界屈曲力的增加，临界屈曲热载荷线性减少；此外，边长比和梯度指数也对微板的稳定性产生一定影响。     

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.     

Thermoelastic buckling behavior of thick functionally graded rectangular plates

Thermoelastic buckling behavior of thick rectangular plate made of functionally graded materials is investigated in this article. The material properties of the plate are assumed to vary continuously through the thickness of the plate according to a power-law distribution. Three types of thermal loading as uniform temperature raise, nonlinear and linear temperature distribution through the thickness of plate are considered. The coupled governing stability equations are derived based on the Reddy’s higher-order shear deformation plate theory using the energy method. The resulted stability equations are decoupled and solved analytically for the functionally graded rectangular plates with two opposite edges simply supported subjected to different types of thermal loading. A comparison of the present results with those available in the literature is carried out to establish the accuracy of the presented analytical method. The influences of power of functionally graded material, plate thickness, aspect ratio, thermal loading conditions and boundary conditions on the critical buckling temperature of aluminum/alumina functionally graded rectangular plates are investigated and discussed in detail. The critical buckling temperatures of thick functionally graded rectangular plates with various boundary conditions are reported for the first time and can be served as benchmark results for researchers to validate their numerical and analytical methods in the future.     

Stress Concentration Factor in Functionally Graded Plates with Circular Holes Subjected to Anti-plane Shear Loading

Stress concentration factors due to the presence of geometrical discontinuities (circular holes) in functionally graded plates are derived. The material property inhomogeneity is assumed to be in the radial direction originating at the center of the plate. Variable separable closed-form solutions are obtained for the stresses and displacements in functionally graded plates (without and with holes) subjected to anti-plane shear loading. The stresses in functionally graded plates without a hole are not homogeneous as it is in the case of homogeneous plates. Either a stress concentration (more than the applied stress) or dilution (less than the applied stress) occurs depending on whether the modulus increases (hardening graded material) or decreases (softening graded material) away from the center of the graded plate without a hole. A novel definition of the stress concentration factor due to the geometrical discontinuity in functionally graded plates is derived. The effect of the circular hole in functionally graded plates is to magnify (compared to homogeneous plates) the stress concentration when the modulus decreases away from the center of the hole (softening material). Beneficial reduction of the stress concentration factor is achieved in hardening functionally graded materials.     

A finite element formulation for analysis of functionally graded plates and shells

Summary A finite element formulation is derived for the thermoelastic analysis of functionally graded (FG) plates and shells. The power-law distribution model is assumed for the composition of the constitutent materials in the thickness direction. The procedure adopted to derive the finite element formulation contains the analytical through-the-thickness integration inherently. Such formulation accounts for the large gradient of the material properties of FG plates and shells through the thickness without using the Gauss points in the thickness direction. The explicit through-the-thickness integration becomes possible due to the proper decomposition of the material properties into the product of a scalar variable and a constant matrix through the thickness. The nonlinear heat-transfer equation is solved for thermal distribution through the thickness by the Rayleigh-Ritz method. According to the results, the formulation accounts for the nonlinear variation in the stress components through the thickness especially for regions with a variation in martial propperties near the free surfaces.     

功能梯度简支矩形板的非线性动力响应

研究了功能梯度简支矩形板在横向简谐激励作用下的非线性动力响应问题。采用幂律分布规律描述功能梯度材料的等效材料参数,基于Galerkin法建立了系统广义坐标的常微分控制方程。利用平均法得到了系统的幅频响应特性,分析了功能梯度矩形薄板的非线性主共振特性。数值算例验证了平均化方法的正确性,揭示了功能梯度平板主共振响应中的多值性和跳跃现象;同时分析发现初始条件会改变功能梯度平板主共振的响应幅值。最后讨论了功能梯度材料的梯度指数对系统幅值响应的影响。     

基于物理中面FGM中厚板自由振动的有限元分析

基于物理中面和一阶剪切变形板理论,研究了不同边界条件下功能梯度材料(FGM)中厚板的自由振动问题．假设功能梯度板的材料性质沿厚度方向按幂函数规律连续变化．根据哈密顿原理建立了FGM板有限元形式的自由振动方程,利用MATLAB软件编写程序进行了计算．通过数值算例,讨论了不同边界条件下FGM中厚板的无量纲频率随材料梯度指数和厚宽比的变化情况,并与经典板理论下的频率进行了比较．