Three-dimensional Elastic Body Research Articles

Modeling a triclinic lattice elastic body based on the linear couple stress theory

Triclinic lattice structures with generalized parallelepiped unit cells are crucial targets in understanding the influence of lattice structure on the mechanical behaviors of lattice elastic bodies. This study models a three-dimensional elastic body with a triclinic lattice microstructure comprising three uniform elastic members rigidly joined at a single point. Based on the linear couple stress theory, a specialized case of the Cosserat theory, this study derives elasticity tensors in the constitutive equations of the triclinic lattice elastic body that depend on the crossing angles between the members. To derive the angular-dependent elasticity tensors, coordinate transformations between oblique and orthogonal coordinate systems are effectively employed in formulating the substitution of deformation and force fields in the continuum for the fields in the lattice structure. Evaluating the elasticity tensors reveals the influences of the crossing angles on the anisotropic tensile and torsional properties, the existence of the angle-covariant lattice number manipulation that cancels out the size effect, and several types of geometrical shape dependence. This study could contribute to the theoretical evaluation of elastic behaviors of both artificially engineered and naturally formed materials with lattice microstructures.

Wave propagation in three-dimensional fractional viscoelastic infinite solid body

In order to analyze the wave propagation in three-dimensional isotropic and viscoelastic body, the Cauchy initial value problem on unbounded domain is considered for the wave equation written as a system of fractional partial differential equations consisting of equation of motion of three-dimensional solid body, equation of strain, as well as of the constitutive equation, obtained by generalizing the classical Hooke’s law of three-dimensional isotropic and elastic body by replacing Lamé coefficients with the relaxation moduli to account for different memory kernels corresponding to the propagation of compressive and shear waves. By the use of the method of integral transforms, namely Laplace and Fourier transforms, the displacement field, as a solution of the Cauchy initial value problem, is expressed through Green’s functions corresponding to both compressive and shear wave propagation. Two approaches are adopted in solving the Cauchy problem: in the first one the displacement field is expressed through the scalar and vector fields, obtained as solutions of wave equations which are consequences of decoupling the wave equation for displacement field, while in the second approach the displacement field is obtained by the action of the resolvent tensor on the initial conditions. Fractional anti-Zener and Zener model I+ID.ID, as well as fractional Burgers model VII, as representatives of models yielding infinite and finite wave propagation speed, are used in numerical calculations and graphical representations of Green’s functions and its short and large time asymptotic behavior.

Finite element study of V-shaped crack-tip fields in a three-dimensional nonlinear strain-limiting elastic body

We study the crack-tip fields in a three-dimensional (3D) plate with a V-shaped crack employing Rajagopal’s theory of elasticity with the strain-limiting effect. The classical linear elastic fracture mechanics (LEFM) theory bears a well-documented inconsistency, i.e., for a static crack problem under the traction-free boundary condition, the LEFM theory predicts singular crack-tip strains. Such results are inconsistent with the constitutive relation derived using the assumption of uniform bound for the strain values. Instead, Rajagopal’s recent theory of elasticity provides implicit constitutive relations between the linearized strain and Cauchy stress in which one can impose an a-priori upper bound for strains. We formulate the strain-limiting theory comprehensively in this article. The linearized strain is expressed as a nonlinear function of the Cauchy stress. The equation for the equilibrium of the linear momentum under a special type of nonlinear constitutive relation reduces to a second-order quasilinear elliptic boundary value problem (BVP). Using a Picard-type linearization and a continuous Galerkin-type finite element technique, we solve several BVPs with tensile and shear displacement loading. The numerical results for the strains, the stresses, the stress intensity factor (SIF), and strain energy density (SED) are obtained for different Poisson values, strain-limiting parametric values, and the V-shaped angles. The growth of the near-tip strains is far less than that of stresses, and the results for both the SIF and SED are also consistent with the results from LEFM. Our work demonstrates that the framework of Rajagopal’s theory of elasticity can provide a basis for developing physically meaningful and mathematically well-posed BVPs to study evolution of cracks, damage, nucleation, or failure in elastic bodies.

Analysis of near-interface cracks in three-dimensional anisotropic multi-materials by efficient BIEM

This paper presents the full analysis of near-interface cracks in three-dimensional linear elastic multimaterial bodies by a weakly singular boundary integral equation method (BIEM). The proposed technique was implemented in a general framework that allowed the treatment of finite cracked bodies with general configurations, material anisotropy, and mode-mixity. A system of integral equations governing unknown data on the boundary, crack surface, and material interface was established using a pair of weakly singular weak-form displacement and traction integral equations together with continuity along the material interface. A symmetric Galerkin boundary element method along with the finite element technique were implemented to solve the governing integral equations. In addition, the special near-front interpolation was employed to enhance the approximation of the relative crack-face displacement in the neighborhood of the crack front. The solved crack-face data were used directly to post-process the stress intensity factors and T-stresses along the crack front. An extensive numerical study was carried out not only to demonstrate the accuracy, convergence, and capability of the proposed technique, but also to explore the influence of material stiffness and distance to the material interface on fracture data along the crack front.

Finite element study of V-shaped crack-tip fields in a three-dimensional nonlinear strain-limiting elastic body

We study the crack-tip fields in a three-dimensional (3D) plate with a V-shaped crack employing Rajagopal’s theory of elasticity with the strain-limiting effect. The classical linear elastic fracture mechanics (LEFM) theory bears a well-documented inconsistency, i.e., for a static crack problem under the traction-free boundary condition, the LEFM theory predicts singular crack-tip strains. Such results are inconsistent with the constitutive relation derived using the assumption of uniform bound for the strain values. Instead, Rajagopal’s recent theory of elasticity provides implicit constitutive relations between the linearized strain and Cauchy stress in which one can impose an a-priori upper bound for strains. We formulate the strain-limiting theory comprehensively in this article. The linearized strain is expressed as a nonlinear function of the Cauchy stress. The equation for the equilibrium of the linear momentum under a special type of nonlinear constitutive relation reduces to a second-order quasilinear elliptic boundary value problem (BVP). Using a Picard-type linearization and a continuous Galerkin-type finite element technique, we solve several BVPs with tensile and shear displacement loading. The numerical results for the strains, the stresses, the stress intensity factor (SIF), and strain energy density (SED) are obtained for different Poisson values, strain-limiting parametric values, and the V-shaped angles. The growth of the near-tip strains is far less than that of stresses, and the results for both the SIF and SED are also consistent with the results from LEFM. Our work demonstrates that the framework of Rajagopal’s theory of elasticity can provide a basis for developing physically meaningful and mathematically well-posed BVPs to study evolution of cracks, damage, nucleation, or failure in elastic bodies.

An estimation of eigenvalues of three-dimensional elastic bodies on the basis of Physics-Informed Neural Networks

An estimation of eigenvalues of three-dimensional elastic bodies on the basis of Physics-Informed Neural Networks

Quadrilateral scaled boundary spectral shell elements for static and dynamic analyses of functionally graded plates and shells

This paper presents new shell elements for static and dynamic analyses of functionally graded plates and shells with the material properties varying through the thickness according to a power-law distribution in terms of the volume fractions of the constituents. In the formulation, a shell element is treated as a three-dimensional linear elastic body and its middle surface is represented with a quadrilateral spectral element. Along the thickness direction, the shell geometry is described by scaling the middle surface while the displacements are approximated by interpolating the displacements on the top, middle and bottom surfaces using quadratic Lagrange shape functions. The assumed natural strain method is applied to eliminate transverse shear locking, membrane locking and curvature thickness locking. The developed shell elements are featured by requiring only the shell mid-surface to be discretized, containing no rotational degree of freedom, adopting the full three-dimensional constitutive law and involving only the in-plane numerical integration. The validity and performance of the formulation are demonstrated through benchmark examples covering plates, cylindrical shells and spherical shells with varying span-to-thickness ratio and material gradation index. Static solutions of two functionally graded shell problems under gravity loads are also provided.

Recovering a bounded elastic body by electromagnetic far-field measurements

<p style='text-indent:20px;'>This paper is concerned with the scattering of a time-harmonic electromagnetic wave by a three-dimensional elastic body. The general transmission conditions are considered to model the interaction between the electromagnetic field and the elastic body on the interface by Voigt's model. The existence of a unique solution is first proved in an appropriate Sobolev space by employing the variational method with the classical Fredholm alternative. The inverse problem is then investigated to recover the elastic body by the scattered wave-field data. It is shown that the shape and location of the body is uniquely determined by the fixed energy magnetic (or electric) far-field measurements corresponding to incident plane waves with all polarizations.</p>

NUMERICAL SOLUTION OF THE DYNAMIC PROBLEM OF AXISYMMETRIC VIBRATIONS OF REINFORCED SHELLS

The reliability of the results obtained in the work is determined by the rigor and correctness of the statements of the initial problems; theoretical substantiation of the finite-difference schemes used; controlled accuracy of numerical calculations; conducting test calculations; compliance of the established regularities with the general properties of oscillations of thin-walled structural elements. The correctness of the formulation of the problems is achieved by using the well-known equations of the theory of shells and rods of the Tymoshenko type, which are approximations of the original equations of the three-dimensional theory of elasticity. When deriving the equations, the equations of oscillations of the multilayer shell in the smooth region and the equations of oscillations of reinforcing ribbed elements (transverse ribs) were obtained. It is not difficult to show that the indicated equations by the classification of equations in partial derivatives are equations of the hyperbolic type, which are an approximation of the oscillating equations of three-dimensional elastic bodies and sufficiently correctly reproduce wave processes in non-homogeneous shell structures, taking into account spatial gaps. Numerical algorithms for approximate solutions of the original equations are based on the use of the integro-interpolation method of constructing difference schemes. When constructing difference schemes, kinematic quantities refer to difference points with integer indices, and the values of deformations and moment forces refer to difference points with half-integer indices. This approximation of the initial kinematic and static values allows the fulfillment of the law of conservation of the total mechanical energy of the elastic structure at the difference level. The numerical algorithm is based on the use of separate finite-difference relations in the smooth domain and on the lines of spatial discontinuities with the second order of accuracy in spatial and temporal coordinates.

Micro-buckling of carbon fibers in shape memory polymer composites under bending in the glass transition temperature region

Abstract The influence of a wide temperature range in the glass transition region of a shape memory polymer (SMP) matrix on micro-buckling of the fiber reinforcements in shape memory polymer composites (SMPC) under large bending deformation is described. Analytical expressions to estimate the strain energy, neutral strain surface, critical buckling surface and half wavelength of the buckled fibers in the SMPC are presented based on the minimum energy method. This study considers the reinforced fibers as three-dimensional elastic bodies and the matrix as a temperature stimulated flat plate. A comprehensive study was performed to understand the dynamic temperature behavior of the micro-buckled fibers and corresponding results were validated by previous works in the literature. The effects of fiber volume fraction and thickness of the SMPC plates on the half wavelength are also discussed along with the simultaneous influence of temperature on the parameters computed in the minimum energy analysis.

Влияние распределенных дислокаций на устойчивость полого нелинейно-упругого шара

The phenomenon of stability loss of a hollow elastic sphere containing distributed dislocations and loaded with external hydrostatic pressure is studied. The study was carried out in the framework of the nonlinear elasticity theory and the continuum theory of continuously distributed dislocations. An exact statement and solution of the stability problem for a three-dimensional elastic body with distributed dislocations are given. The static problem of nonlinear elasticity theory for a body with distributed dislocations is reduced to a system of equations consisting of equilibrium equations, incompatibility equations with a given dislocation density tensor, and constitutive equations of the material. The unperturbed state is caused by external pressure and a spherically symmet-ric distribution of dislocations. For distributed edge dislocations in the framework of a harmonic (semi-linear) mate-rial model, the unperturbed state is defined as an exact spherically symmetric solution to a nonlinear boundary value problem. This solution is valid for any function that characterizes the density of edge dislocations. The perturbed equilibrium state is described by a boundary value problem linearized in the neighborhood of the equilibrium. The analysis of the axisymmetric buckling of the sphere was performed using the bifurcation method. It consists in determining the equilibrium positions of an elastic body, which differ little from the unperturbed state. By solving the linearized problem, the value of the external pressure at which the sphere first loses stability is found. The effect of dislocations on the buckling of thin and thick spherical shells is analyzed.

Stress intensity factors for branching cracks in space structures

One of the mechanisms of destruction of a solid is the formation and growth of “micro” and “macro” cracks. This task is especially acute for space objects. A microcrack in the body of the spacecraft practically does not affect its mechanical strength, but it can lead to depressurization of the body, which can lead to the loss of life support for the crew. Cracks can have different causes of their occurrence, but the mechanism of their growth in modern fracture mechanics is associated with the chosen criterion. For example, in the force criterion, crack growth begins when the stress intensity factors also reach a certain threshold critical value characteristic of a given material. This paper presents the results of calculating stress intensity factors for spatial branching cracks. The calculations were carried out on the basis of the original author's program, which makes it possible to determine the stress and displacement fields in a three-dimensional elastic body weakened by a system of arbitrarily oriented cracks. The program is based on the method of discontinuous displacements, implemented for three-dimensional problems of fracture mechanics in an elastic medium, implemented in C ++ codes. The main purpose of the work is to demonstrate the possibility of determining the stress intensity factors for a spatial branching crack and their dependence on the geometric parameters of the problem.

A new crack-tip singular element for cracks in three-dimensional elastic bodies

Symplectic analytical singular element (SASE) is a special crack-tip element in the framework of finite element method (FEM) for modelling cracks. A group of SASEs, which have been developed by the authors, for various crack problems almost in two-dimensional (2D), but not touched yet more complex ones, i.e., in three-dimensional (3D) domain. The underlying reason lies in the fact that analytical symplectic eigen solution for 3D elastic crack problem is not yet available, a 3D SASE thus cannot be constructed through a simple generalization from the 2D SASEs. This study aims to fill out this 3D gap, we thus use a trail displacement field to construct a 3D SASE. The derived strain and stress fields still possess singularity in the vicinity of crack-tip. Finite element formulation is derived based on the minimum total potential energy principle. It is found that some of key features of a 2D SASE still exist in the developed 3D one, e.g., the SIFs can be calculated accurately without any post-processing. Numerical examples are considered to show the accuracy and performance of the proposed element.

Equilibrium stability of nonlinear elastic sphere with distributed dislocations

An exact formulation and a solution of the stability problem for a three-dimensional elastic body containing distributed dislocations are given. The buckling phenomenon for a hollow nonlinear elastic sphere of the semi-linear (harmonic) material with edge dislocations is studied. The study is carried out within the framework of the continuum theory of continuously distributed dislocations using the bifurcation method of buckling analysis. The bifurcation method is to find the equilibrium positions of an elastic body, which differ little from the subcritical (unperturbed) state. The perturbed equilibrium state is described by a linearized boundary value problem. By solving a homogeneous linear boundary value problem, the minimum critical value of the external pressure at which the sphere loses stability is found. The influence of dislocations on the buckling of both thin and thick spherical shells is analyzed.

Stress concentration around a rectangular cuboid hole in a three-dimensional elastic body under tension loading

Stress concentration caused by holes can be investigated by numerical and analytical methods. Current analytical methods can only solve two-dimensional problems. This paper proposes an analytical study on a three-dimensional stress concentration problem that involves a rectangular cuboid hole in a three-dimensional elastic body under tension loading. Based on the finite element method and U-transformation method, the problem can be expressed as a set of uncoupled equations with cyclic periodicity. Displacements of the three-dimensional elastic body are derived in analytical form to study stress distribution in it. Numerical simulation is conducted using ABAQUS to verify the analytical solution. Stress concentration factors in cases of uniaxial, biaxial, and triaxial tensions and the effect of the side ratio of the hole on them are discussed.

Homogenization of Winkler–Steklov spectral conditions in three-dimensional linear elasticity

We consider a homogenization Winkler–Steklov spectral problem that consists of the elasticity equations for a three-dimensional homogeneous anisotropic elastic body which has a plane part of the surface subject to alternating boundary conditions on small regions periodically placed along the plane. These conditions are of the Dirichlet type and of the Winkler–Steklov type, the latter containing the spectral parameter. The rest of the boundary of the body is fixed, and the period and size of the regions, where the spectral parameter arises, are of order $$\varepsilon $$ . For fixed $$\varepsilon $$ , the problem has a discrete spectrum, and we address the asymptotic behavior of the eigenvalues $$\{\beta _k^\varepsilon \}_{k=1}^{\infty }$$ as $$\varepsilon \rightarrow 0$$ . We show that $$\beta _k^\varepsilon =O(\varepsilon ^{-1})$$ for each fixed k, and we observe a common limit point for all the rescaled eigenvalues $$\varepsilon \beta _k^\varepsilon $$ while we make it evident that, although the periodicity of the structure only affects the boundary conditions, a band-gap structure of the spectrum is inherited asymptotically. Also, we provide the asymptotic behavior for certain “groups” of eigenmodes.

Problems of Control of the Dynamics of Incompletely Defined Three-Dimensional Elastic Bodies. II. Discretely Defined Desired State

Problems of Control of the Dynamics of Incompletely Defined Three-Dimensional Elastic Bodies. II. Discretely Defined Desired State

Problems of Control of the Dynamics of Incompletely Defined Three-Dimensional Elastic Bodies. I. Continuously Defined Desired State

The author solves problems of control of a linearly transformed vector function of displacements of points of a three-dimensional elastic body with the purpose of root-mean-square approximation to its continuously defined values. No constraints are imposed on body’s geometry, and observations of its initial-boundary state are continuously defined. Spatial, surface, and initially distributed external-dynamic perturbations are considered as control factors. The accuracy and uniqueness of the control are evaluated.

Uniqueness and Lipschitz stability of an inverse boundary value problem for time-harmonic elastic waves

We consider the inverse problem of determining the Lamé parameters and the density of a three-dimensional elastic body from the local time-harmonic Dirichlet-to-Neumann map. We prove uniqueness and Lipschitz stability of this inverse problem when the Lamé parameters and the density are assumed to be piecewise constant on a given domain partition.

On Mathematical Models of Dynamics of Three-Dimensional Elastic Bodies. Part I. Bodies with Continuously Observed Initial Boundary State

On Mathematical Models of Dynamics of Three-Dimensional Elastic Bodies. Part I. Bodies with Continuously Observed Initial Boundary State