Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an...Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.展开更多
In order to address the complex uncertainties caused by interfacing between the fuzziness and randomness of the safety problem for embankment engineering projects, and to evaluate the safety of embankment engineering ...In order to address the complex uncertainties caused by interfacing between the fuzziness and randomness of the safety problem for embankment engineering projects, and to evaluate the safety of embankment engineering projects more scientifically and reasonably, this study presents the fuzzy logic modeling of the stochastic finite element method （SFEM） based on the harmonious finite element （HFE） technique using a first-order approximation theorem. Fuzzy mathematical models of safety repertories were introduced into the SFEM to analyze the stability of embankments and foundations in order to describe the fuzzy failure procedure for the random safety performance function. The fuzzy models were developed with membership functions with half depressed gamma distribution, half depressed normal distribution, and half depressed echelon distribution. The fuzzy stochastic mathematical algorithm was used to comprehensively study the local failure mechanism of the main embankment section near Jingnan in the Yangtze River in terms of numerical analysis for the probability integration of reliability on the random field affected by three fuzzy factors. The result shows that the middle region of the embankment is the principal zone of concentrated failure due to local fractures. There is also some local shear failure on the embankment crust. This study provides a referential method for solving complex multi-uncertainty problems in engineering safety analysis.展开更多
To evaluate the water storage and project the future evolution of glaciers, the ice-thickness of glaciers is an essential input. However, direct measurements of ice thickness are labo- rious, not feasible everywhere, ...To evaluate the water storage and project the future evolution of glaciers, the ice-thickness of glaciers is an essential input. However, direct measurements of ice thickness are labo- rious, not feasible everywhere, and necessarily restricted to a small number of glaciers. In this article, we develop a simple method to estimate the ice-thickness along flow-line of mountain glaciers. Different from the traditional method based on shallow ice approximation （SIA）, which gives a relationship be- tween ice thickness, surface slope, and yield stress of glaciers, the improved method considers and pre- sents a simple way to calibrate the influence of valley wall on ice discharge. The required inputs are the glacier surface topography and outlines. This shows the potential of the method for estimating the ice-thickness distribution and volume of glaciers without using of direct thickness measurements.展开更多
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.展开更多
This paper introduces a novel approach for parameter sensitivity evaluation and efficient slope reliability analysis based on quantile-based first-order second-moment method(QFOSM).The core principles of the QFOSM are...This paper introduces a novel approach for parameter sensitivity evaluation and efficient slope reliability analysis based on quantile-based first-order second-moment method(QFOSM).The core principles of the QFOSM are elucidated geometrically from the perspective of expanding ellipsoids.Based on this geometric interpretation,the QFOSM is further extended to estimate sensitivity indices and assess the significance of various uncertain parameters involved in the slope system.The proposed method has the advantage of computational simplicity,akin to the conventional first-order second-moment method(FOSM),while providing estimation accuracy close to that of the first-order reliability method(FORM).Its performance is demonstrated with a numerical example and three slope examples.The results show that the proposed method can efficiently estimate the slope reliability and simultaneously evaluate the sensitivity of the uncertain parameters.The proposed method does not involve complex optimization or iteration required by the FORM.It can provide a valuable complement to the existing approximate reliability analysis methods,offering rapid sensitivity evaluation and slope reliability analysis.展开更多
Multidisciplinary feasible method (MDF) is conventional method to multidisciplinary optimization (MDO) and well-understood by users. It reduces the dimensions of the multidisciplinary optimization problem by using the...Multidisciplinary feasible method (MDF) is conventional method to multidisciplinary optimization (MDO) and well-understood by users. It reduces the dimensions of the multidisciplinary optimization problem by using the design variables as independent optimization variables. However, at each iteration of the conventional optimization procedure, multidisciplinary analysis (MDA) is numerously performed that results in extreme expense and low optimization efficiency. The intrinsic weakness of MDF is due to the times that it loop fixed-point iterations in MDA, which drive us to improve MDF by building inexpensive approximations as surrogates for expensive MDA. An simple example is presented to demonstrate the usefulness of the improved MDF. Results show that a significant reduction in the number of multidisciplinary analysis required for optimization is obtained as compared with original MDF and the efficiency of optimization is increased.展开更多
The shape approximation method has been proven to be rapid and practicable in resolving low-thrust trajectory;however,it still faces the challenges of large deviation from the optimal solution and inability to satisfy...The shape approximation method has been proven to be rapid and practicable in resolving low-thrust trajectory;however,it still faces the challenges of large deviation from the optimal solution and inability to satisfy the specific flight time and fuel mass constraints.In this paper,a modified shape approximation low-thrust model is presented,and a novel constrained optimization algorithm is developed to solve this problem.The proposed method aims at settling the bi-objective optimization orbit involving the twin objectives of minimum flight time and low fuel consumption and enhancing the accuracy of optimized orbit.In particular,a transformed high-order polynomial model based on finite Fourier series is proposed,which can be characterized as a multi-constraint optimization problem.Then,a novel optimization algorithm is specifically developed to optimize the large-scale multi-constraint dynamical equations of shape trajectory.The key performance indicators of the index include minimum flight time,low fuel consumption and bi-objective optimization of the two.Simulation results prove that this approach possesses both the high precision achievable by numerical methods and low computational complexity offered by shape approximation techniques.Besides,the Pareto front of the fuel-time bi-objective optimization orbit is firstly introduced to analyze an intact optimal solution set.Furthermore,we have demonstrated that our proposed approach is appropriate to generate the preliminary orbit for pseudo-spectral method.展开更多
In this paper, we firstly define a decreasing sequence {Pn(S)} by the generation of the Sierpinski gasket where each Pn(S) can be obtained in finite steps. Then we prove that the Hausdorff measure Hs(S) of the Sierpin...In this paper, we firstly define a decreasing sequence {Pn(S)} by the generation of the Sierpinski gasket where each Pn(S) can be obtained in finite steps. Then we prove that the Hausdorff measure Hs(S) of the Sierpinski gasket S can be approximated by {Pn(S)} with Pn(S)/(l + l/2n-3)s≤Hs(S)≤ Pn(S). An algorithm is presented to get Pn(S) for n ≤5. As an application, we obtain the best lower bound of Hs(S) till now: Hs(S)≥0.5631.展开更多
The aim of this paper is to employ fractional order proportional integral derivative(FO-PID)controller and integer order PID controller to control the position of the levitated object in a magnetic levitation system(M...The aim of this paper is to employ fractional order proportional integral derivative(FO-PID)controller and integer order PID controller to control the position of the levitated object in a magnetic levitation system(MLS),which is inherently nonlinear and unstable system.The proposal is to deploy discrete optimal pole-zero approximation method for realization of digital fractional order controller.An approach of phase shaping by slope cancellation of asymptotic phase plots for zeros and poles within given bandwidth is explored.The controller parameters are tuned using dynamic particle swarm optimization(d PSO)technique.Effectiveness of the proposed control scheme is verified by simulation and experimental results.The performance of realized digital FO-PID controller has been compared with that of the integer order PID controllers.It is observed that effort required in fractional order control is smaller as compared with its integer counterpart for obtaining the same system performance.展开更多
An intrinsic extension of Pad′e approximation method, called the generalized Pad′e approximation method, is proposed based on the classic Pad′e approximation theorem. According to the proposed method, the numerator...An intrinsic extension of Pad′e approximation method, called the generalized Pad′e approximation method, is proposed based on the classic Pad′e approximation theorem. According to the proposed method, the numerator and denominator of Pad′e approximant are extended from polynomial functions to a series composed of any kind of function, which means that the generalized Pad′e approximant is not limited to some forms, but can be constructed in different forms in solving different problems. Thus, many existing modifications of Pad′e approximation method can be considered to be the special cases of the proposed method. For solving homoclinic and heteroclinic orbits of strongly nonlinear autonomous oscillators, two novel kinds of generalized Pad′e approximants are constructed. Then, some examples are given to show the validity of the present method. To show the accuracy of the method, all solutions obtained in this paper are compared with those of the Runge–Kutta method.展开更多
The analysis technology of Amplitude Variation with Offset(AVO)is one of the important methods for oil and gas reservoir prediction.Zoeppritz equation and its approximations are the theoretical basis of AVO analysis,w...The analysis technology of Amplitude Variation with Offset(AVO)is one of the important methods for oil and gas reservoir prediction.Zoeppritz equation and its approximations are the theoretical basis of AVO analysis,which assumes that the upper and lower media of a horizontal interface are single-phase media.Limited by this assumption,AVO analysis has limited prediction and identification accuracy for complex porous reservoirs.In view of this,the first-order approximate analytical expressions of oblique elastic wave at an interface of porous media are derived.Firstly,the incident and scattering characteristics of various waves at the interface of porous media are analyzed,and the displacement vectors generated by these elastic waves are described by exponential function.Secondly,the kinematic and dynamic boundary conditions at the interface of porous media are discussed.Thirdly,by substituting the displacement vectors of incident and scattered waves into boundary conditions,the exact analytical equation is derived.Then,considering the symmetry of scattering matrix in the equation,the exact analytical expressions of each scattered wave are obtained.Furthermore,under the assumptions of small incident angle,weak elasticity at an interface of porous media,and ignoring the second-and higherorder terms,the first-order approximate analytical expressions are derived.Establishing a model of sandstone porous media with different porosity in upper and lower media,the correctness of the approximate analytical expressions is verified,and the elastic wave response characteristics of lithology and pore fluids are analyzed.展开更多
Suppose thatλ_(1),λ_(2),λ_(3),λ_(4),λ_(5)are nonzero real numbers,not all of the same sign,andλ_(1)/λ_(2)is irrational and algebraic.Let V be a well-spaced sequence,δ>0.In this paper,it is proved that,for ...Suppose thatλ_(1),λ_(2),λ_(3),λ_(4),λ_(5)are nonzero real numbers,not all of the same sign,andλ_(1)/λ_(2)is irrational and algebraic.Let V be a well-spaced sequence,δ>0.In this paper,it is proved that,for anyε>0,the number of v∈V with v≤N such that the following inequality|λ_(1)p_(1)~2+λ_(2)p_(2)~2+λ_(3)p_(3)~4+λ_(4)p_(4)~4+λ_5p_5~4-v|<v^(-δ)has no solution in prime variables p_(1),p_(2),p_(3),p_(4),p_(5)does not exceed O(N^(29/32+2δ+ε)).展开更多
Let K be a closed convex subset of a real reflexive Banach space E, T：K→K be a nonexpansive mapping, and f：K→K be a fixed weakly contractive （may not be contractive） mapping. Then for any t∈（0, 1）, let x1∈K ...Let K be a closed convex subset of a real reflexive Banach space E, T：K→K be a nonexpansive mapping, and f：K→K be a fixed weakly contractive （may not be contractive） mapping. Then for any t∈（0, 1）, let x1∈K be the unique fixed point of the weak contraction x1→tf（x）＋（1-t）Tx. If T has a fixed point and E admits a weakly sequentially continuous duality mapping from E to E^＊, then it is shown that {xt} converges to a fixed point of T as t→0. The results presented here improve and generalize the corresponding results in （Xu, 2004）.展开更多
In many industrial applications,heat transfer and tangent hyperbolic fluid flow processes have been garnering increasing attention,owing to their immense importance in technology,engineering,and science.These processe...In many industrial applications,heat transfer and tangent hyperbolic fluid flow processes have been garnering increasing attention,owing to their immense importance in technology,engineering,and science.These processes are relevant for polymer solutions,porous industrial materials,ceramic processing,oil recovery,and fluid beds.The present tangent hyperbolic fluid flow and heat transfer model accurately predicts the shear-thinning phenomenon and describes the blood flow characteristics.Therefore,the entropy production analysis of a non-Newtonian tangent hyperbolic material flow through a vertical microchannel with a quadratic density temperature fluctuation(quadratic/nonlinear Boussinesq approximation)is performed in the present study.The impacts of the hydrodynamic flow and Newton’s thermal conditions on the flow,heat transfer,and entropy generation are analyzed.The governing nonlinear equations are solved with the spectral quasi-linearization method(SQLM).The obtained results are compared with those calculated with a finite element method and the bvp4c routine.In addition,the effects of key parameters on the velocity of the hyperbolic tangent material,the entropy generation,the temperature,and the Nusselt number are discussed.The entropy generation increases with the buoyancy force,the pressure gradient factor,the non-linear convection,and the Eckert number.The non-Newtonian fluid factor improves the magnitude of the velocity field.The power-law index of the hyperbolic fluid and the Weissenberg number are found to be favorable for increasing the temperature field.The buoyancy force caused by the nonlinear change in the fluid density versus temperature improves the thermal energy of the system.展开更多
By combining the classical appropriate functions “1, x, x 2” with the method of multiplier enlargement, this paper establishes a theorem to approximate any unbounded continuous functions with modified positive...By combining the classical appropriate functions “1, x, x 2” with the method of multiplier enlargement, this paper establishes a theorem to approximate any unbounded continuous functions with modified positive linear operators. As an example, Hermite Fejér interpolation polynomial operators are analysed and studied, and a general conclusion is obtained.展开更多
A formulation of a differential equation as projection and fixed point pi-Mem alloivs approximations using general piecnvise functions. We prone existence and uniqueness of the up proximate solution* convergence in th...A formulation of a differential equation as projection and fixed point pi-Mem alloivs approximations using general piecnvise functions. We prone existence and uniqueness of the up proximate solution* convergence in the L2 norm and nodal supercnnvergence. These results generalize those obtained earlier by Hulme for continuous piecevjise polynomials and by Delfour-Dubeau for discontinuous pieceuiise polynomials. A duality relationship for the two types of approximations is also given.展开更多
The purpose of this paper is to propose and study local spline approximation methods for singular product integration,for which;i)the precision degree is the highest possible using splint approximation; ii) the nodes ...The purpose of this paper is to propose and study local spline approximation methods for singular product integration,for which;i)the precision degree is the highest possible using splint approximation; ii) the nodes fan be assumed equal to arbitrary points,where the integrand function f is known; iii) the number of the requested evaluations of f at the nodes is low,iv) a satisfactory convergence theory can be proved.展开更多
We use the method of discrete dipole approximation with surface interaction to construct a model in which a plurality of nanoparticles is arranged on the surface of BK7 glass. Nanoparticles are in air medium illuminat...We use the method of discrete dipole approximation with surface interaction to construct a model in which a plurality of nanoparticles is arranged on the surface of BK7 glass. Nanoparticles are in air medium illuminated by evanescent wave generated from total internal reflection. The effects of the wavelength, the polarization of the incident wave, the number of nanoparticles and the spacing of multiple nanoparticles on the field enhancement and extinction efficiency are calculated by our model. Our work could pave the way to improve the field enhancement of multiple nanoparticles systems.展开更多
Since the research of flare slamming prediction is seldom when parametric rolling happens, we present an efficient approximation method for flare slamming analysis of large container ships in parametric rolling condit...Since the research of flare slamming prediction is seldom when parametric rolling happens, we present an efficient approximation method for flare slamming analysis of large container ships in parametric rolling conditions. We adopt a 6-DOF weakly nonlinear time domain model to predict the ship motions of parametric rolling conditions. Unlike previous flare slamming analysis, our proposed method takes roll motion into account to calculate the impact angle and relative vertical velocity between ship sections on the bow flare and wave surface. We use the Wagner model to analyze the slamming impact forces and the slamming occurrence probability. Through numerical simulations, we investigate the maximum flare slamming pressures of a container ship for different speeds and wave conditions. To further clarify the mechanism of flare slamming phenomena in parametric rolling conditions, we also conduct real-time simulations to determine the relationship between slamming pressure and 3-DOF motions, namely roll, pitch, and heave.展开更多
In this work, by choosing an orthonormal basis for the Hilbert space L^2[0, 1], an approximation method for finding approximate solutions of the equation （I ＋ K）x = y is proposed, called Haar wavelet approximation ...In this work, by choosing an orthonormal basis for the Hilbert space L^2[0, 1], an approximation method for finding approximate solutions of the equation （I ＋ K）x = y is proposed, called Haar wavelet approximation method （HWAM）. To prove the applicabifity of the HWAM, a more general applicability theorem on an approximation method （AM） for an operator equation Ax = y is proved first. As an application, applicability of the HWAM is obtained. Fhrthermore, four steps to use the HWAM are listed and three numerical examples are given in order to illustrate the effectiveness of the method.展开更多
基金Supported by the National Natural Science Foundation of China (Grant Nos.52088102 and 51879287)National Key Research and Development Program of China (Grant No.2022YFB2602301)。
文摘Long-term responses of floating structures pose a great concern in their design phase. Existing approaches for addressing long-term extreme responses are extremely cumbersome for adoption. This work aims to develop an approach for the long-term extreme-response analysis of floating structures. A modified gradient-based retrieval algorithm in conjunction with the inverse first-order reliability method(IFORM) is proposed to enable the use of convolution models in long-term extreme analysis of structures with an analytical formula of response amplitude operator(RAO). The proposed algorithm ensures convergence stability and iteration accuracy and exhibits a higher computational efficiency than the traditional backtracking method. However, when the RAO of general offshore structures cannot be analytically expressed, the convolutional integration method fails to function properly. A numerical discretization approach is further proposed for offshore structures in the case when the analytical expression of the RAO is not feasible. Through iterative discretization of environmental contours(ECs) and RAOs, a detailed procedure is proposed to calculate the long-term response extremes of offshore structures. The validity and accuracy of the proposed approach are tested using a floating offshore wind turbine as a numerical example. The long-term extreme heave responses of various return periods are calculated via the IFORM in conjunction with a numerical discretization approach. The environmental data corresponding to N-year structural responses are located inside the ECs, which indicates that the selection of design points directly along the ECs yields conservative design results.
基金supported by the National Natural Science Foundation of China(Grant No.50379046)the Doctoral Fund of the Ministry of Education of China(Grant No.A50221)
文摘In order to address the complex uncertainties caused by interfacing between the fuzziness and randomness of the safety problem for embankment engineering projects, and to evaluate the safety of embankment engineering projects more scientifically and reasonably, this study presents the fuzzy logic modeling of the stochastic finite element method （SFEM） based on the harmonious finite element （HFE） technique using a first-order approximation theorem. Fuzzy mathematical models of safety repertories were introduced into the SFEM to analyze the stability of embankments and foundations in order to describe the fuzzy failure procedure for the random safety performance function. The fuzzy models were developed with membership functions with half depressed gamma distribution, half depressed normal distribution, and half depressed echelon distribution. The fuzzy stochastic mathematical algorithm was used to comprehensively study the local failure mechanism of the main embankment section near Jingnan in the Yangtze River in terms of numerical analysis for the probability integration of reliability on the random field affected by three fuzzy factors. The result shows that the middle region of the embankment is the principal zone of concentrated failure due to local fractures. There is also some local shear failure on the embankment crust. This study provides a referential method for solving complex multi-uncertainty problems in engineering safety analysis.
基金supported by the National Basic Research Program of China (No. 2007CB411501)the Knowledge Innovation Project of the Chinese Academy of Sciences (No. KZCX2-EW-311)+1 种基金the National Natural Science Foundation of China (Nos. 91025012, J0930003/J0109)the Project for Outstanding Young Scientists of the National Natural Science Foundation of China (No. 40121101)
文摘To evaluate the water storage and project the future evolution of glaciers, the ice-thickness of glaciers is an essential input. However, direct measurements of ice thickness are labo- rious, not feasible everywhere, and necessarily restricted to a small number of glaciers. In this article, we develop a simple method to estimate the ice-thickness along flow-line of mountain glaciers. Different from the traditional method based on shallow ice approximation （SIA）, which gives a relationship be- tween ice thickness, surface slope, and yield stress of glaciers, the improved method considers and pre- sents a simple way to calibrate the influence of valley wall on ice discharge. The required inputs are the glacier surface topography and outlines. This shows the potential of the method for estimating the ice-thickness distribution and volume of glaciers without using of direct thickness measurements.
文摘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.
基金supported by the National Natural Science Foundation of China(Grant Nos.52109144,52025094 and 52222905).
文摘This paper introduces a novel approach for parameter sensitivity evaluation and efficient slope reliability analysis based on quantile-based first-order second-moment method(QFOSM).The core principles of the QFOSM are elucidated geometrically from the perspective of expanding ellipsoids.Based on this geometric interpretation,the QFOSM is further extended to estimate sensitivity indices and assess the significance of various uncertain parameters involved in the slope system.The proposed method has the advantage of computational simplicity,akin to the conventional first-order second-moment method(FOSM),while providing estimation accuracy close to that of the first-order reliability method(FORM).Its performance is demonstrated with a numerical example and three slope examples.The results show that the proposed method can efficiently estimate the slope reliability and simultaneously evaluate the sensitivity of the uncertain parameters.The proposed method does not involve complex optimization or iteration required by the FORM.It can provide a valuable complement to the existing approximate reliability analysis methods,offering rapid sensitivity evaluation and slope reliability analysis.
文摘Multidisciplinary feasible method (MDF) is conventional method to multidisciplinary optimization (MDO) and well-understood by users. It reduces the dimensions of the multidisciplinary optimization problem by using the design variables as independent optimization variables. However, at each iteration of the conventional optimization procedure, multidisciplinary analysis (MDA) is numerously performed that results in extreme expense and low optimization efficiency. The intrinsic weakness of MDF is due to the times that it loop fixed-point iterations in MDA, which drive us to improve MDF by building inexpensive approximations as surrogates for expensive MDA. An simple example is presented to demonstrate the usefulness of the improved MDF. Results show that a significant reduction in the number of multidisciplinary analysis required for optimization is obtained as compared with original MDF and the efficiency of optimization is increased.
基金supported by the National Natural Science Foundation of China(Nos.61627810,61790562,61403096).
文摘The shape approximation method has been proven to be rapid and practicable in resolving low-thrust trajectory;however,it still faces the challenges of large deviation from the optimal solution and inability to satisfy the specific flight time and fuel mass constraints.In this paper,a modified shape approximation low-thrust model is presented,and a novel constrained optimization algorithm is developed to solve this problem.The proposed method aims at settling the bi-objective optimization orbit involving the twin objectives of minimum flight time and low fuel consumption and enhancing the accuracy of optimized orbit.In particular,a transformed high-order polynomial model based on finite Fourier series is proposed,which can be characterized as a multi-constraint optimization problem.Then,a novel optimization algorithm is specifically developed to optimize the large-scale multi-constraint dynamical equations of shape trajectory.The key performance indicators of the index include minimum flight time,low fuel consumption and bi-objective optimization of the two.Simulation results prove that this approach possesses both the high precision achievable by numerical methods and low computational complexity offered by shape approximation techniques.Besides,the Pareto front of the fuel-time bi-objective optimization orbit is firstly introduced to analyze an intact optimal solution set.Furthermore,we have demonstrated that our proposed approach is appropriate to generate the preliminary orbit for pseudo-spectral method.
文摘In this paper, we firstly define a decreasing sequence {Pn(S)} by the generation of the Sierpinski gasket where each Pn(S) can be obtained in finite steps. Then we prove that the Hausdorff measure Hs(S) of the Sierpinski gasket S can be approximated by {Pn(S)} with Pn(S)/(l + l/2n-3)s≤Hs(S)≤ Pn(S). An algorithm is presented to get Pn(S) for n ≤5. As an application, we obtain the best lower bound of Hs(S) till now: Hs(S)≥0.5631.
基金supported by the Board of Research in Nuclear Sciences of the Department of Atomic Energy,India(2012/36/69-BRNS/2012)
文摘The aim of this paper is to employ fractional order proportional integral derivative(FO-PID)controller and integer order PID controller to control the position of the levitated object in a magnetic levitation system(MLS),which is inherently nonlinear and unstable system.The proposal is to deploy discrete optimal pole-zero approximation method for realization of digital fractional order controller.An approach of phase shaping by slope cancellation of asymptotic phase plots for zeros and poles within given bandwidth is explored.The controller parameters are tuned using dynamic particle swarm optimization(d PSO)technique.Effectiveness of the proposed control scheme is verified by simulation and experimental results.The performance of realized digital FO-PID controller has been compared with that of the integer order PID controllers.It is observed that effort required in fractional order control is smaller as compared with its integer counterpart for obtaining the same system performance.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11172093 and 11372102)the Hunan Provincial Innovation Foundation for Postgraduate,China(Grant No.CX2012B159)
文摘An intrinsic extension of Pad′e approximation method, called the generalized Pad′e approximation method, is proposed based on the classic Pad′e approximation theorem. According to the proposed method, the numerator and denominator of Pad′e approximant are extended from polynomial functions to a series composed of any kind of function, which means that the generalized Pad′e approximant is not limited to some forms, but can be constructed in different forms in solving different problems. Thus, many existing modifications of Pad′e approximation method can be considered to be the special cases of the proposed method. For solving homoclinic and heteroclinic orbits of strongly nonlinear autonomous oscillators, two novel kinds of generalized Pad′e approximants are constructed. Then, some examples are given to show the validity of the present method. To show the accuracy of the method, all solutions obtained in this paper are compared with those of the Runge–Kutta method.
基金financially supported by the National Natural Science Foundation of China(Grant No.42104131)the Natural Science Foundation of Sichuan Province of China(Grant No.2022NSFSC1140)Open Fund(PLC20211101)of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation
文摘The analysis technology of Amplitude Variation with Offset(AVO)is one of the important methods for oil and gas reservoir prediction.Zoeppritz equation and its approximations are the theoretical basis of AVO analysis,which assumes that the upper and lower media of a horizontal interface are single-phase media.Limited by this assumption,AVO analysis has limited prediction and identification accuracy for complex porous reservoirs.In view of this,the first-order approximate analytical expressions of oblique elastic wave at an interface of porous media are derived.Firstly,the incident and scattering characteristics of various waves at the interface of porous media are analyzed,and the displacement vectors generated by these elastic waves are described by exponential function.Secondly,the kinematic and dynamic boundary conditions at the interface of porous media are discussed.Thirdly,by substituting the displacement vectors of incident and scattered waves into boundary conditions,the exact analytical equation is derived.Then,considering the symmetry of scattering matrix in the equation,the exact analytical expressions of each scattered wave are obtained.Furthermore,under the assumptions of small incident angle,weak elasticity at an interface of porous media,and ignoring the second-and higherorder terms,the first-order approximate analytical expressions are derived.Establishing a model of sandstone porous media with different porosity in upper and lower media,the correctness of the approximate analytical expressions is verified,and the elastic wave response characteristics of lithology and pore fluids are analyzed.
基金Supported by NSFC(Nos.12301006,12471009,12071238,11901566,12001047,11971476)Beijing Natural Science Foundation(No.1242003)。
文摘Suppose thatλ_(1),λ_(2),λ_(3),λ_(4),λ_(5)are nonzero real numbers,not all of the same sign,andλ_(1)/λ_(2)is irrational and algebraic.Let V be a well-spaced sequence,δ>0.In this paper,it is proved that,for anyε>0,the number of v∈V with v≤N such that the following inequality|λ_(1)p_(1)~2+λ_(2)p_(2)~2+λ_(3)p_(3)~4+λ_(4)p_(4)~4+λ_5p_5~4-v|<v^(-δ)has no solution in prime variables p_(1),p_(2),p_(3),p_(4),p_(5)does not exceed O(N^(29/32+2δ+ε)).
文摘Let K be a closed convex subset of a real reflexive Banach space E, T：K→K be a nonexpansive mapping, and f：K→K be a fixed weakly contractive （may not be contractive） mapping. Then for any t∈（0, 1）, let x1∈K be the unique fixed point of the weak contraction x1→tf（x）＋（1-t）Tx. If T has a fixed point and E admits a weakly sequentially continuous duality mapping from E to E^＊, then it is shown that {xt} converges to a fixed point of T as t→0. The results presented here improve and generalize the corresponding results in （Xu, 2004）.
文摘In many industrial applications,heat transfer and tangent hyperbolic fluid flow processes have been garnering increasing attention,owing to their immense importance in technology,engineering,and science.These processes are relevant for polymer solutions,porous industrial materials,ceramic processing,oil recovery,and fluid beds.The present tangent hyperbolic fluid flow and heat transfer model accurately predicts the shear-thinning phenomenon and describes the blood flow characteristics.Therefore,the entropy production analysis of a non-Newtonian tangent hyperbolic material flow through a vertical microchannel with a quadratic density temperature fluctuation(quadratic/nonlinear Boussinesq approximation)is performed in the present study.The impacts of the hydrodynamic flow and Newton’s thermal conditions on the flow,heat transfer,and entropy generation are analyzed.The governing nonlinear equations are solved with the spectral quasi-linearization method(SQLM).The obtained results are compared with those calculated with a finite element method and the bvp4c routine.In addition,the effects of key parameters on the velocity of the hyperbolic tangent material,the entropy generation,the temperature,and the Nusselt number are discussed.The entropy generation increases with the buoyancy force,the pressure gradient factor,the non-linear convection,and the Eckert number.The non-Newtonian fluid factor improves the magnitude of the velocity field.The power-law index of the hyperbolic fluid and the Weissenberg number are found to be favorable for increasing the temperature field.The buoyancy force caused by the nonlinear change in the fluid density versus temperature improves the thermal energy of the system.
文摘By combining the classical appropriate functions “1, x, x 2” with the method of multiplier enlargement, this paper establishes a theorem to approximate any unbounded continuous functions with modified positive linear operators. As an example, Hermite Fejér interpolation polynomial operators are analysed and studied, and a general conclusion is obtained.
基金This research has been supported in part by the Natural Sciences and Engineering Research Council of Canada(Grant OGPIN-336)and by the"Ministere de l'Education du Quebec"(FCAR Grant-ER-0725)
文摘A formulation of a differential equation as projection and fixed point pi-Mem alloivs approximations using general piecnvise functions. We prone existence and uniqueness of the up proximate solution* convergence in the L2 norm and nodal supercnnvergence. These results generalize those obtained earlier by Hulme for continuous piecevjise polynomials and by Delfour-Dubeau for discontinuous pieceuiise polynomials. A duality relationship for the two types of approximations is also given.
基金Work sponsored by"Ministero dell' University"CNR of Italy
文摘The purpose of this paper is to propose and study local spline approximation methods for singular product integration,for which;i)the precision degree is the highest possible using splint approximation; ii) the nodes fan be assumed equal to arbitrary points,where the integrand function f is known; iii) the number of the requested evaluations of f at the nodes is low,iv) a satisfactory convergence theory can be proved.
基金supported by the Zhejiang Provincial Natural Science Foundation of China(No.LGF20C050001)the National Nature Science Foundation of China(No.61805213)。
文摘We use the method of discrete dipole approximation with surface interaction to construct a model in which a plurality of nanoparticles is arranged on the surface of BK7 glass. Nanoparticles are in air medium illuminated by evanescent wave generated from total internal reflection. The effects of the wavelength, the polarization of the incident wave, the number of nanoparticles and the spacing of multiple nanoparticles on the field enhancement and extinction efficiency are calculated by our model. Our work could pave the way to improve the field enhancement of multiple nanoparticles systems.
基金supported by the ChinaMinistry of Education Key Research Project "KSHIP-II Project"(Knowledge-based Ship Design Hyper-Integrated Platform)No.GKZY010004the National Key Basic Research Program of China No.2014CB046804
文摘Since the research of flare slamming prediction is seldom when parametric rolling happens, we present an efficient approximation method for flare slamming analysis of large container ships in parametric rolling conditions. We adopt a 6-DOF weakly nonlinear time domain model to predict the ship motions of parametric rolling conditions. Unlike previous flare slamming analysis, our proposed method takes roll motion into account to calculate the impact angle and relative vertical velocity between ship sections on the bow flare and wave surface. We use the Wagner model to analyze the slamming impact forces and the slamming occurrence probability. Through numerical simulations, we investigate the maximum flare slamming pressures of a container ship for different speeds and wave conditions. To further clarify the mechanism of flare slamming phenomena in parametric rolling conditions, we also conduct real-time simulations to determine the relationship between slamming pressure and 3-DOF motions, namely roll, pitch, and heave.
基金support by the NSFC(11371012,11401359,11471200)the FRF for the Central Universities(GK201301007)the NSRP of Shaanxi Province(2014JQ1010)
文摘In this work, by choosing an orthonormal basis for the Hilbert space L^2[0, 1], an approximation method for finding approximate solutions of the equation （I ＋ K）x = y is proposed, called Haar wavelet approximation method （HWAM）. To prove the applicabifity of the HWAM, a more general applicability theorem on an approximation method （AM） for an operator equation Ax = y is proved first. As an application, applicability of the HWAM is obtained. Fhrthermore, four steps to use the HWAM are listed and three numerical examples are given in order to illustrate the effectiveness of the method.
