To find analytical solutions of nonlinear systems for locating the acoustic emission/microseismic（AE/MS） source without knowing the wave velocity of structures, the sensor location coordinates were simplified as a c...To find analytical solutions of nonlinear systems for locating the acoustic emission/microseismic（AE/MS） source without knowing the wave velocity of structures, the sensor location coordinates were simplified as a cuboid monitoring network. Different locations of sensors on upper and lower surfaces were considered and used to establish nonlinear equations. Based on the proposed functions of time difference of arrivals, the analytical solutions were obtained using five sensors under three networks. The proposed analytical solutions were validated using authentic data of numerical tests and experiments. The results show that located results are consistent with authentic data, and the outstanding characteristics of the new solution are that the solved process is not influenced by the wave velocity knowledge and iterated algorithms.展开更多
To find the analytical solution of the acoustic emission/microseismic(AE/MS) source location coordinates, the sensor location coordinates were optimized and simplified. A cube monitoring network of sensor location was...To find the analytical solution of the acoustic emission/microseismic(AE/MS) source location coordinates, the sensor location coordinates were optimized and simplified. A cube monitoring network of sensor location was selected, and the AE/MS source localization equations were established. A location method with P-wave velocity by analytical solutions (P-VAS) was obtained with these equations. The virtual location tests show that the relocation results of analytical method are fully consistent with the actual coordinates for events both inside and outside the monitoring network; whereas the location error of traditional time difference method is between 0.01 and 0.03 m for events inside the sensor array, and the location errors are larger, which is up to 1080986 m for events outside the sensor array. The broken pencil location tests were carried out in the cross section of 100 mm×98 mm, 350 mm-length granite rock specimen using five AE sensors. Five AE sources were relocated with the conventional method and the P-VAS method. For the four events outside monitoring network, the positioning accuracy by P-VAS method is higher than that by the traditional method, and the location accuracy of the larger one can be increased by 17.61 mm. The results of both virtual and broken pencil location tests show that the proposed analytical solution is effective to improve the positioning accuracy. It can locate the coordinates of AE/MS source only using simple four arithmetic operations, without determining the fitting initial value and iterative calculation, which can be solved by a conventional calculator or Microsoft Excel.展开更多
An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion te...An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion technique. Unlike previous studies, in which the saturated soil deposits were simulated with the single-phase elastic theory, in this paper, they are simulated with Biot＇s dynamic theory for saturated porous media, and the half space is assumed as a single-phase elastic medium. The effects of the dimensionless frequency, the incidence angle of P-wave and the porosity of soil deposits on the surface displacement magnifications of the hemispherical alluvial valley are investigated. Numerical results show that the existence of a saturated hemispherical alluvial valley has much influence on the surface displacement magnifications. It is more reasonable to simulate soil deposits with Biot＇s dynamic theory when evaluating the displacement responses of a hemispherical alluvial valley with an incidence of P-waves.展开更多
Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green func...Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.展开更多
It is always a challenge for a model developer to verify a three-dimensional hydrodynamic model, especially for the baroclinic term over variable topography, due to a lack of observational data sets or suitable analyt...It is always a challenge for a model developer to verify a three-dimensional hydrodynamic model, especially for the baroclinic term over variable topography, due to a lack of observational data sets or suitable analytical solutions. In this paper, exact solutions for the periodic forcing by surface heat flux and wind stress are given by solving the linearized equations of motion neglecting the rotation, advection and horizontal diffusion terms. The temperature at the bottom is set to a prescribed periodic value and a slip condition on flow is enforced at the bottom. The geometry of the quarter annulus, which has been extensively studied for two- and three-dimensional analytical solutions of unstratified water bodies, is used with a general power law variation of the bottom slope in the radial direction and is constant in the azimuthal direction. The analytical solutions are derived in a cylindrical coordinate system, which describes the three-dimensional fluid field in a Cartesian coordinate system. The results presented in this paper should provide a foundation for studying and verifying the baroclinic term over a varied topography in a three-dimensional numerical model.展开更多
Progressive modifications in submarine topography and shorelines drastically affect tidal dynamics in bays.This study examines the influence of topographic slope,bay length,and driving forces on tidal currents followi...Progressive modifications in submarine topography and shorelines drastically affect tidal dynamics in bays.This study examines the influence of topographic slope,bay length,and driving forces on tidal currents following land reclamation.Tidal equa-tions are analytically solved using infinite series,deriving expressions for tidal levels and currents in narrow bays with varying topog-raphy.Tidal levels,influenced by topographic variations,are characterized by amplitude and phase lag of their complex amplitude.These levels demonstrate high sensitivity to longitudinal slope variations but remain relatively stable under lateral slope changes.Un-der constant topographic slopes,even minor changes in bay length can drastically modify amplitude and phase lag,highlighting the sensitivity of tidal dynamics to geometric alterations.Tidal velocity notably increases with steeper longitudinal slopes and modestly rises with elevated lateral slopes.However,changes in longitudinal and lateral slopes do not considerably alter flow patterns.While external forces predominantly regulate tidal velocity with negligible effects on flow patterns,endogenous resistance influences veloci-ty but minimally impacts flow structure.These findings enhance the understanding of tidal responses to geometric and topographic changes,providing valuable guidance for land reclamation projects and coastal management strategies.展开更多
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.展开更多
Underhand cut-and-fill mining has been widely used in underground mining operations,especially when the rock mass or orebody is of poor quality or prone to rockburst due to high stress.In such cases,mining workers sho...Underhand cut-and-fill mining has been widely used in underground mining operations,especially when the rock mass or orebody is of poor quality or prone to rockburst due to high stress.In such cases,mining workers should carry out all production activities under the cemented backfill roof or sill mat instead of a highly fractured and unstable rock roof or a strong rock roof with a high potential of rockburst.Therefore,the stability and required strength of the sill mat are critical issues for mining engineers.In 1991,Mitchell considered that sill mat could fail by caving,sliding,rotation,and flexure.Mitchell also proposed an analytical solution to determine the minimum required strength of the sill mat for each type of failure based on two stiff or immobile rock walls.However,recent publications using numerical modeling and field measurements indicate that the compressive stresses in the sill mat induced by rock wall closure due to a stope excavation beneath the sill mat can be significant.It is thus highly necessary to investigate the required strength of the sill mat by considering rock wall closure.In this study,the crushing failure of sill mat due to rock wall closure generated by underground excavation and a new failure mode called"crushing and caving”is revealed by numerical modeling.An analytical solution corresponding to each failure mode is then developed to estimate the minimum required cohesion(cmin)of the sill mat.A criterion is also proposed to determine if the sill mat fails by crushing or crushing-and-caving failure.The proposed analytical solution does not involve any correction coefficients.The validity of the proposed analytical solution is demonstrated by numerical modeling.The proposed analytical solution can thus be employed to predict the cmin of sill mat subjected to wall closure generated by underlying stope excavation.展开更多
The stress gradient of surrounding rock and reasonable prestress of support are the keys to ensuring the stability of roadways.The elastic-plastic analytical solution for surrounding rock was derived based on unified ...The stress gradient of surrounding rock and reasonable prestress of support are the keys to ensuring the stability of roadways.The elastic-plastic analytical solution for surrounding rock was derived based on unified strength theory.A model for solving the stress gradient of the surrounding rock with the intermediate principal stress parameter b was established.The correctness and applicability of the solution for the stress gradient in the roadway surrounding rock was verified via multiple methods.Furthermore,the laws of stress,displacement,and the plastic zone of the surrounding rock with different b values and prestresses were revealed.As b increases,the stress gradient in the plastic zone increases,and the displacement and plastic zone radius decrease.As the prestress increases,the peak stress shifts toward the sidewalls,and the stress and stress gradient increments decrease.In addition,the displacement increment and plastic zone increment were proposed to characterize the support effect.The balance point of the plastic zone area appears before that of the displacement zone.The relationship between the stress gradient compensation coefficient and the prestress is obtained.This study provides a research method and idea for determining the reasonable prestress of support in roadways.展开更多
In the underhand cut-and-fill mining method,a sill mat(i.e.an artificial horizontal pillar)constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs...In the underhand cut-and-fill mining method,a sill mat(i.e.an artificial horizontal pillar)constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs.A critical issue is to determine the minimum required strength of the sill mat to ensure a safe and cost-effective design.Until now,Mitchell’s analytical solution is the only available option,considering two stiff and immobile rock walls.Unavoidable rock wall closure associated with stope excavation below the sill mat was neglected.This,along with other undefined parameters,explains why Mitchell’s solution is rarely used in sill mat design.A new analytical solution for determining the minimum required strength of the sill mat accounting for wall closure is necessary.In this study,a closed-form analytical solution for estimating rock wall closure generated by stope excavation below a sill mat is developed by using Salamon’s and Flamant’s models.The proposed analytical solution does not contain any coefficients of correction or calibration.Despite several assumptions(or somewhat of oversimplifications)necessary to render a simple analytical solution possible,good agreements are obtained between the rock wall closures predicted by applying the proposed analytical solution and those obtained numerically with FLAC3D for many cases with arbitrarily chosen geometrical and material parameters.The proposed analytical solution is therefore validated and can be used to evaluate the rock wall closure generated by stope excavation below a sill mat.展开更多
In practical engineering construction,multi-layered barriers containing geomembranes are extensively applied to retard the migration of pollutants.However,the associated analytical theory on pollutants diffusion still...In practical engineering construction,multi-layered barriers containing geomembranes are extensively applied to retard the migration of pollutants.However,the associated analytical theory on pollutants diffusion still needs to be further improved.In this work,general analytical solutions are derived for one-dimensional diffusion of degradable organic contaminant(DOC)in the multi-layered media containing geomembranes under a time-varying concentration boundary condition,where the variable substitution and separated variable approaches are employed.These analytical solutions with clear expressions can be used not only to study the diffusion behaviors of DOC in bottom and vertical composite barrier systems,but also to verify other complex numerical models.The proposed general analytical solutions are then fully validated via three comparative analyses,including comparisons with the experimental measurements,an existing analytical solution,and a finite-difference solution.Ultimately,the influences of different factors on the composite cutoff wall’s(CCW,which consists of two soil-bentonite layers and a geomembrane)service performance are investigated through a composite vertical barrier system as the application example.The findings obtained from this investigation can provide scientific guidance for the barrier performance evaluation and the engineering design of CCWs.This application example also exhibits the necessity and effectiveness of the developed analytical solutions.展开更多
The double-beam system is a crucial foundational structure in industry,with extensive application contexts and significant research value.The double-beam system with damping and gyroscopic effects is termed as the dam...The double-beam system is a crucial foundational structure in industry,with extensive application contexts and significant research value.The double-beam system with damping and gyroscopic effects is termed as the damped gyroscopic double-beam system.In such systems,the orthogonality conditions of the undamped double-beam system are no longer applicable,rendering it impossible to decouple them in modal space using the modal superposition method(MSM) to obtain analytical solutions.Based on the complex modal method and state space method,this paper takes the damped pipe-in-pipe(PIP) system as an example to solve this problem.The concepts of the original system and adjoint system are introduced,and the orthogonality conditions of the damped PIP system are given in the state-space.Based on the derived orthogonality conditions,the transient and steady-state response solutions are obtained.In the numerical discussion section,the convergence and accuracy of the solutions are verified.In addition,the dynamic responses of the system under different excitations and initial conditions are studied,and the forward and reverse synchronous vibrations in the PIP system are discussed.Overall,the method presented in this paper provides a convenient way to analyze the dynamics of the damped gyroscopic double-beam system.展开更多
In recent years,magneto-electro-elastic(MEE)cylindrical shells with step-wise thicknesses have shown significant potential in the field of vibration energy harvesting.To aid the design of such energy harvesting device...In recent years,magneto-electro-elastic(MEE)cylindrical shells with step-wise thicknesses have shown significant potential in the field of vibration energy harvesting.To aid the design of such energy harvesting devices,an accurate free vibration analysis of embedded MEE cylindrical shells with step-wise thicknesses is performed within the framework of symplectic mechanics.By using the Legendre transformation,a new known vector is defined to transform the higher-order partial differential governing equations into a set of lower-order ordinary differential equations.Therefore,the original vibration analysis is regarded as an eigen problem in the symplectic space,and analytical solutions can be represented by the symplectic series.In numerical examples,the new analytical solutions are compared with the existing results,and good agreement is observed.Furthermore,the effects of critical design parameters on free vibration characteristics are thoroughly investigated.All numerical results can serve as benchmarks for the development of other approximate or numerical methods.展开更多
The radar radiation source signals hold extremely high reconnaissance value.Accurately positioning these signals constitutes one of the key technologies in safeguarding the security of the electromagnetic space.The po...The radar radiation source signals hold extremely high reconnaissance value.Accurately positioning these signals constitutes one of the key technologies in safeguarding the security of the electromagnetic space.The positioning error in multi-station scenarios is influenced not only by the accuracy of positioning parameter estimation but also by the geometric configuration of the positioning platform.This paper focuses on the direction of arrival(DOA),frequency difference of arrival(FDOA),and time difference of arrival(TDOA)methods,analyzing the optimal configuration,optimal detection area,and optimal position dilution of precision in both elevation-known and elevation-unknown scenarios.Specifically,the paper constructs a signal receiving model,establishes the corresponding positioning equations,and performs dimensional normalization on these equations to derive measurement values in meters.Through differential processing,the position dilution of precision is obtained,which is then used as the optimization function to determine the optimal configuration,optimal detection area,and optimal position dilution of precision.Simulation results validate the accuracy of the proposed formulas.展开更多
An analytical solution for scattering of plane P waves by circular-arc layered alluvial valleys was derived by Fourier-Bessel series expansion technique, and the solution was utilized to analyze the effects of alluvia...An analytical solution for scattering of plane P waves by circular-arc layered alluvial valleys was derived by Fourier-Bessel series expansion technique, and the solution was utilized to analyze the effects of alluvial sequence and their relative stiffness on the scattering of incident waves.展开更多
Using the variable transformation method,the formulae of the axial symmetrical wall temperature distribution during steady heat conduction of a hollow cylinder are derived in this paper.The wall temperature distributi...Using the variable transformation method,the formulae of the axial symmetrical wall temperature distribution during steady heat conduction of a hollow cylinder are derived in this paper.The wall temperature distribution and the wall heat flux distribution in both axial and radial direction can be calculated by the temperature distribution of the liquid medium both inside and outside the cylinder with temperature changing in axial direction.The calculation results are almost consistent with the experience results.The applicative condition of the formulae in this paper consists with most of practice.They can be applied to the engineering calculation of the steady heat conduction.The calculation is simple and accurate.展开更多
In the non-uniform stress field, the surrounding rock plastic zone of the circular roadway shows different shapes under the different confining pressure conditions. Based on the boundary shape characteristics of the p...In the non-uniform stress field, the surrounding rock plastic zone of the circular roadway shows different shapes under the different confining pressure conditions. Based on the boundary shape characteristics of the plastic zone, the characteristic radii of the plastic zone were proposed, namely the horizontal,longitudinal and medial axis radii, which could reflect the plastic zone shapes characteristics and classify the sizes of the key parts. On the theoretical basis of elastic-plastic mechanics, analytical solutions for the characteristic radii were obtained by theoretical deduction, and the relationships between the characteristic radii and key influencing factors were analyzed. Finally, the evaluation criterion of the circular roadway surrounding rock plastic zone shapes, evaluation criterion of the location of potential hazards caused by the roadway surrounding rock and evaluation critical points of roadway dynamic disasters based on characteristic radii were proposed. This work could provide a theoretical basis for stability analysis of the surrounding rock, support design, and guide the prevention and control of dynamic roadway disasters.展开更多
Hollow cylinders are widely used in spacecraft, rockets, weapons, metallurgy, materials, and mechanical manufacturing industries, and so on, hydraulic bulging roll cylinder and hydraulic press work all belong to hollo...Hollow cylinders are widely used in spacecraft, rockets, weapons, metallurgy, materials, and mechanical manufacturing industries, and so on, hydraulic bulging roll cylinder and hydraulic press work all belong to hollow cylinders. However, up till now, the solution of the cylinder subjected to the pressures in the three-dimensional space is still at the stage of the analytical solution to the normal pressure or the approximate solution to the variable pressure by numerical method. The analytical solution to the variable pressure of the cylinder has not yet made any breakthrough in theory and can not meet accurate theoretical analysis and calculation requirements of the cylindrical in Engineering. In view of their importance, the precision calculation and theoretical analysis are required to investigate on engineering. A stress function which meets both the biharmonic equations and boundary conditions is constructed in the three-dimensional space. Furthermore, the analytic solution of a hollow cylinder subjected to exponential function distributed variable pressure on its inner and outer surfaces is deduced. By controlling the pressure subject to exponential function distributed variable pressure in the hydraulic bulging roller without any rolling load, using a static tester to record the strain supported hydraulic bulging roll, and comparing with the theoretical calculation, the experimental test result has a higher degree of agreement with the theoretical calculation. Simultaneously, the famous Lam6 solution can be deduced when given the unlimited length of cylinder along the axis. The analytic solution paves the way for the mathematic building and solution of hollow cylinder with randomly uneven pressure.展开更多
A new unified analytical solution is presented for predicting the range of plastic zone and stress distributions around a deep circular tunnel in a homogeneous isotropic continuous medium. The rock mass, grouting zone...A new unified analytical solution is presented for predicting the range of plastic zone and stress distributions around a deep circular tunnel in a homogeneous isotropic continuous medium. The rock mass, grouting zone and lining are assumed as elastic-perfectly plastic and governed by the unified strength theory(UST). This new solution has made it possible to consider the interaction between seepage pressure, lining, grouting and rock mass, and the intermediate principal stress effect together. Moreover, parametric analysis is carried out to identify the influence of the related parameters on the plastic zone radius. Under the given conditions, the results show that the plastic zone radius decreases with an increasing cohesion, internal friction angle and hydraulic conductivity of lining and unified failure criterion parameter, respectively; whereas the plastic zone radius increases with the growth of elasticity modulus of lining. Comparison of results from the new solution and the other published one shows well agreement and provides confidence in the new solution proposed.展开更多
An analytical solution for scattering of plane P waves by a semi-cylindrical hill was derived by using the wave function expansion method, and convergence of the solution and accuracy of truncation were verified. The ...An analytical solution for scattering of plane P waves by a semi-cylindrical hill was derived by using the wave function expansion method, and convergence of the solution and accuracy of truncation were verified. The effect of incident frequency and incident angle on the surface motion of the hill was discussed, and it was shown that a hill greatly amplifies incident plane P waves, and maximum horizontal displacement amplitudes appear mostly at the inclined incidence of waves, which are located at the half-space; and maximum vertical displacement amplitudes emerge mostly at the vertical incidence of waves, which are situated at the hill.展开更多
基金Projects(11447242,41272304,51209236,51274254)supported by the National Natural Science Foundation of ChinaProject(2015CB060200)supported by the National Basic Research Program of China
文摘To find analytical solutions of nonlinear systems for locating the acoustic emission/microseismic（AE/MS） source without knowing the wave velocity of structures, the sensor location coordinates were simplified as a cuboid monitoring network. Different locations of sensors on upper and lower surfaces were considered and used to establish nonlinear equations. Based on the proposed functions of time difference of arrivals, the analytical solutions were obtained using five sensors under three networks. The proposed analytical solutions were validated using authentic data of numerical tests and experiments. The results show that located results are consistent with authentic data, and the outstanding characteristics of the new solution are that the solved process is not influenced by the wave velocity knowledge and iterated algorithms.
基金Project (10872218) supported by the National Natural Science Foundation of ChinaProject (2010CB732004) supported by the National Basic Research Program of China+1 种基金Project (kjdb2010-6) supported by Doctoral Candidate Innovation Research Support Program of Science & Technology ReviewProject (201105) supported by Scholarship Award for Excellent Doctoral Student of Ministry of Education of China
文摘To find the analytical solution of the acoustic emission/microseismic(AE/MS) source location coordinates, the sensor location coordinates were optimized and simplified. A cube monitoring network of sensor location was selected, and the AE/MS source localization equations were established. A location method with P-wave velocity by analytical solutions (P-VAS) was obtained with these equations. The virtual location tests show that the relocation results of analytical method are fully consistent with the actual coordinates for events both inside and outside the monitoring network; whereas the location error of traditional time difference method is between 0.01 and 0.03 m for events inside the sensor array, and the location errors are larger, which is up to 1080986 m for events outside the sensor array. The broken pencil location tests were carried out in the cross section of 100 mm×98 mm, 350 mm-length granite rock specimen using five AE sensors. Five AE sources were relocated with the conventional method and the P-VAS method. For the four events outside monitoring network, the positioning accuracy by P-VAS method is higher than that by the traditional method, and the location accuracy of the larger one can be increased by 17.61 mm. The results of both virtual and broken pencil location tests show that the proposed analytical solution is effective to improve the positioning accuracy. It can locate the coordinates of AE/MS source only using simple four arithmetic operations, without determining the fitting initial value and iterative calculation, which can be solved by a conventional calculator or Microsoft Excel.
基金Project supported by the National Natural Science Foundation of China (No. 50478062) and Natural Science Foundation of Beijing (No. 8052015).
文摘An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion technique. Unlike previous studies, in which the saturated soil deposits were simulated with the single-phase elastic theory, in this paper, they are simulated with Biot＇s dynamic theory for saturated porous media, and the half space is assumed as a single-phase elastic medium. The effects of the dimensionless frequency, the incidence angle of P-wave and the porosity of soil deposits on the surface displacement magnifications of the hemispherical alluvial valley are investigated. Numerical results show that the existence of a saturated hemispherical alluvial valley has much influence on the surface displacement magnifications. It is more reasonable to simulate soil deposits with Biot＇s dynamic theory when evaluating the displacement responses of a hemispherical alluvial valley with an incidence of P-waves.
基金Project supported by the National Natural Science Foundation of China (No. 50776097)
文摘Analytically solving a three-dimensional （3-D） bioheat transfer problem with phase change during a freezing process is extremely difficult but theoretically important. The moving heat source model and the Green function method are introduced to deal with the cryopreservation process of in vitro biomaterials. Exact solutions for the 3-D temperature transients of tissues under various boundary conditions, such as totally convective cooling, totally fixed temperature cooling and a hybrid between them on tissue surfaces, are obtained. Furthermore, the cryosurgical process in living tissues subject to freezing by a single or multiple cryoprobes is also analytically solved. A closed-form analytical solution to the bioheat phase change process is derived by considering contributions from blood perfusion heat transfer, metabolic heat generation, and heat sink of a cryoprobe. The present method is expected to have significant value for analytically solving complex bioheat transfer problems with phase change.
文摘It is always a challenge for a model developer to verify a three-dimensional hydrodynamic model, especially for the baroclinic term over variable topography, due to a lack of observational data sets or suitable analytical solutions. In this paper, exact solutions for the periodic forcing by surface heat flux and wind stress are given by solving the linearized equations of motion neglecting the rotation, advection and horizontal diffusion terms. The temperature at the bottom is set to a prescribed periodic value and a slip condition on flow is enforced at the bottom. The geometry of the quarter annulus, which has been extensively studied for two- and three-dimensional analytical solutions of unstratified water bodies, is used with a general power law variation of the bottom slope in the radial direction and is constant in the azimuthal direction. The analytical solutions are derived in a cylindrical coordinate system, which describes the three-dimensional fluid field in a Cartesian coordinate system. The results presented in this paper should provide a foundation for studying and verifying the baroclinic term over a varied topography in a three-dimensional numerical model.
基金supported by the National Natural Sci-ence Foundation of China(No.U2106204)the Shanxi Water Conservancy Science and Technology Re-search and Extension Project(No.2023ZF19).
文摘Progressive modifications in submarine topography and shorelines drastically affect tidal dynamics in bays.This study examines the influence of topographic slope,bay length,and driving forces on tidal currents following land reclamation.Tidal equa-tions are analytically solved using infinite series,deriving expressions for tidal levels and currents in narrow bays with varying topog-raphy.Tidal levels,influenced by topographic variations,are characterized by amplitude and phase lag of their complex amplitude.These levels demonstrate high sensitivity to longitudinal slope variations but remain relatively stable under lateral slope changes.Un-der constant topographic slopes,even minor changes in bay length can drastically modify amplitude and phase lag,highlighting the sensitivity of tidal dynamics to geometric alterations.Tidal velocity notably increases with steeper longitudinal slopes and modestly rises with elevated lateral slopes.However,changes in longitudinal and lateral slopes do not considerably alter flow patterns.While external forces predominantly regulate tidal velocity with negligible effects on flow patterns,endogenous resistance influences veloci-ty but minimally impacts flow structure.These findings enhance the understanding of tidal responses to geometric and topographic changes,providing valuable guidance for land reclamation projects and coastal management strategies.
基金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.
基金financial support from the Young Scientist Project of the National Key Research and Development Program of China(Grant No.2021YFC2900600)Beijing Nova Program(Grant No.20220484057)+1 种基金The authors acknowledge the financial support from the Natural Sciences and Engineering Research Council of Canada(Grant No.RGPIN-2018-06902)industrial partners of the Research Institute on Mines and the Environment(RIME UQAT-Polytechnique:https://irme.ca/en/).
文摘Underhand cut-and-fill mining has been widely used in underground mining operations,especially when the rock mass or orebody is of poor quality or prone to rockburst due to high stress.In such cases,mining workers should carry out all production activities under the cemented backfill roof or sill mat instead of a highly fractured and unstable rock roof or a strong rock roof with a high potential of rockburst.Therefore,the stability and required strength of the sill mat are critical issues for mining engineers.In 1991,Mitchell considered that sill mat could fail by caving,sliding,rotation,and flexure.Mitchell also proposed an analytical solution to determine the minimum required strength of the sill mat for each type of failure based on two stiff or immobile rock walls.However,recent publications using numerical modeling and field measurements indicate that the compressive stresses in the sill mat induced by rock wall closure due to a stope excavation beneath the sill mat can be significant.It is thus highly necessary to investigate the required strength of the sill mat by considering rock wall closure.In this study,the crushing failure of sill mat due to rock wall closure generated by underground excavation and a new failure mode called"crushing and caving”is revealed by numerical modeling.An analytical solution corresponding to each failure mode is then developed to estimate the minimum required cohesion(cmin)of the sill mat.A criterion is also proposed to determine if the sill mat fails by crushing or crushing-and-caving failure.The proposed analytical solution does not involve any correction coefficients.The validity of the proposed analytical solution is demonstrated by numerical modeling.The proposed analytical solution can thus be employed to predict the cmin of sill mat subjected to wall closure generated by underlying stope excavation.
基金Project(52274130)supported by the National Natural Science Foundation of ChinaProject(ZR2024ZD22)supported by the Major Basic Research Project of the Shandong Provincial Natural Science Foundation,China+2 种基金Project(2023375)supported by the Guizhou University Research and Innovation Team,ChinaProject(Leading Fund(2023)09)supported by the Natural Science Research Fund of Guizhou University,ChinaProject(JYBSYS2021101)supported by the Open Fund of Key Laboratory of Safe and Effective Coal Mining,Ministry of Education,China。
文摘The stress gradient of surrounding rock and reasonable prestress of support are the keys to ensuring the stability of roadways.The elastic-plastic analytical solution for surrounding rock was derived based on unified strength theory.A model for solving the stress gradient of the surrounding rock with the intermediate principal stress parameter b was established.The correctness and applicability of the solution for the stress gradient in the roadway surrounding rock was verified via multiple methods.Furthermore,the laws of stress,displacement,and the plastic zone of the surrounding rock with different b values and prestresses were revealed.As b increases,the stress gradient in the plastic zone increases,and the displacement and plastic zone radius decrease.As the prestress increases,the peak stress shifts toward the sidewalls,and the stress and stress gradient increments decrease.In addition,the displacement increment and plastic zone increment were proposed to characterize the support effect.The balance point of the plastic zone area appears before that of the displacement zone.The relationship between the stress gradient compensation coefficient and the prestress is obtained.This study provides a research method and idea for determining the reasonable prestress of support in roadways.
基金financial support from the Young Scientist Project of the National Key Research and Development Program of China(Grant No.2021YFC2900600)the Beijing Nova Program(Grant No.20220484057)+1 种基金The authors acknowledge the financial support from the Natural Sciences and Engineering Research Council of Canada(Grant No.RGPIN-2018-06902)industrial partners of the Research Institute on Mines and the Environment(RIME UQAT-Polytechnique:https://irme.ca/en/).
文摘In the underhand cut-and-fill mining method,a sill mat(i.e.an artificial horizontal pillar)constructed by cemented backfill is essential to prevent mine workers from being directly exposed under problematic rock roofs.A critical issue is to determine the minimum required strength of the sill mat to ensure a safe and cost-effective design.Until now,Mitchell’s analytical solution is the only available option,considering two stiff and immobile rock walls.Unavoidable rock wall closure associated with stope excavation below the sill mat was neglected.This,along with other undefined parameters,explains why Mitchell’s solution is rarely used in sill mat design.A new analytical solution for determining the minimum required strength of the sill mat accounting for wall closure is necessary.In this study,a closed-form analytical solution for estimating rock wall closure generated by stope excavation below a sill mat is developed by using Salamon’s and Flamant’s models.The proposed analytical solution does not contain any coefficients of correction or calibration.Despite several assumptions(or somewhat of oversimplifications)necessary to render a simple analytical solution possible,good agreements are obtained between the rock wall closures predicted by applying the proposed analytical solution and those obtained numerically with FLAC3D for many cases with arbitrarily chosen geometrical and material parameters.The proposed analytical solution is therefore validated and can be used to evaluate the rock wall closure generated by stope excavation below a sill mat.
基金Project(2023YFC3707800)supported by the National Key Research and Development Program of China。
文摘In practical engineering construction,multi-layered barriers containing geomembranes are extensively applied to retard the migration of pollutants.However,the associated analytical theory on pollutants diffusion still needs to be further improved.In this work,general analytical solutions are derived for one-dimensional diffusion of degradable organic contaminant(DOC)in the multi-layered media containing geomembranes under a time-varying concentration boundary condition,where the variable substitution and separated variable approaches are employed.These analytical solutions with clear expressions can be used not only to study the diffusion behaviors of DOC in bottom and vertical composite barrier systems,but also to verify other complex numerical models.The proposed general analytical solutions are then fully validated via three comparative analyses,including comparisons with the experimental measurements,an existing analytical solution,and a finite-difference solution.Ultimately,the influences of different factors on the composite cutoff wall’s(CCW,which consists of two soil-bentonite layers and a geomembrane)service performance are investigated through a composite vertical barrier system as the application example.The findings obtained from this investigation can provide scientific guidance for the barrier performance evaluation and the engineering design of CCWs.This application example also exhibits the necessity and effectiveness of the developed analytical solutions.
基金Project supported by the National Natural Science Foundation of China (No. 12272323)。
文摘The double-beam system is a crucial foundational structure in industry,with extensive application contexts and significant research value.The double-beam system with damping and gyroscopic effects is termed as the damped gyroscopic double-beam system.In such systems,the orthogonality conditions of the undamped double-beam system are no longer applicable,rendering it impossible to decouple them in modal space using the modal superposition method(MSM) to obtain analytical solutions.Based on the complex modal method and state space method,this paper takes the damped pipe-in-pipe(PIP) system as an example to solve this problem.The concepts of the original system and adjoint system are introduced,and the orthogonality conditions of the damped PIP system are given in the state-space.Based on the derived orthogonality conditions,the transient and steady-state response solutions are obtained.In the numerical discussion section,the convergence and accuracy of the solutions are verified.In addition,the dynamic responses of the system under different excitations and initial conditions are studied,and the forward and reverse synchronous vibrations in the PIP system are discussed.Overall,the method presented in this paper provides a convenient way to analyze the dynamics of the damped gyroscopic double-beam system.
基金Project supported by the Science and Technology Plan Joint Program of Liaoning Province of China(Natural Science Foundation-Doctoral Research Launch Project)(No.2024-BSLH-027)the Fundamental Research Funds for Undergraduate Universities of Liaoning Province of China(No.LJBKY2024033)+1 种基金the National Natural Science Foundation of China(No.12472064)the Natural Science Foundation of Liaoning Province of China(No.2023-MS-118)。
文摘In recent years,magneto-electro-elastic(MEE)cylindrical shells with step-wise thicknesses have shown significant potential in the field of vibration energy harvesting.To aid the design of such energy harvesting devices,an accurate free vibration analysis of embedded MEE cylindrical shells with step-wise thicknesses is performed within the framework of symplectic mechanics.By using the Legendre transformation,a new known vector is defined to transform the higher-order partial differential governing equations into a set of lower-order ordinary differential equations.Therefore,the original vibration analysis is regarded as an eigen problem in the symplectic space,and analytical solutions can be represented by the symplectic series.In numerical examples,the new analytical solutions are compared with the existing results,and good agreement is observed.Furthermore,the effects of critical design parameters on free vibration characteristics are thoroughly investigated.All numerical results can serve as benchmarks for the development of other approximate or numerical methods.
基金supported by the National Natural Science Foundation of China(Nos.62027801,62301035).
文摘The radar radiation source signals hold extremely high reconnaissance value.Accurately positioning these signals constitutes one of the key technologies in safeguarding the security of the electromagnetic space.The positioning error in multi-station scenarios is influenced not only by the accuracy of positioning parameter estimation but also by the geometric configuration of the positioning platform.This paper focuses on the direction of arrival(DOA),frequency difference of arrival(FDOA),and time difference of arrival(TDOA)methods,analyzing the optimal configuration,optimal detection area,and optimal position dilution of precision in both elevation-known and elevation-unknown scenarios.Specifically,the paper constructs a signal receiving model,establishes the corresponding positioning equations,and performs dimensional normalization on these equations to derive measurement values in meters.Through differential processing,the position dilution of precision is obtained,which is then used as the optimization function to determine the optimal configuration,optimal detection area,and optimal position dilution of precision.Simulation results validate the accuracy of the proposed formulas.
基金State Natural Science Foundation of China (No.59878032).
文摘An analytical solution for scattering of plane P waves by circular-arc layered alluvial valleys was derived by Fourier-Bessel series expansion technique, and the solution was utilized to analyze the effects of alluvial sequence and their relative stiffness on the scattering of incident waves.
文摘Using the variable transformation method,the formulae of the axial symmetrical wall temperature distribution during steady heat conduction of a hollow cylinder are derived in this paper.The wall temperature distribution and the wall heat flux distribution in both axial and radial direction can be calculated by the temperature distribution of the liquid medium both inside and outside the cylinder with temperature changing in axial direction.The calculation results are almost consistent with the experience results.The applicative condition of the formulae in this paper consists with most of practice.They can be applied to the engineering calculation of the steady heat conduction.The calculation is simple and accurate.
基金supported by the National Natural Science Foundation of China (Grant No. 51234006)the National Key Research and Development Program of China (Grant No. 2016YFC0600708)
文摘In the non-uniform stress field, the surrounding rock plastic zone of the circular roadway shows different shapes under the different confining pressure conditions. Based on the boundary shape characteristics of the plastic zone, the characteristic radii of the plastic zone were proposed, namely the horizontal,longitudinal and medial axis radii, which could reflect the plastic zone shapes characteristics and classify the sizes of the key parts. On the theoretical basis of elastic-plastic mechanics, analytical solutions for the characteristic radii were obtained by theoretical deduction, and the relationships between the characteristic radii and key influencing factors were analyzed. Finally, the evaluation criterion of the circular roadway surrounding rock plastic zone shapes, evaluation criterion of the location of potential hazards caused by the roadway surrounding rock and evaluation critical points of roadway dynamic disasters based on characteristic radii were proposed. This work could provide a theoretical basis for stability analysis of the surrounding rock, support design, and guide the prevention and control of dynamic roadway disasters.
基金supported by National Natural Science Foundation of China (Grant No. 50875230)
文摘Hollow cylinders are widely used in spacecraft, rockets, weapons, metallurgy, materials, and mechanical manufacturing industries, and so on, hydraulic bulging roll cylinder and hydraulic press work all belong to hollow cylinders. However, up till now, the solution of the cylinder subjected to the pressures in the three-dimensional space is still at the stage of the analytical solution to the normal pressure or the approximate solution to the variable pressure by numerical method. The analytical solution to the variable pressure of the cylinder has not yet made any breakthrough in theory and can not meet accurate theoretical analysis and calculation requirements of the cylindrical in Engineering. In view of their importance, the precision calculation and theoretical analysis are required to investigate on engineering. A stress function which meets both the biharmonic equations and boundary conditions is constructed in the three-dimensional space. Furthermore, the analytic solution of a hollow cylinder subjected to exponential function distributed variable pressure on its inner and outer surfaces is deduced. By controlling the pressure subject to exponential function distributed variable pressure in the hydraulic bulging roller without any rolling load, using a static tester to record the strain supported hydraulic bulging roll, and comparing with the theoretical calculation, the experimental test result has a higher degree of agreement with the theoretical calculation. Simultaneously, the famous Lam6 solution can be deduced when given the unlimited length of cylinder along the axis. The analytic solution paves the way for the mathematic building and solution of hollow cylinder with randomly uneven pressure.
基金Project(51378309)supported by National Natural Science Foundation of China
文摘A new unified analytical solution is presented for predicting the range of plastic zone and stress distributions around a deep circular tunnel in a homogeneous isotropic continuous medium. The rock mass, grouting zone and lining are assumed as elastic-perfectly plastic and governed by the unified strength theory(UST). This new solution has made it possible to consider the interaction between seepage pressure, lining, grouting and rock mass, and the intermediate principal stress effect together. Moreover, parametric analysis is carried out to identify the influence of the related parameters on the plastic zone radius. Under the given conditions, the results show that the plastic zone radius decreases with an increasing cohesion, internal friction angle and hydraulic conductivity of lining and unified failure criterion parameter, respectively; whereas the plastic zone radius increases with the growth of elasticity modulus of lining. Comparison of results from the new solution and the other published one shows well agreement and provides confidence in the new solution proposed.
基金National Natural Science Foundation of China Under Grant No.50378063Excellent Young Teacher Program of Ministry of Education
文摘An analytical solution for scattering of plane P waves by a semi-cylindrical hill was derived by using the wave function expansion method, and convergence of the solution and accuracy of truncation were verified. The effect of incident frequency and incident angle on the surface motion of the hill was discussed, and it was shown that a hill greatly amplifies incident plane P waves, and maximum horizontal displacement amplitudes appear mostly at the inclined incidence of waves, which are located at the half-space; and maximum vertical displacement amplitudes emerge mostly at the vertical incidence of waves, which are situated at the hill.
