Journal of Elasticity最新文献

A technique for developing convex dual variational principles for the governing PDE of nonlinear elastostatics and elastodynamics is presented. This allows the definition of notions of a variational dual solution and a dual solution corresponding to the PDEs of nonlinear elasticity, even when the latter arise as formal Euler–Lagrange equations corresponding to non-quasiconvex elastic energy functionals whose energy minimizers do not exist. This is demonstrated rigorously in the case of elastostatics for the Saint-Venant Kirchhoff material (in all dimensions), where the existence of variational dual solutions is also proven. The existence of a variational dual solution for the incompressible neo-Hookean material in 2-d is also shown. Stressed and unstressed elastostatic and elastodynamic solutions in 1 space dimension corresponding to a non-convex, double-well energy are computed using the dual methodology. In particular, we show the stability of a dual elastodynamic equilibrium solution for which there are regions of non-vanishing length with negative elastic stiffness, i.e. non-hyperbolic regions, for which the corresponding primal problem is ill-posed and demonstrates an explosive ‘Hadamard instability;’ this appears to have implications for the modeling of physically observed softening behavior in macroscopic mechanical response.

A second-order multiscale theory based on the concept of a Representative Volume Element (RVE) is proposed to link a classical poromechanical model at the RVE scale to a high-order poromechanical model at the macro-scale in the context of finite-strain kinematics. The proposed theory is carefully derived from the Principle of Multiscale Virtual Power, which is a generalization of the Hill-Mandel Principle of Macrohomogeneity. The coupled governing equations of the low-scale and the homogenization rules for the flux and stress-like quantities are obtained by means of standard variational arguments. The main theoretical result is that the minimally constrained space for the pore pressure field allows for non-zero net fluid flow across the RVE boundaries, unlike first-order theories. The direct consequence of this finding is that the present theory can be consistently applied in cases where the low-scale (RVE level) exhibits substantial volume changes (swelling or shrinking) as a consequence of the evolution of the macro-scale kinematics. Details of formulation development and expression for the homogenized tangent operators are presented for those interested in the computational implementation of the model.

Generalizing smooth volumetric growth to the singular case, using de Rham currents and flat chains, we demonstrate how regular boundaries of bodies may evolve to fractals.

The paper considers an axisymmetric stress state of a homogeneous elastic space with a circular disc-shaped crack, one of the edges of which is pressed into a cylindrical circular stamp with static friction. It is assumed that the contact zone is considered under the generalized law of dry friction, i.e. tangential contact stresses are proportional to normal contact pressure, while the proportionality coefficient depends on the radial coordinates of the points of the contacting surfaces and is directly proportional to them. Considering the fact that in this case the Abel images of contact stresses are also related in a similar way, the solution of the problem, with the help of rotation operators and theory of analytical functions, is reduced to an inhomogeneous Riemann problem for two functions and the closed solution in quadratures is constructed. A numerical analysis was carried out and regularities of changes in both normal and shear real contact stresses, as well as rigid displacement of the stamp depending on the physical and geometric parameters were revealed.

Propagation of Rayleigh-type waves is investigated in a half-space composed of nonlocal micropolar thermoelastic material containing void pores. Dispersion relation is derived for a mechanically stress-free and thermally insulated boundary surface of the half-space. The particle motion during the propagation of the waves is found to follow elliptical path. Numerical computations for a specific material are performed to analyze the characteristics of propagating Rayleigh-type waves in detail. Comparison between the phase speed and corresponding attenuation coefficient in some particular cases is also carried out. The effect of various parameters on the characteristics of waves in question is also studied.

We provide (upper and lower) scaling bounds for a singular perturbation model for the cubic-to-tetragonal phase transformation with (partial) displacement boundary data. We illustrate that the order of lamination of the affine displacement data determines the complexity of the microstructure. As in (Rüland and Tribuzio in ESAIM Control Optim. Calc. Var. 29:68, 2023) we heavily exploit careful Fourier space localization methods in distinguishing between the different lamination orders in the data.

For a one dimensional 2-lattice we write down the momentum equation and the equation ruling the shift vector for the dynamic case as a first order system and characterize it in terms of its hyperbolicity. We use similar arguments for problems of increasing difficulty and tackle: one dimensional 3-lattices, three dimensional 2-lattices, three dimensional 3-lattices and finally the most general case of three dimensional ((n+1))-lattices. Our approach is confined to the geometrically and materially linear elastodynamic theory and is valid for generic anisotropic materials that have the multilattice structure. The main finding is that, with the assumptions adopted, the presence of each shift vector is related with zero eigenvalues when the system is written as a first order system.

Recent works have shown that in contrast to classical linear elastic fracture mechanics, endowing crack fronts in a brittle Green-elastic solid with Steigmann-Ogden surface elasticity yields a model that predicts bounded stresses and strains at the crack tips for plane-strain problems. However, singularities persist for anti-plane shear (mode-III fracture) under far-field loading, even when Steigmann-Ogden surface elasticity is incorporated. This work is motivated by obtaining a model of brittle fracture capable of predicting bounded stresses and strains for all modes of loading. We formulate an exact general theory of a three-dimensional solid containing a boundary surface with strain-gradient surface elasticity. For planar reference surfaces parameterized by flat coordinates, the form of surface elasticity reduces to that introduced by Hilgers and Pipkin, and when the surface energy is independent of the surface covariant derivative of the stretching, the theory reduces to that of Steigmann and Ogden. We discuss material symmetry using Murdoch and Cohen’s extension of Noll’s theory. We present a model small-strain surface energy that incorporates resistance to geodesic distortion, satisfies strong ellipticity, and requires the same material constants found in the Steigmann-Ogden theory. Finally, we derive and apply the linearized theory to mode-III fracture in an infinite plate under far-field loading. We prove that there always exists a unique classical solution to the governing integro-differential equation, and in contrast to using Steigmann-Ogden surface elasticity, our model is consistent with the linearization assumption in predicting finite stresses and strains at the crack tips.

A complete closed-form solution is derived to the finite plane elasticity problem associated with an elliptical incompressible liquid inclusion embedded in an infinite compressible hyperelastic matrix of harmonic type subjected to uniform remote in-plane Piola stresses. The internal uniform hydrostatic Piola stresses and deformation gradients within the elliptical incompressible liquid inclusion are obtained. The internal uniform hydrostatic tension and hoop stress on the matrix side along the liquid-solid interface for various shapes of the elliptical liquid inclusion and different elastic constants of the matrix under various remote loadings are calculated.

We study the effective properties of a linear elastic or perfectly viscoplastic composite consisting of a soft matrix reinforced by possibly very stiff fibers. The effective material is characterized by the emergence of a discrepancy between the displacement in the fibers and the overall averaged displacement, giving rise to a concentration strain in the matrix. The effective energy stored in the fibers is a combination of stretching, bending and torsional energies. This work completes and corrects results previously obtained by the author and coworkers, who failed to notice the torsional contribution.