The testing of large structures is limited by high costs and long cycles, making scaling methods an attractive solution. However, the scaling process of elastic rings introduces complexities in multi-parameter geometr...The testing of large structures is limited by high costs and long cycles, making scaling methods an attractive solution. However, the scaling process of elastic rings introduces complexities in multi-parameter geometric distortions, leading to a diminution in the predictive accuracy of the distorted similitude. To address this challenge, this study formulates a novel set of scaling laws, tailored to account for the intricate geometric distortions associated with elastic rings. The proposed scaling laws are formulated based on the intrinsic deformation characteristics of elastic rings, rather than the traditional systemic governing equations. Numerical and experimental cases are conducted to assess the efficacy and precision of the proposed scaling laws, and the obtained results are compared with those achieved by traditional methods. The outcomes demonstrate that the scaling laws put forth by this study significantly enhance the predictive capabilities for deformations of elastic rings.展开更多
Uncertain parameters are widespread in engineering systems.This study investigates the modal analysis of a fluid-conveying pipe subjected to elastic supports with unknown-but-bound parameters.The governing equation fo...Uncertain parameters are widespread in engineering systems.This study investigates the modal analysis of a fluid-conveying pipe subjected to elastic supports with unknown-but-bound parameters.The governing equation for the elastically supported fluid-conveying pipe is transformed into ordinary differential equations using the Galerkin truncation method.The Chebyshev interval approach,integrated with the assumed mode method is then used to investigate the effects of uncertainties of support stiffness,fluid speed,and pipe length on the natural frequencies and mode shapes of the pipe.Additionally,both symmetrical and asymmetrical support stiffnesses are discussed.The accuracy and effectiveness of the Chebyshev interval approach are verified through comparison with the Monte Carlo method.The results reveal that,for the same deviation coefficient,uncertainties in symmetrical support stiffness have a greater impact on the first four natural frequencies than those of the asymmetrical one.There may be significant differences in the sensitivity of natural frequencies and mode shapes of the same order to uncertain parameters.Notably,mode shapes susceptible to uncertain parameters exhibit wider fluctuation intervals near the elastic supports,requiring more attention.展开更多
The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using in...The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.展开更多
This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solv...This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solving a specific numerical problem under the scope of the linear finite element method(LFEM),so the method is termed computational method for analytical solutions with finite elements(CMAS-FE).The primary objective of the CMAS-FE is to construct analytical expressions for displacements and reaction forces at nodes,as well as for strains and stresses at elemental quadrature points,all of which are formulated as infinite series solutions of various orders of Poisson’s ratios.Like the conventional LFEM,the CMAS-FE forms global sparse linear equations,but the Young’s modulus and Poisson’s ratio remain variables(or symbols).By employing a direct inverse method to solve these symbolic linear systems,an analytical expression of the displacement field can be constructed.The CMAS-FE is validated via patch and bending tests,which demonstrate convergence with mesh and term refine-ment.Furthermore,the CMAS-FE is applied to obtain the bending stiffness of a beam structure and to estimate an approximate stress intensity factor for a straight crack within a square-shaped plate.展开更多
Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods rema...Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.展开更多
The two-phase flow in porous media is affected by multiple factors.In the present study,a two-dimensional numerical model of porous media was developed using the actual pore structure of the core sample.The phase fiel...The two-phase flow in porous media is affected by multiple factors.In the present study,a two-dimensional numerical model of porous media was developed using the actual pore structure of the core sample.The phase field method was utilized to simulate the impact of displacement velocity,the water-gas viscosity ratio,and the density ratio on the flow behavior of two-phase fluids in porous media.The effectiveness of displacement was evaluated by analyzing CO_(2)saturation levels.The results indicate that the saturation of CO_(2)in porous media increased as the displacement velocity increased.When the displacement velocity exceeded 0.01 m/s,there was a corresponding increase in CO_(2)saturation.Conversely,when the displacement velocity was below this threshold,the impact on CO_(2)saturation was minimal.An“inflection point,”M3,was present in the viscosity ratio.When the viscosity of CO_(2)is less than 8.937×10^(-5)Pa·s(viscosity ratio below M3),variations in the viscosity of CO_(2)had little impact on its saturation.Conversely,when the viscosity of CO_(2)exceeded 8.937×10^(-5)Pa·s(viscosity ratio greater than M3),saturation increased with an increase in the viscosity ratio.In terms of the density ratio,the saturation of CO_(2)increased monotonically with an increase in the density ratio.Similarly,increasing density ratios resulted in a monotonic increase in CO_(2)saturation,though this trend was less pronounced in numerical simulations.Analysis results of displacement within dead-end pores using pressure and velocity diagrams reveal eddy currents as contributing factors.Finally,the impact of pore throat structure on the formation of dominant channels was examined.展开更多
This paper investigates the active traveling wave vibration control of an elastic supported rotating porous aluminium conical shell(CS)under impact loading.Piezoelectric smart materials in the form of micro fiber comp...This paper investigates the active traveling wave vibration control of an elastic supported rotating porous aluminium conical shell(CS)under impact loading.Piezoelectric smart materials in the form of micro fiber composites(MFCs)are used as actuators and sensors.To this end,a metal pore truncated CS with MFCs attached to its surface is considered.Adding artificial virtual springs at two edges of the truncated CS achieves various elastic supported boundaries by changing the spring stiffness.Based on the first-order shear deformation theory(FSDT),minimum energy principle,and artificial virtual spring technology,the theoretical formulations considering the electromechanical coupling are derived.The comparison of the natural frequency of the present results with the natural frequencies reported in previous literature evaluates the accuracy of the present approach.To study the vibration control,the integral quadrature method in conjunction with the differential quadrature approximation in the length direction is used to discretize the partial differential dynamical system to form a set of ordinary differential equations.With the aid of the velocity negative feedback method,both the time history and the input control voltage on the actuator are demonstrated to present the effects of velocity feedback gain,pore distribution type,semi-vertex angle,impact loading,and rotational angular velocity on the traveling wave vibration control.展开更多
With the development of seismic engineering and seismic exploration of energy, the underground media that westudy are more and more complicated. Conventional anisotropy theory or two-phase isotropy theory is difficult...With the development of seismic engineering and seismic exploration of energy, the underground media that westudy are more and more complicated. Conventional anisotropy theory or two-phase isotropy theory is difficult todescribe anisotropic media containing fluid, such as fractures containing gas, shales containing water Based onBlot theory about two-phase anisotropy, with the use of elastic plane wave equations, we get Christoffel equations.We calculate and analyze the effects of frequency on phase velocity, attenuation, amplitude ratio and polarizationdirection of elastic waves of two-phase, transversely isotropic media. Results show that frequency affects slow Pwave the greatest among the four kinds of waves, i.e., fast P wave, slow P wave, fast S wave and slow S wave.Fluid phase amplitude to solid phase amplitude ratio of fast P wave, fast S wave and slow S wave approaches unitfor large dissipation coefficients. Polarization analysis shows that polarization direction of fluid phase displacement is different from, not parallel to or reverse to, that of solid phase displacement in two-phase anisotropic media.展开更多
A nonlocal study of the vibration responses of functionally graded(FG)beams supported by a viscoelastic Winkler-Pasternak foundation is presented.The damping responses of both the Winkler and Pasternak layers of the f...A nonlocal study of the vibration responses of functionally graded(FG)beams supported by a viscoelastic Winkler-Pasternak foundation is presented.The damping responses of both the Winkler and Pasternak layers of the foundation are considered in the formulation,which were not considered in most literature on this subject,and the bending deformation of the beams and the elastic and damping responses of the foundation as nonlocal by uniting the equivalently differential formulation of well-posed strain-driven(ε-D)and stress-driven(σ-D)two-phase local/nonlocal integral models with constitutive constraints are comprehensively considered,which can address both the stiffness softening and toughing effects due to scale reduction.The generalized differential quadrature method(GDQM)is used to solve the complex eigenvalue problem.After verifying the solution procedure,a series of benchmark results for the vibration frequency of different bounded FG beams supported by the foundation are obtained.Subsequently,the effects of the nonlocality of the foundation on the undamped/damping vibration frequency of the beams are examined.展开更多
Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical sol...Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.展开更多
The bending of rectangular plate is divided into the generalized statically determinate bending and the generalized statically indeterminate bending based on the analysis of the completeness of calculating condition a...The bending of rectangular plate is divided into the generalized statically determinate bending and the generalized statically indeterminate bending based on the analysis of the completeness of calculating condition at the corner point. The former can be solved directly by the equilibrium differential equation and the boundary conditions of four edges of the plate. The latter can be solved by using the superposition principle. Making use of the recommended method, the bending of the plate with all kinds of...展开更多
The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data...The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.展开更多
A new strategy for elastic modulus adjustment is proposed based on the element bearing ratio （EBR）,and the elastic modulus reduction method （EMRM） is proposed for limit load evaluation of frame structures. The EBR...A new strategy for elastic modulus adjustment is proposed based on the element bearing ratio （EBR）,and the elastic modulus reduction method （EMRM） is proposed for limit load evaluation of frame structures. The EBR is defined employing the generalized yield criterion,and the reference EBR is determined by introducing the extrema and the degree of uniformity of EBR in the structure. The elastic modulus in the element with an EBR greater than the reference one is reduced based on the linear elastic finite element analysis and the equilibrium of strain energy. The lower-bound of limit-loads of frame structures are analyzed and the numerical example demonstrates the flexibility,accuracy and effciency of the proposed method.展开更多
With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their conta...With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their contact status will deteriorate under cyclic loading, resulting in ballast degradation. Discrete element method(DEM) was used to research improved performance of ballast bed using elastic sleeper. Clusters were generated by bonding spheres to model real ballasts, while broken bonds were utilized to distinguish breakage. Two kinds of ballast beds with elastic sleeper and conventional sleeper were established, respectively. After applying cyclic loading to the models, differences of mechanical properties between two models were analyzed by contrasting their dynamic behavior indexes, such as particle contact force, sleeper settlement, vibration velocity and acceleration, breakage characteristic. The results illustrate that compared with conventional sleeper, elastic sleeper increases sleeper settlement, while reduces ballast vibration and contact force between particles, which could depress ballast breakage.展开更多
In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squar...In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares （CVMLS） approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin （EFG） method and the ICVMLS approximation, the improved complex variable element-free Galerkin （ICVEFG） method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.展开更多
A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the co...A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.展开更多
Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
It is an important subject to probe the structure in the medium by various kinds of detection methods in the geotechnical engineering. Based on the propagation theory of elastic wave in half-space layered medium, the ...It is an important subject to probe the structure in the medium by various kinds of detection methods in the geotechnical engineering. Based on the propagation theory of elastic wave in half-space layered medium, the propagation characteristics of elastic wave in layered medium with different elastic parameters are discussed using dynamic analysis of finite element method. It is known that the S-wave velocity, density and thickness of layer are related to the properties of the elastic wave including waveform characteristics, spectral characteristics and time-frequency characteristics. We pay special attention to the structure with low velocity interlayer. The impact imaging method is applied to the grouting construction of the immersed tube tunnel. Data acquisition and analytical method are introduced in detail. The grouting effects can be qualitatively evaluated by comparing the characteristics of elastic wave before grouting with those after grouting. Finally, a quantitative evaluation is obtained according to the relationship between energy response of elastic wave and impedance ratio.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.52405095,12272089,and 92360305)the Guangdong Basic and Applied Basic Research Foundation of China(No.2023A1515110557)+4 种基金the Natural Science Foundation of Liaoning Province of China(No.2023-BSBA-102)the Open Fund of National Key Laboratory of Particle Transport and Separation Technology of China(No.WZKF-2024-6)the Open Project of Guangxi Key Laboratory of Automobile Components and Vehicle Technology of China(Nos.2024GKLACVTKF07 and 2024GKLACVTKF06)the Basic Research Projects of Liaoning Provincial Department of Education of China(No.JYTQN2023162)the Fundamental Research Funds for the Central Universities of China(No.N2403022)。
文摘The testing of large structures is limited by high costs and long cycles, making scaling methods an attractive solution. However, the scaling process of elastic rings introduces complexities in multi-parameter geometric distortions, leading to a diminution in the predictive accuracy of the distorted similitude. To address this challenge, this study formulates a novel set of scaling laws, tailored to account for the intricate geometric distortions associated with elastic rings. The proposed scaling laws are formulated based on the intrinsic deformation characteristics of elastic rings, rather than the traditional systemic governing equations. Numerical and experimental cases are conducted to assess the efficacy and precision of the proposed scaling laws, and the obtained results are compared with those achieved by traditional methods. The outcomes demonstrate that the scaling laws put forth by this study significantly enhance the predictive capabilities for deformations of elastic rings.
基金supported by the National Natural Science Foundation of China(Grant Nos.12272211,12072181,and 12121002).
文摘Uncertain parameters are widespread in engineering systems.This study investigates the modal analysis of a fluid-conveying pipe subjected to elastic supports with unknown-but-bound parameters.The governing equation for the elastically supported fluid-conveying pipe is transformed into ordinary differential equations using the Galerkin truncation method.The Chebyshev interval approach,integrated with the assumed mode method is then used to investigate the effects of uncertainties of support stiffness,fluid speed,and pipe length on the natural frequencies and mode shapes of the pipe.Additionally,both symmetrical and asymmetrical support stiffnesses are discussed.The accuracy and effectiveness of the Chebyshev interval approach are verified through comparison with the Monte Carlo method.The results reveal that,for the same deviation coefficient,uncertainties in symmetrical support stiffness have a greater impact on the first four natural frequencies than those of the asymmetrical one.There may be significant differences in the sensitivity of natural frequencies and mode shapes of the same order to uncertain parameters.Notably,mode shapes susceptible to uncertain parameters exhibit wider fluctuation intervals near the elastic supports,requiring more attention.
基金funded by Vice Chancellor of Research at Shiraz University(grant 3GFU2M1820).
文摘The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.
基金supported by the National Natural Science Foundation of China Excellence Research Group Program for“Multiscale Problems in Nonlinear Mechanics”(Grant No.12588201)the National Key R&D Program of China(Grant No.2023YFA1008901)+1 种基金the National Nat-ural Science Foundation of China(Grant No.12172009)supported by“The Fundamental Research Funds for the Central Universities,Peking University”.
文摘This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solving a specific numerical problem under the scope of the linear finite element method(LFEM),so the method is termed computational method for analytical solutions with finite elements(CMAS-FE).The primary objective of the CMAS-FE is to construct analytical expressions for displacements and reaction forces at nodes,as well as for strains and stresses at elemental quadrature points,all of which are formulated as infinite series solutions of various orders of Poisson’s ratios.Like the conventional LFEM,the CMAS-FE forms global sparse linear equations,but the Young’s modulus and Poisson’s ratio remain variables(or symbols).By employing a direct inverse method to solve these symbolic linear systems,an analytical expression of the displacement field can be constructed.The CMAS-FE is validated via patch and bending tests,which demonstrate convergence with mesh and term refine-ment.Furthermore,the CMAS-FE is applied to obtain the bending stiffness of a beam structure and to estimate an approximate stress intensity factor for a straight crack within a square-shaped plate.
基金sponsorship of the National Natural Science Foundation of China(42430809,42274157,42030103,42404132)the Fund of State Key Laboratory of Deep Oil and Gas,China University of Petroleum(East China)(SKLDOG2024-ZYTS-02)+5 种基金the Postdoctoral Fellowship Program of CPSF(GZB20240850)the Postdoctoral Project of Qingdao(QDBSH20240102082)the Fundamental Research Funds for the Central Universities(24CX07004A,24CX06036A)the CNPC Innovation Fund(2024DQ02-0505,2024DQ02-0136)the Innovation fund project for graduate student of China University of Petroleum(East China)the Fundamental Research Funds for the Central Universities(24CX04002A).
文摘Reservoirs with a group of vertical fractures in a vertical transversely isotropic(VTI)background are considered as orthorhombic(ORT)medium.However,fracture detection in ORT medium using seismic inversion methods remains challenging,as it requires the estimation of more than eight parameters.Assuming the reservoir to be a weakly anisotropic ORT medium with small contrasts in the background elastic parameters,a new azimuthal elastic impedance equation was first derived using parameter combinations and mathematical approximations.This equation exhibited almost the same accuracy as the original equation and contained only six model parameters:the compression modulus,anisotropic shear modulus,anisotropic compression modulus,density,normal fracture weakness,and tangential fracture weakness.Subsequently,a stepwise inversion method using second-order derivatives of the elastic impedance was developed to estimate these parameters.Moreover,the Thomsen anisotropy parameter,epsilon,was estimated from the inversion results using the ratio of the anisotropic compression modulus to the compression modulus.Synthetic examples with moderate noise and field data examples confirm the feasibility and effectiveness of the inversion method.The proposed method exhibited accuracy similar to that of previous inversion strategies and could predict richer vertical fracture information.Ultimately,the method was applied to a three-dimensional work area,and the predictions were consistent with logging and geological a priori information,confirming the effectiveness of this method.Summarily,the proposed stepwise inversion method can alleviate the uncertainty of multi-parameter inversion in ORT medium,thereby improving the reliability of fracture detection.
基金National Science Foundation of China,Grant/Award Number:5227090113Shenzhen Science and Technology Program,Grant/Award Numbers:KCXFZ20230731093901003,KCXFZ20211020163816023The Project of Hetao Shenzhen-Hong Kong Science and Technology Innovation Cooperation Zone,Grant/Award Number:HZQB-KCZYB-2020083。
文摘The two-phase flow in porous media is affected by multiple factors.In the present study,a two-dimensional numerical model of porous media was developed using the actual pore structure of the core sample.The phase field method was utilized to simulate the impact of displacement velocity,the water-gas viscosity ratio,and the density ratio on the flow behavior of two-phase fluids in porous media.The effectiveness of displacement was evaluated by analyzing CO_(2)saturation levels.The results indicate that the saturation of CO_(2)in porous media increased as the displacement velocity increased.When the displacement velocity exceeded 0.01 m/s,there was a corresponding increase in CO_(2)saturation.Conversely,when the displacement velocity was below this threshold,the impact on CO_(2)saturation was minimal.An“inflection point,”M3,was present in the viscosity ratio.When the viscosity of CO_(2)is less than 8.937×10^(-5)Pa·s(viscosity ratio below M3),variations in the viscosity of CO_(2)had little impact on its saturation.Conversely,when the viscosity of CO_(2)exceeded 8.937×10^(-5)Pa·s(viscosity ratio greater than M3),saturation increased with an increase in the viscosity ratio.In terms of the density ratio,the saturation of CO_(2)increased monotonically with an increase in the density ratio.Similarly,increasing density ratios resulted in a monotonic increase in CO_(2)saturation,though this trend was less pronounced in numerical simulations.Analysis results of displacement within dead-end pores using pressure and velocity diagrams reveal eddy currents as contributing factors.Finally,the impact of pore throat structure on the formation of dominant channels was examined.
基金Supported by the National Natural Science Foundation of China(Nos.12272056 and 11832002)。
文摘This paper investigates the active traveling wave vibration control of an elastic supported rotating porous aluminium conical shell(CS)under impact loading.Piezoelectric smart materials in the form of micro fiber composites(MFCs)are used as actuators and sensors.To this end,a metal pore truncated CS with MFCs attached to its surface is considered.Adding artificial virtual springs at two edges of the truncated CS achieves various elastic supported boundaries by changing the spring stiffness.Based on the first-order shear deformation theory(FSDT),minimum energy principle,and artificial virtual spring technology,the theoretical formulations considering the electromechanical coupling are derived.The comparison of the natural frequency of the present results with the natural frequencies reported in previous literature evaluates the accuracy of the present approach.To study the vibration control,the integral quadrature method in conjunction with the differential quadrature approximation in the length direction is used to discretize the partial differential dynamical system to form a set of ordinary differential equations.With the aid of the velocity negative feedback method,both the time history and the input control voltage on the actuator are demonstrated to present the effects of velocity feedback gain,pore distribution type,semi-vertex angle,impact loading,and rotational angular velocity on the traveling wave vibration control.
文摘With the development of seismic engineering and seismic exploration of energy, the underground media that westudy are more and more complicated. Conventional anisotropy theory or two-phase isotropy theory is difficult todescribe anisotropic media containing fluid, such as fractures containing gas, shales containing water Based onBlot theory about two-phase anisotropy, with the use of elastic plane wave equations, we get Christoffel equations.We calculate and analyze the effects of frequency on phase velocity, attenuation, amplitude ratio and polarizationdirection of elastic waves of two-phase, transversely isotropic media. Results show that frequency affects slow Pwave the greatest among the four kinds of waves, i.e., fast P wave, slow P wave, fast S wave and slow S wave.Fluid phase amplitude to solid phase amplitude ratio of fast P wave, fast S wave and slow S wave approaches unitfor large dissipation coefficients. Polarization analysis shows that polarization direction of fluid phase displacement is different from, not parallel to or reverse to, that of solid phase displacement in two-phase anisotropic media.
基金the National Natural Science Foundation of China(No.12172169)the China Scholarship Council(CSC)(No.202006830038)the Natural Sciences and Engineering Research Council of Canada(No.RGPIN-2017-03716115112)。
文摘A nonlocal study of the vibration responses of functionally graded(FG)beams supported by a viscoelastic Winkler-Pasternak foundation is presented.The damping responses of both the Winkler and Pasternak layers of the foundation are considered in the formulation,which were not considered in most literature on this subject,and the bending deformation of the beams and the elastic and damping responses of the foundation as nonlocal by uniting the equivalently differential formulation of well-posed strain-driven(ε-D)and stress-driven(σ-D)two-phase local/nonlocal integral models with constitutive constraints are comprehensively considered,which can address both the stiffness softening and toughing effects due to scale reduction.The generalized differential quadrature method(GDQM)is used to solve the complex eigenvalue problem.After verifying the solution procedure,a series of benchmark results for the vibration frequency of different bounded FG beams supported by the foundation are obtained.Subsequently,the effects of the nonlocality of the foundation on the undamped/damping vibration frequency of the beams are examined.
文摘Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.
文摘The bending of rectangular plate is divided into the generalized statically determinate bending and the generalized statically indeterminate bending based on the analysis of the completeness of calculating condition at the corner point. The former can be solved directly by the equilibrium differential equation and the boundary conditions of four edges of the plate. The latter can be solved by using the superposition principle. Making use of the recommended method, the bending of the plate with all kinds of...
基金supported by the National Science and Technology Major Project of China(Grant No. 2011ZX05004-003,2011ZX05014-006-006)the National Key Basic Research Program of China(Grant No. 2013CB228602)the Natural Science Foundation of China(Grant No. 40974066)
文摘The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.
基金supported by the National Natural Science Foundation of China (No. 50768001)the Foundation of New Century Excellent Talents in University (No. NCET-04-0834)the Guangxi Natural Science Foundation (No. 0728026)
文摘A new strategy for elastic modulus adjustment is proposed based on the element bearing ratio （EBR）,and the elastic modulus reduction method （EMRM） is proposed for limit load evaluation of frame structures. The EBR is defined employing the generalized yield criterion,and the reference EBR is determined by introducing the extrema and the degree of uniformity of EBR in the structure. The elastic modulus in the element with an EBR greater than the reference one is reduced based on the linear elastic finite element analysis and the equilibrium of strain energy. The lower-bound of limit-loads of frame structures are analyzed and the numerical example demonstrates the flexibility,accuracy and effciency of the proposed method.
基金Project(U1234211)supported by the National Natural Science Foundation of ChinaProject(2013G009-B)supported by China Railway Corporation
文摘With the development of high-speed and heavy-haul railway in China, problems like insufficient thickness of ballast bed and overlarge track stiffness are obvious. Ballast may break into small particles and their contact status will deteriorate under cyclic loading, resulting in ballast degradation. Discrete element method(DEM) was used to research improved performance of ballast bed using elastic sleeper. Clusters were generated by bonding spheres to model real ballasts, while broken bonds were utilized to distinguish breakage. Two kinds of ballast beds with elastic sleeper and conventional sleeper were established, respectively. After applying cyclic loading to the models, differences of mechanical properties between two models were analyzed by contrasting their dynamic behavior indexes, such as particle contact force, sleeper settlement, vibration velocity and acceleration, breakage characteristic. The results illustrate that compared with conventional sleeper, elastic sleeper increases sleeper settlement, while reduces ballast vibration and contact force between particles, which could depress ballast breakage.
基金supported by the National Natural Science Foundation of China (Grant No.11026223)the Shanghai Leading Academic Discipline Project,China (Grant No.S30106)the Innovation Fund Project for Graduate Student of Shanghai University,China (Grant No.SHUCX112359)
文摘In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares （CVMLS） approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin （EFG） method and the ICVMLS approximation, the improved complex variable element-free Galerkin （ICVEFG） method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.
基金The project supported by the National Natural Science Foundation of China (19772025)
文摘A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
基金the National Basic Research Program (973) of China(No.2011CB013505)
文摘It is an important subject to probe the structure in the medium by various kinds of detection methods in the geotechnical engineering. Based on the propagation theory of elastic wave in half-space layered medium, the propagation characteristics of elastic wave in layered medium with different elastic parameters are discussed using dynamic analysis of finite element method. It is known that the S-wave velocity, density and thickness of layer are related to the properties of the elastic wave including waveform characteristics, spectral characteristics and time-frequency characteristics. We pay special attention to the structure with low velocity interlayer. The impact imaging method is applied to the grouting construction of the immersed tube tunnel. Data acquisition and analytical method are introduced in detail. The grouting effects can be qualitatively evaluated by comparing the characteristics of elastic wave before grouting with those after grouting. Finally, a quantitative evaluation is obtained according to the relationship between energy response of elastic wave and impedance ratio.
