Analytical solution for the effective elastic properties of rocks with the tilted penny-shaped cracks in the transversely isotropic background

Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with the isotropic background,while the explicit model for the cracked rock with the anisotropic background is rarely investigated in spite of such case being often encountered in the earth.Hence,we first studied dependences of the crack opening displacement tensors on the crack dip angle in the coordinate systems formed by symmetry planes of the crack and the background anisotropy,respectively,by forty groups of numerical experiments.Based on the conclusion from the experiments,the analytical solution was derived for the effective elastic properties of the rock with the inclined penny-shaped cracks in the transversely isotropic background.Further,we comprehensively analyzed,according to the developed model,effects of the crack dip angle,background anisotropy,filling fluid and crack density on the effective elastic properties of the cracked rock.The analysis results indicate that the dip angle and background anisotropy can significantly either enhance or weaken the anisotropy degrees of the P-and SH-wave velocities,whereas they have relatively small effects on the SV-wave velocity anisotropy.Moreover,the filling fluid can increase the stiffness coefficients related to the compressional modulus by reducing crack compliance parameters,while its effects on shear coefficients depend on the crack dip angle.The increasing crack density reduces velocities of the dry rock,and decreasing rates of the velocities are affected by the crack dip angle.By comparing with exact numerical results and experimental data,it was demonstrated that the proposed model can achieve high-precision estimations of stiffness coefficients.Moreover,the assumption of the weakly anisotropic background results in the consistency between the proposed model and Hudson's published theory for the orthorhombic rock.

Analyzing the Form-Finding of a Large-Span Transversely Stiffened Suspended Cable System: A Method Considering Construction Processes

The precise control of the shape of transversely stiffened suspended cable systems is crucial. However, existing form-finding methods primarily rely on iterative calculations that treat loads as fixed known conditions. These methods are inefficient and fail to accurately control shape results. In this study, we propose a form-finding method that analyzes the load response of models under different sag and stress levels, taking into account the construction process. To analyze the system, a structural finite element model was established in ANSYS, and geometric nonlinear analysis was conducted using the Newton-Raphson method. The form-finding analysis results demonstrate that the proposed method achieves precise control of shape, with a maximum shape error ranging from 0.33% to 0.98%. Furthermore, the relationships between loads and tension forces are influenced by the deformed shape of the structures, exhibiting significant geometric nonlinear characteristics. Meanwhile, the load response analysis reveals that the stress level of the self-equilibrium state in the transversely stiffened suspended cable system is primarily governed by strength criteria, while shape is predominantly controlled by stiffness criteria. Importantly, by simulating the initial tensioning process as an initial condition, this method solves for a counterweight that satisfies the requirements and achieves a self-equilibrium state with the desired shape. The shape of the self-equilibrium state is precisely controlled by simulating the construction process. Overall, this work presents a new method for analyzing the form-finding process of large-span transversely stiffened suspended cable system, considering the construction process which was often overlooked in previous studies.

Poroelastodynamic responses and elastic moduli of a transversely isotropic porous cylinder under forced deformation test

Existing transversely isotropic poroelastodynamics solutions are limited to infinite domains and without experimental validation. Furthermore, there is a lack of analytical simulations for the elastic moduli dispersion of fluid-saturated porous cylinders. To address these three limitations and investigate the mechanisms of moduli dispersion, we present the analytical solutions of the poromechanical responses and the elastic moduli dispersion of a transversely isotropic, fluid-saturated, finite porous cylinder subjected to a forced deformation test. Through an example, we demonstrate the effects of loading frequency, boundary conditions, and material's anisotropy, dimension, and permeability on the responses of pore pressure,force, displacement, and dynamic elastic moduli of the cylinder. The specimen's responses are significantly influenced by the frequency of the applied load, resulting in a drained state at low frequencies and an undrained state at high frequencies. At high frequencies, the sample behaves identically for an open or a closed lateral boundary, and permeability has insignificant effects. The dynamic elastic moduli are mainly controlled by the loading frequency and the ratio of the sample's radius to its height. Lastly,we show excellent matches between the newly derived analytical solution and laboratory measurements on one clay and two shale samples from Mont Terri.

Solving elastic wave equations in 2D transversely isotropic media by a weighted Runge-Kutta discontinuous Galerkin method

Accurate wave propagation simulation in anisotropic media is important for forward modeling, migration and inversion. In this study, the weighted Runge-Kutta discontinuous Galerkin (RKDG) method is extended to solve the elastic wave equations in 2D transversely isotropic media. The spatial discretization is based on the numerical flux discontinuous Galerkin scheme. An explicit weighted two-step iterative Runge-Kutta method is used as time-stepping algorithm. The weighted RKDG method has good flexibility and applicability of dealing with undulating geometries and boundary conditions. To verify the correctness and effectiveness of this method, several numerical examples are presented for elastic wave propagations in vertical transversely isotropic and tilted transversely isotropic media. The results show that the weighted RKDG method is promising for solving wave propagation problems in complex anisotropic medium.

Effects of 3D deformation and nonlinear stress–strain relationship on the Brazilian test for a transversely isotropic rock

To improve the accuracy of indirect tensile strength for a transversely isotropic rock in the Brazilian test, this study considered the three-dimensional (3D) deformation and the nonlinear stress–strain relationship. A parametric study of a numerical Brazilian test was performed for a general range of elastic constants, revealing that the 3D modeling evaluated the indirect tensile strength up to 40% higher than the plane stress modeling. For the actual Asan gneiss, the 3D model evaluated the indirect tensile strength up to 10% higher and slightly enhanced the accuracy of deformation estimation compared with the plane stress model. The nonlinearity in stress–strain curve of Asan gneiss under uniaxial compression was then considered, such that the evaluated indirect tensile strength was affected by up to 10% and its anisotropy agreed well with the physical intuition. The estimation of deformation was significantly enhanced. The further validation on the nonlinear model is expected as future research.

Analysis of stresses at the center of transversely isotropic Brazilian disk

This article presents the stresses at the center of a Brazilian disk(BD)for transversely isotropic rocks.It is shown that the solution of stresses at the center of an anisotropic disk is a function of the disk radius and the magnitude of applied load,as well as the material orientation with respect to the load axis and two dimensionless ratios with specific physical meanings and limitations.These two dimensionless parameters are the ratios of Young’s modulus and apparent shear modulus,although the ratio of apparent shear modulus will be eliminated if the Saint-Venant assumption is considered.Considerable finite element simulations are carried out to find the stresses at the disk center concerning the material orientation and the two dimensionless parameters.Also,an approximate formula obtained from analytical results,previously proposed in the literature for solving the tensile and compressive stresses at the disk center,is re-written and simplified based on these new definitions.The results of the approximate formula fitted to the analytical results are compared to those obtained from numerical solutions,suggesting a good agreement between the numerical and analytical methods.An approximate equation for the shear stress at the disk center is also formulated based on the numerical results.Finally,the influence of the assumptions for simplification of the proposed formula for the tensile,compressive,and shear stresses at the disk center is discussed,and simple and practical equations are proposed as estimations for the stresses at the center of the BD specimen for low to moderate anisotropic rocks.For highly anisotropic rocks,the reference plots can be used for more accuracy.

Numerical modeling of seismic wavefields in transversely isotropic media with a compact staggered-grid finite difference scheme

To deal with the numerical dispersion problem, by combining the staggeredgrid technology with the compact finite difference scheme, we derive a compact staggered- grid finite difference scheme from the first-order velocity-stress wave equations for the transversely isotropic media. Comparing the principal truncation error terms of the compact staggered-grid finite difference scheme, the staggered-grid finite difference scheme, and the compact finite difference scheme, we analyze the approximation accuracy of these three schemes using Fourier analysis. Finally, seismic wave numerical simulation in transversely isotropic （VTI） media is performed using the three schemes. The results indicate that the compact staggered-grid finite difference scheme has the smallest truncation error, the highest accuracy, and the weakest numerical dispersion among the three schemes. In summary, the numerical modeling shows the validity of the compact staggered-grid finite difference scheme.

The reflection and transmission of elastic waves at a interface between a transversely isotropi celastic soils and saturated soils

According to Biot′s wave equation of transversely isotropic saturated soil, this paper deduces the general equation of the reflection coefficients and transmission coefficients when qP 1 wave goes through from saturated soil to elastic media. The effects of anisotropies and boundary drainage condition on reflection coefficients and transmission coefficients are analyzed by numerical method. The idea of this paper can be applied to the case when qSV wave or qP 2 wave goes through from saturated soils to elastic soils.

Torsional vibration of a pipe pile in transversely isotropic saturated soil

This study considers the torsional vibration of a pipe pile in a transversely isotropic saturated soil layer. Based on Biot＇s poroelastic theory and the constitutive relations of the transversely isotropic medium, the dynamic governing equations of the outer and inner transversely isotropic saturated soil layers are derived. The Laplace transform is used to solve the governing equations of the outer and inner soil layers. The dynamic torsional response of the pipe pile in the frequency domain is derived utilizing 1D elastic theory and the continuous conditions at the interfaces between the pipe pile and the soils. The time domain solution is obtained by Fourier inverse transform. A parametric study is conducted to demonstrate the influence of the anisotropies of the outer and inner soil on the torsional dynamic response of the pipe pile.

Pure bending of simply supported circular plate of transversely isotropic functionally graded material

This paper considers the pure bending problem of simply supported transversely isotropic circular plates with elastic compliance coefficients being arbitrary functions of the thickness coordinate. First, the partial differential equation, which is satisfied by the stress functions for the axisymmetric deformation problem is derived. Then, stress functions are obtained by proper manipulation. The analytical expressions of axial force, bending moment and displacements are then deduced through integration. And then, stress functions are employed to solve problems of transversely isotropic functionally graded circular plate, with the integral constants completely determined from boundary conditions. An elasticity solution for pure bending problem, which coincides with the available solution when degenerated into the elasticity solutions for homogenous circular plate, is thus obtained. A numerical example is finally presented to show the effect of material inhomogeneity on the elastic field in a simply supported circular plate of transversely isotropic functionally graded material (FGM).

THE NON-AXISYMMETRICAL DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED POROELASTIC MEDIA

The Blot's wave equations of transversely isotropic saturated poroelastic media excited hy non-axisymmetrical harmonic source were solved by means of Fourier expansion and Hankel transform. Then the components of total stress in porous media are expressed with the solutions of Biot's wave equations. The method of research on non-axisymmetrical dynamic response of saturated porous media is discussed, and a numerical result is presented.

Diffraction of SH-waves by topographic features in a layered transversely isotropic half-space

Abstract： The scattering of plane SH-waves by topographic features in a layered transversely isotropic （TI） half-space is investigated by using an indirect boundary element method （IBEM）. Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green＇s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green＇s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories.

Dynamic analysis of slab track on multi-layered transversely isotropic saturated soils subjected to train loads

The dynamic responses of a slab track on transversely isotropic saturated soils subjected to moving train loads are investigated by a semi-analytical approach. The track model is described as an upper Euler beam to simulate the rails and a lower Euler beam to model the slab. Rail pads between the rails and slab are represented by a continuous layer of springs and dashpots. A series of point loads are formulated to describe the moving train loads. The governing equations of track-ground systems are solved using the double Fourier transform, and the dynamic responses in the time domain are obtained by the inverse Fourier transform. The results show that a train load with high velocity will generate a larger response in transversely isotropic saturated soil than the lower velocity load, and special attention should be paid on the pore pressure in the vicinity of the ground surface. The anisotropic parameters of a surface soil layer will have greater influence on the displacement and excess pore water pressure than those of the subsoil layer. The traditional design method taking ground soil as homogeneous isotropic soil is unsafe for the case of RE 〈 1 and RG 〈 1, so a transversely isotropic foundation model is of great significance to the design for high train velocities.

Nonlinear amplitude versus angle inversion for transversely isotropic media with vertical symmetry axis using new weak anisotropy approximation equations

In VTI media,the conventional inversion methods based on the existing approximation formulas are difficult to accurately estimate the anisotropic parameters of reservoirs,even more so for unconventional reservoirs with strong seismic anisotropy.Theoretically,the above problems can be solved by utilizing the exact reflection coefficients equations.However,their complicated expression increases the difficulty in calculating the Jacobian matrix when applying them to the Bayesian deterministic inversion.Therefore,the new reduced approximation equations starting from the exact equations are derived here by linearizing the slowness expressions.The relatively simple form and satisfactory calculation accuracy make the reduced equations easy to apply for inversion while ensuring the accuracy of the inversion results.In addition,the blockiness constraint,which follows the differentiable Laplace distribution,is added to the prior model to improve contrasts between layers.Then,the concept of GLI and an iterative reweighted least-squares algorithm is combined to solve the objective function.Lastly,we obtain the iterative solution expression of the elastic parameters and anisotropy parameters and achieve nonlinear AVA inversion based on the reduced equations.The test results of synthetic data and field data show that the proposed method can accurately obtain the VTI parameters from prestack AVA seismic data.

Deformation analysis of transversely isotropic coal-rock mass with porous and cracks

Coal-rock as a typical sedimentary rock has obvious stratification,namely it has transversely isotropic feature.Meanwhile,deformation leads to coal-rock mass having the characteristics of different porous and crack structures as well as local anisotropy.Equivalent axial and circumferential strain' formulas of the pure coal-rock mass specimen with a single crack were derived through the establishment of equivalent mechanical model of standard cylindrical coal-rock specimen,and have been widely used to a variety of media combined different structures containing multiple cracks.The complete stress strain curve of a real coal-rock specimen was obtained by the CTC test.Additionally,according to the comparison with the theoretical value,the theoretical mechanical model could well explain the deformation characteristics of coal-rock mass and verify its validity.Further,following features were analyzed:strain normalized coefficient and elastic modulus(Poisson's ratio) in vertical and parallel direction to the stratification,stratification angle,porosity,pore radius,normal and tangential stiffness of crack,and the relationship of different crack width with different tangential stiffness of crack.Through the analysis above,it substantiate this claim that the theoretical model with better reliability reflects the transversely isotropic nature of the coal-rock and the local anisotropy caused by the porous and cracks.

3-D DYNAMIC RESPONSE OF TRANSVERSELY ISOTROPIC SATURATED SOILS

A study on dynamic response of transversely isotropic saturated poroelastic media under a circular non-axisymmetrical harmonic source has been presented by Huang Yi et al. using the technique of Fourier expansion and Hankel transform. However, the method may not always be valid. The work is extended to the general case being in the rectangular coordinate. The purpose is to study the 3-d dynamic response of transversely isotropic saturated soils under a general source distributing in arbitrary rectangular zoon on the medium surface. Based on Biot＇s theory for fluid- saturated porous media, the 3-d wave motion equations in rectangular coordinate for transversely isotropic saturated poroelastic media were transformed into the two uncoupling governing differential equations of 6-order and 2-order respectively by means of the displacement functions. Then, using the technique of double Fourier transform, the governing differential equations were easily solved. Integral solutions of soil skeleton displacements and pore pressure as well as the total stresses for poroelastic media were obtained. Furthermore, a systematic study on half-space problem in saturated soils was performed. Integral solutions for surface displacements under the general harmonic source distributing on arbitrary surface zone, considering both case of drained surface and undrained surface, were presented.

DYNAMIC CONTACT STIFFNESS OF VIBRATING RIGID SPHERE CONTACTING SEMI-INFINITE TRANSVERSELY ISOTROPIC VISCOELASTIC SOLID

Dynamic contact stiffness at the interface between a vibrating rigid sphere and a semi-infinite transversely isotropic viscoelastic solid is investigated. An oscillating force superimposed onto a static compressive force in the vertical direction excites the vibration of a rigid sphere, which causes variable contact radius and contact pressure distribution in the contact region. The assumption of a sufficiently small oscillating force yields a dynamic contact-pressure distribution of a constant contact radius, which gives dynamic contact stiffness at the interface between the rigid sphere and the semi-infinite solid. Numerical calculations show the influence of vibration frequency of the sphere, and elastic constants of the transversely isotropic solid on dynamic contact stiffness, which benefits quantitative evaluation of elastic constants and orientation of single hexagonal grains by resonance-frequency shifts of the oscillator in resonance ultrasound microscopy.

In-plane dynamic Green's functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic half-space

The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic（TI）half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf（1985）,and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.

A modified failure criterion for transversely isotropic rocks

A modified failure criterion is proposed to determine the strength of transversely isotropic rocks. Me-chanical properties of some metamorphic and sedimentary rocks including gneiss, slate, marble, schist, shale, sandstone and limestone, which show transversely isotropic behavior, were taken into consider-ation. Afterward, introduced triaxial rock strength criterion was modified for transversely isotropic rocks. Through modification process an index was obtained that can be considered as a strength reduction parameter due to rock strength anisotropy. Comparison of the parameter with previous anisotropy in-dexes in literature showed reasonable results for the studied rock samples. The modified criterion was compared to modified Hoek-Brown and Ramamurthy criteria for different transversely isotropic rocks. It can be concluded that the modified failure criterion proposed in this study can be used for predicting the strength of transversely isotropic rocks.

A PENNY-SHAPED CRACK IN A TRANSVERSELY ISOTROPIC PIEZOELECTRIC SOLID:MODES Ⅱ AND Ⅲ PROBLEMS

An exact analysis of the modes Ⅱ and Ⅲ problems of a penny- shaped crack in a transversely isotropic piezoelectric medium is performed in this paper.The potential theory method is employed based on the general solution of three-dimensional piezoelasticity and the four harmonics involved are represented by one complex potential.Previous results in potential theory are then utilized to obtain the exact solution that is expressed in terms of elementary functions.Comparison is made between the current results with those published and good agreement is obtained.