On the basis of introducing the fundamental theory and the basic analysis steps of the sub model method, the strength of the new engine complex assembly structure was analyzed according to the properties of the engin...On the basis of introducing the fundamental theory and the basic analysis steps of the sub model method, the strength of the new engine complex assembly structure was analyzed according to the properties of the engine structures, some of the key parts of the engine were analyzed with refined mesh by sub model method and the error of the FEM solution was estimated by the extrapolation method. The example showed that the sub model can not only analyze the comlex structures without the restriction of the software and hardware of the computers, but get the more precise analysis result also. This method is more suitable for the strength analysis of the complex assembly structure.展开更多
To study the uncertainty quantification of resonant states in open quantum systems,we developed a Bayesian framework by integrating a reduced basis method(RBM)emulator with the Gamow coupled-channel(GCC)approach.The R...To study the uncertainty quantification of resonant states in open quantum systems,we developed a Bayesian framework by integrating a reduced basis method(RBM)emulator with the Gamow coupled-channel(GCC)approach.The RBM,constructed via eigenvector continuation and trained on both bound and resonant configurations,enables the fast and accurate emulation of resonance properties across the parameter space.To identify the physical resonant states from the emulator’s output,we introduce an overlap-based selection technique that effectively isolates true solutions from background artifacts.By applying this framework to unbound nucleus ^(6)Be,we quantified the model uncertainty in the predicted complex energies.The results demonstrate relative errors of 17.48%in the real part and 8.24%in the imaginary part,while achieving a speedup of four orders of magnitude compared with the full GCC calculations.To further investigate the asymptotic behavior of the resonant-state wavefunctions within the RBM framework,we employed a Lippmann–Schwinger(L–S)-based correction scheme.This approach not only improves the consistency between eigenvalues and wavefunctions but also enables a seamless extension from real-space training data to the complex energy plane.By bridging the gap between bound-state and continuum regimes,the L–S correction significantly enhances the emulator’s capability to accurately capture continuum structures in open quantum systems.展开更多
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.展开更多
Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate ...Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate for actual situations.Although soil exhibits material nonlinearity under earthquake excitation,existing research on the TI medium is limited to soil linearity and neglects the nonlinear response of TI sites.A 2D equivalent linear model for a layered TI half-space subjected to seismic waves is derived in the transformed wave number domain using the exact dynamic stiffness matrix of the TI medium.This study introduces a method for determining the effective shear strain of TI sites under oblique wave incidence,and further describes a systematic study on the effects of TI parameters and soil nonlinearity on site responses.Numerical results indicate that seismic responses of the TI medium significantly differ from those of isotropic sites and that the responses are highly dependent on TI parameters,particularly in nonlinear cases,while also being sensitive to incident angle and excitation intensity.Moreover,the differences in peak acceleration and waveform for various TI materials may also be amplified due to the strong nonlinearity.The study provides valuable insights for improving the accuracy of seismic response analysis in engineering applications.展开更多
This paper aims to solve the resonance failure probability and develop an effective method to estimate the effects of variables and failure modes on failure probability of axially functionally graded material(FGM)pipe...This paper aims to solve the resonance failure probability and develop an effective method to estimate the effects of variables and failure modes on failure probability of axially functionally graded material(FGM)pipe conveying fluid.Correspondingly,the natural frequency of axially FGM pipes conveying fluid is calculated using the differential quadrature method(DQM).A variable sensitivity analysis(VSA)is introduced to measure the effect of each random variable,and a mode sensitivity analysis(MSA)is introduced to acquire the importance ranking of failure modes.Then,an active learning Kriging(ALK)method is established to calculate the resonance failure probability and sensitivity indices,which greatly improves the application of resonance reliability analysis for pipelines in engineering practice.Based on the resonance reliability analysis method,the effects of fluid velocity,volume fraction and fluid density of axially FGM pipe conveying fluid on resonance reliability are analyzed.The results demonstrate that the proposed method has great performance in the anti-resonance analysis of pipes conveying fluid.展开更多
In this paper,the homotopy analysis method (HAM) is applied to solve generalized biological populationmodels.The fractional derivatives are described by Caputo's sense.The method introduces a significant improveme...In this paper,the homotopy analysis method (HAM) is applied to solve generalized biological populationmodels.The fractional derivatives are described by Caputo's sense.The method introduces a significant improvementin this field over existing techniques.Results obtained using the scheme presented here agree well with the analyticalsolutions and the numerical results presented in Ref.[6].However,the fundamental solutions of these equations stillexhibit useful scaling properties that make them attractive for applications.展开更多
In order to identify the day and night pollution sources of PM10 in ambient air in Longyan City,the authors analyzed the elemental composition of respirable particulate matters in the day and night ambient air samples...In order to identify the day and night pollution sources of PM10 in ambient air in Longyan City,the authors analyzed the elemental composition of respirable particulate matters in the day and night ambient air samples and various pollution sources which were collected in January 2010 in Longyan with inductivity coupled plasma-mass spectrometry（ICP-MS）.Then chemical mass balance（CMB） model and factor analysis（FA） method were applied to comparatively study the inorganic components in the sources and receptor samples.The results of factor analysis show that the major sources were road dust,waste incineration and mixed sources which contained automobile exhaust,soil dust/secondary dust and coal dust during the daytime in Longyan City,China.There are two major sources of pollution which are soil dust and mixture sources of automobile exhaust and secondary dust during the night in Longyan.The results of CMB show that the major sources are secondary dust,automobile exhaust and road dust during the daytime in Longyan.The major sources are secondary dust,soil dust and automobile exhaust during the night in Longyan.The results of the two methods are similar to each other and the results will guide us to plan to control the PM10 pollution sources in Longyan.展开更多
The validity of correlation analysis between finite element model（FEM） and modal test data is strongly affected by three factors, i.e., quality of excitation and measurement points in modal test,FEM reduction method...The validity of correlation analysis between finite element model（FEM） and modal test data is strongly affected by three factors, i.e., quality of excitation and measurement points in modal test,FEM reduction methods, and correlation check techniques. A new criterion based on modified mode participation（MMP） for choosing the best excitation point is presented. Comparison between this new criterion and mode participation（MP） criterion is made by using Case 1 with a simple printed circuit board（PCB）. The result indicates that this new criterion produces better results. In Case 2, 35 measurement points are selected to perform modal test and correlation analysis while 9 selected in Case 3.System equivalent reduction expansion process（SEREP）, modal assurance criteria（MAC）, coordinate modal assurance criteria（CoMAC）, pseudo orthogonality check（POC） and coordinate orthogonality check（CORTHOG） are used to show the error introduced by modal test in Cases 2 and 3. Case 2 shows that additional errors which cannot be identified by using CoMAC can be found by using CORTHOG.In both Cases 2 and 3, Guyan reduction, improved reduced system（IRS） method, SEREP and Hybrid reduction are compared for accuracy and robustness. The results suggest that the quality of the reduction process is problem dependent. However, the IRS method is an improvement over the Guyan reduction, and the Hybrid reduction is an improvement over the SEREP reduction.展开更多
This paper relates to the deep research on the Splinc Model Method of KED analysis. With the use of cubic B-splinc function as a link’s transverse deflection interpolation function, the principle of virtual displacem...This paper relates to the deep research on the Splinc Model Method of KED analysis. With the use of cubic B-splinc function as a link’s transverse deflection interpolation function, the principle of virtual displacement is presented as a basic theory for the general formulation of the equations of motion, and thus abandoned the kinematic assumption and the instantaneous structure assumption which arc used in the Spline Model Method. In thc same time, the nonlinear terms sue as coupling terms between thc rigid body motion and elastic deformation arc included. New member’s spline models are established. Mass matrix, Coriolis mass matrix, normal and tangential mass matrix, linear stiffness matrix, nonlinear stiffness matrix and rotation matrix arc derived. The kinematic differential equations of a member and system are deduced in the end. The Newmark direct integration method is used as the solution scheme of the kinematic differential equations to get the periodic response.展开更多
In this paper, we consider the giving up smoking model. First, we present the giving up smoking model in fractional order. Then the homotopy analysis method (HAM) is employed to compute an approximate and analytical s...In this paper, we consider the giving up smoking model. First, we present the giving up smoking model in fractional order. Then the homotopy analysis method (HAM) is employed to compute an approximate and analytical solution of the model in fractional order. The obtained results are compaired with those obtained by forth order Runge-Kutta method and nonstandard numerical method in the integer case. Finally, we present some numerical results.展开更多
Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challeng...Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challenge.This study investigates a bacterial meningitis model through deterministic and stochastic versions.Four-compartment population dynamics explain the concept,particularly the susceptible population,carrier,infected,and recovered.The model predicts the nonnegative equilibrium points and reproduction number,i.e.,the Meningitis-Free Equilibrium(MFE),and Meningitis-Existing Equilibrium(MEE).For the stochastic version of the existing deterministicmodel,the twomethodologies studied are transition probabilities and non-parametric perturbations.Also,positivity,boundedness,extinction,and disease persistence are studiedrigorouslywiththe helpofwell-known theorems.Standard and nonstandard techniques such as EulerMaruyama,stochastic Euler,stochastic Runge Kutta,and stochastic nonstandard finite difference in the sense of delay have been presented for computational analysis of the stochastic model.Unfortunately,standard methods fail to restore the biological properties of the model,so the stochastic nonstandard finite difference approximation is offered as an efficient,low-cost,and independent of time step size.In addition,the convergence,local,and global stability around the equilibria of the nonstandard computational method is studied by assuming the perturbation effect is zero.The simulations and comparison of the methods are presented to support the theoretical results and for the best visualization of results.展开更多
A two parameter mathematical model was developed to find the concentration for immobilized enzyme systems in porous spherical particles. This model contains a non-linear term related to reversible Michaelies-Menten ki...A two parameter mathematical model was developed to find the concentration for immobilized enzyme systems in porous spherical particles. This model contains a non-linear term related to reversible Michaelies-Menten kinetics. Analytical expression pertaining to the substrate concentration was reported for all possible values of Thiele module φ and α . In this work, we report the theoretically evaluated steady-state effectiveness factor for immobilized enzyme systems in porous spherical particles. These analytical results were found to be in good agreement with numerical results. Moreover, herein we employ new “Homotopy analysis method” (HAM) to solve non-linear reaction/diffusion equation.展开更多
Due to the size effects of rockfill materials, the settlement difference between numerical simulation and in situ monitoring of rockfill dams is a topic of general concern.The constitutive model parameters obtained fr...Due to the size effects of rockfill materials, the settlement difference between numerical simulation and in situ monitoring of rockfill dams is a topic of general concern.The constitutive model parameters obtained from laboratory triaxial tests often underestimate the deformation of high rockfill dams.Therefore, constitutive model parameters obtained by back analysis were used to calculate and predict the long-term deformation of rockfill dams.Instead of using artificial neural networks (ANNs), the response surface method (RSM) was employed to replace the finite element simulation used in the optimization iteration.Only 27 training samples were required for RSM, improving computational efficiency compared with ANN, which required 300 training samples.RSM can be used to describe the relationship between the constitutive model parameters and dam settlements.The inversion results of the Shuibuya concrete face rockfill dam (CFRD) show that the calculated settlements agree with the measured data, indicating the accuracy and efficiency of RSM.展开更多
Sensitivity analysis （SA） is an effective tool for studying crop models; it is an important link in model localization and plays an important role in crop model calibration and application. The objectives were to （...Sensitivity analysis （SA） is an effective tool for studying crop models; it is an important link in model localization and plays an important role in crop model calibration and application. The objectives were to （i） determine influential and non-influential parameters with respect to above ground biomass （AGB）, canopy cover （CC）, and grain yield of winter wheat in the Beijing area based on the AquaCrop model under different water treatments （rainfall, normal irrigation, and over-irrigation）; and （ii） generate an AquaCrop model that can be used in the Beijing area by setting non-influential parameters to fixed values and adjusting influential parameters according to the SA results. In this study, field experiments were conducted during the 2012-2013,2013-2014, and 2014-2015 winter wheat growing seasons at the National Precision Agriculture Demonstration Research Base in Beijing, China. The extended Fourier amplitude sensitivity test （EFAST） method was used to perform SA of the AquaCrop model using 42 crop parameters, in order to verify the SA results, data from the 2013-2014 growing season were used to calibrate the AquaCrop model, and data from 2012-2013 and 2014-2015 growing seasons were val- idated. For AGB and yield of winter wheat, the total order sensitivity analysis had more sensitive parameters than the first order sensitivity analysis. For the AGB time-series, parameter sensitivity was changed under different water treatments; in comparison with the non-stressful conditions （normal irrigation and over-irrigation）, there were more sensitive parameters under water stress （rainfall）, while root development parameters were more sensitive. For CC with time-series and yield, there were more sensitive parameters under water stress than under no water stress. Two parameters sets were selected to calibrate the AquaCrop model, one group of parameters were under water stress, and the others were under no water stress, there were two more sensitive parameters （growing degree-days （GDD） from sowing to the maximum rooting depth （root） and the maximum effective rooting depth （rtx）） under water stress than under no water stress. The results showed that there was higher accuracy under water stress than under no water stress. This study provides guidelines for AquaCrop model calibration and application in Beijing, China, as well providing guidance to simplify the AquaCrop model and improve its precision, especially when many parameters are used.展开更多
With the fast growth of Chinese economic, more and more capital will be invested in environmental projects. How to select the environmental investment projects （alternatives） for obtaining the best environmental qua...With the fast growth of Chinese economic, more and more capital will be invested in environmental projects. How to select the environmental investment projects （alternatives） for obtaining the best environmental quality and economic benefits is an important problem for the decision makers. The purpose of this paper is to develop a decision-making model to rank a finite number of alternatives with several and sometimes conflicting criteria. A model for ranking the projects of municipal sewage treatment plants is proposed by using exports＇ information and the data of the real projects. And, the ranking result is given based on the PROMETHEE method. Furthermore, by means of the concept of the weight stability intervals （WSI）, the sensitivity of the ranking results to the size of criteria values and the change of weights value of criteria are discussed. The result shows that some criteria, such as “proportion of benefit to project cost”, will influence the ranking result of alternatives very strong while others not. The influence are not only from the value of criterion but also from the changing the weight of criterion. So, some criteria such as “proportion of benefit to project cost” are key critera for ranking the projects. Decision makers must be cautious to them.展开更多
The investigations of physical attributes of oceans,including parameters such as heat flow and bathymetry,have garnered substantial attention and are particularly valuable for examining Earth’s thermal structures and...The investigations of physical attributes of oceans,including parameters such as heat flow and bathymetry,have garnered substantial attention and are particularly valuable for examining Earth’s thermal structures and dynamic processes.Nevertheless,classical plate cooling models exhibit disparities when predicting observed heat flow and seafloor depth for extremely young and old lithospheres.Furthermore,a comprehensive analysis of global heat flow predictions and regional ocean heat flow or bathymetry data with physical models has been lacking.In this study,we employed power-law models derived from the singularity theory of fractal density to meticulously fit the latest ocean heat flow and bathymetry.Notably,power-law models offer distinct advantages over traditional plate cooling models,showcasing robust self-similarity,scale invariance,or scaling properties,and providing a better fit to observed data.The outcomes of our singularity analysis concerning heat flow and bathymetry across diverse oceanic regions exhibit a degree of consistency with the global ocean spreading rate model.In addition,we applied the similarity method to predict a higher resolution(0.1°×0.1°)global heat flow map based on the most recent heat flow data and geological/geophysical observables refined through linear correlation analysis.Regions displaying significant disparities between predicted and observed heat flow are closely linked to hydrothermal vent fields and active structures.Finally,combining the actual bathymetry and predicted heat flow with the power-law models allows for the quantitative and comprehensive detection of anomalous regions of ocean subsidence and heat flow,which deviate from traditional plate cooling models.The anomalous regions of subsidence and heat flow show different degrees of anisotropy,providing new ideas and clues for further analysis of ocean topography or hydrothermal circulation of mid-ocean ridges.展开更多
3D dynamic analysis models of 1000 m deep-ocean mining pipeline, including steel lift pipe, pump, buffer and flexible hose, were established by finite element method (FEM). The coupling effect of steel lift pipe and f...3D dynamic analysis models of 1000 m deep-ocean mining pipeline, including steel lift pipe, pump, buffer and flexible hose, were established by finite element method (FEM). The coupling effect of steel lift pipe and flexible hose, and main external loads of pipeline were considered in the models, such as gravity, buoyancy, hydrodynamic forces, internal and external fluid pressures, concentrated suspension buoyancy on the flexible hose, torsional moment and axial force induced by pump working. Some relevant FEM models and solution techniques were developed, according to various 3D transient behaviors of integrated deep-ocean mining pipeline, including towing motions of track-keeping operation and launch process of pipeline. Meanwhile, an experimental verification system in towing water tank that had similar characteristics of designed mining pipeline was developed to verify the accuracy of the FEM models and dynamic simulation. The experiment results show that the experimental records and simulation results of stress of pipe are coincided. Based on the further simulations of 1 000 m deep-ocean mining pipeline, the simulation results show that, to form configuration of a saddle shape, the total concentrated suspension buoyancy of flexible hose should be 95%?105% of the gravity of flexible hose in water, the first suspension point occupies 1/3 of the total buoyancy, and the second suspension point occupies 2/3 of the total buoyancy. When towing velocity of mining system is less than 0.5 m/s, the towing track of buffer is coincided with the setting route of ship on the whole and the configuration of flexible hose is also kept well.展开更多
Considering that there are some limitations in analyzing the anti-sliding seismic stability of dam-foundation systems with the traditional pseudo-static method and response spectrum method, the dynamic strength reduct...Considering that there are some limitations in analyzing the anti-sliding seismic stability of dam-foundation systems with the traditional pseudo-static method and response spectrum method, the dynamic strength reduction method was used to study the deep anti-sliding stability of a high gravity dam with a complex dam foundation in response to strong earthquake-induced ground action. Based on static anti-sliding stability analysis of the dam foundation undertaken by decreasing the shear strength parameters of the rock mass in equal proportion, the seismic time history analysis was carried out. The proposed instability criterion for the dynamic strength reduction method was that the peak values of dynamic displacements and plastic strain energy change suddenly with the increase of the strength reduction factor. The elasto-plastic behavior of the dam foundation was idealized using the Drucker-Prager yield criterion based on the associated flow rule assumption. The result of elasto-plastic time history analysis of an overflow dam monolith based on the dynamic strength reduction method was compared with that of the dynamic linear elastic analysis, and the reliability of elasto-plastic time history analysis was confirmed. The results also show that the safety factors of the dam-foundation system in the static and dynamic cases are 3.25 and 3.0, respectively, and that the F2 fault has a significant influence on the anti-sliding stability of the high gravity dam. It is also concluded that the proposed instability criterion for the dynamic strength reduction method is feasible.展开更多
The Soil and Water Assessment Tool （SWAT） was implemented in a small forested watershed of the Soan River Basin innorthern Pakistan through application of the sequential uncertainty fitting （SUFI-2） method to inve...The Soil and Water Assessment Tool （SWAT） was implemented in a small forested watershed of the Soan River Basin innorthern Pakistan through application of the sequential uncertainty fitting （SUFI-2） method to investigate the associateduncertainty in runoff and sediment load estimation. The model was calibrated for a 10-year period （1991–2000） with aninitial 4-year warm-up period （1987–1990）, and was validated for the subsequent 10-year period （2001–2010）. Themodel evaluation indices R2 （the coefficient of determination）, NS （the Nash-Sutcliffe efficiency）, and PBIAS （percentbias） for stream flows simulation indicated that there was a good agreement between the measured and simulated flows.To assess the uncertainty in the model outputs, p-factor （a 95% prediction uncertainty, 95PPU） and r-factors （averagewideness width of the 95PPU band divided by the standard deviation of the observed values） were taken into account.The 95PPU band bracketed 72% of the observed data during the calibration and 67% during the validation. The r-factorwas 0.81 during the calibration and 0.68 during the validation. For monthly sediment yield, the model evaluation coefficients（R2 and NS） for the calibration were computed as 0.81 and 0.79, respectively; for validation, they were 0.78and 0.74, respectively. Meanwhile, the 95PPU covered more than 60% of the observed sediment data during calibrationand validation. Moreover, improved model prediction and parameter estimation were observed with the increasednumber of iterations. However, the model performance became worse after the fourth iterations due to an unreasonableparameter estimation. Overall results indicated the applicability of the SWAT model with moderate levels of uncertaintyduring the calibration and high levels during the validation. Thus, this calibrated SWAT model can be used for assessmentof water balance components, climate change studies, and land use management practices.展开更多
文摘On the basis of introducing the fundamental theory and the basic analysis steps of the sub model method, the strength of the new engine complex assembly structure was analyzed according to the properties of the engine structures, some of the key parts of the engine were analyzed with refined mesh by sub model method and the error of the FEM solution was estimated by the extrapolation method. The example showed that the sub model can not only analyze the comlex structures without the restriction of the software and hardware of the computers, but get the more precise analysis result also. This method is more suitable for the strength analysis of the complex assembly structure.
基金supported by the National Key Research and Development Program(MOST 2023YFA1606404 and MOST 2022YFA1602303)the National Natural Science Foundation of China(Nos.12347106,12147101,and 12447122)the China Postdoctoral Science Foundation(No.2024M760489).
文摘To study the uncertainty quantification of resonant states in open quantum systems,we developed a Bayesian framework by integrating a reduced basis method(RBM)emulator with the Gamow coupled-channel(GCC)approach.The RBM,constructed via eigenvector continuation and trained on both bound and resonant configurations,enables the fast and accurate emulation of resonance properties across the parameter space.To identify the physical resonant states from the emulator’s output,we introduce an overlap-based selection technique that effectively isolates true solutions from background artifacts.By applying this framework to unbound nucleus ^(6)Be,we quantified the model uncertainty in the predicted complex energies.The results demonstrate relative errors of 17.48%in the real part and 8.24%in the imaginary part,while achieving a speedup of four orders of magnitude compared with the full GCC calculations.To further investigate the asymptotic behavior of the resonant-state wavefunctions within the RBM framework,we employed a Lippmann–Schwinger(L–S)-based correction scheme.This approach not only improves the consistency between eigenvalues and wavefunctions but also enables a seamless extension from real-space training data to the complex energy plane.By bridging the gap between bound-state and continuum regimes,the L–S correction significantly enhances the emulator’s capability to accurately capture continuum structures in open quantum systems.
基金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.
基金National Natural Science Foundation of China under Grant No.U2139208。
文摘Natural soil generally exhibits significant transverse isotropy(TI)due to weathering and sedimentation,meaning that horizontal moduli differ from their vertical counterpart.The TI mechanical model is more appropriate for actual situations.Although soil exhibits material nonlinearity under earthquake excitation,existing research on the TI medium is limited to soil linearity and neglects the nonlinear response of TI sites.A 2D equivalent linear model for a layered TI half-space subjected to seismic waves is derived in the transformed wave number domain using the exact dynamic stiffness matrix of the TI medium.This study introduces a method for determining the effective shear strain of TI sites under oblique wave incidence,and further describes a systematic study on the effects of TI parameters and soil nonlinearity on site responses.Numerical results indicate that seismic responses of the TI medium significantly differ from those of isotropic sites and that the responses are highly dependent on TI parameters,particularly in nonlinear cases,while also being sensitive to incident angle and excitation intensity.Moreover,the differences in peak acceleration and waveform for various TI materials may also be amplified due to the strong nonlinearity.The study provides valuable insights for improving the accuracy of seismic response analysis in engineering applications.
基金The funding was provided by Laboratory Fund (Grant No.SYJJ200320).
文摘This paper aims to solve the resonance failure probability and develop an effective method to estimate the effects of variables and failure modes on failure probability of axially functionally graded material(FGM)pipe conveying fluid.Correspondingly,the natural frequency of axially FGM pipes conveying fluid is calculated using the differential quadrature method(DQM).A variable sensitivity analysis(VSA)is introduced to measure the effect of each random variable,and a mode sensitivity analysis(MSA)is introduced to acquire the importance ranking of failure modes.Then,an active learning Kriging(ALK)method is established to calculate the resonance failure probability and sensitivity indices,which greatly improves the application of resonance reliability analysis for pipelines in engineering practice.Based on the resonance reliability analysis method,the effects of fluid velocity,volume fraction and fluid density of axially FGM pipe conveying fluid on resonance reliability are analyzed.The results demonstrate that the proposed method has great performance in the anti-resonance analysis of pipes conveying fluid.
文摘In this paper,the homotopy analysis method (HAM) is applied to solve generalized biological populationmodels.The fractional derivatives are described by Caputo's sense.The method introduces a significant improvementin this field over existing techniques.Results obtained using the scheme presented here agree well with the analyticalsolutions and the numerical results presented in Ref.[6].However,the fundamental solutions of these equations stillexhibit useful scaling properties that make them attractive for applications.
基金Supported by the Natural Basic Research Program of China(No.2005CB422207)the Fund of Eco-enviromental Impacts and Protection in Devoloping and Utilizing of Oil-shale Resources(No.OSR-01-06)
文摘In order to identify the day and night pollution sources of PM10 in ambient air in Longyan City,the authors analyzed the elemental composition of respirable particulate matters in the day and night ambient air samples and various pollution sources which were collected in January 2010 in Longyan with inductivity coupled plasma-mass spectrometry（ICP-MS）.Then chemical mass balance（CMB） model and factor analysis（FA） method were applied to comparatively study the inorganic components in the sources and receptor samples.The results of factor analysis show that the major sources were road dust,waste incineration and mixed sources which contained automobile exhaust,soil dust/secondary dust and coal dust during the daytime in Longyan City,China.There are two major sources of pollution which are soil dust and mixture sources of automobile exhaust and secondary dust during the night in Longyan.The results of CMB show that the major sources are secondary dust,automobile exhaust and road dust during the daytime in Longyan.The major sources are secondary dust,soil dust and automobile exhaust during the night in Longyan.The results of the two methods are similar to each other and the results will guide us to plan to control the PM10 pollution sources in Longyan.
基金supported by Science and Technology on Reliability and Environmental Engineering Laboratory,Beihang University
文摘The validity of correlation analysis between finite element model（FEM） and modal test data is strongly affected by three factors, i.e., quality of excitation and measurement points in modal test,FEM reduction methods, and correlation check techniques. A new criterion based on modified mode participation（MMP） for choosing the best excitation point is presented. Comparison between this new criterion and mode participation（MP） criterion is made by using Case 1 with a simple printed circuit board（PCB）. The result indicates that this new criterion produces better results. In Case 2, 35 measurement points are selected to perform modal test and correlation analysis while 9 selected in Case 3.System equivalent reduction expansion process（SEREP）, modal assurance criteria（MAC）, coordinate modal assurance criteria（CoMAC）, pseudo orthogonality check（POC） and coordinate orthogonality check（CORTHOG） are used to show the error introduced by modal test in Cases 2 and 3. Case 2 shows that additional errors which cannot be identified by using CoMAC can be found by using CORTHOG.In both Cases 2 and 3, Guyan reduction, improved reduced system（IRS） method, SEREP and Hybrid reduction are compared for accuracy and robustness. The results suggest that the quality of the reduction process is problem dependent. However, the IRS method is an improvement over the Guyan reduction, and the Hybrid reduction is an improvement over the SEREP reduction.
文摘This paper relates to the deep research on the Splinc Model Method of KED analysis. With the use of cubic B-splinc function as a link’s transverse deflection interpolation function, the principle of virtual displacement is presented as a basic theory for the general formulation of the equations of motion, and thus abandoned the kinematic assumption and the instantaneous structure assumption which arc used in the Spline Model Method. In thc same time, the nonlinear terms sue as coupling terms between thc rigid body motion and elastic deformation arc included. New member’s spline models are established. Mass matrix, Coriolis mass matrix, normal and tangential mass matrix, linear stiffness matrix, nonlinear stiffness matrix and rotation matrix arc derived. The kinematic differential equations of a member and system are deduced in the end. The Newmark direct integration method is used as the solution scheme of the kinematic differential equations to get the periodic response.
文摘In this paper, we consider the giving up smoking model. First, we present the giving up smoking model in fractional order. Then the homotopy analysis method (HAM) is employed to compute an approximate and analytical solution of the model in fractional order. The obtained results are compaired with those obtained by forth order Runge-Kutta method and nonstandard numerical method in the integer case. Finally, we present some numerical results.
基金Deanship of Research and Graduate Studies at King Khalid University for funding this work through large Research Project under Grant Number RGP2/302/45supported by the Deanship of Scientific Research,Vice Presidency forGraduate Studies and Scientific Research,King Faisal University,Saudi Arabia(Grant Number A426).
文摘Based on theWorld Health Organization(WHO),Meningitis is a severe infection of the meninges,the membranes covering the brain and spinal cord.It is a devastating disease and remains a significant public health challenge.This study investigates a bacterial meningitis model through deterministic and stochastic versions.Four-compartment population dynamics explain the concept,particularly the susceptible population,carrier,infected,and recovered.The model predicts the nonnegative equilibrium points and reproduction number,i.e.,the Meningitis-Free Equilibrium(MFE),and Meningitis-Existing Equilibrium(MEE).For the stochastic version of the existing deterministicmodel,the twomethodologies studied are transition probabilities and non-parametric perturbations.Also,positivity,boundedness,extinction,and disease persistence are studiedrigorouslywiththe helpofwell-known theorems.Standard and nonstandard techniques such as EulerMaruyama,stochastic Euler,stochastic Runge Kutta,and stochastic nonstandard finite difference in the sense of delay have been presented for computational analysis of the stochastic model.Unfortunately,standard methods fail to restore the biological properties of the model,so the stochastic nonstandard finite difference approximation is offered as an efficient,low-cost,and independent of time step size.In addition,the convergence,local,and global stability around the equilibria of the nonstandard computational method is studied by assuming the perturbation effect is zero.The simulations and comparison of the methods are presented to support the theoretical results and for the best visualization of results.
文摘A two parameter mathematical model was developed to find the concentration for immobilized enzyme systems in porous spherical particles. This model contains a non-linear term related to reversible Michaelies-Menten kinetics. Analytical expression pertaining to the substrate concentration was reported for all possible values of Thiele module φ and α . In this work, we report the theoretically evaluated steady-state effectiveness factor for immobilized enzyme systems in porous spherical particles. These analytical results were found to be in good agreement with numerical results. Moreover, herein we employ new “Homotopy analysis method” (HAM) to solve non-linear reaction/diffusion equation.
基金supported by the National Natural Science Foundation of China(Grant No.51579193)the Science and Technology Planning Project of Guizhou Province(Grant No.[2016]1154)
文摘Due to the size effects of rockfill materials, the settlement difference between numerical simulation and in situ monitoring of rockfill dams is a topic of general concern.The constitutive model parameters obtained from laboratory triaxial tests often underestimate the deformation of high rockfill dams.Therefore, constitutive model parameters obtained by back analysis were used to calculate and predict the long-term deformation of rockfill dams.Instead of using artificial neural networks (ANNs), the response surface method (RSM) was employed to replace the finite element simulation used in the optimization iteration.Only 27 training samples were required for RSM, improving computational efficiency compared with ANN, which required 300 training samples.RSM can be used to describe the relationship between the constitutive model parameters and dam settlements.The inversion results of the Shuibuya concrete face rockfill dam (CFRD) show that the calculated settlements agree with the measured data, indicating the accuracy and efficiency of RSM.
基金supported by the National Natural Science Foundation of China(41571416)the Natural Science Foundation of Beijing,China(4152019)the Beijing Academy of Agricultural and Forestry Sciences Innovation Capacity Construction Specific Projects,China(KJCX20150409)
文摘Sensitivity analysis （SA） is an effective tool for studying crop models; it is an important link in model localization and plays an important role in crop model calibration and application. The objectives were to （i） determine influential and non-influential parameters with respect to above ground biomass （AGB）, canopy cover （CC）, and grain yield of winter wheat in the Beijing area based on the AquaCrop model under different water treatments （rainfall, normal irrigation, and over-irrigation）; and （ii） generate an AquaCrop model that can be used in the Beijing area by setting non-influential parameters to fixed values and adjusting influential parameters according to the SA results. In this study, field experiments were conducted during the 2012-2013,2013-2014, and 2014-2015 winter wheat growing seasons at the National Precision Agriculture Demonstration Research Base in Beijing, China. The extended Fourier amplitude sensitivity test （EFAST） method was used to perform SA of the AquaCrop model using 42 crop parameters, in order to verify the SA results, data from the 2013-2014 growing season were used to calibrate the AquaCrop model, and data from 2012-2013 and 2014-2015 growing seasons were val- idated. For AGB and yield of winter wheat, the total order sensitivity analysis had more sensitive parameters than the first order sensitivity analysis. For the AGB time-series, parameter sensitivity was changed under different water treatments; in comparison with the non-stressful conditions （normal irrigation and over-irrigation）, there were more sensitive parameters under water stress （rainfall）, while root development parameters were more sensitive. For CC with time-series and yield, there were more sensitive parameters under water stress than under no water stress. Two parameters sets were selected to calibrate the AquaCrop model, one group of parameters were under water stress, and the others were under no water stress, there were two more sensitive parameters （growing degree-days （GDD） from sowing to the maximum rooting depth （root） and the maximum effective rooting depth （rtx）） under water stress than under no water stress. The results showed that there was higher accuracy under water stress than under no water stress. This study provides guidelines for AquaCrop model calibration and application in Beijing, China, as well providing guidance to simplify the AquaCrop model and improve its precision, especially when many parameters are used.
基金Shanghai Leading Academic Discipline Project (T0502)Shanghai Municipal Educational Commission Project (05EZ32).
文摘With the fast growth of Chinese economic, more and more capital will be invested in environmental projects. How to select the environmental investment projects （alternatives） for obtaining the best environmental quality and economic benefits is an important problem for the decision makers. The purpose of this paper is to develop a decision-making model to rank a finite number of alternatives with several and sometimes conflicting criteria. A model for ranking the projects of municipal sewage treatment plants is proposed by using exports＇ information and the data of the real projects. And, the ranking result is given based on the PROMETHEE method. Furthermore, by means of the concept of the weight stability intervals （WSI）, the sensitivity of the ranking results to the size of criteria values and the change of weights value of criteria are discussed. The result shows that some criteria, such as “proportion of benefit to project cost”, will influence the ranking result of alternatives very strong while others not. The influence are not only from the value of criterion but also from the changing the weight of criterion. So, some criteria such as “proportion of benefit to project cost” are key critera for ranking the projects. Decision makers must be cautious to them.
基金supported by the Guangdong Province Introduced Innovative R&D Team of Big Data-Mathematical Earth Sciences and Extreme Geological Events Team(grant number 2021ZT09H399)the National Natural Science Foundation of China(grant number 42430111,42050103).
文摘The investigations of physical attributes of oceans,including parameters such as heat flow and bathymetry,have garnered substantial attention and are particularly valuable for examining Earth’s thermal structures and dynamic processes.Nevertheless,classical plate cooling models exhibit disparities when predicting observed heat flow and seafloor depth for extremely young and old lithospheres.Furthermore,a comprehensive analysis of global heat flow predictions and regional ocean heat flow or bathymetry data with physical models has been lacking.In this study,we employed power-law models derived from the singularity theory of fractal density to meticulously fit the latest ocean heat flow and bathymetry.Notably,power-law models offer distinct advantages over traditional plate cooling models,showcasing robust self-similarity,scale invariance,or scaling properties,and providing a better fit to observed data.The outcomes of our singularity analysis concerning heat flow and bathymetry across diverse oceanic regions exhibit a degree of consistency with the global ocean spreading rate model.In addition,we applied the similarity method to predict a higher resolution(0.1°×0.1°)global heat flow map based on the most recent heat flow data and geological/geophysical observables refined through linear correlation analysis.Regions displaying significant disparities between predicted and observed heat flow are closely linked to hydrothermal vent fields and active structures.Finally,combining the actual bathymetry and predicted heat flow with the power-law models allows for the quantitative and comprehensive detection of anomalous regions of ocean subsidence and heat flow,which deviate from traditional plate cooling models.The anomalous regions of subsidence and heat flow show different degrees of anisotropy,providing new ideas and clues for further analysis of ocean topography or hydrothermal circulation of mid-ocean ridges.
基金Project(DY105-3-2-2) supported by China Ocean Mineral Resources Research and Development Association(COMRA)Project(50675226) supported by the National Natural Science Foundation of China
文摘3D dynamic analysis models of 1000 m deep-ocean mining pipeline, including steel lift pipe, pump, buffer and flexible hose, were established by finite element method (FEM). The coupling effect of steel lift pipe and flexible hose, and main external loads of pipeline were considered in the models, such as gravity, buoyancy, hydrodynamic forces, internal and external fluid pressures, concentrated suspension buoyancy on the flexible hose, torsional moment and axial force induced by pump working. Some relevant FEM models and solution techniques were developed, according to various 3D transient behaviors of integrated deep-ocean mining pipeline, including towing motions of track-keeping operation and launch process of pipeline. Meanwhile, an experimental verification system in towing water tank that had similar characteristics of designed mining pipeline was developed to verify the accuracy of the FEM models and dynamic simulation. The experiment results show that the experimental records and simulation results of stress of pipe are coincided. Based on the further simulations of 1 000 m deep-ocean mining pipeline, the simulation results show that, to form configuration of a saddle shape, the total concentrated suspension buoyancy of flexible hose should be 95%?105% of the gravity of flexible hose in water, the first suspension point occupies 1/3 of the total buoyancy, and the second suspension point occupies 2/3 of the total buoyancy. When towing velocity of mining system is less than 0.5 m/s, the towing track of buffer is coincided with the setting route of ship on the whole and the configuration of flexible hose is also kept well.
基金supported by the National Basic Research Program of China (973 Program,Grant No.2007CB714104)the National Natural Science Foundation of China (Grant No. 50779011)the Innovative Project for Graduate Students of Jiangsu Province (Grant No. CX09B_155Z)
文摘Considering that there are some limitations in analyzing the anti-sliding seismic stability of dam-foundation systems with the traditional pseudo-static method and response spectrum method, the dynamic strength reduction method was used to study the deep anti-sliding stability of a high gravity dam with a complex dam foundation in response to strong earthquake-induced ground action. Based on static anti-sliding stability analysis of the dam foundation undertaken by decreasing the shear strength parameters of the rock mass in equal proportion, the seismic time history analysis was carried out. The proposed instability criterion for the dynamic strength reduction method was that the peak values of dynamic displacements and plastic strain energy change suddenly with the increase of the strength reduction factor. The elasto-plastic behavior of the dam foundation was idealized using the Drucker-Prager yield criterion based on the associated flow rule assumption. The result of elasto-plastic time history analysis of an overflow dam monolith based on the dynamic strength reduction method was compared with that of the dynamic linear elastic analysis, and the reliability of elasto-plastic time history analysis was confirmed. The results also show that the safety factors of the dam-foundation system in the static and dynamic cases are 3.25 and 3.0, respectively, and that the F2 fault has a significant influence on the anti-sliding stability of the high gravity dam. It is also concluded that the proposed instability criterion for the dynamic strength reduction method is feasible.
基金supported by the Centre of Excellence in Water Resources Engineering, University of Engineering and Technology Lahore, and local authorities in Pakistan
文摘The Soil and Water Assessment Tool （SWAT） was implemented in a small forested watershed of the Soan River Basin innorthern Pakistan through application of the sequential uncertainty fitting （SUFI-2） method to investigate the associateduncertainty in runoff and sediment load estimation. The model was calibrated for a 10-year period （1991–2000） with aninitial 4-year warm-up period （1987–1990）, and was validated for the subsequent 10-year period （2001–2010）. Themodel evaluation indices R2 （the coefficient of determination）, NS （the Nash-Sutcliffe efficiency）, and PBIAS （percentbias） for stream flows simulation indicated that there was a good agreement between the measured and simulated flows.To assess the uncertainty in the model outputs, p-factor （a 95% prediction uncertainty, 95PPU） and r-factors （averagewideness width of the 95PPU band divided by the standard deviation of the observed values） were taken into account.The 95PPU band bracketed 72% of the observed data during the calibration and 67% during the validation. The r-factorwas 0.81 during the calibration and 0.68 during the validation. For monthly sediment yield, the model evaluation coefficients（R2 and NS） for the calibration were computed as 0.81 and 0.79, respectively; for validation, they were 0.78and 0.74, respectively. Meanwhile, the 95PPU covered more than 60% of the observed sediment data during calibrationand validation. Moreover, improved model prediction and parameter estimation were observed with the increasednumber of iterations. However, the model performance became worse after the fourth iterations due to an unreasonableparameter estimation. Overall results indicated the applicability of the SWAT model with moderate levels of uncertaintyduring the calibration and high levels during the validation. Thus, this calibrated SWAT model can be used for assessmentof water balance components, climate change studies, and land use management practices.
