A Hamiltonian system is derived for the plane elasticity problem of two-dimensional dodecagonal quasicrystals by introducing the simple state function. By using symplectic elasticity approach, the analytic solutions o...A Hamiltonian system is derived for the plane elasticity problem of two-dimensional dodecagonal quasicrystals by introducing the simple state function. By using symplectic elasticity approach, the analytic solutions of the phonon and phason displacements are obtained further for the quasicrystal plates. In addition, the effectiveness of the approach is verified by comparison with the data of the finite integral transformation method.展开更多
Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with th...Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with the isotropic background,while the explicit model for the cracked rock with the anisotropic background is rarely investigated in spite of such case being often encountered in the earth.Hence,we first studied dependences of the crack opening displacement tensors on the crack dip angle in the coordinate systems formed by symmetry planes of the crack and the background anisotropy,respectively,by forty groups of numerical experiments.Based on the conclusion from the experiments,the analytical solution was derived for the effective elastic properties of the rock with the inclined penny-shaped cracks in the transversely isotropic background.Further,we comprehensively analyzed,according to the developed model,effects of the crack dip angle,background anisotropy,filling fluid and crack density on the effective elastic properties of the cracked rock.The analysis results indicate that the dip angle and background anisotropy can significantly either enhance or weaken the anisotropy degrees of the P-and SH-wave velocities,whereas they have relatively small effects on the SV-wave velocity anisotropy.Moreover,the filling fluid can increase the stiffness coefficients related to the compressional modulus by reducing crack compliance parameters,while its effects on shear coefficients depend on the crack dip angle.The increasing crack density reduces velocities of the dry rock,and decreasing rates of the velocities are affected by the crack dip angle.By comparing with exact numerical results and experimental data,it was demonstrated that the proposed model can achieve high-precision estimations of stiffness coefficients.Moreover,the assumption of the weakly anisotropic background results in the consistency between the proposed model and Hudson's published theory for the orthorhombic rock.展开更多
Based on the Landau-Lifshitz-Gilbert(LLG)equation,the precession relaxation of magnetization is studied when the external field H is parallel to the uniaxial anisotropic field H_(k).The evolution of three-component ma...Based on the Landau-Lifshitz-Gilbert(LLG)equation,the precession relaxation of magnetization is studied when the external field H is parallel to the uniaxial anisotropic field H_(k).The evolution of three-component magnetization is solved analytically under the condition of H=nH_(k)(n=3,1 and 0).It is found that with an increase of H or a decrease of the initial polar angle of magnetization,the relaxation time decreases and the angular frequency of magnetization increases.For comparison,the analytical solution for H_(k)=0 is also given.When the magnetization becomes stable,the angular frequency is proportional to the total effective field acting on the magnetization.The analytical solutions are not only conducive to the understanding of the precession relaxation of magnetization,but also can be used as a standard model to test the numerical calculation of LLG equation.展开更多
Classical localization methods use Cartesian or Polar coordinates, which require a priori range information to determine whether to estimate position or to only find bearings. The modified polar representation (MPR) u...Classical localization methods use Cartesian or Polar coordinates, which require a priori range information to determine whether to estimate position or to only find bearings. The modified polar representation (MPR) unifies near-field and farfield models, alleviating the thresholding effect. Current localization methods in MPR based on the angle of arrival (AOA) and time difference of arrival (TDOA) measurements resort to semidefinite relaxation (SDR) and Gauss-Newton iteration, which are computationally complex and face the possible diverge problem. This paper formulates a pseudo linear equation between the measurements and the unknown MPR position,which leads to a closed-form solution for the hybrid TDOA-AOA localization problem, namely hybrid constrained optimization(HCO). HCO attains Cramér-Rao bound (CRB)-level accuracy for mild Gaussian noise. Compared with the existing closed-form solutions for the hybrid TDOA-AOA case, HCO provides comparable performance to the hybrid generalized trust region subproblem (HGTRS) solution and is better than the hybrid successive unconstrained minimization (HSUM) solution in large noise region. Its computational complexity is lower than that of HGTRS. Simulations validate the performance of HCO achieves the CRB that the maximum likelihood estimator (MLE) attains if the noise is small, but the MLE deviates from CRB earlier.展开更多
We search for analytical wave solutions of an electronically and biologically important model named as the Fitzhugh–Nagumo model with truncated M-fractional derivative, in which the expafunction and extended sinh-Gor...We search for analytical wave solutions of an electronically and biologically important model named as the Fitzhugh–Nagumo model with truncated M-fractional derivative, in which the expafunction and extended sinh-Gordon equation expansion(ESh GEE) schemes are utilized. The solutions obtained include dark, bright, dark-bright, periodic and other kinds of solitons. These analytical wave solutions are gained and verified with the use of Mathematica software. These solutions do not exist in literature. Some of the solutions are demonstrated by 2D, 3D and contour graphs. This model is mostly used in circuit theory, transmission of nerve impulses, and population genetics. Finally, both the schemes are more applicable, reliable and significant to deal with the fractional nonlinear partial differential equations.展开更多
Multiple objectives to be optimized simultaneously are prevalent in real-life problems. This paper develops a new Pareto Method for bi-objective optimization which yields analytical solutions. The Pareto optimal front...Multiple objectives to be optimized simultaneously are prevalent in real-life problems. This paper develops a new Pareto Method for bi-objective optimization which yields analytical solutions. The Pareto optimal front is obtained in closed-form, enabling the derivation of various solutions in a convenient and efficient way. The advantage of analytical solution is the possibility of deriving accurate, exact and well-understood solutions, which is especially useful for policy analysis. An extension of the method to include multiple objectives is provided with the objectives being classified into two types. Such an extension expands the applicability of the developed techniques.展开更多
This paper investigates the validity and shortcomings of the existing analytical solution for the ultimate bearing capacity of a pile embedded in a rock mass using the modified HoekeBrown failure criterion.Although th...This paper investigates the validity and shortcomings of the existing analytical solution for the ultimate bearing capacity of a pile embedded in a rock mass using the modified HoekeBrown failure criterion.Although this criterion is considered a reference value for empirical and numerical calculations,some limitations of its basic simplifications have not been clarified yet.This research compares the analytical results obtained from the novel discontinuity layout optimization(DLO)method and the numerical solutions from the finite difference method(FDM).The limitations of the analytical solution are considered by comparing different DLO failure modes,thus allowing for the first time a critical evaluation of its scope and conditioning for implementation.Errors of up to 40%in the bearing capacity and unrealistic failure modes are the main issues in the analytical solution.The main aspects of the DLO method are also analyzed with an emphasis on the linearization of the rock failure criterion and the accuracy resulting from the discretization size.The analysis demonstrates DLO as a very efficient and accurate tool to address the pile tip bearing capacity,presenting considerable advantages over other methods.展开更多
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.展开更多
In regard to unconventional oil reservoirs,the transient dual-porosity and triple-porosity models have been adopted to describe the fluid flow in the complex fracture network.It has been proven to cause inaccurate pro...In regard to unconventional oil reservoirs,the transient dual-porosity and triple-porosity models have been adopted to describe the fluid flow in the complex fracture network.It has been proven to cause inaccurate production evaluations because of the absence of matrix-macrofracture communication.In addition,most of the existing models are solved analytically based on Laplace transform and numerical inversion.Hence,an approximate analytical solution is derived directly in real-time space considering variable matrix blocks and simultaneous matrix depletion.To simplify the derivation,the simultaneous matrix depletion is divided into two parts:one part feeding the macrofractures and the other part feeding the microfractures.Then,a series of partial differential equations(PDEs)describing the transient flow and boundary conditions are constructed and solved analytically by integration.Finally,a relationship between oil rate and production time in real-time space is obtained.The new model is verified against classical analytical models.When the microfracture system and matrix-macrofracture communication is neglected,the result of the new model agrees with those obtained with the dual-porosity and triple-porosity model,respectively.Certainly,the new model also has an excellent agreement with the numerical model.The model is then applied to two actual tight oil wells completed in western Canada sedimentary basin.After identifying the flow regime,the solution suitably matches the field production data,and the model parameters are determined.Through these output parameters,we can accurately forecast the production and even estimate the petrophysical properties.展开更多
We report on the magnetohydrodynamic impact on the axisymmetric flow of Al_(2)O_(3)/Cu nanoparticles suspended in H_(2)O past a stretched/shrinked sheet.With the use of partial differential equations and the correspon...We report on the magnetohydrodynamic impact on the axisymmetric flow of Al_(2)O_(3)/Cu nanoparticles suspended in H_(2)O past a stretched/shrinked sheet.With the use of partial differential equations and the corresponding thermophysical characteristics of nanoparticles,the physical flow process is illustrated.The resultant nonlinear system of partial differential equations is converted into a system of ordinary differential equations using the suitable similarity transformations.The transformed differential equations are solved analytically.Impacts of the magnetic parameter,solid volume fraction and stretching/shrinking parameter on momentum and temperature distribution have been analyzed and interpreted graphically.The skin friction and Nusselt number were also evaluated.In addition,existence of dual solution was deduced for the shrinking sheet and unique solution for the stretching one.Further,Al_(2)O_(3)/H_(2)O nanofluid flow has better thermal conductivity on comparing with Cu/H_(2)O nanofluid.Furthermore,it was found that the first solutions of the stream are stable and physically realizable,whereas those of the second ones are unstable.展开更多
Second-order axially moving systems are common models in the field of dynamics, such as axially moving strings, cables, and belts. In the traditional research work, it is difficult to obtain closed-form solutions for ...Second-order axially moving systems are common models in the field of dynamics, such as axially moving strings, cables, and belts. In the traditional research work, it is difficult to obtain closed-form solutions for the forced vibration when the damping effect and the coupling effect of multiple second-order models are considered.In this paper, Green's function method based on the Laplace transform is used to obtain closed-form solutions for the forced vibration of second-order axially moving systems. By taking the axially moving damping string system and multi-string system connected by springs as examples, the detailed solution methods and the analytical Green's functions of these second-order systems are given. The mode functions and frequency equations are also obtained by the obtained Green's functions. The reliability and convenience of the results are verified by several examples. This paper provides a systematic analytical method for the dynamic analysis of second-order axially moving systems, and the obtained Green's functions are applicable to different second-order systems rather than just string systems. In addition, the work of this paper also has positive significance for the study on the forced vibration of high-order systems.展开更多
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.展开更多
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.展开更多
This paper analyses the bending of rectangular orthotropic plates on a Winkler elastic foundation.Appropriate definition of symplectic inner product and symplectic space formed by generalized displacements establish d...This paper analyses the bending of rectangular orthotropic plates on a Winkler elastic foundation.Appropriate definition of symplectic inner product and symplectic space formed by generalized displacements establish dual variables and dual equations in the symplectic space.The operator matrix of the equation set is proven to be a Hamilton operator matrix.Separation of variables and eigenfunction expansion creates a basis for analyzing the bending of rectangular orthotropic plates on Winkler elastic foundation and obtaining solutions for plates having any boundary condition.There is discussion of symplectic eigenvalue problems of orthotropic plates under two typical boundary conditions,with opposite sides simply supported and opposite sides clamped.Transcendental equations of eigenvalues and symplectic eigenvectors in analytical form given.Analytical solutions using two examples are presented to show the use of the new methods described in this paper.To verify the accuracy and convergence,a fully simply supported plate that is fully and simply supported under uniformly distributed load is used to compare the classical Navier method,the Levy method and the new method.Results show that the new technique has good accuracy and better convergence speed than other methods,especially in relation to internal forces.A fully clamped rectangular plate on Winkler foundation is solved to validate application of the new methods,with solutions compared to those produced by the Galerkin method.展开更多
The theory of an eddy viscosity model is applied to the study of the flow in a compound channel which is partially vegetated. The governing equation is constituted by analyzing the longitudinal forces acting on the un...The theory of an eddy viscosity model is applied to the study of the flow in a compound channel which is partially vegetated. The governing equation is constituted by analyzing the longitudinal forces acting on the unit volume where the effect of the vegetation on the flow is considered as a drag force item, The compound channel is divided into 3 sub-regions in the transverse direction, and the coefficients in every region's differential equations were solved simultaneously. Thus, the analytical solution of the transverse distribution of the depth-averaged velocity for uniform flow in a partially vegetated compound channel was obtained. The results can be used to predict the transverse distribution of bed shear stress, which has an important effect on the transportation of sediment. By comparing the analytical results with the measured data, the analytical solution in this paper is shown to be sufficiently accurate to predict most hydraulic features for engineering design purposes.展开更多
This paper presents an analytical solution of the one-dimensional consolidation in unsaturated soil with a finite thickness under vertical loading and confinements in the lateral directions. The boundary contains the ...This paper presents an analytical solution of the one-dimensional consolidation in unsaturated soil with a finite thickness under vertical loading and confinements in the lateral directions. The boundary contains the top surface permeable to water and air and the bottom impermeable to water and air. The analytical solution is for Fredlund's one-dimensional consolidation equation in unsaturated soils. The transfer relationship between the state vectors at top surface and any depth is obtained by using the Laplace transform and Cayley-Hamilton mathematical methods to the governing equations of water and air, Darcy's law and Fick's law. Excess pore-air pressure, excess pore-water pressure and settlement in the Laplace-transformed domain are obtained by using the Laplace transform with the initial conditions and boundary conditions. By performing inverse Laplace transforms, the analytical solutions are obtained in the time domain. A typical example illustrates the consolidation characteristics of unsaturated soil from an- alytical results. Finally, comparisons between the analytical solutions and results of the finite difference method indicate that the analytical solution is correct.展开更多
The stress concentration problem in structures with a circular or elliptic hole can be investigated by analytical methods. For the problem with a rectangular hole, only approximate results are derived. This paper dedu...The stress concentration problem in structures with a circular or elliptic hole can be investigated by analytical methods. For the problem with a rectangular hole, only approximate results are derived. This paper deduces the analytical solutions to the stress concentration problem in plates with a rectangular hole under biaxial tensions. By using the U-transformation technique and the finite element method, the analytical displacement solutions of the finite element equations are derived in the series form. Therefore, the stress concentration can then be discussed easily and conveniently. For plate problem the bilinear rectangular element with four nodes is taken as an example to demonstrate the applicability of the proposed method. The stress concentration factors for various ratios of height to width of the hole are obtained.展开更多
Based on Huang's accurate tri-sectional nonlin- ear kinematic equation (1997), a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeabili...Based on Huang's accurate tri-sectional nonlin- ear kinematic equation (1997), a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction: The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation (PDE) sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient (TPG) has a great effect on the external moving boundary but has little effect on the internal moving boundary.展开更多
According to the physiological and anatomical characteristics of small intestine,neglecting the effect of its motility on the distribution and absorption of drug and nutrient,Y.Miyamoto et al.proposed a model of two-d...According to the physiological and anatomical characteristics of small intestine,neglecting the effect of its motility on the distribution and absorption of drug and nutrient,Y.Miyamoto et al.proposed a model of two-dimensional laminar flow in a circular porous tube with permeable wall and calculated the concentration profile of drugby numerical analysis.In this paper,we give a steady-state analytical solution of the above model including deactivationterm.The obtained results are in agreement with the results of their numerical analysis. Moreover the analytical solution presented in this paper reveals the relation among the physiological parameters of the model and describes the basic absorption rule of drug and nutrient through the intestinal wall and hence pro- vides a theoretical basis for determining the permeability and reflection coefficient through in situ experiments.展开更多
基金Project supported by the National Natural Science Foundation of China (Grant Nos.12261064 and 11861048)the Natural Science Foundation of Inner Mongolia,China (Grant Nos.2021MS01004 and 2022QN01008)the High-level Talents Scientific Research Start-up Foundation of Inner Mongolia University (Grant No.10000-21311201/165)。
文摘A Hamiltonian system is derived for the plane elasticity problem of two-dimensional dodecagonal quasicrystals by introducing the simple state function. By using symplectic elasticity approach, the analytic solutions of the phonon and phason displacements are obtained further for the quasicrystal plates. In addition, the effectiveness of the approach is verified by comparison with the data of the finite integral transformation method.
基金We would like to acknowledge all the reviewers and editors and the sponsorship of National Natural Science Foundation of China(42030103)the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(2021QNLM020001-6)the Laoshan National Laboratory of Science and Technology Foundation(LSKJ202203400).
文摘Seismic prediction of cracks is of great significance in many disciplines,for which the rock physics model is indispensable.However,up to now,multitudinous analytical models focus primarily on the cracked rock with the isotropic background,while the explicit model for the cracked rock with the anisotropic background is rarely investigated in spite of such case being often encountered in the earth.Hence,we first studied dependences of the crack opening displacement tensors on the crack dip angle in the coordinate systems formed by symmetry planes of the crack and the background anisotropy,respectively,by forty groups of numerical experiments.Based on the conclusion from the experiments,the analytical solution was derived for the effective elastic properties of the rock with the inclined penny-shaped cracks in the transversely isotropic background.Further,we comprehensively analyzed,according to the developed model,effects of the crack dip angle,background anisotropy,filling fluid and crack density on the effective elastic properties of the cracked rock.The analysis results indicate that the dip angle and background anisotropy can significantly either enhance or weaken the anisotropy degrees of the P-and SH-wave velocities,whereas they have relatively small effects on the SV-wave velocity anisotropy.Moreover,the filling fluid can increase the stiffness coefficients related to the compressional modulus by reducing crack compliance parameters,while its effects on shear coefficients depend on the crack dip angle.The increasing crack density reduces velocities of the dry rock,and decreasing rates of the velocities are affected by the crack dip angle.By comparing with exact numerical results and experimental data,it was demonstrated that the proposed model can achieve high-precision estimations of stiffness coefficients.Moreover,the assumption of the weakly anisotropic background results in the consistency between the proposed model and Hudson's published theory for the orthorhombic rock.
基金Project supported by the National Key R&D Program of China (Grant No.2021YFB3501300)the National Natural Science Foundation of China (Grant Nos.91963201 and 12174163)the 111 Project (Grant No.B20063)。
文摘Based on the Landau-Lifshitz-Gilbert(LLG)equation,the precession relaxation of magnetization is studied when the external field H is parallel to the uniaxial anisotropic field H_(k).The evolution of three-component magnetization is solved analytically under the condition of H=nH_(k)(n=3,1 and 0).It is found that with an increase of H or a decrease of the initial polar angle of magnetization,the relaxation time decreases and the angular frequency of magnetization increases.For comparison,the analytical solution for H_(k)=0 is also given.When the magnetization becomes stable,the angular frequency is proportional to the total effective field acting on the magnetization.The analytical solutions are not only conducive to the understanding of the precession relaxation of magnetization,but also can be used as a standard model to test the numerical calculation of LLG equation.
基金supported by the National Natural Science Foundation of China (62101359)Sichuan University and Yibin Municipal People’s Government University and City Strategic Cooperation Special Fund Project (2020CDYB-29)+1 种基金the Science and Technology Plan Transfer Payment Project of Sichuan Province (2021ZYSF007)the Key Research and Development Program of Science and Technology Department of Sichuan Province (2020YFS0575,2021KJT0012-2 021YFS-0067)。
文摘Classical localization methods use Cartesian or Polar coordinates, which require a priori range information to determine whether to estimate position or to only find bearings. The modified polar representation (MPR) unifies near-field and farfield models, alleviating the thresholding effect. Current localization methods in MPR based on the angle of arrival (AOA) and time difference of arrival (TDOA) measurements resort to semidefinite relaxation (SDR) and Gauss-Newton iteration, which are computationally complex and face the possible diverge problem. This paper formulates a pseudo linear equation between the measurements and the unknown MPR position,which leads to a closed-form solution for the hybrid TDOA-AOA localization problem, namely hybrid constrained optimization(HCO). HCO attains Cramér-Rao bound (CRB)-level accuracy for mild Gaussian noise. Compared with the existing closed-form solutions for the hybrid TDOA-AOA case, HCO provides comparable performance to the hybrid generalized trust region subproblem (HGTRS) solution and is better than the hybrid successive unconstrained minimization (HSUM) solution in large noise region. Its computational complexity is lower than that of HGTRS. Simulations validate the performance of HCO achieves the CRB that the maximum likelihood estimator (MLE) attains if the noise is small, but the MLE deviates from CRB earlier.
文摘We search for analytical wave solutions of an electronically and biologically important model named as the Fitzhugh–Nagumo model with truncated M-fractional derivative, in which the expafunction and extended sinh-Gordon equation expansion(ESh GEE) schemes are utilized. The solutions obtained include dark, bright, dark-bright, periodic and other kinds of solitons. These analytical wave solutions are gained and verified with the use of Mathematica software. These solutions do not exist in literature. Some of the solutions are demonstrated by 2D, 3D and contour graphs. This model is mostly used in circuit theory, transmission of nerve impulses, and population genetics. Finally, both the schemes are more applicable, reliable and significant to deal with the fractional nonlinear partial differential equations.
文摘Multiple objectives to be optimized simultaneously are prevalent in real-life problems. This paper develops a new Pareto Method for bi-objective optimization which yields analytical solutions. The Pareto optimal front is obtained in closed-form, enabling the derivation of various solutions in a convenient and efficient way. The advantage of analytical solution is the possibility of deriving accurate, exact and well-understood solutions, which is especially useful for policy analysis. An extension of the method to include multiple objectives is provided with the objectives being classified into two types. Such an extension expands the applicability of the developed techniques.
文摘This paper investigates the validity and shortcomings of the existing analytical solution for the ultimate bearing capacity of a pile embedded in a rock mass using the modified HoekeBrown failure criterion.Although this criterion is considered a reference value for empirical and numerical calculations,some limitations of its basic simplifications have not been clarified yet.This research compares the analytical results obtained from the novel discontinuity layout optimization(DLO)method and the numerical solutions from the finite difference method(FDM).The limitations of the analytical solution are considered by comparing different DLO failure modes,thus allowing for the first time a critical evaluation of its scope and conditioning for implementation.Errors of up to 40%in the bearing capacity and unrealistic failure modes are the main issues in the analytical solution.The main aspects of the DLO method are also analyzed with an emphasis on the linearization of the rock failure criterion and the accuracy resulting from the discretization size.The analysis demonstrates DLO as a very efficient and accurate tool to address the pile tip bearing capacity,presenting considerable advantages over other methods.
基金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.
基金This study was supported by Basic Research Project from Jiangmen Science and Technology Bureau(Grant No.2220002000356)China University of Petroleum(Beijing)(Grand No.2462023BJRC007)The Guangdong Basic and Applied Basic Research Foundation(No.2022A1515110376).
文摘In regard to unconventional oil reservoirs,the transient dual-porosity and triple-porosity models have been adopted to describe the fluid flow in the complex fracture network.It has been proven to cause inaccurate production evaluations because of the absence of matrix-macrofracture communication.In addition,most of the existing models are solved analytically based on Laplace transform and numerical inversion.Hence,an approximate analytical solution is derived directly in real-time space considering variable matrix blocks and simultaneous matrix depletion.To simplify the derivation,the simultaneous matrix depletion is divided into two parts:one part feeding the macrofractures and the other part feeding the microfractures.Then,a series of partial differential equations(PDEs)describing the transient flow and boundary conditions are constructed and solved analytically by integration.Finally,a relationship between oil rate and production time in real-time space is obtained.The new model is verified against classical analytical models.When the microfracture system and matrix-macrofracture communication is neglected,the result of the new model agrees with those obtained with the dual-porosity and triple-porosity model,respectively.Certainly,the new model also has an excellent agreement with the numerical model.The model is then applied to two actual tight oil wells completed in western Canada sedimentary basin.After identifying the flow regime,the solution suitably matches the field production data,and the model parameters are determined.Through these output parameters,we can accurately forecast the production and even estimate the petrophysical properties.
基金LMP acknowledges financial support from ANID through Convocatoria Nacional Subvención a Instalación en la Academia Convocatoria Año 2021,Grant SA77210040。
文摘We report on the magnetohydrodynamic impact on the axisymmetric flow of Al_(2)O_(3)/Cu nanoparticles suspended in H_(2)O past a stretched/shrinked sheet.With the use of partial differential equations and the corresponding thermophysical characteristics of nanoparticles,the physical flow process is illustrated.The resultant nonlinear system of partial differential equations is converted into a system of ordinary differential equations using the suitable similarity transformations.The transformed differential equations are solved analytically.Impacts of the magnetic parameter,solid volume fraction and stretching/shrinking parameter on momentum and temperature distribution have been analyzed and interpreted graphically.The skin friction and Nusselt number were also evaluated.In addition,existence of dual solution was deduced for the shrinking sheet and unique solution for the stretching one.Further,Al_(2)O_(3)/H_(2)O nanofluid flow has better thermal conductivity on comparing with Cu/H_(2)O nanofluid.Furthermore,it was found that the first solutions of the stream are stable and physically realizable,whereas those of the second ones are unstable.
基金Project supported by the National Natural Science Foundation of China (No. 12272323)。
文摘Second-order axially moving systems are common models in the field of dynamics, such as axially moving strings, cables, and belts. In the traditional research work, it is difficult to obtain closed-form solutions for the forced vibration when the damping effect and the coupling effect of multiple second-order models are considered.In this paper, Green's function method based on the Laplace transform is used to obtain closed-form solutions for the forced vibration of second-order axially moving systems. By taking the axially moving damping string system and multi-string system connected by springs as examples, the detailed solution methods and the analytical Green's functions of these second-order systems are given. The mode functions and frequency equations are also obtained by the obtained Green's functions. The reliability and convenience of the results are verified by several examples. This paper provides a systematic analytical method for the dynamic analysis of second-order axially moving systems, and the obtained Green's functions are applicable to different second-order systems rather than just string systems. In addition, the work of this paper also has positive significance for the study on the forced vibration of high-order systems.
基金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.
基金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.
基金supported by the National Natural Science Foundation of China (10772039 and 10632030)the National Basic Research Program of China (973 Program) (2010CB832704)
文摘This paper analyses the bending of rectangular orthotropic plates on a Winkler elastic foundation.Appropriate definition of symplectic inner product and symplectic space formed by generalized displacements establish dual variables and dual equations in the symplectic space.The operator matrix of the equation set is proven to be a Hamilton operator matrix.Separation of variables and eigenfunction expansion creates a basis for analyzing the bending of rectangular orthotropic plates on Winkler elastic foundation and obtaining solutions for plates having any boundary condition.There is discussion of symplectic eigenvalue problems of orthotropic plates under two typical boundary conditions,with opposite sides simply supported and opposite sides clamped.Transcendental equations of eigenvalues and symplectic eigenvectors in analytical form given.Analytical solutions using two examples are presented to show the use of the new methods described in this paper.To verify the accuracy and convergence,a fully simply supported plate that is fully and simply supported under uniformly distributed load is used to compare the classical Navier method,the Levy method and the new method.Results show that the new technique has good accuracy and better convergence speed than other methods,especially in relation to internal forces.A fully clamped rectangular plate on Winkler foundation is solved to validate application of the new methods,with solutions compared to those produced by the Galerkin method.
基金the National Natural Science Foundation of China(Nos.50679061,50709025and50749031)
文摘The theory of an eddy viscosity model is applied to the study of the flow in a compound channel which is partially vegetated. The governing equation is constituted by analyzing the longitudinal forces acting on the unit volume where the effect of the vegetation on the flow is considered as a drag force item, The compound channel is divided into 3 sub-regions in the transverse direction, and the coefficients in every region's differential equations were solved simultaneously. Thus, the analytical solution of the transverse distribution of the depth-averaged velocity for uniform flow in a partially vegetated compound channel was obtained. The results can be used to predict the transverse distribution of bed shear stress, which has an important effect on the transportation of sediment. By comparing the analytical results with the measured data, the analytical solution in this paper is shown to be sufficiently accurate to predict most hydraulic features for engineering design purposes.
文摘This paper presents an analytical solution of the one-dimensional consolidation in unsaturated soil with a finite thickness under vertical loading and confinements in the lateral directions. The boundary contains the top surface permeable to water and air and the bottom impermeable to water and air. The analytical solution is for Fredlund's one-dimensional consolidation equation in unsaturated soils. The transfer relationship between the state vectors at top surface and any depth is obtained by using the Laplace transform and Cayley-Hamilton mathematical methods to the governing equations of water and air, Darcy's law and Fick's law. Excess pore-air pressure, excess pore-water pressure and settlement in the Laplace-transformed domain are obtained by using the Laplace transform with the initial conditions and boundary conditions. By performing inverse Laplace transforms, the analytical solutions are obtained in the time domain. A typical example illustrates the consolidation characteristics of unsaturated soil from an- alytical results. Finally, comparisons between the analytical solutions and results of the finite difference method indicate that the analytical solution is correct.
基金supported by the National Natural Science Foundation of China (No.10772202)the Chinese PostdoctoralScience Foundation (No.20060400757).
文摘The stress concentration problem in structures with a circular or elliptic hole can be investigated by analytical methods. For the problem with a rectangular hole, only approximate results are derived. This paper deduces the analytical solutions to the stress concentration problem in plates with a rectangular hole under biaxial tensions. By using the U-transformation technique and the finite element method, the analytical displacement solutions of the finite element equations are derived in the series form. Therefore, the stress concentration can then be discussed easily and conveniently. For plate problem the bilinear rectangular element with four nodes is taken as an example to demonstrate the applicability of the proposed method. The stress concentration factors for various ratios of height to width of the hole are obtained.
基金supported by the National Natural Science Foundation of China(11102237)Program for Changjiang Scholars and Innovative Research Team in University(IRT1294)+1 种基金Specialized Research Fund for the Doctoral Program of Higher Education(20110133120012)China Scholarship Council(CSC)
文摘Based on Huang's accurate tri-sectional nonlin- ear kinematic equation (1997), a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction: The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation (PDE) sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient (TPG) has a great effect on the external moving boundary but has little effect on the internal moving boundary.
基金The project supported by NSF of Shandong Province
文摘According to the physiological and anatomical characteristics of small intestine,neglecting the effect of its motility on the distribution and absorption of drug and nutrient,Y.Miyamoto et al.proposed a model of two-dimensional laminar flow in a circular porous tube with permeable wall and calculated the concentration profile of drugby numerical analysis.In this paper,we give a steady-state analytical solution of the above model including deactivationterm.The obtained results are in agreement with the results of their numerical analysis. Moreover the analytical solution presented in this paper reveals the relation among the physiological parameters of the model and describes the basic absorption rule of drug and nutrient through the intestinal wall and hence pro- vides a theoretical basis for determining the permeability and reflection coefficient through in situ experiments.