Three-dimensional Linear Elasticity Research Articles

Gradient preserving Operator Inference: Data-driven reduced-order models for equations with gradient structure

Hamiltonian Operator Inference has been introduced in Sharma et al. (2022) to learn structure-preserving reduced-order models (ROMs) for Hamiltonian systems. This approach constructs a low-dimensional model using only data and knowledge of the Hamiltonian function. Such ROMs can keep the intrinsic structure of the system, allowing them to capture the physics described by the governing equations. In this work, we extend this approach to more general systems that are either conservative or dissipative in energy, and which possess a gradient structure. We derive the optimization problems for inferring structure-preserving ROMs that preserve the gradient structure. We further derive an a priori error estimate for the reduced-order approximation. To test the algorithms, we consider semi-discretized partial differential equations with gradient structure, such as the parameterized wave and Korteweg–de-Vries equations, and equations of three-dimensional linear elasticity in the conservative case and the one- and two-dimensional Allen-Cahn equations in the dissipative case. The numerical results illustrate the accuracy, structure-preservation properties, and predictive capabilities of the gradient-preserving Operator Inference ROMs.

A family of three-dimensional virtual elements for Hellinger-Reissner elasticity problems

We present a family of Virtual Element Methods for three-dimensional linear elasticity problems based on the Hellinger-Reissner variational principle. A convergence and stability analysis is developed. Moreover, using the hybridization technique and exploiting the information derived from this procedure, we show how to compute a better approximation for the displacement field. The numerical experiments confirm the theoretical predictions.

Development of quadratic enhanced assumed strain elements for three-dimensional linear elasticity

Enhanced assumed strain (EAS) elements are widely used in the literature to address the issue of locking associated with conventional elements. However, existing literature in the context of three-dimensional (3D) elasticity problems is predominantly restricted to the eight-node linear EAS elements. Thus, existing 3D EAS elements do not exploit the superior performance offered by quadratic elements over linear elements. In the current work, we propose a novel twenty-seven node quadratic EAS element, which, to the best of our knowledge, is the first such attempt in the literature. Additionally, the manuscript also presents a six-node wedge and an eighteen-node wedge EAS element. The proposed elements are derived methodically by investigating the interrelation between the two-field Hellinger–Reissner (HR) and three-field Veubeke–Hu–Washizu (VHW) variational formulations. The robustness and performance of the proposed EAS elements is demonstrated through numerous examples.

A multiple-data-based direct method for inverse problem in three-dimensional linear elasticity

Estimating the unknown Young’s modulus and boundary conditions of an elastic object from locally observed boundary data is an inverse problem which is usually solved by iterative methods. In this paper, a multiple-data-based direct (MD) method is proposed to solve this inverse problem without iterations. The proposed method applies a global normalized objective function and a global regularization coefficient τg for energy-like regularization to ensure the convergence of the unknown displacement boundary conditions. A coarse-mesh based preprocessing algorithm is proposed to speed up the determination of the regularization coefficient. A multi-force-position method is proposed to obtain observation data, which improves accuracy of the MD method when the number of experiments is limited. Both numerical and physical experiments were conducted. Compared with the single-data-based direct (SD) methods proposed in our previous work, the MD method yields more accurate estimation of Young’s modulus and boundary conditions in both numerical and physical experiments.

Computation of Three-Dimensional Linear Elasticity Problem using Mesh and Mesh-free Particle Modelling Techniques

Computation of Three-Dimensional Linear Elasticity Problem using Mesh and Mesh-free Particle Modelling Techniques

Computation of Three-Dimensional Linear Elasticity Problem using Mesh and Mesh-free Particle Modelling Techniques

Computation of Three-Dimensional Linear Elasticity Problem using Mesh and Mesh-free Particle Modelling Techniques

Damage as a Material Phase Transition

AbstractWe propose paradigmatic examples to show how material damage phenomena can be efficiently described as a solid-solid phase transition. Starting from the pioneering work of J.L. Ericksen (J. Elast. 5(3):191–201, 1975) and the extensions of R.L. Fosdick and other authors to three-dimensional non linear elasticity, we describe the insurgence of damage as a hard → soft transition between two material states (damage and undamaged) characterized by two different energy wells. We consider the two separate constitutive assumptions of a simple Neo-Hookean type damageable material and a more complex microstructure inspired damageable Gent type material with variable limit threshold of the first invariant. In both cases we study two different deformation shear classes, one homogeneous and the other one inhomogeneous and obtain fully analytic description of the system damage response under cyclic loading. The considered constitutive assumptions and deformation classes are aimed at attaining fully analytic descriptions. On the other hand, we remark that the proposed, Griffith type, variational approach of damage, based on two different energy density functions for the damaged and undamaged material phases, and a resulting non (rank-one) convex energy, can be extended to systems with more complex energy functions, possibly with a larger number of wells representing an increasing degree of damage.

Inverse of initial stress in carbon fiber reinforced polymer laminates using lamb waves and deep neural network

The prediction of the initial stress in composites is essential for the non-destructive testing (NDT) and structural health monitoring (SHM) of carbon fibre reinforced polymer (CFRP). This paper examines the potential of Lamb waves in the inverse of initial stress by calculating the influence of initial stress on the dispersion characteristics of Lamb waves propagating in multilayered CFRP laminates. By introducing the mechanics of incremental deformation into the linear three-dimensional elasticity theory, the Legendre orthogonal polynomial expansion (LOPE) method is used to mathematically model the Lamb wave propagating in multilayered CFRP laminates subjected to horizontal and vertical homogeneous initial stresses. Then, a three-hidden-layers Feed Forward Deep Neural Network (DNN) with Back Propagation (BP) algorithm is constructed to invert the magnitude and direction of the initial stresses. The input features are the phase velocities of fundamental Lamb wave A0 mode at five different frequencies. Both training and testing samples are obtained by LOPE forward calculation. An ablation experiment is presented to compare the two different activation functions. Finally, the accuracy of the inverse is verified by comparing with the available outcomes of LOPE forward calculation.

Varvec{J}-, M- and varvec{L}-integrals of line charges and line forces

In this work, the varvec{J}-, M- and varvec{L}-integrals of two line charges and two line forces in generalized plane strain are derived in the framework of three-dimensional (3D) electrostatics and three-dimensional compatible linear elasticity, respectively, in order to study the interaction between them. The key point in this derivation is achieved by expressing the varvec{J}-, M- and varvec{L}-integrals of point charges and point forces in three dimensions in terms of the corresponding three-dimensional Green functions, that is, the three-dimensional Green function of the Laplace operator and the three-dimensional Green tensor of the Navier operator, respectively. The major mathematical tool used in deriving the varvec{J}-, M- and L_3-integrals of line charges and line forces from the corresponding varvec{J}-, M- and varvec{L}-integrals of point charges and point forces is the method of embedding or method of descent of Green functions to two dimensions (2D) from the corresponding Green functions in 3D. The analytical expressions of varvec{J}-, M- and L_3-integrals of line charges and line forces in antiplane and plane strain are derived and discussed. The varvec{J}-integral is the electrostatic part of the Lorentz force (electrostatic interaction force) between two line charges in electrostatics and the Cherepanov force between two line forces in elasticity. The M-integral of two line sources (charges or forces) equals half the electrostatic interaction energy in electrostatics and half the elastic interaction energy in elasticity of these two line sources, respectively. The L_3-integral of two line sources (charges or forces) is the z-component of the configurational vector moment or the rotational moment representing the total torque about the z-axis caused by the interaction of the two line sources. An important outcome is that the varvec{J}-integral is twice the negative gradient of the M-integral, leading to the result that the {varvec{J}}-integral for line charges and line forces is a conservative force (with the corresponding interaction energy playing the role of the potential energy) and consequently an irrotational vector field. Finally, the obtained varvec{J}-, M- and L_3-integrals being functions of the distance and of the angle, are able to describe the physical behaviour of the interaction of two line sources (charges and forces), since they represent the fundamental and necessary quantities, that is, the interaction force, interaction energy and total torque produced by the interaction of these two line sources.

Numerical Study on the Buckling Behavior of FG Porous Spherical Caps Reinforced by Graphene Platelets

The buckling response of functionally graded (FG) porous spherical caps reinforced by graphene platelets (GPLs) is assessed here, including both symmetric and uniform porosity patterns in the metal matrix, together with five different GPL distributions. The Halpin-Tsai model is here applied, together with an extended rule of mixture to determine the elastic properties and mass density of the selected shells, respectively. The equilibrium equations of the pre-buckling state are here determined according to a linear three-dimensional (3D) elasticity basics and principle of virtual work, whose solution is determined from classical finite elements. The buckling load is, thus, obtained based on the nonlinear Green strain field and generalized geometric stiffness concept. A large parametric investigation studies the sensitivity of the natural frequencies of FG porous spherical caps reinforced by GPLs to different parameters, namely, the porosity coefficients and distributions, together with different polar angles and stiffness coefficients of the elastic foundation, but also different GPL patterns and weight fractions of graphene nanofillers. Results denote that the maximum and minimum buckling loads are reached for GPL-X and GPL-O distributions, respectively. Additionally, the difference between the maximum and minimum critical buckling loads for different porosity distributions is approximately equal to 90%, which belong to symmetric distributions. It is also found that a high weight fraction of GPLs and a high porosity coefficient yield the highest and lowest effects of the structure on the buckling loads of the structure for an amount of 100% and 12.5%, respectively.

Derivation of an elastic shell model of order as tends to infinity

Starting from three-dimensional linear elasticity and performing the dimensional reduction by integration over the thickness, we derive a general form of the areal strain energy density for elastic shells. To obtain the new constitutive model, we do not approximate the deformation fields as polynomials in the thickness coordinate, but rather we keep all terms in the thickness-wise series expansions. As a result, we deduce the explicit form of the shell strain energy density in which the constitutive coefficients are expressed as integrals depending on the thickness h and on the initial curvature. Then, to obtain the shell model of order , we expand the integral coefficients in the strain energy function as power series of h and truncate them to the power . In the case , we recover the classical strain energy density for combined bending and stretching of linear shells, which leads to the Koiter model. Finally, we prove that the proposed shell strain energy function is coercive for any , as well as for the general case .

Spectral and Evolution Analysis of Composite Elastic Plates with High Contrast

We analyse the behaviour of thin composite plates whose material properties vary periodically in-plane and possess a high degree of contrast between the individual components. Starting from the equations of three-dimensional linear elasticity that describe soft inclusions embedded in a relatively stiff thin-plate matrix, we derive the corresponding asymptotically equivalent two-dimensional plate equations. Our approach is based on recent results concerning decomposition of deformations with bounded scaled symmetrised gradients. Using an operator-theoretic approach, we calculate the limit resolvent and analyse the associated limit spectrum and effective evolution equations. We obtain our results under various asymptotic relations between the size of the soft inclusions (equivalently, the period) and the plate thickness as well as under various scaling combinations between the contrast, spectrum, and time. In particular, we demonstrate significant qualitative differences between the asymptotic models obtained in different regimes.

A linear two-dimensional mathematical model for thin two-layer plates with partial shear interaction, with a view towards application to laminated glass

This work presents a consistent derivation, from three-dimensional linear elasticity, of a two-dimensional mathematical model describing the bending and in-plane stretching behaviours, under a general system of quasi-static distributed loads and prescribed support displacements, of thin two-layer plates with partial shear interaction. This layerwise model is specifically tailored for the requirements posed by the analysis of laminated glass plates commonly used in building structures (consisting of two thin glass layers bonded together by an adhesive interlayer). Our approach, based on Podio-Guidugli’s method of internal constraints, avoids mutually contradictory assumptions (not uncommon in the literature on structural mechanics) and yields a complete two-dimensional characterisation of displacement, strain and stress fields that exactly satisfy the field equations of three-dimensional linear elasticity and the boundary conditions at the end faces. The choice of generalised variables is designed to bring to light the following fundamental conclusion and physical insight: the resulting two-dimensional boundary value problem is a combination of the equations of Kirchhoff and Mindlin plates (with specified rigidities).Two examples illustrate the application of the proposed model: (i) the cylindrical bending of plate strips and (ii) a family of problems with Navier-type analytical solution. The solutions exhibit continuity across the whole range of zero, partial and full interaction between the layers. Moreover, the Navier-type solutions are consistently in close agreement with the results of three-dimensional finite element analyses. On the contrary, analogous results previously reported in the literature exhibit considerable deviations. An explanation for these discrepancies is discussed in detail.

A FEM-based direct method for identification of Young’s modulus and boundary conditions in three-dimensional linear elasticity from local observation

Simultaneously identifying unknown Young’s modulus and boundary conditions of a linear elastic material with the measurements only on a part of its boundary is an inverse problem, which is usually solved by iterative methods in most studies. In this paper, we present a direct inverse method based on the finite element method (FEM) to recover these unknown parameters without iteration. A novel energy-like regularization term is developed to ensure the convergence of the estimated parameters with noisy input. In order to increase the accuracy of estimation, a novel regularization coefficient and its related preprocessing steps are designed based on the Morozov’s discrepancy principle to obtain the optimal regularization parameter, which is independent to the observation noise. In numerical experiments, several three-dimensional objects with different shapes and distribution of observed and unknown regions were used to verify the proposed method under different observation noises. Physical experiments were also conducted on a silicone phantom. The results of both numerical and physical experiments successfully recovered the Young’s modulus and displacement boundary conditions with expected accuracies close to their corresponding observation noise levels.

Parallel Scalability of Three-Level FROSch Preconditioners to 220000 Cores using the Theta Supercomputer

The parallel performance of the three-level fast and robust overlapping Schwarz (FROSch) preconditioners is investigated for linear elasticity. The FROSch framework is part of the Trilinos software library and contains a parallel implementation of different preconditioners with energy minimizing coarse spaces of generalized Dryja--Smith--Widlund type. The three-level extension is constructed by a recursive application of the FROSch preconditioner to the coarse problem. In this paper, the additional steps in the implementation in order to apply the FROSch preconditioner recursively are described in detail. Furthermore, it is shown that no explicit geometric information is needed in the recursive application of the preconditioner. In particular, the rigid body modes, including the rotations, can be interpolated on the coarse level without additional geometric information. Parallel results for a three-dimensional linear elasticity problem obtained on the Theta supercomputer (Argonne Leadership Computing Facility, Argonne, IL) using up to 220 000 cores are discussed and compared to results obtained on the SuperMUC-NG supercomputer (Leibniz Supercomputing Centre, Garching, Germany). Notably, it can be observed that a hierarchical communication operation in FROSch related to the coarse operator starts to dominate the computing time on Theta, which has a dragonfly interconnect, for 100 000 message passing interface (MPI) ranks or more. The same operation, however, scales well and stays within the order of a second in all experiments performed on SuperMUC-NG, which uses a fat tree network. Using hybrid MPI/OpenMP parallelization, the onset of the MPI communication problem on Theta can be delayed. Further analysis of the performance of FROSch on large supercomputers with dragonfly interconnects will be necessary.

Three dimensional free vibration analysis of functionally graded graphene reinforced composite laminated cylindrical panel

In this paper, natural frequency analysis of functionally graded graphene platelet (GPL) reinforced composite cylindrical panel is investigated. Halpin-Tsai equations are employed to obtain the mechanical properties of the structure. The linear three-dimensional elasticity theory based on the Hamilton principle and the numerical finite element method are used for obtaining the governing equations of motion. Four patterns of GPL distributions such as: FG-X, FG-V, FG-O and UD are assumed through the thickness of shell. The effects of various parameters such as different distributions and weight fractions of GPL, zigzag and armchair lay-up, geometry and different boundary conditions on the natural frequencies and mode shapes of shell have been investigated. The results indicate that the maximum and minimum natural frequencies belong to FG-X pattern and FG-O pattern, respectively.

Effect of extreme in-plane and transverse stiffness ratios on frequencies and load transfer between face sheets and core of rectangular sandwich plates

Challenging issues for sandwich structures include accurately finding their frequencies and how the load transfers from face sheets to the core when their material stiffness ratio (FCSR) in either the axial or the transverse directions varies from 1 to 107. Here we use both the equivalent single layer and the layer wise shear and normal deformable plate theories to analyze problems for rectangular plates of different thickness/length ratios. The Ritz method with Jacobi polynomials as basis functions is employed to numerically solve three-dimensional linear elasticity theory equations. Transverse stresses are calculated using a one-step stress recovery scheme. It is found that the FCSR in the transverse direction determines the minimum degree of complete polynomials in the thickness coordinate needed to accurately compute the lowest six frequencies and the transverse stresses. New results include the following: (i) the FCSR in the transverse direction has a greater influence than that in the axial direction on the computed frequencies and transverse stresses, and (ii) whereas both transverse normal and shear stresses in the core contribute to transmitting loads between the two face sheets for FCSR < 103, however, only the former is effective for FCSR ≥ 105.

Analysis of Longitudinal Guided Wave Propagation in the Functionally Graded Hollow Cylinder Using State-Vector Formalism and Legendre Polynomial Hybrid Approach

A numerical approach for analyzing the guided wave propagation behavior in a functionally graded material (FGM) hollow cylinder is presented. Based on the three-dimensional linear elasticity theory, the state matrix form of displacement and stress component is introduced into the differential governing equation. Combining with boundary conditions, the dispersion characteristic equations can be constructed. Meanwhile, taking the orthogonal completeness of Legendre series into consideration, the exhausting root findings of transcendental equation can be avoided, and it transforms the redundant integral operators into analytical expressions by means of its recursive properties. Moreover, it is not only unnecessary to stratify FGMs, and introducing the univariate nonlinear regression to modelling the arbitrary gradient distribution, and transforming the redundant integral operator into an analytical expression by means of recursive property. It should be pointed out that the coefficient matrix consists only zero and constant terms, which are related to the elastic constants and normalized wave number. And it provides a fast and effectively approach to extracting the dispersion curves, displacement distribution and stress profile, simultaneously. Meanwhile, the proposed method can get rid of the exhausted root-finding algorithm, and accordingly distinguish the different modes and the dispersion curves simultaneously. Finally, the typical non-stratified computation of dispersion curves of longitudinal guided wave for FGM hollow cylinder is realized. The accuracy of the proposed method is further demonstrated through comparison with the available data from Legendre polynomial method. Several numerical cases about iron based alumina FGM cylinder are studied, and the convergence of the present polynomial approach is discussed. Then the effect of the circumferential wave number, gradient variation, diameter-thickness ratio on the dispersion characteristics are illustrated. Moreover, the distribution characteristics of displacement and stress profiles of F(1,1) mode at a given frequency are discussed in detail. It is shown that the state-vector formalism and Legendre polynomial hybrid approach is a powerful tool in solving wave motions in FGM elastic cylinder, which lays a theoretical foundation for the non-destructive evaluation and quantitative characterization of the gradient structure characteristics of functionally graded hollow cylinder.

A hybrid mimetic spectral element method for three-dimensional linear elasticity problems

We introduce a domain decomposition structure-preserving method based on a hybrid mimetic spectral element method for three-dimensional linear elasticity problems in curvilinear conforming structured meshes. The method is an equilibrium method which satisfies pointwise equilibrium of forces. The domain decomposition is established through hybridization which first allows for an inter-element normal stress discontinuity and then enforces the normal stress continuity using a Lagrange multiplier which turns out to be the displacement in the trace space. Dual basis functions are employed to simplify the discretization and to obtain a higher sparsity. Numerical tests supporting the method are presented.

A displacement potential function using complex variables for numerical computations of three-dimensional elasticity problems

This paper shows the development of a displacement potential function based on the Galerkin potential using complex variables. The displacement potential function results in a more suitable method for numerical calculations since it avoids the strenuous integration process associated with stress potential methods. Completeness of the displacement potential function is demonstrated. The displacement potential function was applied to the solution of the first fundamental problem of elasticity over a three-dimensional domain with known boundary conditions. It’s application for numerical calculations is demonstrated by solving the pure shear problem over a three-dimensional unit hexahedral cell. Finally, the obtained numerical results are compared against finite element results, proving the validity of the displacement potential function in solving three-dimensional linear elasticity problems.