Kirchhoff integral migration imaging is widely used in industrial production due to its advantages of not being limited by observation systems,good target imaging effects,and high computational efficiency.Vertical tra...Kirchhoff integral migration imaging is widely used in industrial production due to its advantages of not being limited by observation systems,good target imaging effects,and high computational efficiency.Vertical transversely isotropic(VTI)media,as a typical anisotropic media,has always been a primary focus of anisotropic migration imaging research.We focus on the problems of low accuracy and efficiency in travel-time calculations associated with conventional Kirchhoff integral migration for VTI media.A travel-time calculation method based on physical-informed neural network(PINN)for VTI media is introduced into the process of Kirchhoff integral migration imaging.Model experiments and field data processing have shown that the travel-time calculation based on PINN can significantly improve both the accuracy and efficiency when compared to traditional finite difference algorithms,there-by enabling high-precision Kirchhoff integral migration imaging for VTI media.展开更多
The main purpose of this paper is to use the Chelyshkov-collocation spectral method for solving nonlinear Quadratic integral equations of Volterra type.The method is based on the approximate solutions in terms of Chel...The main purpose of this paper is to use the Chelyshkov-collocation spectral method for solving nonlinear Quadratic integral equations of Volterra type.The method is based on the approximate solutions in terms of Chelyshkov polynomials with unknown coefficients.The Chelyshkov polynomials and their properties are employed to derive the operational matrices of integral and product.The application of these operational matrices for solving the mentioned problem is explained.The error analysis of the proposed method is investigated.Finally,some numerical examples are provided to demonstrate the efficiency of the method.展开更多
Fatigue analysis of engine turbine blade is an essential issue.Due to various uncertainties during the manufacture and operation,the fatigue damage and life of turbine blade present randomness.In this study,the random...Fatigue analysis of engine turbine blade is an essential issue.Due to various uncertainties during the manufacture and operation,the fatigue damage and life of turbine blade present randomness.In this study,the randomness of structural parameters,working condition and vibration environment are considered for fatigue life predication and reliability assessment.First,the lowcycle fatigue problem is modelled as stochastic static system with random parameters,while the high-cycle fatigue problem is considered as stochastic dynamic system under random excitations.Then,to deal with the two failure modes,the novel Direct Probability Integral Method(DPIM)is proposed,which is efficient and accurate for solving stochastic static and dynamic systems.The probability density functions of accumulated damage and fatigue life of turbine blade for low-cycle and high-cycle fatigue problems are achieved,respectively.Furthermore,the time–frequency hybrid method is advanced to enhance the computational efficiency for governing equation of system.Finally,the results of typical examples demonstrate high accuracy and efficiency of the proposed method by comparison with Monte Carlo simulation and other methods.It is indicated that the DPIM is a unified method for predication of random fatigue life for low-cycle and highcycle fatigue problems.The rotational speed,density,fatigue strength coefficient,and fatigue plasticity index have a high sensitivity to fatigue reliability of engine turbine blade.展开更多
This paper extends the one-dimensional(1D)nonlocal strain gradient integral model(NStraGIM)to the two-dimensional(2D)Kirchhoff axisymmetric nanoplates,based on nonlocal strain gradient integral relations formulated al...This paper extends the one-dimensional(1D)nonlocal strain gradient integral model(NStraGIM)to the two-dimensional(2D)Kirchhoff axisymmetric nanoplates,based on nonlocal strain gradient integral relations formulated along both the radial and circumferential directions.By transforming the proposed integral constitutive equations into the equivalent differential forms,complemented by the corresponding constitutive boundary conditions(CBCs),a well-posed mathematical formulation is established for analyzing the axisymmetric bending and buckling of annular/circular functionally graded(FG)sandwich nanoplates.The boundary conditions at the inner edge of a solid nanoplate are derived by L'H?spital's rule.The numerical solution is obtained by the generalized differential quadrature method(GDQM).The accuracy of the proposed model is validated through comparison with the data from the existing literature.A parameter study is conducted to demonstrate the effects of FG sandwich parameters,size parameters,and nonlocal gradient parameters.展开更多
Due to the difficulties in obtaining large deformation mining subsidence using differential Interferometric Synthetic Aperture Radar （D-InSAR） alone, a new algorithm was proposed to extract large deformation mining ...Due to the difficulties in obtaining large deformation mining subsidence using differential Interferometric Synthetic Aperture Radar （D-InSAR） alone, a new algorithm was proposed to extract large deformation mining subsidence using D-InSAR technique and probability integral method. The details of the algorithm are as follows：the control points set, containing correct phase unwrapping points on the subsidence basin edge generated by D-InSAR and several observation points （near the maximum subsidence and inflection points）, was established at first; genetic algorithm （GA） was then used to optimize the parameters of probability integral method; at last, the surface subsidence was deduced according to the optimum parameters. The results of the experiment in Huaibei mining area, China, show that the presented method can generate the correct mining subsidence basin with a few surface observations, and the relative error of maximum subsidence point is about 8.3%, which is much better than that of conventional D-InSAR （relative error is 68.0%）.展开更多
In this paper, Kirchhoff formula has been transformed from surface integral form into a line integral form. The new form of the formula can be applied to separate geometrical optical field from diffraction field, and ...In this paper, Kirchhoff formula has been transformed from surface integral form into a line integral form. The new form of the formula can be applied to separate geometrical optical field from diffraction field, and reduce the time of numerical computation greatly. Based on the new form, an analytical formula of diffraction field in the far zone has been presented for the polygonal aperture illuminated by a uniform plane wave.展开更多
In this paper,a deep collocation method(DCM)for thin plate bending problems is proposed.This method takes advantage of computational graphs and backpropagation algorithms involved in deep learning.Besides,the proposed...In this paper,a deep collocation method(DCM)for thin plate bending problems is proposed.This method takes advantage of computational graphs and backpropagation algorithms involved in deep learning.Besides,the proposed DCM is based on a feedforward deep neural network(DNN)and differs from most previous applications of deep learning for mechanical problems.First,batches of randomly distributed collocation points are initially generated inside the domain and along the boundaries.A loss function is built with the aim that the governing partial differential equations(PDEs)of Kirchhoff plate bending problems,and the boundary/initial conditions are minimised at those collocation points.A combination of optimizers is adopted in the backpropagation process to minimize the loss function so as to obtain the optimal hyperparameters.In Kirchhoff plate bending problems,the C^1 continuity requirement poses significant difficulties in traditional mesh-based methods.This can be solved by the proposed DCM,which uses a deep neural network to approximate the continuous transversal deflection,and is proved to be suitable to the bending analysis of Kirchhoff plate of various geometries.展开更多
In order to improve rolled strip quality, precise plate shape control theory should be established. Roll flat- tening theory is an important part of the plate shape theory. To improve the accuracy of roll flattening c...In order to improve rolled strip quality, precise plate shape control theory should be established. Roll flat- tening theory is an important part of the plate shape theory. To improve the accuracy of roll flattening calculation based on semi infinite body model, especially near the two roll barrel edges, a new and more accurate roll flattening model is proposed. Based on boundary integral equation method, an analytical model for solving a finite length semi infinite body is established. The lateral surface displacement field of the finite length semi-infinite body is simulated by finite element method （FEM） and lateral surface displacement decay functions are established. Based on the boundary integral equation method, the numerical solution of the finite length semi-infinite body under the distribu ted force is obtained and an accurate roll flattening model is established. Different from the traditional semi-infinite body model, the matrix form of the new roll flattening model is established through the mathematical derivation. The result from the new model is more consistent with that by FEM especially near the edges.展开更多
The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the bas...The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.展开更多
An adaptive cell-based domain integration method（CDIM） is proposed for the treatment of domain integrals in 3D boundary element method（BEM）. The domain integrals are computed in background cells rather than volume...An adaptive cell-based domain integration method（CDIM） is proposed for the treatment of domain integrals in 3D boundary element method（BEM）. The domain integrals are computed in background cells rather than volume elements. The cells are created from the boundary elements based on an adaptive oct-tree structure and no other discretization is needed. Cells containing the boundary elements are subdivided into smaller sub-cells adaptively according to the sizes and levels of the boundary elements; and the sub-cells outside the domain are deleted to obtain the desired accuracy. The method is applied in the 3D potential and elasticity problems in this paper.展开更多
For deep tunnel projects,selecting an appropriate initial support distance is critical to improving the self-supporting capacity of surrounding rock.In this work,an intuitive method for determining the tunnel’s initi...For deep tunnel projects,selecting an appropriate initial support distance is critical to improving the self-supporting capacity of surrounding rock.In this work,an intuitive method for determining the tunnel’s initial support distance was proposed.First,based on the convergence-confinement method,a three-dimensional analytical model was constructed by combining an analytical solution of a non-circular tunnel with the Tecplot software.Then,according to the integral failure criteria of rock,the failure tendency coefficients of hard surrounding rock were computed and the spatial distribution plots of that were constructed.On this basis,the tunnel’s key failure positions were identified,and the relationship between the failure tendency coefficient at key failure positions and their distances from the working face was established.Finally,the distance from the working face that corresponds to the critical failure tendency coefficient was taken as the optimal support distance.A practical project was used as an example,and a reasonable initial support distance was successfully determined by applying the developed method.Moreover,it is found that the stability of hard surrounding rock decreases rapidly within the range of 1.0D(D is the tunnel diameter)from the working face,and tends to be stable outside the range of 1.0D.展开更多
On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the...On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software ＂Auto mapping system of submarine topography and geomorphology MBChart＂. This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China＇ s next Five-year Plan.展开更多
This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic So...This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.展开更多
The Z component and X component profiles of seismic waves extracted with the prestack Kirchhoff integral migration could approximate to the primary wave （P wave） and converted shear wave （PS wave） profiles under c...The Z component and X component profiles of seismic waves extracted with the prestack Kirchhoff integral migration could approximate to the primary wave （P wave） and converted shear wave （PS wave） profiles under certain conditions. The relative change of their reflection amplitude reflects the formation stress anomaly and subsurface media anisotropy. The principle and method for extracting amplitude ratios were studied and the application of amplitude ratio profiles was also examined when processing and interpreting actual seismic data. The amplitude ratio profile is an effective supplementary means of identifying the stratigraphic boundary and lithology.展开更多
For the longitudinal seismic response analysis of a tunnel structure under asynchronous earthquake excitations,a longitudinal integral response deformation method classified as a practical approach is proposed in this...For the longitudinal seismic response analysis of a tunnel structure under asynchronous earthquake excitations,a longitudinal integral response deformation method classified as a practical approach is proposed in this paper.The determinations of the structural critical moments when maximal deformations and internal forces in the longitudinal direction occur are deduced as well.When applying the proposed method,the static analysis of the free-field computation model subjected to the least favorable free-field deformation at the tunnel buried depth is performed first to calculate the equivalent input seismic loads.Then,the equivalent input seismic loads are imposed on the integral tunnel-foundation computation model to conduct the static calculation.Afterwards,the critical longitudinal seismic responses of the tunnel are obtained.The applicability of the new method is verified by comparing the seismic responses of a shield tunnel structure in Beijing,determined by the proposed procedure and by a dynamic time-history analysis under a series of obliquely incident out-of-plane and in-plane waves.The results show that the proposed method has a clear concept with high accuracy and simple progress.Meanwhile,this method provides a feasible way to determine the critical moments of the longitudinal seismic responses of a tunnel structure.Therefore,the proposed method can be effectively applied to analyze the seismic response of a long-line underground structure subjected to non-uniform excitations.展开更多
In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] a...In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] are used and are utilized as a basis in Galerkin method to approximate the solution of integral equations. Then, in some examples the mentioned wavelets are compared with each other.展开更多
A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous H...A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.展开更多
In this article, by introducing characteristic singular integral operator and associate singular integral equations （SIEs）, the authors discuss the direct method of solution for a class of singular integral equation...In this article, by introducing characteristic singular integral operator and associate singular integral equations （SIEs）, the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.展开更多
In this paper,the classical composite middle rectangle rule for the computation of Cauchy principal value integral(the singular kernel 1=(x-s))is discussed.With the density function approximated only while the singula...In this paper,the classical composite middle rectangle rule for the computation of Cauchy principal value integral(the singular kernel 1=(x-s))is discussed.With the density function approximated only while the singular kernel is calculated analysis,then the error functional of asymptotic expansion is obtained.We construct a series to approach the singular point.An extrapolation algorithm is presented and the convergence rate of extrapolation algorithm is proved.At last,some numerical results are presented to confirm the theoretical results and show the efficiency of the algorithms.展开更多
基金funded by the Open Fund of No.1 Institute of Geology and Mineral Resources of Shan-dong Province(No.2022DY03)the National Natural Sci-ence Foundation of China(Nos.42074138 and 42206195)the Wenhai Program of the ST Fund of Laoshan Labo-ratory(No.LSKJ202204803).
文摘Kirchhoff integral migration imaging is widely used in industrial production due to its advantages of not being limited by observation systems,good target imaging effects,and high computational efficiency.Vertical transversely isotropic(VTI)media,as a typical anisotropic media,has always been a primary focus of anisotropic migration imaging research.We focus on the problems of low accuracy and efficiency in travel-time calculations associated with conventional Kirchhoff integral migration for VTI media.A travel-time calculation method based on physical-informed neural network(PINN)for VTI media is introduced into the process of Kirchhoff integral migration imaging.Model experiments and field data processing have shown that the travel-time calculation based on PINN can significantly improve both the accuracy and efficiency when compared to traditional finite difference algorithms,there-by enabling high-precision Kirchhoff integral migration imaging for VTI media.
文摘The main purpose of this paper is to use the Chelyshkov-collocation spectral method for solving nonlinear Quadratic integral equations of Volterra type.The method is based on the approximate solutions in terms of Chelyshkov polynomials with unknown coefficients.The Chelyshkov polynomials and their properties are employed to derive the operational matrices of integral and product.The application of these operational matrices for solving the mentioned problem is explained.The error analysis of the proposed method is investigated.Finally,some numerical examples are provided to demonstrate the efficiency of the method.
基金supports of the National Natural Science Foundation of China(Nos.12032008,12102080)the Fundamental Research Funds for the Central Universities,China(No.DUT23RC(3)038)are much appreciated。
文摘Fatigue analysis of engine turbine blade is an essential issue.Due to various uncertainties during the manufacture and operation,the fatigue damage and life of turbine blade present randomness.In this study,the randomness of structural parameters,working condition and vibration environment are considered for fatigue life predication and reliability assessment.First,the lowcycle fatigue problem is modelled as stochastic static system with random parameters,while the high-cycle fatigue problem is considered as stochastic dynamic system under random excitations.Then,to deal with the two failure modes,the novel Direct Probability Integral Method(DPIM)is proposed,which is efficient and accurate for solving stochastic static and dynamic systems.The probability density functions of accumulated damage and fatigue life of turbine blade for low-cycle and high-cycle fatigue problems are achieved,respectively.Furthermore,the time–frequency hybrid method is advanced to enhance the computational efficiency for governing equation of system.Finally,the results of typical examples demonstrate high accuracy and efficiency of the proposed method by comparison with Monte Carlo simulation and other methods.It is indicated that the DPIM is a unified method for predication of random fatigue life for low-cycle and highcycle fatigue problems.The rotational speed,density,fatigue strength coefficient,and fatigue plasticity index have a high sensitivity to fatigue reliability of engine turbine blade.
基金Project supported by the National Natural Science Foundation of China(No.12172169)the Priority Academic Program Development of Jiangsu Higher Education Institutions。
文摘This paper extends the one-dimensional(1D)nonlocal strain gradient integral model(NStraGIM)to the two-dimensional(2D)Kirchhoff axisymmetric nanoplates,based on nonlocal strain gradient integral relations formulated along both the radial and circumferential directions.By transforming the proposed integral constitutive equations into the equivalent differential forms,complemented by the corresponding constitutive boundary conditions(CBCs),a well-posed mathematical formulation is established for analyzing the axisymmetric bending and buckling of annular/circular functionally graded(FG)sandwich nanoplates.The boundary conditions at the inner edge of a solid nanoplate are derived by L'H?spital's rule.The numerical solution is obtained by the generalized differential quadrature method(GDQM).The accuracy of the proposed model is validated through comparison with the data from the existing literature.A parameter study is conducted to demonstrate the effects of FG sandwich parameters,size parameters,and nonlocal gradient parameters.
基金Project (BK20130174) supported by the Basic Research Project of Jiangsu Province (Natural Science Foundation) Project (1101109C) supported by Jiangsu Planned Projects for Postdoctoral Research Funds,China+1 种基金Project (201325) supported by the Key Laboratory of Geo-informatics of State Bureau of Surveying and Mapping,ChinaProject (SZBF2011-6-B35) supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions,China
文摘Due to the difficulties in obtaining large deformation mining subsidence using differential Interferometric Synthetic Aperture Radar （D-InSAR） alone, a new algorithm was proposed to extract large deformation mining subsidence using D-InSAR technique and probability integral method. The details of the algorithm are as follows：the control points set, containing correct phase unwrapping points on the subsidence basin edge generated by D-InSAR and several observation points （near the maximum subsidence and inflection points）, was established at first; genetic algorithm （GA） was then used to optimize the parameters of probability integral method; at last, the surface subsidence was deduced according to the optimum parameters. The results of the experiment in Huaibei mining area, China, show that the presented method can generate the correct mining subsidence basin with a few surface observations, and the relative error of maximum subsidence point is about 8.3%, which is much better than that of conventional D-InSAR （relative error is 68.0%）.
文摘In this paper, Kirchhoff formula has been transformed from surface integral form into a line integral form. The new form of the formula can be applied to separate geometrical optical field from diffraction field, and reduce the time of numerical computation greatly. Based on the new form, an analytical formula of diffraction field in the far zone has been presented for the polygonal aperture illuminated by a uniform plane wave.
文摘In this paper,a deep collocation method(DCM)for thin plate bending problems is proposed.This method takes advantage of computational graphs and backpropagation algorithms involved in deep learning.Besides,the proposed DCM is based on a feedforward deep neural network(DNN)and differs from most previous applications of deep learning for mechanical problems.First,batches of randomly distributed collocation points are initially generated inside the domain and along the boundaries.A loss function is built with the aim that the governing partial differential equations(PDEs)of Kirchhoff plate bending problems,and the boundary/initial conditions are minimised at those collocation points.A combination of optimizers is adopted in the backpropagation process to minimize the loss function so as to obtain the optimal hyperparameters.In Kirchhoff plate bending problems,the C^1 continuity requirement poses significant difficulties in traditional mesh-based methods.This can be solved by the proposed DCM,which uses a deep neural network to approximate the continuous transversal deflection,and is proved to be suitable to the bending analysis of Kirchhoff plate of various geometries.
基金Item Sponsored by National Natural Science Foundation of China(51075353)
文摘In order to improve rolled strip quality, precise plate shape control theory should be established. Roll flat- tening theory is an important part of the plate shape theory. To improve the accuracy of roll flattening calculation based on semi infinite body model, especially near the two roll barrel edges, a new and more accurate roll flattening model is proposed. Based on boundary integral equation method, an analytical model for solving a finite length semi infinite body is established. The lateral surface displacement field of the finite length semi-infinite body is simulated by finite element method （FEM） and lateral surface displacement decay functions are established. Based on the boundary integral equation method, the numerical solution of the finite length semi-infinite body under the distribu ted force is obtained and an accurate roll flattening model is established. Different from the traditional semi-infinite body model, the matrix form of the new roll flattening model is established through the mathematical derivation. The result from the new model is more consistent with that by FEM especially near the edges.
文摘The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.
基金Financial support for the project from the National Natural Science Foundation of China(No.51609181)
文摘An adaptive cell-based domain integration method（CDIM） is proposed for the treatment of domain integrals in 3D boundary element method（BEM）. The domain integrals are computed in background cells rather than volume elements. The cells are created from the boundary elements based on an adaptive oct-tree structure and no other discretization is needed. Cells containing the boundary elements are subdivided into smaller sub-cells adaptively according to the sizes and levels of the boundary elements; and the sub-cells outside the domain are deleted to obtain the desired accuracy. The method is applied in the 3D potential and elasticity problems in this paper.
基金Project(2021JLM-49) supported by Natural Science Basic Research Program of Shaanxi-Joint Fund of Hanjiang to Weihe River Valley Water Diversion Project,ChinaProject(42077248) supported by the National Natural Science Foundation of China
文摘For deep tunnel projects,selecting an appropriate initial support distance is critical to improving the self-supporting capacity of surrounding rock.In this work,an intuitive method for determining the tunnel’s initial support distance was proposed.First,based on the convergence-confinement method,a three-dimensional analytical model was constructed by combining an analytical solution of a non-circular tunnel with the Tecplot software.Then,according to the integral failure criteria of rock,the failure tendency coefficients of hard surrounding rock were computed and the spatial distribution plots of that were constructed.On this basis,the tunnel’s key failure positions were identified,and the relationship between the failure tendency coefficient at key failure positions and their distances from the working face was established.Finally,the distance from the working face that corresponds to the critical failure tendency coefficient was taken as the optimal support distance.A practical project was used as an example,and a reasonable initial support distance was successfully determined by applying the developed method.Moreover,it is found that the stability of hard surrounding rock decreases rapidly within the range of 1.0D(D is the tunnel diameter)from the working face,and tends to be stable outside the range of 1.0D.
基金This study was supported by Projects under contract Nos DY105 China's 0cean-03-01-01 and DY105-03-01-07the National Natural Science Foundation of China under contract No.40506017the Youth Foundation of Marine High-tech Project of China under contract No.2002AA616010.
文摘On the basis of three geological models and several orebody boundaries, a method of grid subdivision and integral has been proposed to calculate and evaluate the resources of cobalt-rich crusts on the seamounts in the central Pacific Ocean. The formulas of this method are deduced and the interface of program module is designed. The method is carried out in the software ＂Auto mapping system of submarine topography and geomorphology MBChart＂. This method and program will possibly become a potential tool to calculate the resources of seamounts and determine the target diggings for China＇ s next Five-year Plan.
基金Project supported by the Natural Science Foundation of China(10371009)Research Fund for the Doctoral Program Higher Education
文摘This paper determines the exact error order on optimization of adaptive direct methods of approximate solution of the class of Fredholm integral equations of the second kind with kernel belonging to the anisotropic Sobolev classes, and also gives an optimal algorithm.
文摘The Z component and X component profiles of seismic waves extracted with the prestack Kirchhoff integral migration could approximate to the primary wave （P wave） and converted shear wave （PS wave） profiles under certain conditions. The relative change of their reflection amplitude reflects the formation stress anomaly and subsurface media anisotropy. The principle and method for extracting amplitude ratios were studied and the application of amplitude ratio profiles was also examined when processing and interpreting actual seismic data. The amplitude ratio profile is an effective supplementary means of identifying the stratigraphic boundary and lithology.
基金National Natural Science Foundation of China under Grant No.51478247。
文摘For the longitudinal seismic response analysis of a tunnel structure under asynchronous earthquake excitations,a longitudinal integral response deformation method classified as a practical approach is proposed in this paper.The determinations of the structural critical moments when maximal deformations and internal forces in the longitudinal direction occur are deduced as well.When applying the proposed method,the static analysis of the free-field computation model subjected to the least favorable free-field deformation at the tunnel buried depth is performed first to calculate the equivalent input seismic loads.Then,the equivalent input seismic loads are imposed on the integral tunnel-foundation computation model to conduct the static calculation.Afterwards,the critical longitudinal seismic responses of the tunnel are obtained.The applicability of the new method is verified by comparing the seismic responses of a shield tunnel structure in Beijing,determined by the proposed procedure and by a dynamic time-history analysis under a series of obliquely incident out-of-plane and in-plane waves.The results show that the proposed method has a clear concept with high accuracy and simple progress.Meanwhile,this method provides a feasible way to determine the critical moments of the longitudinal seismic responses of a tunnel structure.Therefore,the proposed method can be effectively applied to analyze the seismic response of a long-line underground structure subjected to non-uniform excitations.
文摘In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] are used and are utilized as a basis in Galerkin method to approximate the solution of integral equations. Then, in some examples the mentioned wavelets are compared with each other.
文摘A direct method to find the first integral for two-dimensional autonomous system in polar coordinates is suggested. It is shown that if the equation of motion expressed by differential 1-forms for a given autonomous Hamiltonian system is multiplied by a set of multiplicative functions, then the general expression of the first integral can be obtained, An example is given to illustrate the application of the results.
基金Project was supported by RFDP of Higher Education and NNSF of China, SF of Wuhan University
文摘In this article, by introducing characteristic singular integral operator and associate singular integral equations （SIEs）, the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.
基金The work of Jin Li was supported by National Natural Science Foundation of China(Grant No.11471195)China Postdoctoral Science Foundation(Grant No.2015T80703)+1 种基金Shan-dong Provincial Natural Science Foundation of China(Grant No.ZR2016JL006)Na-tional Natural Science Foundation of China(Grant No.11771398).
文摘In this paper,the classical composite middle rectangle rule for the computation of Cauchy principal value integral(the singular kernel 1=(x-s))is discussed.With the density function approximated only while the singular kernel is calculated analysis,then the error functional of asymptotic expansion is obtained.We construct a series to approach the singular point.An extrapolation algorithm is presented and the convergence rate of extrapolation algorithm is proved.At last,some numerical results are presented to confirm the theoretical results and show the efficiency of the algorithms.
