The mechanical properties of bedding rock in cold regions are frequently affected by freeze-thaw(F-T)cycles and ani-sotropy.Research on the mechanical characteristics of rock damage under the combined action of F-T an...The mechanical properties of bedding rock in cold regions are frequently affected by freeze-thaw(F-T)cycles and ani-sotropy.Research on the mechanical characteristics of rock damage under the combined action of F-T and bedding angles is limited,and most traditional rock damage models cannot accurately capture these characteristics.We performed axial compression tests to ex-plore the strength characteristics of bedding slates at the bedding angles ofβ=0°,30°,45°,60°,and 90°under various F-T cycles.The experimental findings suggest that the elastic modulus and uniaxial compressive strength of the slate declined exponentially as the number of F-T cycles increased.Axial compressive strength was characterized by a U-shaped tendency with the bedding angle.This study proposes a damage model for the uniaxial compressive strength of transversely isotropic rock,which integrates the F-T effect,utilizing the enhanced anisotropic Hoek-Brown strength criterion and a statistical damage model.This model was validated using experimental data.This statistical damage model can precisely capture the dual attributes of rock mass strength reduction with F-T cy-cles and variations arising from the loading direction.展开更多
Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential...Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential. This paper presents a comprehensive analysis of infinite transversely isotropic poroelasticity under a fluid source, based on Biot's theory, aiming to uncover new and previously unexplored insights in the literature. We begin our study by deriving a general solution for fluid-saturated, transversely isotropic poroelastic materials in terms of harmonic functions that satisfy sixth-order homogeneous partial differential equations, using potential theory and Almansi's theorem. Based on these general solutions and potential functions, we construct a Green's function for a point fluid source, introducing three new harmonic functions with undetermined constants. These constants are determined by enforcing continuity and equilibrium conditions. Substituting these into the general solution yields fundamental solutions for poroelasticity that provide crucial support for a wide range of project problems. Numerical results and comparisons with existing literature are provided to illustrate physical mechanisms through contour plots. Our observations reveal that all components tend to zero in the far field and become singular at the concentrated source. Additionally, the contours exhibit rapid changes near the point fluid source but display gradual variations at a distance from it. These findings highlight the intricate behavior of the system under point liquid loading, offering valuable insights for further research and practical applications.展开更多
Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate ...Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate for actual situations.Although soil exhibits material nonlinearity under earthquake excitation,existing research on the TI medium is limited to soil linearity and neglects the nonlinear response of TI sites.A 2D equivalent linear model for a layered TI half-space subjected to seismic waves is derived in the transformed wave number domain using the exact dynamic stiffness matrix of the TI medium.This study introduces a method for determining the effective shear strain of TI sites under oblique wave incidence,and further describes a systematic study on the effects of TI parameters and soil nonlinearity on site responses.Numerical results indicate that seismic responses of the TI medium significantly differ from those of isotropic sites and that the responses are highly dependent on TI parameters,particularly in nonlinear cases,while also being sensitive to incident angle and excitation intensity.Moreover,the differences in peak acceleration and waveform for various TI materials may also be amplified due to the strong nonlinearity.The study provides valuable insights for improving the accuracy of seismic response analysis in engineering applications.展开更多
We present an assumed enhanced strain finite element framework for the simulation of tensile fracturing processes in transversely isotropic rocks.Fractures along the weak bedding planes and through the anisotropic roc...We present an assumed enhanced strain finite element framework for the simulation of tensile fracturing processes in transversely isotropic rocks.Fractures along the weak bedding planes and through the anisotropic rock matrix are treated with distinct enrichment,and a recently proposed dualmechanism tensile failure criterion for transversely isotropic rocks is adopted to determine crack initiation for the two failure modes.The cohesive crack model is adopted to characterize the response of embedded cracks.As for the numerical implementation of the proposed framework,both algorithms for the update of local history variables at Gauss points and of the global finite element system are derived.Four boundary-value problem simulations are carried out with the proposed framework,including uniaxial tension tests of Argillite,pre-notched square loaded in tension,three-point bending tests on Longmaxi shale,and simulations of tensile cracks induced by a strip load around a tunnel in transversely isotropic rocks.Simulation results reveal that the proposed framework can properly capture the tensile strength anisotropy and the anisotropic evolution of tensile cracks in transversely isotropic rocks.展开更多
Despite the significance of the high flexibility exhibited by short DNAs,there remains an incomplete understanding of their anomalous persistence length.In this study,we propose to model each basic unit of gene sequen...Despite the significance of the high flexibility exhibited by short DNAs,there remains an incomplete understanding of their anomalous persistence length.In this study,we propose to model each basic unit of gene sequences within short DNAs as a transversely isotropic ring.Our comprehensive model analysis successfully replicates the observed high flexibility of short DNAs and also displays the impact of sequence dependence,aligning with experimental findings.Furthermore,our analysis suggests that the bending behavior of short DNAs can be effectively described by the Timoshenko beam theory with the consideration of shear.展开更多
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 th...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.展开更多
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 vel...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 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...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.展开更多
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 satur...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.展开更多
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, th...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.展开更多
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 d...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 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 ...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.展开更多
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 ...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.展开更多
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 si...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.展开更多
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 wit...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.展开更多
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 st...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.展开更多
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-s...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 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 ...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 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 fo...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.展开更多
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, sand...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.展开更多
基金supported by the Qingdao Postdoctoral Science Foundation(No.862205040054)the International Research Fellowship from the Japan Society for the Promotion of Science(Postdoctoral Fellowships for Research in Japan(Standard))the National Natural Science Foundation of China(No.52078093).
文摘The mechanical properties of bedding rock in cold regions are frequently affected by freeze-thaw(F-T)cycles and ani-sotropy.Research on the mechanical characteristics of rock damage under the combined action of F-T and bedding angles is limited,and most traditional rock damage models cannot accurately capture these characteristics.We performed axial compression tests to ex-plore the strength characteristics of bedding slates at the bedding angles ofβ=0°,30°,45°,60°,and 90°under various F-T cycles.The experimental findings suggest that the elastic modulus and uniaxial compressive strength of the slate declined exponentially as the number of F-T cycles increased.Axial compressive strength was characterized by a U-shaped tendency with the bedding angle.This study proposes a damage model for the uniaxial compressive strength of transversely isotropic rock,which integrates the F-T effect,utilizing the enhanced anisotropic Hoek-Brown strength criterion and a statistical damage model.This model was validated using experimental data.This statistical damage model can precisely capture the dual attributes of rock mass strength reduction with F-T cy-cles and variations arising from the loading direction.
基金supported by the National Natural Science Foundation of China (Grant Nos. 12272269, 11972257,11832014 and 11472193)the Shanghai Pilot Program for Basic Researchthe Shanghai Gaofeng Project for University Academic Program Development。
文摘Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential. This paper presents a comprehensive analysis of infinite transversely isotropic poroelasticity under a fluid source, based on Biot's theory, aiming to uncover new and previously unexplored insights in the literature. We begin our study by deriving a general solution for fluid-saturated, transversely isotropic poroelastic materials in terms of harmonic functions that satisfy sixth-order homogeneous partial differential equations, using potential theory and Almansi's theorem. Based on these general solutions and potential functions, we construct a Green's function for a point fluid source, introducing three new harmonic functions with undetermined constants. These constants are determined by enforcing continuity and equilibrium conditions. Substituting these into the general solution yields fundamental solutions for poroelasticity that provide crucial support for a wide range of project problems. Numerical results and comparisons with existing literature are provided to illustrate physical mechanisms through contour plots. Our observations reveal that all components tend to zero in the far field and become singular at the concentrated source. Additionally, the contours exhibit rapid changes near the point fluid source but display gradual variations at a distance from it. These findings highlight the intricate behavior of the system under point liquid loading, offering valuable insights for further research and practical applications.
基金National Natural Science Foundation of China under Grant No.U2139208。
文摘Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate for actual situations.Although soil exhibits material nonlinearity under earthquake excitation,existing research on the TI medium is limited to soil linearity and neglects the nonlinear response of TI sites.A 2D equivalent linear model for a layered TI half-space subjected to seismic waves is derived in the transformed wave number domain using the exact dynamic stiffness matrix of the TI medium.This study introduces a method for determining the effective shear strain of TI sites under oblique wave incidence,and further describes a systematic study on the effects of TI parameters and soil nonlinearity on site responses.Numerical results indicate that seismic responses of the TI medium significantly differ from those of isotropic sites and that the responses are highly dependent on TI parameters,particularly in nonlinear cases,while also being sensitive to incident angle and excitation intensity.Moreover,the differences in peak acceleration and waveform for various TI materials may also be amplified due to the strong nonlinearity.The study provides valuable insights for improving the accuracy of seismic response analysis in engineering applications.
基金supported by the National Natural Science Foundation of China(Grant Nos.52038005 and 52201326)the fellowship of China Postdoctoral Science Foundation(Grant No.2022M721883)Tsinghua University Initiative Scientific Research Program.
文摘We present an assumed enhanced strain finite element framework for the simulation of tensile fracturing processes in transversely isotropic rocks.Fractures along the weak bedding planes and through the anisotropic rock matrix are treated with distinct enrichment,and a recently proposed dualmechanism tensile failure criterion for transversely isotropic rocks is adopted to determine crack initiation for the two failure modes.The cohesive crack model is adopted to characterize the response of embedded cracks.As for the numerical implementation of the proposed framework,both algorithms for the update of local history variables at Gauss points and of the global finite element system are derived.Four boundary-value problem simulations are carried out with the proposed framework,including uniaxial tension tests of Argillite,pre-notched square loaded in tension,three-point bending tests on Longmaxi shale,and simulations of tensile cracks induced by a strip load around a tunnel in transversely isotropic rocks.Simulation results reveal that the proposed framework can properly capture the tensile strength anisotropy and the anisotropic evolution of tensile cracks in transversely isotropic rocks.
基金supported by Zhejiang Provincial Natural Science Foundation of China(Grant No.LZ23A020004)the National Natural Science Foundation of China(Grant Nos.12372318 and 11872334).
文摘Despite the significance of the high flexibility exhibited by short DNAs,there remains an incomplete understanding of their anomalous persistence length.In this study,we propose to model each basic unit of gene sequences within short DNAs as a transversely isotropic ring.Our comprehensive model analysis successfully replicates the observed high flexibility of short DNAs and also displays the impact of sequence dependence,aligning with experimental findings.Furthermore,our analysis suggests that the bending behavior of short DNAs can be effectively described by the Timoshenko beam theory with the consideration of shear.
基金We would like to acknowledge all the reviewers and editors and the sponsorship of National Natural Science Foundation of China(42030103)the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(2021QNLM020001-6)the Laoshan National Laboratory of Science and Technology Foundation(LSKJ202203400).
文摘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.
基金supported by the National High-Tech Research and Development Program of China(Grant No.2006AA06Z202)the Open Fund of the Key Laboratory of Geophysical Exploration of CNPC(Grant No.GPKL0802)+1 种基金the Graduate Student Innovation Fund of China University of Petroleum(East China)(Grant No.S2008-1)the Program for New Century Excellent Talents in University(Grant No.NCET-07-0845)
文摘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 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.
文摘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.
基金The 111 Project under Grant No.B13024the National Natural Science Foundation of China under Grant Nos.U1134207 and 51378177the Program for New Century Excellent Talents in University under Grant No.NCET-12-0843
文摘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.
基金Project (Nos. 10472102 and 10432030) supported by the NationalNatural Science Foundation of China
文摘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 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.
基金National Natural Science Foundation of China under Grant Nos.51578373 and 51578372
文摘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.
基金the National Basic Research Program of China under Grant No.2013CB036405the Key Research Program of the Chinese Academy of Sciences under Grant No.KZZD-EW-05the Natural Science Foundation of China under Grant Nos.41402317,51209201 and 51279198
文摘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.
文摘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.
基金supported by the State Key Basic Research Project of China(No.2011CB201201)the National Natural Science Foundation of China(Nos.51134018 and 11172318)the Key Technologies R&D Program of China(No.2008BAB36B07)
文摘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.
基金National Natural Science Foundation of China under grant No.51578373 and 51578372the Natural Science Foundation of Tianjin Municipality under Grant No.16JCYBJC21600
文摘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 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.
基金supported by the National Natural Science Foundation of China(No.10602053)the Research Grants from Institute of Crustal Dynamics(Nos.ZDJ2007-2 and ZDJ2007-28).
文摘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.
文摘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.
