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.展开更多
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.展开更多
Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
Convection driven by a spatially non-uniform internal heat source between two horizontal isothermal walls is studied by theoretical analysis and numerical simulation,in order to explore the bounds of the temperature a...Convection driven by a spatially non-uniform internal heat source between two horizontal isothermal walls is studied by theoretical analysis and numerical simulation,in order to explore the bounds of the temperature and the vertical heat flux.Specifically,the rigorous lower bound of the weighted average temperature<QT>is derived analytically,by decomposing the temperature field into a background profile and a fluctuation part.This bound obtained for the first time to consider non-uniform heat sources is found to be compatible with the existing bound obtained in uniform internal heat convection.Of physical importance,an analytical relationship is derived as an inequality connecting<QT>and the average vertical heat flux<wT>,by employing the average heat flux on the bottom wall(qb)as an intermediary variable.It clarifies the intrinsic relation between the lower bound of<QT>and the upper bound of<wT>,namely,these two bounds are essentially equivalent providing an easy way to obtain one from another.Furthermore,the analytical bounds are extensively demonstrated through a comprehensive series of direct numerical simulations.展开更多
A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illus...A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.展开更多
Two new versions of accelerated first-order methods for minimizing convex composite functions are proposed. In this paper, we first present an accelerated first-order method which chooses the step size 1/ Lk to be 1/ ...Two new versions of accelerated first-order methods for minimizing convex composite functions are proposed. In this paper, we first present an accelerated first-order method which chooses the step size 1/ Lk to be 1/ L0 at the beginning of each iteration and preserves the computational simplicity of the fast iterative shrinkage-thresholding algorithm. The first proposed algorithm is a non-monotone algorithm. To avoid this behavior, we present another accelerated monotone first-order method. The proposed two accelerated first-order methods are proved to have a better convergence rate for minimizing convex composite functions. Numerical results demonstrate the efficiency of the proposed two accelerated first-order methods.展开更多
The surrounding rock is prone to large-scale loosening and failure after the excavation of shallow large-span caverns because of the thin overlying strata and large cross-section span.The rational design of bolt suppo...The surrounding rock is prone to large-scale loosening and failure after the excavation of shallow large-span caverns because of the thin overlying strata and large cross-section span.The rational design of bolt support is very important to the safety control of surrounding rock as a common support means.The control mechanism and design method of bolt support for shallow-buried large-span caverns is carried out.The calculation method of bolt prestress and length based on arched failure and collapsed failure mode is established.The influence mechanism of different influencing factors on the bolt prestress and length is clarified.At the same time,the constant resistance energy-absorbing bolt with high strength and high toughness is developed,and the comparative test of mechanical properties is carried out.On this basis,the design method of high prestressed bolt support for shallow-buried large-span caverns is put forward,and the field test is carried out in Qingdao metro station in China.The monitoring results show that the maximum roof settlement is 6.8 mm after the new design method is adopted,and the effective control of the shallow-buried large-span caverns is realized.The research results can provide theoretical and technical support for the safety control of shallow-buried large-span caverns.展开更多
In this article, we use the spin coherent state transformation and the ground state variational method to theoretically calculate the ground function. In order to consider the influence of the atom-atom interaction on...In this article, we use the spin coherent state transformation and the ground state variational method to theoretically calculate the ground function. In order to consider the influence of the atom-atom interaction on the extended Dicke model's ground state properties, the mean photon number, the scaled atomic population and the average ground energy are displayed. Using the self-consistent field theory to solve the atom-atom interaction, we discover the system undergoes a first-order quantum phase transition from the normal phase to the superradiant phase, but a famous Dicke-type second-order quantum phase transition without the atom-atom interaction. Meanwhile, the atom-atom interaction makes the phase transition point shift to the lower atom-photon collective coupling strength.展开更多
Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic s...Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.展开更多
The prediction of central bursting defects in the rod extrusion process through conical dies using the upper bound analysisis investigated. A kinematically admissible velocity field, including the radial and angular v...The prediction of central bursting defects in the rod extrusion process through conical dies using the upper bound analysisis investigated. A kinematically admissible velocity field, including the radial and angular velocity components, is proposed. A newcriterion is presented to predict the occurrence of the central bursting defects. Parameter bobt, which represents the risk probability ofcracking, is proposed. It is calculated using the shape of the boundary at the entrance by minimizing the total power dissipationduring the extrusion process. When bobt is equal to or greater than bcr, central bursting occurs. Furthermore, the quantitativerelationships between central bursting defects and process parameters （semi die angle, reduction in area and frictional factor） arestudied. The results show that the central bursting defects are affected primarily by the reduction in area and the friction factor. Thepresented criterion is verified by comparing with the FEM simulation data and the results of the published paper.展开更多
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.展开更多
In this paper,we formulate two new families of fourth-order explicit exponential Runge–Kutta(ERK)methods with four stages for solving first-order differential systems y'(t)+M y(t)=f(y(t)).The order conditions of ...In this paper,we formulate two new families of fourth-order explicit exponential Runge–Kutta(ERK)methods with four stages for solving first-order differential systems y'(t)+M y(t)=f(y(t)).The order conditions of these ERK methods are derived by comparing the Taylor series of the exact solution,which are exactly identical to the order conditions of explicit Runge–Kutta methods,and these ERK methods reduce to classical Runge–Kutta methods once M→0.Moreover,we analyze the stability properties and the convergence of these new methods.Several numerical examples are implemented to illustrate the accuracy and efficiency of these ERK methods by comparison with standard exponential integrators.展开更多
At present, associated flow rule of traditional plastic theory is adopted in the slip line field theory and upper bound method of geotechnical materials. So the stress characteristic line conforms to the velocity line...At present, associated flow rule of traditional plastic theory is adopted in the slip line field theory and upper bound method of geotechnical materials. So the stress characteristic line conforms to the velocity line. It is proved that geotechnical materials do not abide by the associated flow rule. It is impossible for the stress characteristic line to conform to the velocity line. Generalized plastic mechanics theoretically proved that plastic potential surface intersects the Mohr-Coulomb yield surface with an angle, so that the velocity line must be studied by non-associated flow rule. According to limit analysis theory, the theory of slip line field is put forward in this paper, and then the ultimate beating capacity of strip footing is obtained based on the associated flow rule and the non-associated flow nile individually. These two results are identical since the ultimate bearing capacity is independent of flow role. On the contrary, the velocity fields of associated and non-associated flow rules are different which shows the velocity field based on the associat- ed flow rule is incorrect.展开更多
In this paper, we combine Graeffe matrices with the classical numerical method of Dandelin-Graeffe to estimate bounds for the moduli of the zeros of polynomials. Furthermore, we give some examples showing significant ...In this paper, we combine Graeffe matrices with the classical numerical method of Dandelin-Graeffe to estimate bounds for the moduli of the zeros of polynomials. Furthermore, we give some examples showing significant gain for the convergence towards the polynomials dominant zeros moduli.展开更多
Two non-probabilistic, set-theoretical methods for determining the maximum and minimum impulsive responses of structures to uncertain-but-bounded impulses are presented. They are, respectively, based on the theories o...Two non-probabilistic, set-theoretical methods for determining the maximum and minimum impulsive responses of structures to uncertain-but-bounded impulses are presented. They are, respectively, based on the theories of interval mathematics and convex models. The uncertain-but-bounded impulses are assumed to be a convex set, hyper-rectangle or ellipsoid. For the two non-probabilistic methods, less prior information is required about the uncertain nature of impulses than the probabilistic model. Comparisons between the interval analysis method and the convex model, which are developed as an anti-optimization problem of finding the least favorable impulsive response and the most favorable impulsive response, are made through mathematical analyses and numerical calculations. The results of this study indicate that under the condition of the interval vector being determined from an ellipsoid containing the uncertain impulses, the width of the impulsive responses predicted by the interval analysis method is larger than that by the convex model; under the condition of the ellipsoid being determined from an interval vector containing the uncertain impulses, the width of the interval impulsive responses obtained by the interval analysis method is smaller than that by the convex model.展开更多
In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial deriv...In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.展开更多
In this paper we report a sparse truncated Newton algorithm for handling large-scale simple bound nonlinear constrained minimixation problem. The truncated Newton method is used to update the variables with indices ou...In this paper we report a sparse truncated Newton algorithm for handling large-scale simple bound nonlinear constrained minimixation problem. The truncated Newton method is used to update the variables with indices outside of the active set, while the projected gradient method is used to update the active variables. At each iterative level, the search direction consists of three parts, one of which is a subspace truncated Newton direction, the other two are subspace gradient and modified gradient directions. The subspace truncated Newton direction is obtained by solving a sparse system of linear equations. The global convergence and quadratic convergence rate of the algorithm are proved and some numerical tests are given.展开更多
In this paper, we obtain the Berry-Esseen bound for identically distributed random variables by Stein method. The results obtained generalize the results of Shao and Su （2006） and Stein （1986）.
The interaction and its variation between magnetic grains in two kinds of magnetic recording tapes are investigated by first-order reversal curves （FORC） and the 5M method. The composition and microstructure of the ...The interaction and its variation between magnetic grains in two kinds of magnetic recording tapes are investigated by first-order reversal curves （FORC） and the 5M method. The composition and microstructure of the samples are characterised by x-ray diffraction and scanning electron microscope. The FORC diagram can provide more accurate information of the interaction and its variation, but the 5M curves cannot. The positive interaction field and the large variation of the interaction field have opposite effects on the δM curve.展开更多
基金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 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.
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
基金supported by the National Natural Science Foundation of China(Grant Nos.92252202,92152301,12293000,12293002,12302320,and 12388101)the Fundamental Research Funds for the Central Universities.
文摘Convection driven by a spatially non-uniform internal heat source between two horizontal isothermal walls is studied by theoretical analysis and numerical simulation,in order to explore the bounds of the temperature and the vertical heat flux.Specifically,the rigorous lower bound of the weighted average temperature<QT>is derived analytically,by decomposing the temperature field into a background profile and a fluctuation part.This bound obtained for the first time to consider non-uniform heat sources is found to be compatible with the existing bound obtained in uniform internal heat convection.Of physical importance,an analytical relationship is derived as an inequality connecting<QT>and the average vertical heat flux<wT>,by employing the average heat flux on the bottom wall(qb)as an intermediary variable.It clarifies the intrinsic relation between the lower bound of<QT>and the upper bound of<wT>,namely,these two bounds are essentially equivalent providing an easy way to obtain one from another.Furthermore,the analytical bounds are extensively demonstrated through a comprehensive series of direct numerical simulations.
文摘A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.
基金Sponsored by the National Natural Science Foundation of China(Grant No.11461021)the Natural Science Basic Research Plan in Shaanxi Province of China(Grant No.2017JM1014)
文摘Two new versions of accelerated first-order methods for minimizing convex composite functions are proposed. In this paper, we first present an accelerated first-order method which chooses the step size 1/ Lk to be 1/ L0 at the beginning of each iteration and preserves the computational simplicity of the fast iterative shrinkage-thresholding algorithm. The first proposed algorithm is a non-monotone algorithm. To avoid this behavior, we present another accelerated monotone first-order method. The proposed two accelerated first-order methods are proved to have a better convergence rate for minimizing convex composite functions. Numerical results demonstrate the efficiency of the proposed two accelerated first-order methods.
基金Project(2023YFC3805700) supported by the National Key Research and Development Program of ChinaProjects(42477166,42277174) supported by the National Natural Science Foundation of China+2 种基金Project(2024JCCXSB01) supported by the Fundamental Research Funds for the Central Universities,ChinaProject(KFJJ24-01M) supported by the State Key Laboratory of Explosion Science and Safety Protection,Beijing Institute of Technology,ChinaProject(HLCX-2024-04) supported by the Open Foundation of Collaborative Innovation Center of Green Development and Ecological Restoration of Mineral Resources,China。
文摘The surrounding rock is prone to large-scale loosening and failure after the excavation of shallow large-span caverns because of the thin overlying strata and large cross-section span.The rational design of bolt support is very important to the safety control of surrounding rock as a common support means.The control mechanism and design method of bolt support for shallow-buried large-span caverns is carried out.The calculation method of bolt prestress and length based on arched failure and collapsed failure mode is established.The influence mechanism of different influencing factors on the bolt prestress and length is clarified.At the same time,the constant resistance energy-absorbing bolt with high strength and high toughness is developed,and the comparative test of mechanical properties is carried out.On this basis,the design method of high prestressed bolt support for shallow-buried large-span caverns is put forward,and the field test is carried out in Qingdao metro station in China.The monitoring results show that the maximum roof settlement is 6.8 mm after the new design method is adopted,and the effective control of the shallow-buried large-span caverns is realized.The research results can provide theoretical and technical support for the safety control of shallow-buried large-span caverns.
基金Supported by the National Natural Science Foundation of China under Grant Nos.11275118,11404198,91430109,61505100,51502189the Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi Province(STIP)under Grant No.2014102+2 种基金the Launch of the Scientific Research of Shanxi University under Grant No.011151801004the National Fundamental Fund of Personnel Training under Grant No.J1103210The Natural Science Foundation of Shanxi Province under Grant No.2015011008
文摘In this article, we use the spin coherent state transformation and the ground state variational method to theoretically calculate the ground function. In order to consider the influence of the atom-atom interaction on the extended Dicke model's ground state properties, the mean photon number, the scaled atomic population and the average ground energy are displayed. Using the self-consistent field theory to solve the atom-atom interaction, we discover the system undergoes a first-order quantum phase transition from the normal phase to the superradiant phase, but a famous Dicke-type second-order quantum phase transition without the atom-atom interaction. Meanwhile, the atom-atom interaction makes the phase transition point shift to the lower atom-photon collective coupling strength.
文摘Based on the lower bound theorem of limit analysis, a solution procedure for limit analysis of three_dimensional elastoplastic structures was established using conventional boundary element method (BEM). The elastic stress field for lower bound limit analysis was computed directly by three_dimensional boundary element method (3_D BEM). The self_equilibrium stress field was constructed by the linear combination of several self_equilibrium “basis vectors” which can be computed by elastic_plastic incremental iteration of 3_D BEM analysis. The lower bound limit analysis problem was finally reduced to a series of nonlinear programming sub_problems with relatively few optimal variables. The complex method was used to solve the nonlinear programming sub_problems. The numerical results show that the present solution procedure has good accuracy and high efficiency.
文摘The prediction of central bursting defects in the rod extrusion process through conical dies using the upper bound analysisis investigated. A kinematically admissible velocity field, including the radial and angular velocity components, is proposed. A newcriterion is presented to predict the occurrence of the central bursting defects. Parameter bobt, which represents the risk probability ofcracking, is proposed. It is calculated using the shape of the boundary at the entrance by minimizing the total power dissipationduring the extrusion process. When bobt is equal to or greater than bcr, central bursting occurs. Furthermore, the quantitativerelationships between central bursting defects and process parameters （semi die angle, reduction in area and frictional factor） arestudied. The results show that the central bursting defects are affected primarily by the reduction in area and the friction factor. Thepresented criterion is verified by comparing with the FEM simulation data and the results of the published paper.
基金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 NSFC(Grant No.12071419)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2024MA056)。
文摘In this paper,we formulate two new families of fourth-order explicit exponential Runge–Kutta(ERK)methods with four stages for solving first-order differential systems y'(t)+M y(t)=f(y(t)).The order conditions of these ERK methods are derived by comparing the Taylor series of the exact solution,which are exactly identical to the order conditions of explicit Runge–Kutta methods,and these ERK methods reduce to classical Runge–Kutta methods once M→0.Moreover,we analyze the stability properties and the convergence of these new methods.Several numerical examples are implemented to illustrate the accuracy and efficiency of these ERK methods by comparison with standard exponential integrators.
文摘At present, associated flow rule of traditional plastic theory is adopted in the slip line field theory and upper bound method of geotechnical materials. So the stress characteristic line conforms to the velocity line. It is proved that geotechnical materials do not abide by the associated flow rule. It is impossible for the stress characteristic line to conform to the velocity line. Generalized plastic mechanics theoretically proved that plastic potential surface intersects the Mohr-Coulomb yield surface with an angle, so that the velocity line must be studied by non-associated flow rule. According to limit analysis theory, the theory of slip line field is put forward in this paper, and then the ultimate beating capacity of strip footing is obtained based on the associated flow rule and the non-associated flow nile individually. These two results are identical since the ultimate bearing capacity is independent of flow role. On the contrary, the velocity fields of associated and non-associated flow rules are different which shows the velocity field based on the associat- ed flow rule is incorrect.
文摘In this paper, we combine Graeffe matrices with the classical numerical method of Dandelin-Graeffe to estimate bounds for the moduli of the zeros of polynomials. Furthermore, we give some examples showing significant gain for the convergence towards the polynomials dominant zeros moduli.
基金The project supported by the National Outstanding Youth Science Foundation of China (10425208)the National Natural Science Foundation of ChinaInstitute of Engineering Physics of China (10376002) The English text was polished by Keren Wang
文摘Two non-probabilistic, set-theoretical methods for determining the maximum and minimum impulsive responses of structures to uncertain-but-bounded impulses are presented. They are, respectively, based on the theories of interval mathematics and convex models. The uncertain-but-bounded impulses are assumed to be a convex set, hyper-rectangle or ellipsoid. For the two non-probabilistic methods, less prior information is required about the uncertain nature of impulses than the probabilistic model. Comparisons between the interval analysis method and the convex model, which are developed as an anti-optimization problem of finding the least favorable impulsive response and the most favorable impulsive response, are made through mathematical analyses and numerical calculations. The results of this study indicate that under the condition of the interval vector being determined from an ellipsoid containing the uncertain impulses, the width of the impulsive responses predicted by the interval analysis method is larger than that by the convex model; under the condition of the ellipsoid being determined from an interval vector containing the uncertain impulses, the width of the interval impulsive responses obtained by the interval analysis method is smaller than that by the convex model.
文摘In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.
基金The research was supported by the State Education Grant for Retumed Scholars
文摘In this paper we report a sparse truncated Newton algorithm for handling large-scale simple bound nonlinear constrained minimixation problem. The truncated Newton method is used to update the variables with indices outside of the active set, while the projected gradient method is used to update the active variables. At each iterative level, the search direction consists of three parts, one of which is a subspace truncated Newton direction, the other two are subspace gradient and modified gradient directions. The subspace truncated Newton direction is obtained by solving a sparse system of linear equations. The global convergence and quadratic convergence rate of the algorithm are proved and some numerical tests are given.
基金Supported by the National Natural Science Foundation of China (11101364)the Zhejiang Natural Science Foundation of China (Y6110110)
文摘In this paper, we obtain the Berry-Esseen bound for identically distributed random variables by Stein method. The results obtained generalize the results of Shao and Su （2006） and Stein （1986）.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 50672008 and 50971023)Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20090006120019)
文摘The interaction and its variation between magnetic grains in two kinds of magnetic recording tapes are investigated by first-order reversal curves （FORC） and the 5M method. The composition and microstructure of the samples are characterised by x-ray diffraction and scanning electron microscope. The FORC diagram can provide more accurate information of the interaction and its variation, but the 5M curves cannot. The positive interaction field and the large variation of the interaction field have opposite effects on the δM curve.
