Monitoring the change in horizontal stress with multi-wave time-lapse seismic response based on nonlinear elasticity theory

Monitoring the change in horizontal stress from the geophysical data is a tough challenge, and it has a crucial impact on broad practical scenarios which involve reservoir exploration and development, carbon dioxide (CO_(2)) injection and storage, shallow surface prospecting and deep-earth structure description. The change in in-situ stress induced by hydrocarbon production and localized tectonic movements causes the changes in rock mechanic properties (e.g. wave velocities, density and anisotropy) and further causes the changes in seismic amplitudes, phases and travel times. In this study, the nonlinear elasticity theory that regards the rock skeleton (solid phase) and pore fluid as an effective whole is used to characterize the effect of horizontal principal stress on rock overall elastic properties and the stress-dependent anisotropy parameters are therefore formulated. Then the approximate P-wave, SV-wave and SH-wave angle-dependent reflection coefficient equations for the horizontal-stress-induced anisotropic media are proposed. It is shown that, on the different reflectors, the stress-induced relative changes in reflectivities (i.e., relative difference) of elastic parameters (i.e., P- and S-wave velocities and density) are much less than the changes in contrasts of anisotropy parameters. Therefore, the effects of stress change on the reflectivities of three elastic parameters are reasonably neglected to further propose an AVO inversion approach incorporating P-, SH- and SV-wave information to estimate the change in horizontal principal stress from the corresponding time-lapse seismic data. Compared with the existing methods, our method eliminates the need for man-made rock-physical or fitting parameters, providing more stable predictive power. 1D test illustrates that the estimated result from time-lapse P-wave reflection data shows the most reasonable agreement with the real model, while the estimated result from SH-wave reflection data shows the largest bias. 2D test illustrates the feasibility of the proposed inversion method for estimating the change in horizontal stress from P-wave time-lapse seismic data.

NONLINEAR RESPONSES OF A FLUID-CONVEYING PIPE EMBEDDED IN NONLINEAR ELASTIC FOUNDATIONS

The nonlinear responses of planar motions of a fluid-conveying pipe embedded in nonlinear elastic foundations are investigated via the differential quadrature method discretization （DQMD） of the governing partial differential equation. For the analytical model, the effect of the nonlinear elastic foundation is modeled by a nonlinear restraining force. By using an iterative algorithm, a set of ordinary differential dynamical equations derived from the equation of motion of the system are solved numerically and then the bifurcations are analyzed. The numerical results, in which the existence of chaos is demonstrated, are presented in the form of phase portraits of the oscillations. The intermittency transition to chaos has been found to arise.

Implementing the Node Based Smoothed Finite Element Method as User Element in Abaqus for Linear and Nonlinear Elasticity

In this paper,the node based smoothed-strain Abaqus user element(UEL)in the framework of finite element method is introduced.The basic idea behind of the node based smoothed finite element(NSFEM)is that finite element cells are divided into subcells and subcells construct the smoothing domain associated with each node of a finite element cell[Liu,Dai and Nguyen-Thoi(2007)].Therefore,the numerical integration is globally performed over smoothing domains.It is demonstrated that the proposed UEL retains all the advantages of the NSFEM,i.e.,upper bound solution,overly soft stiffness and free from locking in compressible and nearly-incompressible media.In this work,the constant strain triangular(CST)elements are used to construct node based smoothing domains,since any complex two dimensional domains can be discretized using CST elements.This additional challenge is successfully addressed in this paper.The efficacy and robustness of the proposed work is obtained by several benchmark problems in both linear and nonlinear elasticity.The developed UEL and the associated files can be downloaded from https://github.com/nsundar/NSFEM.

PERTURBATION ANALYSIS FOR WAVE EQUATION OF NONLINEAR ELASTIC ROD

The longitudinal oscillation of a nonlinear elastic rod with lateral inertia was studied. Based on the far field and simple wave theory, a nonlinear SchrSdinger （NLS） equation was established under the assumption of small amplitude and long wavelength. It is found that there are NLS envelop solitons in this system. Finally the soliton solution of the NLS equation was presented.

Symmetry solutions of a nonlinear elastic wave equation with third-order anharmonic corrections

Lie symmetry method is applied to analyze a nonlinear elastic wave equation for longitudinal deformations with third-order anharmonic corrections to the elastic energy. Symmetry algebra is found and reductions to second-order ordinary differential equations （ODEs） are obtained through invariance under different symmetries. The reduced ODEs are further analyzed to obtain several exact solutions in an explicit form. It was observed in the literature that anharmonic corrections generally lead to solutions with time-dependent singularities in finite times singularities, we also obtain solutions which Along with solutions with time-dependent do not exhibit time-dependent singularities.

Tunable three-dimensional nonreciprocal transmission in a layered nonlinear elastic wave metamaterial by initial stresses

In this work,the three-dimensional(3 D)propagation behaviors in the nonlinear phononic crystal and elastic wave metamaterial with initial stresses are investigated.The analytical solutions of the fundamental wave and second harmonic with the quasilongitudinal(qP)and quasi-shear(qS_(1) and qS_(2))modes are derived.Based on the transfer and stiffness matrices,band gaps with initial stresses are obtained by the Bloch theorem.The transmission coefficients are calculated to support the band gap property,and the tunability of the nonreciprocal transmission by the initial stress is discussed.This work is expected to provide a way to tune the nonreciprocal transmission with vector characteristics.

A MICROMECHANICAL ANALYSIS OF THE NONLINEAR ELASTIC AND VISCOELASTIC CONSTITUTIVE RELATION OF A POLYMER FILLED WITH RIGID PARTICLES

Micromechanical theory is applied to study the nonlinear elastic and viscoelastic constitutive relations of polymeric matrix filled with high rigidity solid particles. It is shown that Eshelby's method can be extended to the case of nonlinear matrix and Eshelby's tensor still exists provided that Poisson's ratio of the nonlinear matrix assumes constant value in deforming process and the rigidity of elastic filling particles is much higher than that of the matrix. A new method for averaging process is proposed to overcome the difficulty that occured in applying the ordinary equivalent inclusion method or the seff-consistant method to nonlinear matrices. A rather simple constitutive equation is obtained finally and the strengthening effect of solid particles to composites is investigated.

Collision enhanced hyper-damping in nonlinear elastic metamaterial

Nonlinear elastic metamaterial,a topic which has attracted extensive attention in recent years,can enable broadband vibration reduction under relatively large amplitude.The combination of damping and strong nonlinearity in metamaterials may entail extraordinary effects and offer the capability for low-frequency and broadband vibration reduction.However,there exists a clear lack of proper design methods as well as the deficiency in understanding properties arising from this concept.To tackle this problem,this paper numerically demonstrates that the nonlinear elastic metamaterials,consisting of sandwich damping layers and collision resonators,can generate very robust hyper-damping effect,conducive to efficient and broadband vibration suppression.The collision-enhanced hyper damping is persistently presented in a large parameter space,ranging from small to large amplitudes,and for small and large damping coefficients.The achieved robust effects greatly enlarge the application scope of nonlinear metamaterials.We report the design concept,properties and mechanisms of the hyper-damping and its effect on vibration transmission.This paper reveals new properties offered by nonlinear elastic metamaterials,and offers a robust method for achieving efficient low-frequency and broadband vibration suppression.

A RATE TYPE METHOD FOR LARGE DEFORMATION PROBLEMS OF NONLINEAR ELASTICITY

In this paper. we obtain the rate-tvpe constitutive expressions of the nonlinearisotropic elastieity by using the Jaumann, Truesdell and Green-Naghdi stress rgterespectively,Through analysing the simple shear deformation for Mooney-Rivlin material.three kinds of rate-type constitutive equations are verified to be equivalent to the original equation.Rate-type variational prineiples are also presented, and the Ritzmethod is used to obtain the numerical solution of a reetangular rubber membraneunder uniaxial streteh.

Seismic fluid identification using a nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo method

Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo （MCMC） method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient （PCG） algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.

Analytic Static Deflection Solutions of Uniform Cantilever Beams Resting on Nonlinear Elastic Rotational Boundary

In the paper,the analytic static deflection solutions of uniform cantilever beams resting on nonlinear elastic rotational boundary are developed by the Modified Adomian Decomposition Method(MADM).If the applied force function is an analytic function,then the deflection function can be derived and expressed in Maclaurin series.A recurrence relation for the coefficients of the Maclaurin series is derived.It is shown that the proposed solution method is accurate and efficient.The solution method can be successfully applied to the uniform cantilever beam and non-linear elastic rotational boundary problem.

Nonlinear phenomena in vibrations of embedded carbon nanotubes conveying viscous fluid

Various nonlinear phenomena such as bifurcations and chaos in the responses of carbon nanotubes(CNTs)are recognized as being major contributors to the inaccuracy and instability of nanoscale mechanical systems.Therefore,the main purpose of this paper is to predict the nonlinear dynamic behavior of a CNT conveying viscousfluid and supported on a nonlinear elastic foundation.The proposed model is based on nonlocal Euler–Bernoulli beam theory.The Galerkin method and perturbation analysis are used to discretize the partial differential equation of motion and obtain the frequency-response equation,respectively.A detailed parametric study is reported into how the nonlocal parameter,foundation coefficients,fluid viscosity,and amplitude and frequency of the external force influence the nonlinear dynamics of the system.Subharmonic,quasi-periodic,and chaotic behaviors and hardening nonlinearity are revealed by means of the vibration time histories,frequency-response curves,bifurcation diagrams,phase portraits,power spectra,and Poincarémaps.Also,the results show that it is possible to eliminate irregular motion in the whole range of external force amplitude by selecting appropriate parameters.

Solution of nonlinear wave equation of elastic rod

The longitudinal oscillation of a nonlinear elastic rod with lateral inertia are studied. A nonlinear wave equation is derived. The equation is solved by the method of full approximation.

Derivation of the Orthotropic Nonlinear Elastic Material Law Driven by Low-Cost Data(DDONE)

Orthotropic nonlinear elastic materials are common in nature and widely used by various industries.However,there are only limited constitutive models available in today's commercial software(e.g.,ABAQUS,ANSYS,etc.)that adequately describe their mechanical behavior.Moreover,the material parameters in these constitutive models are also difficult to calibrate through low-cost,widely available experimental setups.Therefore,it is paramount to develop new ways to model orthotropic nonlinear elastic materials.In this work,a data-driven orthotropic nonlinear elastic(DDONE)approach is proposed,which builds the constitutive response from stress–strain data sets obtained from three designed uniaxial tensile experiments.The DDONE approach is then embedded into a finite element(FE)analysis framework to solve boundary-value problems(BVPs).Illustrative examples(e.g.,structures with an orthotropic nonlinear elastic material)are presented,which agree well with the simulation results based on the reference material model.The DDONE approach generally makes accurate predictions,but it may lose accuracy when certain stress–strain states that appear in the engineering structure depart significantly from those covered in the data sets.Our DDONE approach is thus further strengthened by a mapping function,which is verified by additional numerical examples that demonstrate the effectiveness of our modified approach.Moreover,artificial neural networks(ANNs)are employed to further improve the computational efficiency and stability of the proposed DDONE approach.

The Time-Dependent Von Kármán Shell Equation as a Limit of Three-Dimensional Nonlinear Elasticity

The asymptotic behaviour of solutions of three-dimensional nonlinear elastodynamics in a thin shell is considered,as the thickness h of the shell tends to zero.Given the appropriate scalings of the applied force and of the initial data in terms of h,it’s verified that three-dimesional solutions of the nonlinear elastodynamic equations converge to solutions of the time-dependent von Kármán equations or dynamic linear equations for shell of arbitrary geometry.

Multiple bifurcations and local energy minimizers in thermoelastic martensitic transformations

Thermoelastic martensitic transformations in shape memory alloys can be modeled on the basis of nonlinear elastic theory.Microstructures of fine phase mixtures are local energy minimizers of the total energy.Using a one-dimensional effective model,we have shown that such microstructures are inhomogeneous solutions of the nonlinear Euler-Lagrange equation and can appear upon loading or unloading to certain critical conditions,the bifurcation conditions.A hybrid numerical method is utilized to calculate the inhomogeneous solutions with a large number of interfaces.The characteristics of the solutions are clarified by three parameters：the number of interfaces,the interface thickness,and the oscillating amplitude.Approximated analytical expressions are obtained for the interface and inhomogeneity energies through the numerical solutions.

Multi-crack imaging using nonclassical nonlinear acoustic method

Solid materials with cracks exhibit the nonclassical nonlinear acoustical behavior. The micro-defects in solid materials can be detected by nonlinear elastic wave spectroscopy （NEWS） method with a time-reversal （TR） mirror. While defects lie in viscoelastic solid material with different distances from one another, the nonlinear and hysteretic stress-strain relation is established with Preisach-Mayergoyz （PM） model in crack zone. Pulse inversion （PI） and TR methods are used in numerical simulation and defect locations can be determined from images obtained by the maximum value. Since false-positive defects might appear and degrade the imaging when the defects are located quite closely, the maximum value imaging with a time window is introduced to analyze how defects affect each other and how the fake one occurs. Furthermore, NEWS-TR- NEWS method is put forward to improve NEWS-TR scheme, with another forward propagation （NEWS） added to the existing phases （NEWS and TR）. In the added phase, scanner locations are determined by locations of all defects imaged in previous phases, so that whether an imaged defect is real can be deduced. NEWS-TR-NEWS method is proved to be effective to distinguish real defects from the false-positive ones. Moreover, it is also helpful to detect the crack that is weaker than others during imaging procedure.

Nonlinear dynamic characteristics of piles embedded in rock

The nonlinear dynamic characteristics of a pile embedded in a rock were investigated. Suppose that both the materials of the pile and the soil around the pile obey nonlinear elastic and linear viscoelastic constitutive relations. The nonlinear partial differential equation governing the dynamic characteristics of the pile was first derived. The Galerkin method was used to simplify the equation and to obtain a nonlinear ordinary differential equation. The methods in nonlinear dynamics were employed to solve the simplified dynamical system, and the time-path curves, phase-trajectory diagrams, power spectrum, Poincare sections and bifurcation and chaos diagrams of the motion of the pile were obtained. The effects of parameters on the dynamic characteristics of the system were also considered in detail.

Nonlinear transverse vibrations of a slightly curved beam carrying a concentrated mass

In this study, a slightly curved Euler Bernoulli beam carrying a concentrated mass was handled. The beam was resting on an elastic foundation and simply supported at both ends. Effects of the concentrated mass on nonlin- ear vibrations were investigated. Sinusoidal and parabolic type functions were used as curvature functions. Equations of motion have cubic nonlinearities because of elongations during vibrations. Damping and harmonic excitation terms were added to the equations of motion. Method of mul- tiple scales, a perturbation technique, was used for solving integro-differential equation analytically. Natural frequen- cies were calculated exactly for different mass ratios, mass locations, curvature functions, and linear elastic foundation coefficients. Amplitude-phase modulation equations were found by considering primary resonance case. Effects of nonlinear terms on natural frequencies were calculated. Frequency-amplitude and frequency-response graphs were plotted. Finally effects of concentrated mass and chosen curvature function on nonlinear vibrations were investigated.

ON OUT-OF-PLANE BUCKLING OF ANISOTROPIC ELASTIC PLATE SUBJECTED TO A SIMPLE SHEAR

Out-of-plane buckling of anisotropic elastic plate subjected to asimple shear is investigated. From exact 3-D equilibrium conditionsof anisotropic elastic body with a plane of elastic symmetry atcritical Configuration, the equation for buckling direction (bucklingwave direction) parameter is derived and the Shape functions ofpossible buckling modes are obtained. The traction free boundaryconditions which must Hold on the upper and lower surfaces of platelead to a linear eigenvalue problem whose nontrivial solutions Arejust the possible buckling modes for the plate.