Analysis of 2D non-axisymmetric elasticity and thermoelasticity problems for radially inhomogeneous hollow cylinders

An analytical approach to solve plane static non-axisymmetric elasticity and thermoelasticity problems for radially inhomogeneous hollow cylinders is presented. This approach is based upon the direct integration method proposed by Vihak (Vigak). The essence of the method mentioned is in the integration of the original differential equilibrium equations, which are independent of the stress–strain relations. This gives the opportunity to deduce the relations, which are invariant with respect to various properties of the material, for the stress-tensor components. From these relations each of the stress-tensor components have been expressed in terms of the governing one. A solution of the equation for the governing stress in the form of Fourier series is presented. To determine the Fourier coefficients, an integral Volterra-type equation is derived and solved by a simple iteration method with rapid convergence. Other stress-tensor components are expressed through the obtained governing stress in the form of an explicit functional dependence on force and thermal loadings.     

Green's matrices for the 2-D Lame system in regions of irregular shape

A method is proposed for the construction of Green's matrices for mixed boundary value problems in regions of irregular shape for the displacement formulation of the plane problem in theory of elasticity. The method is based on the boundary integral equation approach where a kernel matrix B satisfies the 2-D homogeneous Lame system inside the region. This leads to a regular boundary integral equation where the compensating load is applied to the boundary. The Green's matrix is consequently expressed in terms of the kernel matrix B, the fundamental solution matrix of the homogeneous Lame system and a kernel matrix of the inverse regular integral operator. To calculate the stress components, the kernel matrices are differentiated under the integral sign. The proposed method appears highly effective in computing both displacements and stresses.     

HYBRID-TREFFTZ FINITE ELEMENT FORMULATION FOR SIMULATION OF SINGULAR STRESS FIELDS

An equilibrium hybrid-Trefftz formulation based on the direct approximation of the stress and boundary displacement fields is presented. The general solution of the governing differential equations is used to approximate the stress field and the boundary displacements are represented by polynomial functions. When singular solutions are implemented to model local high stress gradients due to concentrated loads or to the presence of wedges or cracks, rational functions are used to approximate the boundary displacements in the neighbourhood of such singular stress points. The equilibrium conditions and the kinematic boundary conditions are locally satisfied. The remaining fundamental relations—the compatibility conditions, the static boundary conditions and the constitutive relations—are enforced in a weighted residual form so designed as to preserve the duality and constitutive reciprocity. The resulting governing system is symmetric and all intervening structural operators have boundary integral expressions. Numerical applications are presented to illustrate the performance of the formulation.     

压电热弹性材料四边简支层合板的精确解

根据压电热弹性材料的控制方程和热传导关系,重构压电热弹性材料的本构关系,通过新本构关系并结合压电热弹性材料热平衡方程,得出压电热弹性材料机-电-热耦合问题的齐次状态方程。应用精细积分法,状态方程可独立求解。此方法在分析压电热弹性体耦合问题时,避免了求解关于热传导方程和热平衡方程的二阶微分方程,大大减少了数值计算的工作量。     

Green's functions for Kirchhoff plates of irregular shape

A method is proposed for the construction of Green's functions for the Sophie Germain equation in regions of irregular shape with mixed boundary conditions imposed. The method is based on the boundary integral equation approach where a kernel vector function B satisfies the biharmonic equation inside the region. This leads to a regular boundary integral equation where the compensating loads and moments are applied to the boundary. Green's function is consequently expressed in terms of the kernel vector function B, the fundamental solution function of the biharmonic equation, and kernel functions of the inverse regular integral operators. To compute moments and forces, the kernel functions are differentiated under the integral sign. The proposed method appears highly effective in computing both displacements and stress components.     

Penny-shaped crack in a sphere embedded in an infinite medium

We consider the problem of determining the stress intensity factor and the crack energy in an Isotropie, homogeneous elastic sphere embedded in an infinite Isotropie, homogeneous elastic medium when there is a diametrical crack in the sphere. We assume that the crack is opened by an internal pressure and the sphere is bonded to the surrounding material. The problem is reduced to the solution of a Fredholm integral equation of the second kind in the auxiliary function φ($\text{t}$). Expressions for the stress intensity factor and the crack energy are obtained in terms of φ($\text{t}$). The integral equation is solved numerically and the numerical values of the stress intensity factor and the crack energy are graphed.     

STRESS- AND ROTATION-BASED HIERARCHIC MODELS FOR LAMINATED COMPOSITES

A hierarchic sequence of equilibrium models in terms of stresses assumed to be not a priori symmetric is derived for cylindrical bending of laminated composites, using first-order stress functions. The stress field of each hierarchic model satisfies a priori (i) the translational equilibrium equations and the stress boundary conditions of two-dimensional elasticity, and (ii) the continuity requirement for the transverse shear and normal stresses at the lamina interfaces. The levels of hierarchy correspond to the degree to which the two first-order compatibility equations and the rotational equilibrium equation of two-dimensional elasticity are satisfied. The numerical solution is based on Fraeijs de Veubeke's dual mixed variational principle, employing the p-version of the finite element method. The number of degrees of freedom is independent of the number of the layers in the laminate. Results are obtained directly for the stresses and rotations; the displacement field is obtained in the post-processing phase by integration. Numerical results with comparisons show the capability of the mathematical and numerical models proposed.     

On integral equation approaches to the mechanics of fibre-crack interaction

The paper examines the problem of a penny-shaped crack which is formed by the development of a crack in both the fibre and the matrix of a composite consisting of an isolated elastic fibre located in an elastic matrix of infinite extent. The composite region is subjected to a uniform strain field in the direction of the fibre. The paper presents two integral-equation based approaches for the analysis of the problem. The first approach considers the formulation of the complete integral equations governing the associated elasticity problem for a two material region. The second approach considers the boundary integral equation formulation of the problem. Both methods entail the numerical solution of the governing integral equations. The solutions to these integral equations are used to evaluate the stress intensity factor at the boundary of the penny-shaped crack.     

Investigation of singularity orders and eigen-angular functions for V-notches in radially inhomogeneous materials

The singularities of V-notches in a material whose modulus of elasticity varies as a power function of the r-coordinate are investigated. The eigen-equations are derived by substituting the asymptotic expansions of displacement fields around a notch vertex into the elasticity equilibrium equations. The radial boundary conditions are also presented by a combination of stress singularity orders and eigen-angular functions. The singularity orders and eigen-angular functions can be calculated simultaneously by solving a set of eigen ordinary differential equations with variable coefficients, which are solved by the interpolating matrix method developed by some of the authors before without the iterative process. The accuracy of the proposed method is verified by comparing the present results with the reference ones when the inhomogeneous material is degraded into a homogeneous one. Then, the stress singularities of V-notches in the radially inhomogeneous material under the plane and anti-plane loadings are investigated, respectively. The results show that the stress singularity of a V-notch in the radially inhomogeneous material under the plane loading is more serious than the one under the anti-plane loading. The plane V-notch under the clamped–clamped boundary condition presents the stress singularity at a smaller notch angle α than the one under the free–free boundary condition. The radially inhomogeneous bulgy V-notch even presents singularity. In addition, the stress singularity becomes stronger with the increase of the exponent c in the variation function of the elasticity modulus.     

On the thermo-elastostatics of heterogeneous materials: I. General integral equation

We consider a linearly thermoelastic composite medium, which consists of a homogeneous matrix containing a statistically inhomogeneous random set of inclusions, when the concentration of the inclusions is a function of the coordinates (so-called functionally graded materials). The composite medium is subjected to essentially inhomogeneous loading by the fields of the stresses, temperature, and body forces (e.g. for a centrifugal load). The general integral equations connecting the stress and strain fields in the point being considered and the surrounding points are obtained for the random fields of inclusions. The method is based on a centering procedure of subtraction from both sides of a known initial integral equation their statistical averages obtained without any auxiliary assumptions, such as effective field hypothesis implicitly exploited in the known centering methods. In so doing the size of a region including the inclusions acting on a separate one is finite, i.e. the locality principle takes place.     

An internal-variable rod model for stress-induced phase transitions in a slender SMA layer. I: Asymptotic equations and a two-phase solution

In part I of this series of two papers, an internal-variable rod model is proposed to study the stress-induced phase transitions in a slender shape memory alloy (SMA) layer. To study the mechanical responses of SMAs, two independent energy functions are adopted: the Helmholtz free energy and the rate of mechanical dissipation. A phase state variable is introduced to describe the phase transition process. Starting from the 2-D governing system and by using the coupled series-asymptotic expansion method, one single equilibrium equation is derived, which involves the leading order term of the axial strain and the phase state functions. Further by using the phase transition criteria, the evolution laws of the phase state functions corresponding to the outer loop of the stress–strain response are derived. As a result, the governing ODEs for the purely loading and purely unloading processes are obtained, which are called the asymptotic rod equations. The two-phase solution of the asymptotic rod equations in an infinitely long layer is then constructed. An explicit solution for the phase volume fraction of the corresponding inhomogeneous deformation is deduced, which appears to be the first analytical expression for this important quantity in stress-induced phase transitions. The key parameters in terms of original material constants for the phase transitions and deformations are also identified.     

New Method for the Solution of Dynamic Problems of the Theory of Elasticity and Fracture Mechanics

We propose a new approach to the solution of dynamic problems of the theory of elasticity and fracture mechanics based on the application of the finite-difference method only with respect to time. In this case, the equations of motion are split into homogeneous and inhomogeneous systems of differential equations for the determination of displacements at time nodes. For trivial initial conditions, only the homogeneous system of differential equations is preserved (of the same type as in the dynamic problem in Laplace transforms). Its efficient solution can be obtained by the methods of boundary integral equations or boundary elements.     

Large deflection analysis of beams with variable stiffness

Summary. In this paper, the Analog Equation Method (AEM), a BEM-based method, is employed to the nonlinear analysis of a Bernoulli-Euler beam with variable stiffness undergoing large deflections, under general boundary conditions which maybe nonlinear. As the cross-sectional properties of the beam vary along its axis, the coefficients of the differential equations governing the equilibrium of the beam are variable. The formulation is in terms of the displacements. The governing equations are derived in both deformed and undeformed configuration and the deviations of the two approaches are studied. Using the concept of the analog equation, the two coupled nonlinear differential equations with variable coefficients are replaced by two uncoupled linear ones pertaining to the axial and transverse deformation of a substitute beam with unit axial and bending stiffness, respectively, under fictitious load distributions. Besides the effectiveness and accuracy of the developed method, a significant advantage is that the displacements as well as the stress resultants are computed at any cross-section of the beam using the respective integral representations as mathematical formulae. Several beams are analyzed under various boundary conditions and loadings to illustrate the merits of the method as well as its applicability, efficiency and accuracy.     

Multiple 3D inhomogeneous inclusions in a half space under contact loading

A semi-analytic solution is given for multiple three-dimensional inhomogeneous inclusions of arbitrary shape in an isotropic half space under contact loading. The solution takes into account interactions between all the inhomogeneous inclusions as well as the interaction between the inhomogeneous inclusions and the loading indenter. In formulating the governing equations for the inhomogeneous inclusion problem, the inhomogeneous inclusions are treated as homogenous inclusions with initial eigenstrains plus unknown equivalent eigenstrains, according to Eshelby’s equivalent inclusion method. Such a treatment converts the original contact problem concerning an inhomogeneous half space into a homogeneous half-space contact problem, for which governing equations with unknown contact load distribution can be conveniently formulated. All the governing equations are solved iteratively using the Conjugate Gradient Method. The iterative process is performed until the convergence of the half-space surface displacements, which are the sum of the displacements due to the contact load and the inhomogeneous inclusions, is achieved. Finally, the obtained solution is applied to two example cases: a single inhomogeneity in a half space subjected to indentation and a stringer of inhomogeneities in an indented half-space. The validation of the solution is done by modeling a layer of film as an inhomogeneity and comparing the present solution with the analytic solution for elastic indentation of thin films. This general solution is expected to have wide applications in addressing engineering problems concerning inelastic deformation and material dissimilarity as well as contact loading.     

Doubly periodic volume?surface integral equation formulation for modelling metamaterials

A generalised volume-surface integral equation is extended by way of the periodic Green's function to model arbitrarily complex designs of metamaterials consisting of high-contrast inhomogeneous anisotropic material regions as well as metallic inclusions. The unique aspect of the formulation is the integration of boundary and volume integral equations to increase modelling efficiency and capability. Specifically, the boundary integral approach with equivalent surface currents is adopted over regions consisting of piecewise homogeneous materials as well as metallic perfect electric/magnetic conductor inclusions, whereas the volume integral equation is employed only in inhomogeneous and/or anisotropic material regions. Because the periodic Green's function only needs to be evaluated for the equivalent surface currents enclosing an inhomogeneous and/or anisotropic region, matrix fill time is much less as compared to using a volume formulation. Furthermore, the incorporation of curvilinear finite elements allows for greater geometrical modelling flexibility for arbitrarily shaped high-contrast regions found in typical designs of engineered metamaterials     

A magneto-electro-elastic material with a penny-shaped crack subjected to temperature loading

Summary The analysis of intensity factors for a penny-shaped crack under thermal, mechanical, electrical and magnetic boundary conditions becomes a very important topic in fracture mechanics. An exact solution is derived for the problem of a penny-shaped crack in a magneto-electro-thermo-elastic material in a temperature field. The problem is analyzed within the framework of the theory of linear magneto-electro-thermo-elasticity. The coupling features of transversely isotropic magneto-electro-thermo-elastic solids are governed by a system of partial differential equations with respect to the elastic displacements, the electric potential, the magnetic potential and the temperature field. The heat conduction equation and equilibrium equations for an infinite magneto-electro-thermo-elastic media are solved by means of the Hankel integral transform. The mathematical formulations for the crack conditions are derived as a set of dual integral equations, which, in turn, are reduced to Abel's integral equation. Solution of Abel's integral equation is applied to derive the elastic, electric and magnetic fields as well as field intensity factors. The intensity factors of thermal stress, electric displacement and magnetic induction are derived explicitly for approximate (impermeable or permeable) and exact (a notch of finite thickness crack) conditions. Due to its explicitness, the solution is remarkable and should be of great interest in the magneto-electro-thermo-elastic material analysis and design.     

离散多层圆筒在热冲击载荷下的弹性动力响应

离散多层圆筒由薄内筒和倾角错绕的钢带层组成,具有制造简便、成本低等优点。预测筒体在热冲击载荷下的热应力对强度设计和安全操作具有重要的应用价值。该文首次研究了离散多层圆筒在热冲击载荷作用下的热弹性动态响应。将内筒和钢带层的径向位移分别分解为满足给定应力边界条件的准静态解和满足初始条件的动态解,准静态解通过齐次线性方法确定,热弹性动态解通过有限Hankel积分变换和Laplace变换确定。根据内外层界面处位移连续条件,得到层间压力关于时间的第二类Volterra积分方程,利用Hermit二次三项式插值方法可求得该层间应力。最后将离散多层圆筒的热弹性动力响应与单层厚壁圆筒的响应进行了比较,并分析了钢带缠绕倾角和材料参数对热弹性动力响应的影响。     

The singular function boundary integral method for a two-dimensional fracture problem

The singular function boundary integral method (SFBIM) originally developed for Laplacian problems with boundary singularities is extended for solving two-dimensional fracture problems formulated in terms of the Airy stress function. Our goal is the accurate, direct computation of the associated stress intensity factors, which appear as coefficients in the asymptotic expansion of the solution near the crack tip. In the SFBIM, the leading terms of the asymptotic solution are used to approximate the solution and to weight the governing biharmonic equation in the Galerkin sense. The discretized equations are reduced to boundary integrals by means of Green's theorem and the Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers. The numerical results on a model problem show that the method converges extremely fast and yields accurate estimates of the leading stress intensity factors.     

Three-dimensional analytical solution for functionally graded magneto–electro-elastic circular plates subjected to uniform load

The problem of a functionally graded, transversely isotropic, magneto–electro-elastic circular plate acted on by a uniform load is considered. The displacements and electric potential are represented by appropriate polynomials in the radial coordinate, of which the coefficients depends on the thickness coordinate, and are called the generalized displacement functions. The governing equations as well as the boundary conditions for these generalized displacement functions are derived from the original equations of equilibrium for axisymmetric problems and the boundary conditions on the upper and lower surfaces of the plate. Explicit expressions are then obtained through a step-by-step integration scheme, with five integral constants determinable from the boundary conditions at the cylindrical surface in the Saint Venant’s sense. The analytical solution is suited to arbitrary variations of material properties along the thickness direction, and can be readily degenerated into those for homogeneous plates. A particular circular plate, with some material constants being the exponential functions of the thickness coordinate, is finally considered for illustration.     

On the dynamic crack propagation in an interface with spatially varying elastic properties

This article provides a comprehensive theoretical treatment of a finite crack propagating in an interfacial layer with spatially varying elastic properties under antiplane loading condition. The theoretical formulations governing the steady state solution are based upon the use of an integral transform technique. The resulting dynamic stress intensity factor of the propagating cracks is obtained by solving the appropriate singular integral equations, using Chebyshev polynomials, for different inhomogeneous materials. Numerical examples are provided to verify the technique and to show the effect of the thickness of the interfacial layer and the material properties upon the dynamic stress intensity factor of the crack and the associated singularity transition.