Effect of material non-homogeneity on the inhomogeneous shearing deformation of a Gent slab subjected to a temperature gradient

It is well known that most rubber-like materials are non-homogeneous due to either imperfect manufacturing conditions or the action of severe thermo-oxidative environments in many practical applications. In this study, within the context of finite thermoelasticity, we theoretically analyze the inhomogeneous shearing deformation of a non-homogeneous rubber-like slab subjected to a thermal gradient across its thickness. The major objective of this study is to investigate the effect of the material non-homogeneity, which is the material-coordinate dependence of the material response functions, on the stress-strain fields for a given temperature gradient. First, we show the existence of a simple shearing deformation from which the generalized shear modulus and the generalized thermal conductivity of the slab could be obtained. Based on this information, the Gent material model is generalized to take the material non-homogeneity and the temperature dependence of the stress into account. To analyze the inhomogeneous shearing deformation of the non-homogeneous slab, deformation and temperature fields are postulated; then the decoupled temperature field is obtained analytically by solving the local energy balance equation. Finally, the static equilibrium equations are solved considering the linear temperature field. Our results show that the spatial pattern and the degree of the material non-homogeneity have profound effects on the stress-strain fields. The shear strain becomes nearly homogeneous and the stresses are relatively small for a certain spatial variation of the material non-homogeneity. This result suggests the possibility of designing a novel class of materials: functionally graded rubber-elastic materials (FGREMs).     

非均匀弹性材料中反平面的运动裂纹

用解析方法研究了非均匀弹性材料中反平面运动裂纹问题。首先采用余弦变换求解非均匀材料的基本方程，然后根据混合边值条件建立裂纹运动的对偶积分方程，再把对偶积分方程化为第二类Fredholm积分方程。给出了数值算例，计算结果表明材料的非均匀性对动应力强度因子有较大的影响。     

Effect of thermal gradient on vibration of non-homogeneous visco-elastic elliptic plate of variable thickness

An analysis and numerical results are presented for free transverse vibrations of non-homogeneous visco-elastic elliptic plate whose temperature and thickness spatial variations both are parabolic along a line through plate centre. The variation in density is assumed as parabolic along a line through plate centre, which raise non-homogeneity of the plate materials and make problem interesting as introducing variation in non-homogeneity of the material mass density reduce the problem practical importance in comparison to homogenous plate. For visco-elastic, the basic elastic and viscous elements are combined. We have taken Kelvin model for visco-elasticity that is the combination of the elastic and viscous elements in parallel. Here the elastic element means the spring and the viscous element means the dashpot. The governing differential equation of motion has been solved by Galerkin’s technique. Deflection, time period and logarithmic decrement corresponding to the first two modes of vibrations of a clamped non-homogeneous visco-elastic elliptic plate for various values of taper constant, thermal constants, non-homogeneity constant and aspect ratio are obtained and shown graphically.     

Corner stress singularities in an FGM thin plate

Describing the behaviors of stress singularities correctly is essential for obtaining accurate numerical solutions of complicated problems with stress singularities. This analysis derives asymptotic solutions for functionally graded material (FGM) thin plates with geometrically induced stress singularities. The classical thin plate theory is used to establish the equilibrium equations for FGM thin plates. It is assumed that the Young’s modulus varies along the thickness and Poisson’s ratio is constant. The eigenfunction expansion method is employed to the equilibrium equations in terms of displacement components for an asymptotic analysis in the vicinity of a sharp corner. The characteristic equations for determining the stress singularity order at the corner vertex and the corresponding corner functions are explicitly given for different combinations of boundary conditions along the radial edges forming the sharp corner. The non-homogeneous elasticity properties are present only in the characteristic equations corresponding to boundary conditions involving simple support. Finally, the effects of material non-homogeneity following a power law on the stress singularity orders are thoroughly examined by showing the minimum real values of the roots of the characteristic equations varying with the material properties and vertex angle.     

The bending of elastic circular ring of non-homogeneous and variable cross section under the actions of arbitrary loads

On the basis of the stepped reduction method suggested in [1], we investigate the problem of the bending of elastic circular ring of non-homogeneous and variable cross section under the actions of arbitrary loads. The general solution of this problem is obtained so that it can be used for the calculations of strength and rigidity of practical problems such as arch, tunnel etc. In order to examine results of this paper and explain the application of this new method, an example is brought out at the end of this paper. Circular ring and arch are commonly used structures in engineering. Timoshenko, S.^[2], Barber, J. R.^[3], Tsumura Rimitsu^[4] et al. have studied these problems of bending, but, so far as we know, it has been solely restricted to the general solution of homogeneous uniform cross section ring. The only known solution for the problems with variable cross section ones has been solely restricted to the solution of special case of flexural rigidity in linear function of coordinates. On account of fundamental equations of the non-homogeneous variable cross section problem being variable coefficients, it is very difficult to solve them. In this paper, we use the stepped reduction method suggested in [1] to transform the variable coefficient differential equation into equivalent constant coefficient one. After introducing virtual internal forces, we obtain general solution of an elastic circular ring with non-homogeneity and variable cross section under the actions of arbitrary loads.     

Finite Griffith crack propagating in a non-homogeneous medium

The problem of a Griffith crack of constant length propagating at a uniform speed in a plane non-homogeneous medium under uniform load is investigated. The equilibrium equations for the non-homogeneous medium are solved by using the Fourier transforms and then the problem is reduced to the solution of dual integral equations. Solving the dual integral equations we obtain the expression for the dynamic stress intensity factor at the edge of the crack. Finally the numerical results for the stress intensity factor are obtained which are displayed graphically to show the effect of the material non-homogeneity on the stress intensity factor.     

Dynamic fracture mechanics analysis for composite material with material non-homogeneity in thickness direction

The problem considered here is the response of a non-homogeneous composite material containing some cracks subjected to dynamic loading. It is assumed that the composite material is orthotropic and all the material properties depend only on the coordinatey (along the thickness direction). In the analysis, the elastic region is divided into a number of plies of infinite length. The material properties are taken to be constants for each ply. By utilizing the Laplace transform and Fourier transform technique, the general solutions for plies are derived. The singular integral equations of the entire elastic region are obtained and solved by the virtual displacement principle. Attention is focused on the time-dependent full field solutions of stress intensity factor(SIF) and strain energy release rate. As a numerical illustration, the dynamic stress intensity factor of a substrate/functionally graded film structure with two cracks under suddenly applied forces on cracks face are presented for various material non-homogeneity parameters.     

The non-homogeneous biharmonic plate equation: fundamental solutions

Real materials and structural components are often non-homogeneous, either by design or because of the physical composition and imperfections in the underlying material. Thus, analytical solutions for non-homogeneous materials under mechanical loads are of considerable interest to engineers and have widespread applications, given the prevalence of these materials in fields as diverse as aerospace, construction, electronics, etc. More precisely, those are essentially composites with carefully manufactured properties that yield desirable mechanical characteristics and properties, such as optimal arrangement of the material, minimum weight, etc. To this end, the displacement fundamental solution (or Green’s function) corresponding to a point force for the non-homogenous biharmonic equation in two dimensions are derived in this work by employing a conformal mapping technique in conjunction with the Radon transformation. These functions, besides being useful in their own right, can also be used within the context of integral equation formulations for the solution of boundary-value problems. Finally, a series of numerical examples that deal with the non-homogeneous plate on elastic foundation problem serve to illustrate the present method.     

Thermal response of a partially insulated interface crack in a graded coating-substrate structure under thermo-mechanical disturbance: Energy release and density

Consider the thermal fracture problem of a functionally graded coating-substrate structure of finite thickness with a partially insulated interface crack subjected to thermal-mechanical supply. A new model is proposed that the heat conduction through the crack region occurs and the temperature drop across the crack surfaces is the result of the thermal resistance. For the first time, real fundamental solutions are derived for the fracture analysis of functionally graded materials. The complicated mixed boundary problems of equations of heat conduction and elasticity are converted analytically into singular integral equations, which are solved numerically. The asymptotic expressions with higher order terms for the singular integral kernels are considered to improve the accuracy and efficiency of the numerical integration. Explicit expressions of various failure modes including stress intensity factors, energy release rate and strain energy density, are provided. Numerical results are presented to illustrate the effects of non-homogeneity parameters and the dimensionless thermal resistance on the temperature distribution along the crack surfaces and extended crack line, the thermal stress intensity factors and minimum strain energy density.     

Dynamic stress intensity factors for homogeneous and smoothly heterogeneous materials using the interaction integral method

Dynamic stress intensity factors (DSIFs) are important fracture parameters in understanding and predicting dynamic fracture behavior of a cracked body. To evaluate DSIFs for both homogeneous and non-homogeneous materials, the interaction integral (conservation integral) originally proposed to evaluate SIFs for a static homogeneous medium is extended to incorporate dynamic effects and material non-homogeneity, and is implemented in conjunction with the finite element method (FEM). The technique is implemented and verified using benchmark problems. Then, various homogeneous and non-homogeneous cracked bodies under dynamic loading are employed to investigate dynamic fracture behavior such as the variation of DSIFs for different material property profiles, the relation between initiation time and the domain size (for integral evaluation), and the contribution of each distinct term in the interaction integral.     

Convective cooling/heating induced thermal stresses in a fluid saturated porous medium undergoing local thermal non-equilibrium

Local thermal non-equilibrium (LTNE) may have profound effects on the pore pressure and thermal stresses in fluid saturated porous media under transient thermal loads. This work investigates the temperature, pore pressure, and thermal stress distributions in a porous medium subjected to convective cooling/heating on its boundary. The LTNE thermo-poroelasticity equations are solved by means of Laplace transform for two fundamental problems in petroleum engineering and nuclear waste storage applications, i.e., an infinite porous medium containing a cylindrical hole or a spherical cavity subjected to symmetrical thermo-mechanical loads on the cavity boundary. Numerical examples are presented to examine the effects of LTNE under convective cooling/heating conditions on the temperature, pore pressure and thermal stresses around the cavities. The results show that the LTNE effects become more pronounced when the convective heat transfer boundary conditions are employed. For the cylindrical hole problem of a sandstone formation, the thermally induced pore pressure and the magnitude of thermal stresses are significantly higher than the corresponding values in the classical poroelasticity, which is particularly true under convective cooling with moderate Biot numbers. For the spherical cavity problem of a clay medium, the LTNE effect may become significant depending on the boundary conditions employed in the classical theory.     

热对流条件下考虑球形夹杂分布的材料热弹接触摩擦热影响分析

机械传动关键活动零部件接触副往往受到力载荷和摩擦热载荷的耦合作用,使得接触界面间的接触力学行为的分析变得极其复杂. 利用基于等效夹杂方法建立的考虑热对流非均质材料热弹接触力学分析模型研究不同摩擦系数、夹杂位置和材料属性等参数对材料表面及内部温升及热应力分布影响规律. 此外,进一步分析了接触副材料中含分布球形夹杂时摩擦热造成的影响. 结果表明:接触副表面温升梯度受热对流系数的影响较大;下表面温升和热应力随摩擦系数增大而增大;分布夹杂则将接触副材料下表面温升及热应力分布变得更为复杂.      

Nonlinear delamination buckling and expansion of functionally graded laminated piezoelectric composite shells

In this paper, an analytical method is presented to investigate the nonlinear buckling and expansion behaviors of local delaminations near the surface of functionally graded laminated piezoelectric composite shells subjected to the thermal, electrical and mechanical loads, where the mid-plane nonlinear geometrical relation of delaminations is considered. In examples, the effects of thermal loading, electric field strength, the stacking patterns of functionally graded laminated piezoelectric composite shells and the patterns of delaminations on the critical axial loading of locally delaminated buckling are described and discussed. Finally, the possible growth directions of local buckling for delaminated sub-shells are described by calculating the expanding forces along the length and short axis of the delaminated sub-shells.     

Interface crack problems for mode II of double dissimilar orthotropic composite materials

The fracture problems near the similar orthotropic composite materials are interface crack tip for mode Ⅱ of double disstudied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized hi-harmonic equations, the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions, a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about bimaterial engineering parameters. According to the uniqueness theorem of limit, both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same, the stress singularity exponents, stress intensity factors and stresses for mode II crack of the orthotropic single material are obtained.     

The stress state of a plate subjected to a piecewise homogeneous load

We consider the stress-strain state of a plate having a doubly connected domain S bounded from the outside by a circle of radius R and from the inside by an ellipse with two rectilinear cuts. The cuts lie symmetrically on the x-axis. The plate is subjected to various forces: the hole contour (the ellipse) is under the action of uniformly distributed forces of intensity q, and the cut shores are free of loads; at the points ±ib of the imaginary axis, the plate is under the action of a lumped force P.The solution of the problem is reduced to determining two analytic functions φ(z) and ψ(z) satisfying certain boundary conditions (depending on the type of the acting loads).We use the Kolosov-Muskhelishvili method to reduce the problem to a system of linear algebraic equations for the coefficients in the expansions of the functions φ(z) and ψ(z). The solution thus obtained is illustrated by numerical examples.     

High order behavior of sandwich plates with free edges––edge effects

The edge effects of a sandwich plate with a “soft” core and free edges, i.e. the plate is supported only at the lower face-sheet, and the upper face-sheet and the core are free of stresses at their edges, using the high order approach (HSAPT), are presented. The two-dimensional analysis consists of a mathematical formulation that uses the classical thin plate theory for the face-sheets and a three-dimensional elasticity theory for the core. The governing equations and the required boundary conditions are derived explicitly through variational principals, yielding a system of eight partial differential equations. The non-homogeneous differential equations system is numerically solved using a modification of the extended Kantorovich method (MEKM). The model presented enables a two-dimensional solution of the stress and displacement fields when subjected to a general scheme of loads. It is applicable to any type of boundary conditions that can be applied separately on each face-sheet and on the core. A numerical study is presented, and it examines the behavior and the two-dimensional stress field of a sandwich plate with free edges, at the upper face-sheet and core, subjected to thermal and uniformly distributed loads, for various boundary conditions at the lower face-sheet. For completeness, the MEKM solution of the two-dimensional high order model is verified through comparison with a three-dimensional Finite Element model revealing good correlation. Furthermore, the problems involved in the construction of an appropriate three-dimensional FE model of a full scale sandwich plate that require large computer resources are discussed.The numerical study yields that the peeling (normal) stresses, which reach their maximum values at the edges of the sandwich plate, using a one-dimensional analysis, varies also in the transverse direction from a maximum value in the middle of the edge, descending towards the corners. Moreover, the nature of variation along the boundaries strongly depends on the type of loading and the transverse boundary conditions. The substantial variation of the stress field in the transverse direction clearly shows the necessity of a two-dimensional analysis and the inefficiencies of the one-dimensional model.     

Bending analysis of a ceramic-metal arched bridge using a mixed first-order theory

In this research, the bending analysis of an arched bridge is presented based on a mixed first-order thick beam one dimensional plate theory. The present arched bridge is considered as a beam with boundary conditions at its edges, which may be simply-supported, and between these two edges, the beam may have quadratic thickness variation. The bridge consists of two layers; the upper flat one is made from an isotropic homogeneous material such as ceramic, and the lower arched layer is made from an isotropic non-homogeneous functionally graded ceramic-metal material. The upper-surface of the arched layer, which represents the interface between the two layers, is ceramic-rich material while the lower-surface of the arched layer is metal-rich material. This structure eliminates interface problem of the arched bridge and thus the stress distributions are smooth. A closed form solution is developed for the static response of such bridge subjected to different distributed loads. The effects of many parameters on the displacements and stresses are investigated. The sample numerical examples presented herein for bending response of the present arched bridge should serve as references for future comparisons.     

Mixed-mode fracture analysis of orthotropic functionally graded materials under mechanical and thermal loads

Mixed-mode fracture problems of orthotropic functionally graded materials (FGMs) are examined under mechanical and thermal loading conditions. In the case of mechanical loading, an embedded crack in an orthotropic FGM layer is considered. The crack is assumed to be loaded by arbitrary normal and shear tractions that are applied to its surfaces. An analytical solution based on the singular integral equations and a numerical approach based on the enriched finite elements are developed to evaluate the mixed-mode stress intensity factors and the energy release rate under the given mechanical loading conditions. The use of this dual approach methodology allowed the verifications of both methods leading to a highly accurate numerical predictive capability to assess the effects of material orthotropy and nonhomogeneity constants on the crack tip parameters. In the case of thermal loading, the response of periodic cracks in an orthotropic FGM layer subjected to transient thermal stresses is examined by means of the developed enriched finite element method. The results presented for the thermally loaded layer illustrate the influences of the material property gradation profiles and crack periodicity on the transient fracture mechanics parameters.     

Mixed Boundary Value Problem in a Non-homogeneous Medium Under Steady Distribution of Temperature

This paper is concerned with a mixed boundary value problem of a non-homogeneous medium under steady distribution of temperature whose elastic constants are exponential functions of y. By using Fourier cosine transforms the mixed boundary value problem of heat conduction is reduced to a Fredholm integral equation of the second kind. Then the elastic problem of the non-homogeneous semi-infinite half-plane under distribution of load over a plane face is solved. The influence of the non-homogeneity of the material on the thermal stress distribution is illustrated graphically.