A cascading failure of landslide dams caused by strong earthquakes or torrential rains in mountainous river valleys can pose great threats to people’s lives,properties,and infrastructures.In this study,based on the t...A cascading failure of landslide dams caused by strong earthquakes or torrential rains in mountainous river valleys can pose great threats to people’s lives,properties,and infrastructures.In this study,based on the three-dimensional Reynoldsaveraged Navier-Stokes equations(RANS),the renormalization group(RNG)k-εturbulence model,suspended and bed load transport equations,and the instability discriminant formula of dam breach side slope,and the explicit finite volume method(FVM),a detailed numerical simulation model for calculating the hydro-morphodynamic characteristics of cascading dam breach process has been developed.The developed numerical model can simulate the breach hydrograph and the dam breach morphology evolution during the cascading failure process of landslide dams.A model test of the breaches of two cascading landslide dams has been used as the validation case.The comparison of the calculated and measured results indicates that the breach hydrograph and the breach morphology evolution process of the upstream and downstream dams are generally consistent with each other,and the relative errors of the key breaching parameters,i.e.,the peak breach flow and the time to peak of each dam,are less than±5%.Further,the comparison of the breach hydrographs of the upstream and downstream dams shows that there is an amplification effect of the breach flood on the cascading landslide dam failures.Three key parameters,i.e.,the distance between the upstream and the downstream dams,the river channel slope,and the downstream dam height,have been used to study the flood amplification effect.The parameter sensitivity analyses show that the peak breach flow at the downstream dam decreases with increasing distance between the upstream and the downstream dams,and the downstream dam height.Further,the peak breach flow at the downstream dam first increases and then decreases with steepening of the river channel slope.When the flood caused by the upstream dam failure flows to the downstream dam,it can produce a surge wave that overtops and erodes the dam crest,resulting in a lowering of the dam crest elevation.This has an impact on the failure occurrence time and the peak breach flow of the downstream dam.The influence of the surge wave on the downstream dam failure process is related to the volume of water that overtops the dam crest and the erosion characteristics of dam material.Moreover,the cascading failure case of the Xiaogangjian and Lower Xiaogangjian landslide dams has also been used as the representative case for validating the model.In comparisons of the calculated and measured breach hydrographs and final breach morphologies,the relative errors of the key dam breaching parameters are all within±10%,which verify the rationality of the model is applicable to real-world cases.Overall,the numerical model developed in this study can provide important technical support for the risk assessment and emergency treatment of failures of cascading landslide dams.展开更多
This work highlights the unparalleled efficiency of the “n<sup>th</sup>-Order Function/ Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (n<sup>th</sup>-FASAM-N) by con...This work highlights the unparalleled efficiency of the “n<sup>th</sup>-Order Function/ Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (n<sup>th</sup>-FASAM-N) by considering the well-known Nordheim-Fuchs reactor dynamics/safety model. This model describes a short-time self-limiting power excursion in a nuclear reactor system having a negative temperature coefficient in which a large amount of reactivity is suddenly inserted, either intentionally or by accident. This nonlinear paradigm model is sufficiently complex to model realistically self-limiting power excursions for short times yet admits closed-form exact expressions for the time-dependent neutron flux, temperature distribution and energy released during the transient power burst. The n<sup>th</sup>-FASAM-N methodology is compared to the extant “n<sup>th</sup>-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (n<sup>th</sup>-CASAM-N) showing that: (i) the 1<sup>st</sup>-FASAM-N and the 1<sup>st</sup>-CASAM-N methodologies are equally efficient for computing the first-order sensitivities;each methodology requires a single large-scale computation for solving the “First-Level Adjoint Sensitivity System” (1<sup>st</sup>-LASS);(ii) the 2<sup>nd</sup>-FASAM-N methodology is considerably more efficient than the 2<sup>nd</sup>-CASAM-N methodology for computing the second-order sensitivities since the number of feature-functions is much smaller than the number of primary parameters;specifically for the Nordheim-Fuchs model, the 2<sup>nd</sup>-FASAM-N methodology requires 2 large-scale computations to obtain all of the exact expressions of the 28 distinct second-order response sensitivities with respect to the model parameters while the 2<sup>nd</sup>-CASAM-N methodology requires 7 large-scale computations for obtaining these 28 second-order sensitivities;(iii) the 3<sup>rd</sup>-FASAM-N methodology is even more efficient than the 3<sup>rd</sup>-CASAM-N methodology: only 2 large-scale computations are needed to obtain the exact expressions of the 84 distinct third-order response sensitivities with respect to the Nordheim-Fuchs model’s parameters when applying the 3<sup>rd</sup>-FASAM-N methodology, while the application of the 3<sup>rd</sup>-CASAM-N methodology requires at least 22 large-scale computations for computing the same 84 distinct third-order sensitivities. Together, the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N methodologies are the most practical methodologies for computing response sensitivities of any order comprehensively and accurately, overcoming the curse of dimensionality in sensitivity analysis.展开更多
This work presents the “n<sup>th</sup>-Order Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (abbreviated as “n<sup>th</sup>-FASAM-N”), which will be shown to be the...This work presents the “n<sup>th</sup>-Order Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (abbreviated as “n<sup>th</sup>-FASAM-N”), which will be shown to be the most efficient methodology for computing exact expressions of sensitivities, of any order, of model responses with respect to features of model parameters and, subsequently, with respect to the model’s uncertain parameters, boundaries, and internal interfaces. The unparalleled efficiency and accuracy of the n<sup>th</sup>-FASAM-N methodology stems from the maximal reduction of the number of adjoint computations (which are considered to be “large-scale” computations) for computing high-order sensitivities. When applying the n<sup>th</sup>-FASAM-N methodology to compute the second- and higher-order sensitivities, the number of large-scale computations is proportional to the number of “model features” as opposed to being proportional to the number of model parameters (which are considerably more than the number of features).When a model has no “feature” functions of parameters, but only comprises primary parameters, the n<sup>th</sup>-FASAM-N methodology becomes identical to the extant n<sup>th</sup> CASAM-N (“n<sup>th</sup>-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems”) methodology. Both the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N methodologies are formulated in linearly increasing higher-dimensional Hilbert spaces as opposed to exponentially increasing parameter-dimensional spaces thus overcoming the curse of dimensionality in sensitivity analysis of nonlinear systems. Both the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N are incomparably more efficient and more accurate than any other methods (statistical, finite differences, etc.) for computing exact expressions of response sensitivities of any order with respect to the model’s features and/or primary uncertain parameters, boundaries, and internal interfaces.展开更多
In order to understand the effect of hardening ductility parameters and softening ductility parameters of the concrete damage plastic model in LS-DYNA,a sensitivity and reliability analysis of these parameters through...In order to understand the effect of hardening ductility parameters and softening ductility parameters of the concrete damage plastic model in LS-DYNA,a sensitivity and reliability analysis of these parameters through a convenient cube unit test was conducted. The results showed that the peak strength strain was independent of the hardening ductility parameter DH,but affected by AH,BH,and CH. The softening ductility was mainly related to the softening ductility parameter AS,but not affected by the damage ductility exponent BS. In case that the model with default parameters failed to match the AS-controlled damage softening phase,an optimized model with an AS correction was developed. The corrected model with the AS value of 2 matched well with the code model,and exhibited good feasibility in predicting the stress-strain curve of different grades of concrete. Moreover,the practicability of the corrected model was further validated by the conventional triaxial test. The simulated curve exhibited favorable consistence with the trial curve. Therefore,the model with parameter correction could provide a prospective reference for predicting the mechanical properties of concrete.展开更多
Input torque is the main power to maintain bipedal walking of robot, and can be calculated from trajectory planning and dynamic modeling on biped robot. During bipedal walking, the input torque is usually required to ...Input torque is the main power to maintain bipedal walking of robot, and can be calculated from trajectory planning and dynamic modeling on biped robot. During bipedal walking, the input torque is usually required to be adjusted due to some uncertain parameters arising from objective or subjective factors in the dynamical model to maintain the pre-planned stable trajectory. Here, a planar 5-link biped robot is used as an illustrating example to investigate the effects of uncertain parameters on the input torques. Kine-matic equations of the biped robot are firstly established by the third-order spline curves based on the trajectory planning method, and the dynamic modeling is accomplished by taking both the certain and uncertain parameters into account. Next, several evaluation indices on input torques are intro-duced to perform sensitivity analysis of the input torque with respect to the uncertain parameters. Finally, based on the Monte Carlo simulation, the values of evaluation indices on input torques are presented, from which all the robot param-eters are classified into three categories, i.e., strongly sensi-tive, sensitive and almost insensitive parameters.展开更多
This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>...This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.展开更多
Population viability analysis(PVA)is a tool to evaluate the risk of extinction for endangered species and aid conservation decision-making.The quality of PVA output is dependent on parameters related to population dyn...Population viability analysis(PVA)is a tool to evaluate the risk of extinction for endangered species and aid conservation decision-making.The quality of PVA output is dependent on parameters related to population dynamics and life-history;however,it has been difficult to collect this information for the giant panda(Aliuropoda melanoleuca),a rare and endangered mammal native to China,confined to some 30 fragmented habitat patches.Since giant pandas are long-lived,mature late,have lower reproductive rates,and show little sexual dimorphism,obtaining data to perform adequate PVA has been difficult.Here,we develop a parameter sensitivity index by modeling the dynamics of six giant panda populations in the Minshan Mountains,in order to determine the parameters most influential to giant panda populations.Our data shows that the giant panda populations are most sensitive to changes in four female parameters:initial breeding age,reproductive rate,mortality rate between age 0 and 1,and mortality rate of adults.The parameter sensitivity index strongly correlated with initial population size,as smaller populations were more sensitive to changes in these four variables.This model suggests that demographic parameters of females have more influence on the results of PVA,indicating that females may play a more important role in giant panda population dynamics than males.Consequently,reintroduction of female individuals to a small giant panda population should be a high priority for conservation efforts.Our findings form a technical basis for the coming program of giant panda reintroduction,and inform which parameters are crucial to successfully and feasibly monitoring wild giant panda populations.展开更多
In this paper, the axial-flux permanent magnet driver is modeledand analyzed in a simple and novel way under three-dimensional cylindricalcoordinates. The inherent three-dimensional characteristics of the deviceare co...In this paper, the axial-flux permanent magnet driver is modeledand analyzed in a simple and novel way under three-dimensional cylindricalcoordinates. The inherent three-dimensional characteristics of the deviceare comprehensively considered, and the governing equations are solved bysimplifying the boundary conditions. The axial magnetization of the sectorshapedpermanent magnets is accurately described in an algebraic form bythe parameters, which makes the physical meaning more explicit than thepurely mathematical expression in general series forms. The parameters of theBessel function are determined simply and the magnetic field distribution ofpermanent magnets and the air-gap is solved. Furthermore, the field solutionsare completely analytical, which provides convenience and satisfactoryaccuracy for modeling a series of electromagnetic performance parameters,such as the axial electromagnetic force density, axial electromagnetic force,and electromagnetic torque. The correctness and accuracy of the analyticalmodels are fully verified by three-dimensional finite element simulations and a15 kW prototype and the results of calculations, simulations, and experimentsunder three methods are highly consistent. The influence of several designparameters on magnetic field distribution and performance is studied and discussed.The results indicate that the modeling method proposed in this papercan calculate the magnetic field distribution and performance accurately andrapidly, which affords an important reference for the design and optimizationof axial-flux permanent magnet drivers.展开更多
A challenge for the development of Land Surface Models(LSMs) is improving transpiration of water exchange and photosynthesis of carbon exchange between terrestrial plants and the atmosphere, both of which are governed...A challenge for the development of Land Surface Models(LSMs) is improving transpiration of water exchange and photosynthesis of carbon exchange between terrestrial plants and the atmosphere, both of which are governed by stoma in leaves. In the photosynthesis module of these LSMs, variations of parameters arising from diversity in plant functional types(PFTs) and climate remain unclear. Identifying sensitive parameters among all photosynthetic parameters before parameter estimation can not only reduce operation cost, but also improve the usability of photosynthesis models worldwide. Here, we analyzed 13 parameters of a biochemically-based photosynthesis model(FvCB), implemented in many LSMs, using two sensitivity analysis(SA) methods(i.e., the Sobol’ method and the Morris method) for setting up the parameter ensemble. Three different model performance metrics, i.e.,Root Mean Squared Error(RMSE), Nash Sutcliffe efficiency(NSE), and Standard Deviation(STDEV) were introduced for model assessment and sensitive parameters identification. The results showed that among all photosynthetic parameters only a small portion of parameters were sensitive, and the sensitive parameters were different across plant functional types: maximum rate of Rubisco activity(Vcmax25), maximum electron transport rate(Jmax25), triose phosphate use rate(TPU) and dark respiration in light(Rd) were sensitive in broad leafevergreen trees(BET), broad leaf-deciduous trees(BDT) and needle leaf-evergreen trees(NET), while only Vcmax25and TPU are sensitive in short vegetation(SV), dwarf trees and shrubs(DTS), and agriculture and grassland(AG). The two sensitivity analysis methods suggested a strong SA coherence;in contrast, different model performance metrics led to different SA results. This misfit suggests that more accurate values of sensitive parameters, specifically, species specific and seasonal variable parameters, are required to improve the performance of the FvCB model.展开更多
To enhance the output torque and minimize the torque ripple of coaxial magnetic gear(CMG),a novel auxiliary flux modulator CMG with unequal magnetic poles is proposed.This design incorporates an inner rotor with an as...To enhance the output torque and minimize the torque ripple of coaxial magnetic gear(CMG),a novel auxiliary flux modulator CMG with unequal magnetic poles is proposed.This design incorporates an inner rotor with an asymmetric sector and a trapezoidal combined N-S pole structure,featuring Halbach arrays for the arrangement of permanent magnets(PMs).The outer rotor PMs adopt a Spoke-type configuration.To optimize the CMG for high output torque and low torque ripple,a sensitivity analysis is conducted to identify key size parameters that significantly influence the optimization objectives.Based on the sensitivity hierarchy of these parameters,a multi-objective optimization analysis is performed using a genetic algorithm(GA)to determine the optimal structural parameter values of the CMG.In addition,a coaxial magnetic gear(CMG)topology with 4 inner and 17 outer pole pairs is adopted,and the parametric model is established.Finally,the electromagnetic properties of the CMG are evaluated using the finite element method.The results indicate a remarkable reduction in torque ripple,specifically by 46.15%.展开更多
During bipedal walking,it is critical to detect and adjust the robot postures by feedback control to maintain its normal state amidst multi-source random disturbances arising from some unavoidable uncertain factors.Th...During bipedal walking,it is critical to detect and adjust the robot postures by feedback control to maintain its normal state amidst multi-source random disturbances arising from some unavoidable uncertain factors.The radical basis function(RBF)neural network model of a five-link biped robot is established,and two certain disturbances and a randomly uncertain disturbance are then mixed with the optimal torques in the network model to study the performance of the biped robot by several evaluation indices and a specific Poincar′e map.In contrast with the simulations,the response varies as desired under optimal inputting while the output is fluctuating in the situation of disturbance driving.Simulation results from noise inputting also show that the dynamics of the robot is less sensitive to the disturbance of knee joint input of the swing leg than those of the other three joints,the response errors of the biped will be increasing with higher disturbance levels,and especially there are larger output fluctuations in the knee and hip joints of the swing leg.展开更多
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.展开更多
To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale co...To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale content,porosity,fluid type and saturation.Then with analyzing the change rules of elastic parameters along with the above influence factors and the cross-plots among elastic parameters,the sensitive elastic parameters of tight sandstone reservoir are determined,and the rock physics template of sweet spot is constructed to guide pre-stack seismic inversion.The results show that velocity ratio and Poisson impedance are the most sensitive elastic parameters to indicate the lithologic and gas-bearing properties of sweet spot in tight sandstone reservoir.The high-quality sweet spot is characterized by the lower velocity ratio and Poisson impedance.Finally,the actual seismic data are selected to predict the sweet spots in tight sandstone gas reservoirs,so as to verify the validity of the rock physical simulation results.The significant consistency between the relative logging curves and inversion results in different wells implies that the utilization of well-bore rock physical simulation can guide the prediction of sweet spot in tight sandstone gas reservoirs.展开更多
The randomness of rock joint development is an important factor in the uncertainty of geotechnical engineering stability.In this study,a method is proposed to evaluate the reliability of intermittent jointed rock slop...The randomness of rock joint development is an important factor in the uncertainty of geotechnical engineering stability.In this study,a method is proposed to evaluate the reliability of intermittent jointed rock slope.The least squares support vector machine(LSSVM)evolved by a bacterial foraging optimization algorithm(BFOA)is used to establish a response surface model to express the mapping relationship between the intermittent joint parameters and the slope safety factor.The training samples are obtained from the numerical calculation based on the joint finite element method during this process.Considering the randomness of the intermittent joint parameters in the actual project,each parameter is evaluated at different locations on the site,and its distribution characteristics are counted.According to these statistical results,a large number of parameter combinations are obtained through Monte Carlo sampling.The trained machine learning mapping model is used to obtain the slope safety factor corresponding to each group,and these results are then used to obtain the slope reliability.When the research results were applied to slope disaster treatment along the Yalu River in China’s Jilin Province,it was found that the joint length and joint inclination angle both play key roles in rock slope stability,which should receive more attention in the slope treatment.In summary,this study establishes a method for evaluating the reliability of intermittent jointed rock slope based on an evolutionary SVM model,and its feasibility is verified by engineering application.展开更多
Stress Joint (SJ) plays a key role in the Top Tensioned Riser (TTR) system for deep water engineering. A preliminary design method of tapered SJ is proposed in the paper, which could help designers obtain accurate...Stress Joint (SJ) plays a key role in the Top Tensioned Riser (TTR) system for deep water engineering. A preliminary design method of tapered SJ is proposed in the paper, which could help designers obtain accurate design data. After a further sensitive analysis is carried out, the related parameters choice and control methods are recommended in the engineering practice. By taking the extreme environment conditions into consideration, the effects of bending stress reduction and curve control are analyzed, and the 3-D FE models are established by ABQOUS for numerical evaluation to verify the correctness of design results. At last, dynamic analysis and fatigue analysis, based on actual project, are carried out with designed stress joint. The analysis results prove the feasibility and guidance of this method in the practical engineering applications.展开更多
In order to develop a coupled basin scale model of ocean circulation and biogeochemical cycling,we present a biogeochemical model including 12 components to study the ecosystem in the China coastal seas(CCS).The for...In order to develop a coupled basin scale model of ocean circulation and biogeochemical cycling,we present a biogeochemical model including 12 components to study the ecosystem in the China coastal seas(CCS).The formulation of phytoplankton mortality and zooplankton growth are modified according to biological characteristics of CCS.The four sensitivity biological parameters,zooplankton assimilation efficiency rate(ZooAE_N),zooplankton basal metabolism rate(ZooBM),maximum specific growth rate of zooplankton(μ_(20)) and maximum chlorophyll to carbon ratio(Chl2C_m) are obtained in sensitivity experiments for the phytoplankton,and experiments about the parameter μ_(20'),half-saturation for phytoplankton NO_3 uptake(K_(NO_3)) and remineralization rate of small detritusN(SDeRRN) are conducted.The results demonstrate that the biogeochemical model is quite sensitive to the zooplankton grazing parameter when it ranges from 0.1 to 1.2 d^(-1).The K_(NO_3) and SDeRRN also play an important role in determining the nitrogen cycle within certain ranges.The sensitive interval of KNO_3 is from 0.1 to 1.5(mmol/m^3)^(-1),and interval of SEdRRN is from 0.01 and 0.1 d^(-1).The observational data from September 1998 to July 2000 obtained at SEATS station are used to validate the performance of biological model after parameters optimization.The results show that the modified model has a good capacity to reveal the biological process features,and the sensitivity analysis can save computational resources greatly during the model simulation.展开更多
Numerical simulation approaches have been widely applied to study mining induced subsidence,and they are potential methods to study the flooding induced uplift for abandoned mines.This paper gives an overview about di...Numerical simulation approaches have been widely applied to study mining induced subsidence,and they are potential methods to study the flooding induced uplift for abandoned mines.This paper gives an overview about different numerical approaches to simulate uplift induced by flooding abandoned underground mines,including three different hydraulic conditions,considering both unconfined and confined water conditions.Four basic simulation schemes using 1-dimensional rock column models verified by analytical solutions demonstrate these procedures.The results reveal that flooding induced uplift is mainly related to the pore pressure in the mine goaf.The parameter study documents that height and stiffness of the mine goaf have the strongest influence on maximum surface uplift.展开更多
Microwave remote sensing is one of the most useful methods for observing the ocean parameters. The Doppler frequency or interferometric phase of the radar echoes can be used for an ocean surface current speed retrieva...Microwave remote sensing is one of the most useful methods for observing the ocean parameters. The Doppler frequency or interferometric phase of the radar echoes can be used for an ocean surface current speed retrieval,which is widely used in spaceborne and airborne radars. While the effect of the ocean currents and waves is interactional. It is impossible to retrieve the ocean surface current speed from Doppler frequency shift directly. In order to study the relationship between the ocean surface current speed and the Doppler frequency shift, a numerical ocean surface Doppler spectrum model is established and validated with a reference. The input parameters of ocean Doppler spectrum include an ocean wave elevation model, a directional distribution function, and wind speed and direction. The suitable ocean wave elevation spectrum and the directional distribution function are selected by comparing the ocean Doppler spectrum in C band with an empirical geophysical model function(CDOP). What is more, the error sensitivities of ocean surface current speed to the wind speed and direction are analyzed. All these simulations are in Ku band. The simulation results show that the ocean surface current speed error is sensitive to the wind speed and direction errors. With VV polarization, the ocean surface current speed error is about 0.15 m/s when the wind speed error is 2 m/s, and the ocean surface current speed error is smaller than 0.3 m/s when the wind direction error is within 20° in the cross wind direction.展开更多
It is stressed in this paper that the seismic hazard analysis is a complex system with multi parameters and its parameter sensitivity should be studied as a whole. We have inspected the effects of parameters continui...It is stressed in this paper that the seismic hazard analysis is a complex system with multi parameters and its parameter sensitivity should be studied as a whole. We have inspected the effects of parameters continuing variation on the result and the effects of exceeding probability and upper bound magnitude of potential source on the parameter sensitivity. Some useful cognitions were got. We have overall calculated the interactions among the parameters including spatial distribution function of earthquake and suggested the concept of relative interaction. The relative interaction is not great than 10% with annual exceeding probability 2×10 -3 . The relative interaction between upper bound magnitude and annual mean rate and the one between upper bound magnitude and space distribution function are 17% and 22% respectively for potantial sources higher upper bound magnitudes in annual exceeding probability 10 -4 . It is convenient to judge whether considering the interactions in actual application. The importance and necessity of calculating the result of the probabilistic distribution are indicated. The effect of the result distribution by different parameter level weights is analyzed. The whole picture of multi parameter sensitivity is revealed. It makes necessary preparation for the reasonable uncertainty correction.展开更多
The mathematical model of single degree of freedom(DOF)nonlinear autonomous bearing system under constant flow supporting model is deduced.The single DOF nonlinear autonomous bearing system is transformed with the met...The mathematical model of single degree of freedom(DOF)nonlinear autonomous bearing system under constant flow supporting model is deduced.The single DOF nonlinear autonomous bearing system is transformed with the method of linear and nonlinear treatment,the mathematical expression and parameters sensitivity of relative error of stiffness and damping are presented.Finally,the main factors of magnetic-liquid double suspension bearing(MLDSB)are analyzed,and the influence on bearing performance indicators of single DOF nonlinear autonomous bearing system of main factors is revealed.The results show that linear stiffness/damping is the first part of equivalent stiffness/damping,and the second and third parts are high order minor term of Tayor series transform.The film thickness,the magnetic-liquid proportionality coefficient,the mass of rotor are the major influence factor of the bearing performance.The research can provide the theoretical reference for the design and nonlinear analysis of MLDSB.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant Nos.U22A20602,U2040221).
文摘A cascading failure of landslide dams caused by strong earthquakes or torrential rains in mountainous river valleys can pose great threats to people’s lives,properties,and infrastructures.In this study,based on the three-dimensional Reynoldsaveraged Navier-Stokes equations(RANS),the renormalization group(RNG)k-εturbulence model,suspended and bed load transport equations,and the instability discriminant formula of dam breach side slope,and the explicit finite volume method(FVM),a detailed numerical simulation model for calculating the hydro-morphodynamic characteristics of cascading dam breach process has been developed.The developed numerical model can simulate the breach hydrograph and the dam breach morphology evolution during the cascading failure process of landslide dams.A model test of the breaches of two cascading landslide dams has been used as the validation case.The comparison of the calculated and measured results indicates that the breach hydrograph and the breach morphology evolution process of the upstream and downstream dams are generally consistent with each other,and the relative errors of the key breaching parameters,i.e.,the peak breach flow and the time to peak of each dam,are less than±5%.Further,the comparison of the breach hydrographs of the upstream and downstream dams shows that there is an amplification effect of the breach flood on the cascading landslide dam failures.Three key parameters,i.e.,the distance between the upstream and the downstream dams,the river channel slope,and the downstream dam height,have been used to study the flood amplification effect.The parameter sensitivity analyses show that the peak breach flow at the downstream dam decreases with increasing distance between the upstream and the downstream dams,and the downstream dam height.Further,the peak breach flow at the downstream dam first increases and then decreases with steepening of the river channel slope.When the flood caused by the upstream dam failure flows to the downstream dam,it can produce a surge wave that overtops and erodes the dam crest,resulting in a lowering of the dam crest elevation.This has an impact on the failure occurrence time and the peak breach flow of the downstream dam.The influence of the surge wave on the downstream dam failure process is related to the volume of water that overtops the dam crest and the erosion characteristics of dam material.Moreover,the cascading failure case of the Xiaogangjian and Lower Xiaogangjian landslide dams has also been used as the representative case for validating the model.In comparisons of the calculated and measured breach hydrographs and final breach morphologies,the relative errors of the key dam breaching parameters are all within±10%,which verify the rationality of the model is applicable to real-world cases.Overall,the numerical model developed in this study can provide important technical support for the risk assessment and emergency treatment of failures of cascading landslide dams.
文摘This work highlights the unparalleled efficiency of the “n<sup>th</sup>-Order Function/ Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (n<sup>th</sup>-FASAM-N) by considering the well-known Nordheim-Fuchs reactor dynamics/safety model. This model describes a short-time self-limiting power excursion in a nuclear reactor system having a negative temperature coefficient in which a large amount of reactivity is suddenly inserted, either intentionally or by accident. This nonlinear paradigm model is sufficiently complex to model realistically self-limiting power excursions for short times yet admits closed-form exact expressions for the time-dependent neutron flux, temperature distribution and energy released during the transient power burst. The n<sup>th</sup>-FASAM-N methodology is compared to the extant “n<sup>th</sup>-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (n<sup>th</sup>-CASAM-N) showing that: (i) the 1<sup>st</sup>-FASAM-N and the 1<sup>st</sup>-CASAM-N methodologies are equally efficient for computing the first-order sensitivities;each methodology requires a single large-scale computation for solving the “First-Level Adjoint Sensitivity System” (1<sup>st</sup>-LASS);(ii) the 2<sup>nd</sup>-FASAM-N methodology is considerably more efficient than the 2<sup>nd</sup>-CASAM-N methodology for computing the second-order sensitivities since the number of feature-functions is much smaller than the number of primary parameters;specifically for the Nordheim-Fuchs model, the 2<sup>nd</sup>-FASAM-N methodology requires 2 large-scale computations to obtain all of the exact expressions of the 28 distinct second-order response sensitivities with respect to the model parameters while the 2<sup>nd</sup>-CASAM-N methodology requires 7 large-scale computations for obtaining these 28 second-order sensitivities;(iii) the 3<sup>rd</sup>-FASAM-N methodology is even more efficient than the 3<sup>rd</sup>-CASAM-N methodology: only 2 large-scale computations are needed to obtain the exact expressions of the 84 distinct third-order response sensitivities with respect to the Nordheim-Fuchs model’s parameters when applying the 3<sup>rd</sup>-FASAM-N methodology, while the application of the 3<sup>rd</sup>-CASAM-N methodology requires at least 22 large-scale computations for computing the same 84 distinct third-order sensitivities. Together, the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N methodologies are the most practical methodologies for computing response sensitivities of any order comprehensively and accurately, overcoming the curse of dimensionality in sensitivity analysis.
文摘This work presents the “n<sup>th</sup>-Order Feature Adjoint Sensitivity Analysis Methodology for Nonlinear Systems” (abbreviated as “n<sup>th</sup>-FASAM-N”), which will be shown to be the most efficient methodology for computing exact expressions of sensitivities, of any order, of model responses with respect to features of model parameters and, subsequently, with respect to the model’s uncertain parameters, boundaries, and internal interfaces. The unparalleled efficiency and accuracy of the n<sup>th</sup>-FASAM-N methodology stems from the maximal reduction of the number of adjoint computations (which are considered to be “large-scale” computations) for computing high-order sensitivities. When applying the n<sup>th</sup>-FASAM-N methodology to compute the second- and higher-order sensitivities, the number of large-scale computations is proportional to the number of “model features” as opposed to being proportional to the number of model parameters (which are considerably more than the number of features).When a model has no “feature” functions of parameters, but only comprises primary parameters, the n<sup>th</sup>-FASAM-N methodology becomes identical to the extant n<sup>th</sup> CASAM-N (“n<sup>th</sup>-Order Comprehensive Adjoint Sensitivity Analysis Methodology for Nonlinear Systems”) methodology. Both the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N methodologies are formulated in linearly increasing higher-dimensional Hilbert spaces as opposed to exponentially increasing parameter-dimensional spaces thus overcoming the curse of dimensionality in sensitivity analysis of nonlinear systems. Both the n<sup>th</sup>-FASAM-N and the n<sup>th</sup>-CASAM-N are incomparably more efficient and more accurate than any other methods (statistical, finite differences, etc.) for computing exact expressions of response sensitivities of any order with respect to the model’s features and/or primary uncertain parameters, boundaries, and internal interfaces.
基金Supported by the National Natural Science Foundation of China(10272109)
文摘In order to understand the effect of hardening ductility parameters and softening ductility parameters of the concrete damage plastic model in LS-DYNA,a sensitivity and reliability analysis of these parameters through a convenient cube unit test was conducted. The results showed that the peak strength strain was independent of the hardening ductility parameter DH,but affected by AH,BH,and CH. The softening ductility was mainly related to the softening ductility parameter AS,but not affected by the damage ductility exponent BS. In case that the model with default parameters failed to match the AS-controlled damage softening phase,an optimized model with an AS correction was developed. The corrected model with the AS value of 2 matched well with the code model,and exhibited good feasibility in predicting the stress-strain curve of different grades of concrete. Moreover,the practicability of the corrected model was further validated by the conventional triaxial test. The simulated curve exhibited favorable consistence with the trial curve. Therefore,the model with parameter correction could provide a prospective reference for predicting the mechanical properties of concrete.
基金supported by the National Natural Science Foundation of China (11142013, 11172260 and 11072214)the Doctoral Fund of Ministry of Education of China (20110101110016)the Fundamental Research Funds for the Central Universities of China(2011QNA4001)
文摘Input torque is the main power to maintain bipedal walking of robot, and can be calculated from trajectory planning and dynamic modeling on biped robot. During bipedal walking, the input torque is usually required to be adjusted due to some uncertain parameters arising from objective or subjective factors in the dynamical model to maintain the pre-planned stable trajectory. Here, a planar 5-link biped robot is used as an illustrating example to investigate the effects of uncertain parameters on the input torques. Kine-matic equations of the biped robot are firstly established by the third-order spline curves based on the trajectory planning method, and the dynamic modeling is accomplished by taking both the certain and uncertain parameters into account. Next, several evaluation indices on input torques are intro-duced to perform sensitivity analysis of the input torque with respect to the uncertain parameters. Finally, based on the Monte Carlo simulation, the values of evaluation indices on input torques are presented, from which all the robot param-eters are classified into three categories, i.e., strongly sensi-tive, sensitive and almost insensitive parameters.
文摘This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.
文摘Population viability analysis(PVA)is a tool to evaluate the risk of extinction for endangered species and aid conservation decision-making.The quality of PVA output is dependent on parameters related to population dynamics and life-history;however,it has been difficult to collect this information for the giant panda(Aliuropoda melanoleuca),a rare and endangered mammal native to China,confined to some 30 fragmented habitat patches.Since giant pandas are long-lived,mature late,have lower reproductive rates,and show little sexual dimorphism,obtaining data to perform adequate PVA has been difficult.Here,we develop a parameter sensitivity index by modeling the dynamics of six giant panda populations in the Minshan Mountains,in order to determine the parameters most influential to giant panda populations.Our data shows that the giant panda populations are most sensitive to changes in four female parameters:initial breeding age,reproductive rate,mortality rate between age 0 and 1,and mortality rate of adults.The parameter sensitivity index strongly correlated with initial population size,as smaller populations were more sensitive to changes in these four variables.This model suggests that demographic parameters of females have more influence on the results of PVA,indicating that females may play a more important role in giant panda population dynamics than males.Consequently,reintroduction of female individuals to a small giant panda population should be a high priority for conservation efforts.Our findings form a technical basis for the coming program of giant panda reintroduction,and inform which parameters are crucial to successfully and feasibly monitoring wild giant panda populations.
基金supported by the National Natural Science Foundation of China under Grant[52077027]Liaoning Province Science and Technology Major Project[No.2020JH1/10100020].
文摘In this paper, the axial-flux permanent magnet driver is modeledand analyzed in a simple and novel way under three-dimensional cylindricalcoordinates. The inherent three-dimensional characteristics of the deviceare comprehensively considered, and the governing equations are solved bysimplifying the boundary conditions. The axial magnetization of the sectorshapedpermanent magnets is accurately described in an algebraic form bythe parameters, which makes the physical meaning more explicit than thepurely mathematical expression in general series forms. The parameters of theBessel function are determined simply and the magnetic field distribution ofpermanent magnets and the air-gap is solved. Furthermore, the field solutionsare completely analytical, which provides convenience and satisfactoryaccuracy for modeling a series of electromagnetic performance parameters,such as the axial electromagnetic force density, axial electromagnetic force,and electromagnetic torque. The correctness and accuracy of the analyticalmodels are fully verified by three-dimensional finite element simulations and a15 kW prototype and the results of calculations, simulations, and experimentsunder three methods are highly consistent. The influence of several designparameters on magnetic field distribution and performance is studied and discussed.The results indicate that the modeling method proposed in this papercan calculate the magnetic field distribution and performance accurately andrapidly, which affords an important reference for the design and optimizationof axial-flux permanent magnet drivers.
基金supported by the CAS"Light of West China"Program (No.[2020]82)Key technology projects of Inner Mongolia Autonomous Region (Grant No.2020GG0306)+1 种基金Science and Technology Plan Projects of Alxa League (Grant No.AMY2020-18)Natural Science Foundation of Gansu Province (No.21JR7RA038).
文摘A challenge for the development of Land Surface Models(LSMs) is improving transpiration of water exchange and photosynthesis of carbon exchange between terrestrial plants and the atmosphere, both of which are governed by stoma in leaves. In the photosynthesis module of these LSMs, variations of parameters arising from diversity in plant functional types(PFTs) and climate remain unclear. Identifying sensitive parameters among all photosynthetic parameters before parameter estimation can not only reduce operation cost, but also improve the usability of photosynthesis models worldwide. Here, we analyzed 13 parameters of a biochemically-based photosynthesis model(FvCB), implemented in many LSMs, using two sensitivity analysis(SA) methods(i.e., the Sobol’ method and the Morris method) for setting up the parameter ensemble. Three different model performance metrics, i.e.,Root Mean Squared Error(RMSE), Nash Sutcliffe efficiency(NSE), and Standard Deviation(STDEV) were introduced for model assessment and sensitive parameters identification. The results showed that among all photosynthetic parameters only a small portion of parameters were sensitive, and the sensitive parameters were different across plant functional types: maximum rate of Rubisco activity(Vcmax25), maximum electron transport rate(Jmax25), triose phosphate use rate(TPU) and dark respiration in light(Rd) were sensitive in broad leafevergreen trees(BET), broad leaf-deciduous trees(BDT) and needle leaf-evergreen trees(NET), while only Vcmax25and TPU are sensitive in short vegetation(SV), dwarf trees and shrubs(DTS), and agriculture and grassland(AG). The two sensitivity analysis methods suggested a strong SA coherence;in contrast, different model performance metrics led to different SA results. This misfit suggests that more accurate values of sensitive parameters, specifically, species specific and seasonal variable parameters, are required to improve the performance of the FvCB model.
文摘To enhance the output torque and minimize the torque ripple of coaxial magnetic gear(CMG),a novel auxiliary flux modulator CMG with unequal magnetic poles is proposed.This design incorporates an inner rotor with an asymmetric sector and a trapezoidal combined N-S pole structure,featuring Halbach arrays for the arrangement of permanent magnets(PMs).The outer rotor PMs adopt a Spoke-type configuration.To optimize the CMG for high output torque and low torque ripple,a sensitivity analysis is conducted to identify key size parameters that significantly influence the optimization objectives.Based on the sensitivity hierarchy of these parameters,a multi-objective optimization analysis is performed using a genetic algorithm(GA)to determine the optimal structural parameter values of the CMG.In addition,a coaxial magnetic gear(CMG)topology with 4 inner and 17 outer pole pairs is adopted,and the parametric model is established.Finally,the electromagnetic properties of the CMG are evaluated using the finite element method.The results indicate a remarkable reduction in torque ripple,specifically by 46.15%.
基金supported by the Science Fund for Creative Research Groups of National Natural Science Foundation of China(51221004)the National Natural Science Foundation of China(11172260,11372270,and 51375434)+2 种基金the Higher School Specialized Research Fund for the Doctoral Program(20110101110016)the Science and technology project of Zhejiang Province(2013C31086)the Fundamental Research Funds forthe Central Universities of China(2013XZZX005)
文摘During bipedal walking,it is critical to detect and adjust the robot postures by feedback control to maintain its normal state amidst multi-source random disturbances arising from some unavoidable uncertain factors.The radical basis function(RBF)neural network model of a five-link biped robot is established,and two certain disturbances and a randomly uncertain disturbance are then mixed with the optimal torques in the network model to study the performance of the biped robot by several evaluation indices and a specific Poincar′e map.In contrast with the simulations,the response varies as desired under optimal inputting while the output is fluctuating in the situation of disturbance driving.Simulation results from noise inputting also show that the dynamics of the robot is less sensitive to the disturbance of knee joint input of the swing leg than those of the other three joints,the response errors of the biped will be increasing with higher disturbance levels,and especially there are larger output fluctuations in the knee and hip joints of the swing leg.
基金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.
基金supported by the National Key R&D Program of China(Grant No.2018YFC1405900)the Major Projects of National Science and Technology(Grant Nos.2016ZX05011-002,2016ZX05027-002-005)+3 种基金the National Natural Science Foundation of China(Grant No.41806073)the Natural Science Foundation of Shandong Province(Grant No.ZR2017BD014)Shandong Provincial Key Laboratory of Depositional Mineralization and Sedimentary Minerals,Shandong University of Science and Technology(Grant No.DMSM2017042)the Fundamental Research Funds for the Central Universities(Grant No.201964016)
文摘To establish the relationship among reservoir characteristics and rock physical parameters,we construct the well-bore rock physical models firstly,considering the influence factors,such as mineral composition,shale content,porosity,fluid type and saturation.Then with analyzing the change rules of elastic parameters along with the above influence factors and the cross-plots among elastic parameters,the sensitive elastic parameters of tight sandstone reservoir are determined,and the rock physics template of sweet spot is constructed to guide pre-stack seismic inversion.The results show that velocity ratio and Poisson impedance are the most sensitive elastic parameters to indicate the lithologic and gas-bearing properties of sweet spot in tight sandstone reservoir.The high-quality sweet spot is characterized by the lower velocity ratio and Poisson impedance.Finally,the actual seismic data are selected to predict the sweet spots in tight sandstone gas reservoirs,so as to verify the validity of the rock physical simulation results.The significant consistency between the relative logging curves and inversion results in different wells implies that the utilization of well-bore rock physical simulation can guide the prediction of sweet spot in tight sandstone gas reservoirs.
基金The authors sincerely appreciate the support from the National Natural Science Foundation of China[Grant Nos.51678101,52078093]Liaoning Revitalization Talents Program[Grant No.XLYC1905015]the Doctoral innovation Program of Dalian Maritime University[Grant No.BSCXXM016].
文摘The randomness of rock joint development is an important factor in the uncertainty of geotechnical engineering stability.In this study,a method is proposed to evaluate the reliability of intermittent jointed rock slope.The least squares support vector machine(LSSVM)evolved by a bacterial foraging optimization algorithm(BFOA)is used to establish a response surface model to express the mapping relationship between the intermittent joint parameters and the slope safety factor.The training samples are obtained from the numerical calculation based on the joint finite element method during this process.Considering the randomness of the intermittent joint parameters in the actual project,each parameter is evaluated at different locations on the site,and its distribution characteristics are counted.According to these statistical results,a large number of parameter combinations are obtained through Monte Carlo sampling.The trained machine learning mapping model is used to obtain the slope safety factor corresponding to each group,and these results are then used to obtain the slope reliability.When the research results were applied to slope disaster treatment along the Yalu River in China’s Jilin Province,it was found that the joint length and joint inclination angle both play key roles in rock slope stability,which should receive more attention in the slope treatment.In summary,this study establishes a method for evaluating the reliability of intermittent jointed rock slope based on an evolutionary SVM model,and its feasibility is verified by engineering application.
基金supported by the National High Technology Research and Development Program of China (863 Program,Grant No. 2008AA09A105-04)
文摘Stress Joint (SJ) plays a key role in the Top Tensioned Riser (TTR) system for deep water engineering. A preliminary design method of tapered SJ is proposed in the paper, which could help designers obtain accurate design data. After a further sensitive analysis is carried out, the related parameters choice and control methods are recommended in the engineering practice. By taking the extreme environment conditions into consideration, the effects of bending stress reduction and curve control are analyzed, and the 3-D FE models are established by ABQOUS for numerical evaluation to verify the correctness of design results. At last, dynamic analysis and fatigue analysis, based on actual project, are carried out with designed stress joint. The analysis results prove the feasibility and guidance of this method in the practical engineering applications.
基金The National Natural Science Foundation of China under contract Nos 41206023,41222038 and 41076011the National Basic Research Project(973 Program)of China under contract No.2011CB403606+2 种基金the China-Korea Joint Ocean Research Center"Cooperation on the Development of Basic Technologies for the Yellow Sea and East China Sea Operational Oceanographic System(YOOS)"the Public Science and Technology Research Funds Projects of Ocean under contrcat No.201205018the"Strategic Priority Research Program"of the Chinese Academy of Sciences,under contract No.XDA01020304
文摘In order to develop a coupled basin scale model of ocean circulation and biogeochemical cycling,we present a biogeochemical model including 12 components to study the ecosystem in the China coastal seas(CCS).The formulation of phytoplankton mortality and zooplankton growth are modified according to biological characteristics of CCS.The four sensitivity biological parameters,zooplankton assimilation efficiency rate(ZooAE_N),zooplankton basal metabolism rate(ZooBM),maximum specific growth rate of zooplankton(μ_(20)) and maximum chlorophyll to carbon ratio(Chl2C_m) are obtained in sensitivity experiments for the phytoplankton,and experiments about the parameter μ_(20'),half-saturation for phytoplankton NO_3 uptake(K_(NO_3)) and remineralization rate of small detritusN(SDeRRN) are conducted.The results demonstrate that the biogeochemical model is quite sensitive to the zooplankton grazing parameter when it ranges from 0.1 to 1.2 d^(-1).The K_(NO_3) and SDeRRN also play an important role in determining the nitrogen cycle within certain ranges.The sensitive interval of KNO_3 is from 0.1 to 1.5(mmol/m^3)^(-1),and interval of SEdRRN is from 0.01 and 0.1 d^(-1).The observational data from September 1998 to July 2000 obtained at SEATS station are used to validate the performance of biological model after parameters optimization.The results show that the modified model has a good capacity to reveal the biological process features,and the sensitivity analysis can save computational resources greatly during the model simulation.
基金The first author would like to thank the China Scholarship Council(CSC No.201806430001)for the financial support for his study in Germany.The authors are grateful to the anonymous reviewer for the recommendations to improve the paper。
文摘Numerical simulation approaches have been widely applied to study mining induced subsidence,and they are potential methods to study the flooding induced uplift for abandoned mines.This paper gives an overview about different numerical approaches to simulate uplift induced by flooding abandoned underground mines,including three different hydraulic conditions,considering both unconfined and confined water conditions.Four basic simulation schemes using 1-dimensional rock column models verified by analytical solutions demonstrate these procedures.The results reveal that flooding induced uplift is mainly related to the pore pressure in the mine goaf.The parameter study documents that height and stiffness of the mine goaf have the strongest influence on maximum surface uplift.
基金The National Natural Science Foundation of China under contract No.41606202the National Key Research and Development Program of China under contract No.2016YFC1401002the Open Fund of Key Laboratory of State Oceanic Administration(SOA) for Space Ocean Remote Sensing and Application under contract No.201601001
文摘Microwave remote sensing is one of the most useful methods for observing the ocean parameters. The Doppler frequency or interferometric phase of the radar echoes can be used for an ocean surface current speed retrieval,which is widely used in spaceborne and airborne radars. While the effect of the ocean currents and waves is interactional. It is impossible to retrieve the ocean surface current speed from Doppler frequency shift directly. In order to study the relationship between the ocean surface current speed and the Doppler frequency shift, a numerical ocean surface Doppler spectrum model is established and validated with a reference. The input parameters of ocean Doppler spectrum include an ocean wave elevation model, a directional distribution function, and wind speed and direction. The suitable ocean wave elevation spectrum and the directional distribution function are selected by comparing the ocean Doppler spectrum in C band with an empirical geophysical model function(CDOP). What is more, the error sensitivities of ocean surface current speed to the wind speed and direction are analyzed. All these simulations are in Ku band. The simulation results show that the ocean surface current speed error is sensitive to the wind speed and direction errors. With VV polarization, the ocean surface current speed error is about 0.15 m/s when the wind speed error is 2 m/s, and the ocean surface current speed error is smaller than 0.3 m/s when the wind direction error is within 20° in the cross wind direction.
文摘It is stressed in this paper that the seismic hazard analysis is a complex system with multi parameters and its parameter sensitivity should be studied as a whole. We have inspected the effects of parameters continuing variation on the result and the effects of exceeding probability and upper bound magnitude of potential source on the parameter sensitivity. Some useful cognitions were got. We have overall calculated the interactions among the parameters including spatial distribution function of earthquake and suggested the concept of relative interaction. The relative interaction is not great than 10% with annual exceeding probability 2×10 -3 . The relative interaction between upper bound magnitude and annual mean rate and the one between upper bound magnitude and space distribution function are 17% and 22% respectively for potantial sources higher upper bound magnitudes in annual exceeding probability 10 -4 . It is convenient to judge whether considering the interactions in actual application. The importance and necessity of calculating the result of the probabilistic distribution are indicated. The effect of the result distribution by different parameter level weights is analyzed. The whole picture of multi parameter sensitivity is revealed. It makes necessary preparation for the reasonable uncertainty correction.
基金the National Natural Science Foundation of China(No.51705445)General Project of Natural Science Foundation of Hebei Province(No.E2016203324)Youth Fund Project of Scientific Research Project of Hebei University(No.QN202013).
文摘The mathematical model of single degree of freedom(DOF)nonlinear autonomous bearing system under constant flow supporting model is deduced.The single DOF nonlinear autonomous bearing system is transformed with the method of linear and nonlinear treatment,the mathematical expression and parameters sensitivity of relative error of stiffness and damping are presented.Finally,the main factors of magnetic-liquid double suspension bearing(MLDSB)are analyzed,and the influence on bearing performance indicators of single DOF nonlinear autonomous bearing system of main factors is revealed.The results show that linear stiffness/damping is the first part of equivalent stiffness/damping,and the second and third parts are high order minor term of Tayor series transform.The film thickness,the magnetic-liquid proportionality coefficient,the mass of rotor are the major influence factor of the bearing performance.The research can provide the theoretical reference for the design and nonlinear analysis of MLDSB.