Elastic Equation Research Articles

An improved stochastic inversion method for 3D elastic impedance under the prior constraints of random medium parameters

Three-dimensional (3D) elastic impedance inversion is widely used in oil and gas exploration with its excellent identification of reservoirs and noise immunity. However, modifying the prior information directly in inversion will undoubtedly destroy the spatial autocorrelation structure of the elastic impedance. Maintaining the elastic impedance spatial autocorrelation structure in seismic inversion is of great significance for guiding the reservoir identification and prediction. The 3D prior information construction method is proposed under the constraint of random medium parameters. This method establishes an intrinsic connection between random medium parameters and elastic impedance based on seismic data and well data, which provides a theoretical basis for simulating prior information satisfying the spatial correlation of the subsurface medium. To ensure the invariance of the spatial structure term in the inversion process, the random term and the spatial structure term in the simulation process are separated by the gradual deformation method (GDM), and the priori constrained objective function is further constructed to reduce the inversion uncertainty. Numerical example shows that the relative error of the proposed simulation method is smaller than that of conventional sequential Gaussian simulation, and the details are more accurately depicted based on the priori inversion results in this study. In addition, the elastic impedance equation based on the p-wave modulus is derived by combining with the reservoir sensitive parameters in field data. The field examples document that this method can effectively perform high-resolution geophysical prediction of underground thin-layer reservoirs.

Quantitative osteoporosis diagnosis of porous cancellous bone using poroelastodynamic modal analysis.

Osteoporosis is a skeletal disease characterized by reduced bone mass and microarchitectural deterioration, leading to increased fragility. This study presents a novel three-dimensional poroelastodynamic model for analyzing cancellous bone free vibration responses. The model incorporates the Navier-Stokes equations of linear elasticity and the Biot theory of porous media, allowing the investigation of osteoporosis-related changes. The analysis considers parameters like porosity, density, elasticity, Poisson ratio, and viscosity of bone marrow within the porous medium. Our findings indicate that natural frequencies of cancellous bone play a crucial role in osteoporosis prediction. By incorporating experimental data from 12 mouse femurs, we unveil insights into osteoporosis prediction. Increased porosity reduces bone stiffness, lowering natural frequencies. However, it also increases bone mass loss relative to stiffness, leading to higher frequencies. Therefore, the natural frequencies of osteoporotic bone are always higher than the natural frequencies of normal bone. Additionally, an increase in bone marrow within the pores, while increasing damping effects, also increases natural frequencies, which is another indication of osteoporosis growth in bone. The presence of bone marrow within the pores further influences natural frequencies, providing additional insights into osteoporosis growth. Thinner and smaller bones are found to be more susceptible to osteoporosis compared to larger and bigger bones due to their higher natural frequencies at equivalent porosity levels.

Small-amplitude waves in a floating poroelastic plate forcing by vertical pitching plate

The linear two-dimensional problem of flexural-gravity waves generated by an oscillating rigid plate build-in a floating poroelastic plate is studied. The problem is coupled. The plate deflections and the hydrodynamics loads are determined at the same time. The liquid under the poroelastic plate is inviscid and incompressible. Dynamics of the floating plate is described by a thin elastic plate equation. Porosity of the floating plate is taken into account only through the liquid flux into the plate. The velocity of the inflow is assumed to be governed by Darcy's law being proportional to the hydrodynamic pressure at the plate/liquid interface. Two cases of the oscillating rigid plate with and without its part in the liquid are considered. The problems are solved by the Fourier transform method for non-zero porosity and by the vertical mode method for elastic plates with zero porosity. The deflection and strain distributions are analyzed depending on the excitation frequency and the porosity. Two models of floating plate porosity, where the hydrostatic pressure is included into Darcy's law (Zavyalova's model) and excluded (Meylan's model), are compared. Plate porosity induces damping to the system. It is shown that the damping rate is non-monotonic with respect to the plate porosity.

Analysis‐aware defeaturing of complex geometries with Neumann features

AbstractLocal modifications of a computational domain are often performed in order to simplify the meshing process and to reduce computational costs and memory requirements. However, removing geometrical features of a domain often introduces a non‐negligible error in the solution of a differential problem in which it is defined. In this work, we extend the results from Buffa et al. (Math Models Methods Appl Sci. 2022; 32(02):359–402.) by studying the case of domains containing an arbitrary number of distinct Neumann features, and by performing an analysis on Poisson's, linear elasticity, and Stokes' equations. We introduce a simple, computationally cheap, reliable, and efficient a posteriori estimator of the geometrical defeaturing error. Moreover, we also introduce a geometric refinement strategy that accounts for the defeaturing error: starting from a fully defeatured geometry, the algorithm determines at each iteration step which features need to be added to the geometrical model to reduce the defeaturing error. These important features are then added to the (partially) defeatured geometrical model at the next iteration, until the solution attains a prescribed accuracy. A wide range of two‐ and three‐dimensional numerical experiments are finally reported to illustrate this work.

Analysis of Stress Concentration in Functionally Graded Plates with Linearly Increasing Young’s Modulus

In this article, the strain and stress analyses of functionally graded plates with circular holes that are subject to a uniaxial far-field traction load are analytically considered. The Young’s modulus is assumed to vary linearly along the radial direction around the hole. The adoption of such a type of inhomogeneity variation can be justified as follows. Firstly, and among all the possible variations of stiffness, the linear one is indeed the simplest inhomogeneity distribution. Surprisingly however, according to our knowledge extent, the associated elastic fields were not yet addressed in the literature. Secondly, a linearly varying stiffness could reasonably imply a remarkable advantage from a technological point of view. In fact, unlike nonlinearly varying stiffness plates, manufacturing routes are only required to handle constant variations throughout the radial domain. After recalling the basic equations for plane stress elasticity, the displacement, strain, and stress fields around the hole were numerically tackled and discussed for different stiffness ratios. A comparison was also carried out with other Young’s modulus distributions that have been commonly employed in the literature.

Numerical model of nonlinear elastic bulk wave propagation in solids for non-destructive evaluation

Nonlinear ultrasonic techniques can be difficult and non-intuitive to understand due to the range of wave mixing combinations available between similar or different wave modes. To overcome this, a numerical model that uses a Finite Difference Time Domain (FDTD) scheme, without a staggered grid system, to solve the nonlinear elastic bulk wave equations in two dimensions is proposed in this paper, with the purpose of better understanding nonlinear ultrasonic techniques. Both material and geometrical nonlinearities are considered and a stress-type boundary condition is used to model the excitation. The energy transfer between frequency bands is shown and the amplitude trend of the harmonics are validated against available theories. It is then used to simulate the nonlinear ultrasonic field generated by a finite-size element of an array in a solid to better understand its behaviour. An array-element directivity at the second harmonic frequency is shown. Since the FDTD model does not account for attenuation, a correction using a ray-based approach is applied to the simulated result for direct comparison against experimental measurements. The field obtained from a typical array used in non-destructive testing is then simulated to characterise its behaviour. Experimental comparison of the nonlinear displacement fields for different focal positions showed good agreement, further validating the FDTD model. This FDTD model opens opportunities to virtually experiment and design appropriate nonlinear ultrasonic array inspections whilst isolating and understanding contribution from material nonlinearity.

Accuracy verification of fast adaptive galerkin finite volume solver for thermal and mechanical load crack propagation in triangular meshes

This study presents the details of the Galerkin finite volume method for triangular unstructured meshes and its application for the fracture analysis of 2D elasticity problems under simultaneous thermal and mechanical loading. The heat diffusion and elasticity equations are discretized using the Galerkin finite volume method (which, due to omitting matrix manipulation calculations, computations are considerably faster than FEM and XFEM Solvers). The maximum tangential stress criterion is presented after briefly describing the interaction integral formulation used to calculate the stress intensity factors in thermo-mechanical problems. The accuracy of the computed results of the present strategy under thermo-mechanical loads is demonstrated using some benchmark test cases. First, a thermal stress analysis of a plate with an inclined central crack under thermal boundary conditions is performed. Second, the development of a crack under mechanical stress boundary conditions is modeled. Third, crack propagation under both thermal and mechanical boundary conditions is simulated. The present modeling strategy's results are compared with those reported in previous numerical works to verify the accuracy. The stress intensity factors and predicted crack trajectories are utilized to assess the accuracy of computed results and investigate the quality of crack simulation by the proposed numerical method.

Isogeometric Tearing and Interconnecting solvers for linear elasticity in multi-patch Isogeometric Analysis with theory for two dimensional domains

We consider the linear elasticity equation, discretized with multi-patch Isogeometric Analysis. A standard discretization error analysis is based on Korn’s inequality, which degrades for certain geometries, such as long and thin cantilevers. This phenomenon is known as geometry locking. We observe that high-order methods, like Isogeometric Analysis is beneficial in such a setting. The main focus of this paper is the construction and analysis of a domain decomposition solver, namely an Isogeometric Tearing and Interconnecting (IETI) solver, where we prove that the convergence behavior does not depend on the constant of Korn’s inequality for the overall domain, but only on the corresponding constants for the individual patches. Moreover, our analysis is explicit in the choice of the spline degree. Numerical experiments are provided which demonstrates the efficiency of the proposed solver.

Deformation of an elastic second gradient spherical body under equatorial line density of dead forces

We consider deformations of an elastic body having initially a spherical shape. Assumed deformation energy depends on the first and second gradient of displacements. We apply an equatorial line density of dead loads, that are forces per unit line length directed in radial direction and applied along the equator of the sphere. We restrict ourselves our analysis to the case of linearized second strain gradient isotropic elasticity (for which the more general energy was determined by Mindlin) with only one characteristic length. Differently to what happens in first gradient continua, i.e. in classic linear elasticity, we show that for the particular class second gradient continua considered here these forces do not determine infinite displacements in the direction of applied dead line forces. Instead, using a series method for the solution of the considered boundary-value problem, we demonstrate that the displacements are finite. So in the deformed configuration there is not the formation of an edge at the material points where the forces are applied. Further investigations are therefore needed for establishing if this elastic-regime edge formation is made possible: (I) either in the case of more general linear elastic constitutive equations or (II) only when large deformations are considered or (III) if non-elastic phenomena are involved.

A Splitting Nearly-Analytic Symplectic Partitioned Runge–Kutta Method for Solving 2D Elastic Wave Equations

In this paper, we present a new splitting nearly-analytic symplectic partitioned Runge–Kutta (SNSPRK) method for the two-dimensional (2D) elastic wave equations. It is an extension to elastic wave equation of our recent work on the locally one-dimensional nearly-analytic symplectic partitioned Runge–Kutta (LOD-NSPRK) method for the 2D acoustic wave equations. The method is based on the spatial differential operator-split technique, in which the resulting spatial discrete matrices are symmetric unlike the conventional nearly-analytic symplectic partitioned Runge–Kutta (NSPRK) method. The stability condition is given, which has more relax restriction for the time step than the NSPRK schemes. To show the performance of the new method, numerical experiments are given. Numerical results illustrate that the SNSPRK method has better long-term calculation capability as time proceeds than the NSPRK method.

Randomized vector iterative linear solvers of high precision for large dense system

Abstract In this paper we suggest randomized linear solvers with a focus on refinement issue to achieve a high precision while maintaining all the advantages of the Monte Carlo method for solving systems of large dimension with dense matrices. It is shown that each iterative refinement step reduces the error by one order of magnitude. The crucial point of the suggested method is, in contrast to the standard Monte Carlo method, that the randomized vector algorithm computes the entire solution column at once, rather than a single component. This makes it possible to efficiently construct the iterative refinement method. We apply the developed method for solving a system of elasticity equations.

Short-time existence of a quasi-stationary fluid–structure interaction problem for plaque growth

Abstract We address a quasi-stationary fluid–structure interaction problem coupled with cell reactions and growth, which comes from the plaque formation during the stage of the atherosclerotic lesion in human arteries. The blood is modeled by the incompressible Navier-Stokes equation, while the motion of vessels is captured by a quasi-stationary equation of nonlinear elasticity. The growth happens when both cells in fluid and solid react, diffuse and transport across the interface, resulting in the accumulation of foam cells, which are exactly seen as the plaques. Via a fixed-point argument, we derive the local well-posedness of the nonlinear system, which is sustained by the analysis of decoupled linear systems.

Simple and multiple traveling waves in a reaction-diffusion-mechanics model

In this paper, we consider a reaction-diffusion system with mechanical deformation of medium. This system consists of an excitable system bidirectionally coupled with an elasticity equation. The main content consists of two parts. First, for γ>0 sufficiently small, simple pulse of homoclinic type exists. We prove that the traveling pulse is linearly stable. Specifically, there is at most one nontrivial eigenvalue near the origin and it is negative. Second, for γ>0 large, we show existence of double twisted front-back wave loop, indicating bifurcations of various complicated traveling waves, including N-front, N-back and wave train. Then we prove that N-fronts are linearly stable. Our arguments are mainly based on geometric singular perturbation theory, exponential dichotomy, heteroclinic bifurcation and the Melnikov method.

Численное решение двумерной задачи определения скорости распространения сейсмических волн в неоднородных cредах с памятью

The numerical method for two-dimensional inverse dynamic seismic problem for a viscoelastic isotropic medium is presented. The system of differential equations of elasticity for isotropic medium of memory type is considered as a mathematical model. The unknown values are the displacement, the memory function of the medium (the kernel of the integral term) and the propagation velocity of elastic waves in a weakly horizontally inhomogeneous medium. Additional information for the inverse problem is the response displacement measured on the surface. The method is based on reducing the inverse problem to a system of Volterra-type integral equations and their sequential numerical implementation. The results of the study are analyzed and compared with the analytical solution. It is shown that the results are in satisfactory agreement.

Higher-order trigonometric series-based analytical solution to free transverse vibration of suspended laminated composite slabs

Suspended floors can be designed as effective pendulum seismic isolators for high-rise buildings. However, research studies on the free vibration of suspended laminated composite slabs are relatively limited. This paper presents a comprehensive analytical solution to the free vibration problem of suspended laminated composite slabs. The equations of elasticity are used to establish the equation of motion. The critical factors influencing free vibration are explicitly considered, including the material anisotropy, number of layers, elastic properties, rotatory inertia, and transverse shear of plates. The higher-order trigonometric series are utilized to solve the dynamic equations. The developed analytical solution is a complete and realistic form of the existing analytical solutions in literature and has the capability of converging fast. The proposed analytical solution does not have a dependency on the shape function as well as separate-of-variable forms unlike the Finite Element Method (FEM). Moreover, the FEM requires an in-depth mesh convergence analysis, particularly for higher-magnitude natural frequencies, which is not the case for the proposed analytical solution for any infinite range of eigenfrequencies. The proposed solution procedure is first verified by the simplified problems available in the literature. For more complex problems, the finite element analysis results obtained from Abaqus are employed to validate the analytical solution. The comparison demonstrates that the proposed analytical solution is accurate and reliable for a wide range of case-study examples. It is found that the natural frequencies of suspended laminated composite slabs are significantly influenced by the crucial factors under investigation.

Numerical simulations for largely deformed beams and rings adopting a nontensile smoothed particle hydrodynamics algorithm

Three typical elastic problems, including beam bending, truss extension and compression, and two-rings collision are simulated with smoothed particle hydrodynamics (SPH) using Lagrangian and Eulerian algorithms. A contact-force model for elastic collisions and equation of state for pressure arising in colliding elastic bodies are also analytically derived. Numerical validations, on using the corresponding theoretical models, are carried out for the beam bending, truss extension and compression simulations. Numerical instabilities caused by largely deformed particle configurations in finite/large elastic deformations are analysed. The numerical experiments show that the algorithms handle small deformations well, but only the Lagrangian algorithm can handle large elastic deformations. The numerical results obtained from the Lagrangian algorithm also show a good agreement with the theoretical values. doi: 10.1017/S1446181123000160

A Hybrid Reproducing Kernel Particle Method for Three-Dimensional Elasticity Problems

This study presents a fast meshless method called the hybrid reproducing kernel particle method (HRKPM) for the solution of three-dimensional (3D) elasticity problems. The equilibrium equations of 3D elasticity are divided into three groups of equations, and two equilibrium equations are contained in each group. By coupling the discrete equations for solving two arbitrary groups of equations, the complete solution of 3D elasticity can be obtained. For an arbitrary group of equations, the 3D elasticity problem is transformed into a series of associated two-dimensional (2D) ones, which is solved by the RKPM to derive the discrete formulae. The discrete equations of 2D problems are combined using the difference method in dimension splitting direction. Then, arbitrarily choosing another group of equilibrium equations, the discrete equation of another group of 2D problems can be obtained similarly. By combining the discrete equations for these two groups of 2D problems, the solution to an original 3D problem will be reached. The numerical results show that the HRKPM performs better than RKPM in solution efficiency.

Microseismic wavefield modelling in anisotropic elastic media using integral equation method

AbstractIn this paper, we present a frequency‐domain volume integral method to model the microseismic wavefield in heterogeneous anisotropic‐elastic media. The elastic wave equation is written as an integral equation of the Lippmann–Schwinger type, and the seismic source is represented as a general moment tensor. The actual medium is split into a background medium and a scattered medium. The background part of the displacement field is computed analytically, but the scattered part requires a numerical solution. The existing matrix‐based implementation of the integral equation is computationally inefficient to model the wavefield in three‐dimensional earth. An integral equation for the particle displacement is, hence, formulated in a matrix‐free manner through the application of the Fourier transform. The biconjugate gradient stabilized method is used to iteratively obtain the solution of this equation. The integral equation method is naturally target oriented, and it is not necessary to fully discretize the model. This is very helpful in the microseismic wavefield computation at receivers in the borehole in many cases; say, for example, we want to focus only on the fluid injection zone in the reservoir–overburden system and not on the whole subsurface region. Additionally, the integral equation system matrix has a low condition number. This provides us flexibility in the selection of the grid size, especially at low frequencies for given wave velocities. Considering all these factors, we apply the numerical scheme to three different models in order of increasing geological complexity. We obtain the elastic displacement fields corresponding to the different types of moment tensor sources, which prove the utility of this method in microseismic. The generated synthetic data are intended to be used in inversion for the microseismic source and model parameters.

Multiple solutions of a nonlinear elastic beam equation supported at two ends

In this paper, several new existence theorems on pairs of solutions are obtained for the fourth order boundary value problem ${u^{(4)}}(t) = h(t,u(t))$ subject to $u(0) = u(1) = u''(0) = u''(1) = 0$. Further, if $h(t,u) = \mu f(t,u)$, we prove that there exists $ \mu ^ * > 0$ such that there is no solution to the above boundary value problem for $\mu > {\mu ^ * }$, and that there are multiple solutions of the above boundary value problem for $0 < \mu < {\mu^ * }$.

Variable-stiffness curved laminated-beams by curvilinear fibers with arbitrarily layup – Vibrational features by sine-based higher-order beam model with renewed-constitutive relations and improved-kinematics

The vibration characteristics of variable-stiffness based composite curved beams with general layups are examined using a higher-order beam theory that accounts for three-dimensional structural response features. The variable stiffness is spatially considered in the laminated beam by introducing the curvilinear fibres-based layers. The structural model is here represented by a sine-based shear deformable model considering zig-zag function for improving in-plane response and higher order polynomial for capturing through-thickness stretching effect. The appropriate constitutive relationships for laminated composite curved beam with arbitrary layup producing stiffness coupling effects are deduced from three-dimensional elasticity equations. The beam governing equations are then formed applying Hamilton’s principle and the free vibration features are predicted employing the eigenvalue analysis. The effectiveness of the curved laminated beam element developed here is tested against the available analytical solutions in the literature. A detailed numerical analysis by opting for different structural parameters like curved beam angle, thickness, ply-angle, layup, and beam end condition are made to investigate the vibration characteristics of general layup constant stiffness laminated curved beams and by including additional parameter, curvilinear fibre angle variation from centre-to-edge, for variable stiffness laminates.