The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critica...The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.展开更多
This paper proposes a non-intrusive computational method for mechanical dynamic systems involving a large-scale of interval uncertain parameters,aiming to reduce the computational costs and improve accuracy in determi...This paper proposes a non-intrusive computational method for mechanical dynamic systems involving a large-scale of interval uncertain parameters,aiming to reduce the computational costs and improve accuracy in determining bounds of system response.The screening method is firstly used to reduce the scale of active uncertain parameters.The sequential high-order polynomials surrogate models are then used to approximate the dynamic system’s response at each time step.To reduce the sampling cost of constructing surrogate model,the interaction effect among uncertain parameters is gradually added to the surrogate model by sequentially incorporating samples from a candidate set,which is composed of vertices and inner grid points.Finally,the points that may produce the bounds of the system response at each time step are searched using the surrogate models.The optimization algorithm is used to locate extreme points,which contribute to determining the inner points producing system response bounds.Additionally,all vertices are also checked using the surrogate models.A vehicle nonlinear dynamic model with 72 uncertain parameters is presented to demonstrate the accuracy and efficiency of the proposed uncertain computational method.展开更多
After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir develo...After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.展开更多
For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for s...For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for solving the dynamical equations of the constrained Birkhoffian system is provided.First the differential equations of motion for the constrained Birkhoffian system as well as for the corresponding free Birkhoffian system are established.Secondly a system of auxiliary equations is constructed and the general solution of the equations is found.Finally by varying the parameters and utilizing the properties of the generalized canonical transformation of the Birkhoffian system the solution of the problem can be obtained.The proposed method reveals the inherent relationship between the solution of a free Birkhoffian system and that of a constrained Birkhoffian system. The research results are of universal significance which can be further used in a variety of constrained mechanical systems such as non-conservative systems and nonholonomic systems etc.展开更多
Modal parameters can accurately characterize the structural dynamic properties and assess the physical state of the structure.Therefore,it is particularly significant to identify the structural modal parameters accordi...Modal parameters can accurately characterize the structural dynamic properties and assess the physical state of the structure.Therefore,it is particularly significant to identify the structural modal parameters according to the monitoring data information in the structural health monitoring(SHM)system,so as to provide a scientific basis for structural damage identification and dynamic model modification.In view of this,this paper reviews methods for identifying structural modal parameters under environmental excitation and briefly describes how to identify structural damages based on the derived modal parameters.The paper primarily introduces data-driven modal parameter recognition methods(e.g.,time-domain,frequency-domain,and time-frequency-domain methods,etc.),briefly describes damage identification methods based on the variations of modal parameters(e.g.,natural frequency,modal shapes,and curvature modal shapes,etc.)and modal validation methods(e.g.,Stability Diagram and Modal Assurance Criterion,etc.).The current status of the application of artificial intelligence(AI)methods in the direction of modal parameter recognition and damage identification is further discussed.Based on the pre-vious analysis,the main development trends of structural modal parameter recognition and damage identification methods are given to provide scientific references for the optimized design and functional upgrading of SHM systems.展开更多
Ion optics are crucial components of ion thrusters and the study of the different ion extraction solutions used since the beginning of the electric propulsion era is essential to understand the evolution of ion engine...Ion optics are crucial components of ion thrusters and the study of the different ion extraction solutions used since the beginning of the electric propulsion era is essential to understand the evolution of ion engines. This work describes ion engine grids＇ main functions, parameters and issues related to thermal expansion and sputter erosion, and then introduces a review of ion optics used for significant launched and tested ion thrusters since 1970. Configurations, geometries, materials and fabrication methods are analyzed to understand when typical ion thrusters use two or three grids, what are the thicknesses and aperture sizes of the screen, accelerator and decelerator grids, why molybdenum and carbon-based materials such as pyrolytic graphite and carbon–carbon are the best available options for ion optics and what is the manufacturing method for each material.展开更多
In the traditional strength reduction method,the cohesion and the friction angle adopt the same reduction parameter,resulting in equivalent proportional reduction.This method does not consider the different effects of...In the traditional strength reduction method,the cohesion and the friction angle adopt the same reduction parameter,resulting in equivalent proportional reduction.This method does not consider the different effects of the cohesion and friction angle on the stability of the same slope and is defective to some extent.Regarding this defect,a strength reduction method based on double reduction parameters,which adopts different reduction parameters,is proposed.The core of the double-parameter reduction method is the matching reduction principle of the slope with different angles.This principle is represented by the ratio of the reduction parameter of the cohesion to that of the friction angle,described as η.With the increase in the slopeangle,ηincreases; in particular,when the slope angle is 45°,tηis 1.0.Through the matching reduction principle,different safety margin factors can be calculated for the cohesion and friction angle.In combination with these two safety margin factors,a formula for calculating the overall safety factor of the slope is proposed,reflecting the different contributions of the cohesion and friction angle to the slope stability.Finally,it is shown that the strength reduction method based on double reduction parameters acquires a larger safety factor than the classic limit equilibrium method,but the calculation results are very close to those obtained by the limit equilibrium method.展开更多
This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to impro...This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.展开更多
A new numerical technique named as fuzzy finite difference method is proposed to solve the heat conduction problems with fuzzy uncertainties in both the phys- ical parameters and initial/boundary conditions. In virtue...A new numerical technique named as fuzzy finite difference method is proposed to solve the heat conduction problems with fuzzy uncertainties in both the phys- ical parameters and initial/boundary conditions. In virtue of the level-cut method, the difference discrete equations with fuzzy parameters are equivalently transformed into groups of interval equations. New stability analysis theory suited to fuzzy difference schemes is developed. Based on the parameter perturbation method, the interval ranges of the uncertain temperature field can be approximately predicted. Subsequently, fuzzy solutions to the original difference equations are obtained by the fuzzy resolution theorem. Two numerical examples are given to demonstrate the feasibility and efficiency of the presented method for solving both steady-state and transient heat conduction problems.展开更多
The inverse heat conduction method is one of methods to identify the casting simulation parameters. A new inverse method was presented according to the Tikhonov regularization theory. One appropriate regularized funct...The inverse heat conduction method is one of methods to identify the casting simulation parameters. A new inverse method was presented according to the Tikhonov regularization theory. One appropriate regularized functional was established, and the functional was solved by the sensitivity coefficient and Newtonaphson iteration method. Moreover, the orthogonal experimental design was used to estimate the appropriate initial value and variation domain of each variable to decrease the number of iteration and improve the identification accuracy and efficiency. It illustrated a detailed case of AlSiTMg sand mold casting and the temperature measurement experiment was done. The physical properties of sand mold and the interracial heat transfer coefficient were identified at the meantime. The results indicated that the new regularization method was efficient in overcoming the ill-posedness of the inverse heat conduction problem and improving the stability and accuracy of the solutions.展开更多
A hybrid numerical-experimental approach to identify elastic modulus of a textile composite panel using vibration test data is proposed and investi- gated. Homogenization method is adopted to predict the initial value...A hybrid numerical-experimental approach to identify elastic modulus of a textile composite panel using vibration test data is proposed and investi- gated. Homogenization method is adopted to predict the initial values of elastic parameters of the composite, and parameter identification is transformed to an optimization problem in which the objective function is the minimization of the discrepancies between the experimental and numerical modal data. Case study is conducted employing a woven fabric reinforced composite panel. Three parameters （Ell, E22, G12） with higher sensitivities are selected to be identified. It is shown that the elastic parameters can be accurately identified from experimental modal data.展开更多
The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their ...The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time.These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation.Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters calculation without any need of adjoint flux. The most straightforward scheme is known as the ‘‘prompt method'' and has been used widely in literature. The main objective of this article is dedicated to point kinetic parameters calculation in Tehran research reactor(TRR) using deterministic as well as probabilistic techniques. WIMS-D5B and CITATION codes have been used in deterministic calculation of forward and adjoint fluxes in the TRR core. On the other hand, the MCNP Monte Carlo code has been employed in the ‘‘prompt method''scheme for effective delayed neutron fraction evaluation.Deterministic results have been cross-checked with probabilistic ones and validated with SAR and experimental data. In comparison with experimental results, the relativedifferences of deterministic as well as probabilistic methods are 7.6 and 3.2%, respectively. These quantities are10.7 and 6.4%, respectively, in comparison with SAR report.展开更多
This paper focuses on studying the integration method of a generalized Birkhoffian system.The method of variation on parameters for the dynamical equations of a generalized Birkhoffian system is presented.The procedur...This paper focuses on studying the integration method of a generalized Birkhoffian system.The method of variation on parameters for the dynamical equations of a generalized Birkhoffian system is presented.The procedure for solving the problem can be divided into two steps：the first step,a system of auxiliary equations is constructed and its general solution is given;the second step,the parameters are varied,and the solution of the problem is obtained by using the properties of generalized canonical transformation.The method of variation on parameters for the generalized Birkhoffian system is of universal significance,and we take a nonholonomic system and a nonconservative system as examples to illustrate the application of the results of this paper.展开更多
The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows ...The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows through an aquifer,its composition and temperature may variation dependent on the aquifer condition through which it flows.Thus,hydrologic investigations can also provide useful information about the subsurface geology of a region.But because such studies investigate processes that follow under the Earth's shallow,obtaining the information necessary to answer these questions is not continuously easy.Springs,which discharge groundwater table directly,afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity(T)and storativity(S)are vital for the evaluation of groundwater resources.There are several methods to estimate the accurate aquifer parameters(i.e.hydrograph analysis,pumping test,etc.).In initial days,these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory.The simultaneous information on the hydraulic behavior of the well(borehole)that provides on this method,the reservoir and the reservoir boundaries,are important for efficient aquifer and well data management and analysis.The most common in-situ test is pumping test performed on wells,which involves the measurement of the fall and increase of groundwater level with respect to time.The alteration in groundwater level(drawdown/recovery)is caused due to pumping of water from the well.Theis(1935)was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer.It is essential to know the transmissivity(T=Kb,where b is the aquifer thickness;pumping flow rate,Q=TW(dh/dl)flow through an aquifer)and storativity(confined aquifer:S=bS_s,unconfined:S=S_y),for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer.The determination of aquifer's parameters is an important basis for groundwater resources evaluation,numerical simulation,development and protection as well as scientific management.For determining aquifer's parameters,pumping test is a main method.A case study shows that these techniques have been fast speed and high correctness.The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.展开更多
The Hansen solubility parameters(HSP)are frequently used for solvent selection and characterization of polymers,and are directly related to the suspension behavior of pigments in solvent mixtures.The performance of cu...The Hansen solubility parameters(HSP)are frequently used for solvent selection and characterization of polymers,and are directly related to the suspension behavior of pigments in solvent mixtures.The performance of currently available group contribution(GC)methods for HSP were evaluated and found to be insufficient for computer-aided product design(CAPD)of paints and coatings.A revised and,for this purpose,improved GC method is presented for estimating HSP of organic compounds,intended for organic pigments.Due to the significant limitations of GC methods,an uncertainty analysis and parameter confidence intervals are provided in order to better quantify the estimation accuracy of the proposed approach.Compared to other applicable GC methods,the prediction error is reduced significantly with average absolute errors of 0.45 MPa^(1/2),1.35 MPa^(1/2),and 1.09 MPa^(1/2) for the partial dispersion(δD),polar(δP)and hydrogen-bonding(δH)solubility parameters respectively for a database of 1106 compounds.The performance for organic pigments is comparable to the overall method performance,with higher average errors forδD and lower average errors forδP andδH.展开更多
Among the methods used for evaluating the potential hydraulic erodibility of rock,the most common methods are those based on the correlation between the force of flowing water and the capacity of a rock to resist eros...Among the methods used for evaluating the potential hydraulic erodibility of rock,the most common methods are those based on the correlation between the force of flowing water and the capacity of a rock to resist erosion,such as Annandale’s and Pells’methods.The capacity of a rock to resist erosion is evaluated based on erodibility indices that are determined from specific geomechanical parameters of a rock mass.These indices include unconfined compressive strength(UCS)of rock,rock block size,joint shear strength,a block’s shape and orientation relative to the direction of flow,joint openings,and the nature of the surface to be potentially eroded.However,it is difficult to determine the relevant geomechanical parameters for evaluating the hydraulic erodibility of rock.The assessment of eroded unlined spillways of dams has shown that the capacity of a rock to resist erosion is not accurately evaluated.Using more than 100 case studies,we develop a method to determine the relevant geomechanical parameters for evaluating the hydraulic erodibility of rock in unlined spillways.The UCS of rock is found not to be a relevant parameter for evaluating the hydraulic erodibility of rock.On the other hand,we find that the use of three-dimensional(3D)block volume measurements,instead of the block size factor used in Annandale’s method,improves the rock block size estimation.Furthermore,the parameter representing the effect of a rock block’s shape and orientation relative to the direction of flow,as considered in Pells’method,is more accurate than the parameter adopted by Annandale’s method.展开更多
Particle swarm optimization （PSO）, like other evolutionary algorithms is a population-based stochastic algorithm inspired from the metaphor of social interaction in birds, insects, wasps, etc. It has been used for f...Particle swarm optimization （PSO）, like other evolutionary algorithms is a population-based stochastic algorithm inspired from the metaphor of social interaction in birds, insects, wasps, etc. It has been used for finding promising solutions in complex search space through the interaction of particles in a swarm. It is a well recognized fact that the performance of evolutionary algorithms to a great extent depends on the choice of appropriate strategy/operating parameters like population size, crossover rate, mutation rate, crossover operator, etc. Generally, these parameters are selected through hit and trial process, which is very unsystematic and requires rigorous experimentation. This paper proposes a systematic based on Taguchi method reasoning scheme for rapidly identifying the strategy parameters for the PSO algorithm. The Taguchi method is a robust design approach using fractional factorial design to study a large number of parameters with small number of experiments. Computer simulations have been performed on two benchmark functionsiRosenbrock function and Griewank functionito validate the approach.展开更多
A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variable...A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.展开更多
This paper discusses a kind of implicit iterative methods with some variable parameters, which are called control parameters, for solving ill-posed operator equations. The theoretical results show that the new methods...This paper discusses a kind of implicit iterative methods with some variable parameters, which are called control parameters, for solving ill-posed operator equations. The theoretical results show that the new methods always lead to optimal convergence rates and have some other important features, especially the methods can be implemented parallelly.展开更多
During the service lives of ballasted tracks,the ballast experiences degradation,such as breakage and fragmentation,which reduces the stability of the tracks.Ballast breakage is directly related to the shapes,sizes,an...During the service lives of ballasted tracks,the ballast experiences degradation,such as breakage and fragmentation,which reduces the stability of the tracks.Ballast breakage is directly related to the shapes,sizes,and stress states of particles.The key to solving this problem is to determine the breakage parameters of railway ballast.In this study,through uniaxial compression and uniaxial shear tests,the factors associated with ballast breakage for a simple stress state were obtained.Then,a refined discrete element model of railway ballast was established.The Box-Behnken method was used to design a simulation test,and a response surface method was used to obtain the optimal ballast breakage parameters.Lastly,the results of a direct shear test were compared with the simulation results to verify the correctness of the parameters and to study the shear resistance of the ballast aggregates.The results showed that for ballast particles with sizes of 22.4–63.0 mm,the characteristic stresses of compression and shear were not significantly affected by the size of the particles,but were greatly affected by their shape.In particular,long particles were more likely to break.During the direct shear test,the cumulative ballast breakage ratio gradually increased to 13.97%.The optimal breakage parameters of the ballast determined by this approach have high application value in the management of ballasted tracks.展开更多
基金National Key R&D Program of China(No.2017YFB1304000)Fundamental Research Funds for the Central Universities,China(No.2232023G-05-1)。
文摘The high-speed winding spindle employs a flexible support system incorporating rubber O-rings.By precisely configuring the structural parameters and the number of the O-rings,the spindle can stably surpass its critical speed points and maintain operational stability across the entire working speed range.However,the support stiffness and damping of rubber O-rings exhibit significant nonlinear frequency dependence.Conventional experimental methods for deriving equivalent stiffness and damping,based on the principle of the forced non-resonance method,require fabricating custom setups for each O-ring specification and conducting vibration tests at varying frequencies,resulting in low efficiency and high costs.This study proposes a hybrid simulation-experimental method for dynamic parameter identification.Firstly,the frequency-dependent dynamic parameters of a specific O-ring support system are experimentally obtained.Subsequently,a corresponding parametric finite element model is established to simulate and solve the equivalent elastic modulus and equivalent stiffness-damping coefficient of this O-ring support system.Ultimately,after iterative simulation,the simulated and experimental results achieve a 99.7%agreement.The parametric finite element model developed herein can directly simulate and inversely estimate frequency-dependent dynamic parameters for O-rings of different specifications but identical elastic modulus.
基金supported by the National Natural Science Foundation of China(Grant No.12272142)Fundamental Research Funds for the Central Universities(Grant No.2172021XXJS048)。
文摘This paper proposes a non-intrusive computational method for mechanical dynamic systems involving a large-scale of interval uncertain parameters,aiming to reduce the computational costs and improve accuracy in determining bounds of system response.The screening method is firstly used to reduce the scale of active uncertain parameters.The sequential high-order polynomials surrogate models are then used to approximate the dynamic system’s response at each time step.To reduce the sampling cost of constructing surrogate model,the interaction effect among uncertain parameters is gradually added to the surrogate model by sequentially incorporating samples from a candidate set,which is composed of vertices and inner grid points.Finally,the points that may produce the bounds of the system response at each time step are searched using the surrogate models.The optimization algorithm is used to locate extreme points,which contribute to determining the inner points producing system response bounds.Additionally,all vertices are also checked using the surrogate models.A vehicle nonlinear dynamic model with 72 uncertain parameters is presented to demonstrate the accuracy and efficiency of the proposed uncertain computational method.
基金funded by the 14th Five-Year Plan Major Science and Technology Project of CNOOC project number KJGG2021-0506.
文摘After a long period of water flooding development,the oilfield has entered the middle and high water cut stage.The physical properties of reservoirs are changed by water erosion,which directly impacts reservoir development.Conventional numerical reservoir simulation methodologies typically employ static assumptions for model construction,presuming invariant reservoir geological parameters throughout the development process while neglecting the reservoir’s temporal evolution characteristics.Although such simplifications reduce computational complexity,they introduce substantial descriptive inaccuracies.Therefore,this paper proposes a meshless numerical simulation method for reservoirs that considers time-varying characteristics.This method avoids the meshing in traditional numerical simulation methods.From the fluid flow perspective,the reservoir’s computational domain is discretized into a series of connection units.An influence domain with a certain radius centered on the nodes is selected,and one-dimensional connection units are established between the nodes to achieve the characterization of the flow topology structure of the reservoir.In order to reflect the dynamic evolution of the reservoir’s physical properties during the water injection development process,the time-varying characteristics are incorporated into the formula of the seepage characteristic parameters in the meshless calculation.The change relationship of the permeability under different surface fluxes is considered to update the calculated connection conductivity in real time.By combining with the seepage control equation for solution,a time-varying meshless numerical simulation method is formed.The results show that compared with the numerical simulationmethod of the connection elementmethod(CEM)that only considers static parameters,this method has higher simulation accuracy and can better simulate the real migration and distribution of oil and water in the reservoir.Thismethod improves the accuracy of reservoir numerical simulation and the development effect of oilfields,providing a scientific basis for optimizing the water injection strategy,adjusting the production plan,and extending the effective production cycle of the oilfield.
基金The National Natural Science Foundation of China(No.10972151,11272227)
文摘For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for solving the dynamical equations of the constrained Birkhoffian system is provided.First the differential equations of motion for the constrained Birkhoffian system as well as for the corresponding free Birkhoffian system are established.Secondly a system of auxiliary equations is constructed and the general solution of the equations is found.Finally by varying the parameters and utilizing the properties of the generalized canonical transformation of the Birkhoffian system the solution of the problem can be obtained.The proposed method reveals the inherent relationship between the solution of a free Birkhoffian system and that of a constrained Birkhoffian system. The research results are of universal significance which can be further used in a variety of constrained mechanical systems such as non-conservative systems and nonholonomic systems etc.
基金supported by the Innovation Foundation of Provincial Education Department of Gansu(2024B-005)the Gansu Province National Science Foundation(22YF7GA182)the Fundamental Research Funds for the Central Universities(No.lzujbky2022-kb01)。
文摘Modal parameters can accurately characterize the structural dynamic properties and assess the physical state of the structure.Therefore,it is particularly significant to identify the structural modal parameters according to the monitoring data information in the structural health monitoring(SHM)system,so as to provide a scientific basis for structural damage identification and dynamic model modification.In view of this,this paper reviews methods for identifying structural modal parameters under environmental excitation and briefly describes how to identify structural damages based on the derived modal parameters.The paper primarily introduces data-driven modal parameter recognition methods(e.g.,time-domain,frequency-domain,and time-frequency-domain methods,etc.),briefly describes damage identification methods based on the variations of modal parameters(e.g.,natural frequency,modal shapes,and curvature modal shapes,etc.)and modal validation methods(e.g.,Stability Diagram and Modal Assurance Criterion,etc.).The current status of the application of artificial intelligence(AI)methods in the direction of modal parameter recognition and damage identification is further discussed.Based on the pre-vious analysis,the main development trends of structural modal parameter recognition and damage identification methods are given to provide scientific references for the optimized design and functional upgrading of SHM systems.
文摘Ion optics are crucial components of ion thrusters and the study of the different ion extraction solutions used since the beginning of the electric propulsion era is essential to understand the evolution of ion engines. This work describes ion engine grids＇ main functions, parameters and issues related to thermal expansion and sputter erosion, and then introduces a review of ion optics used for significant launched and tested ion thrusters since 1970. Configurations, geometries, materials and fabrication methods are analyzed to understand when typical ion thrusters use two or three grids, what are the thicknesses and aperture sizes of the screen, accelerator and decelerator grids, why molybdenum and carbon-based materials such as pyrolytic graphite and carbon–carbon are the best available options for ion optics and what is the manufacturing method for each material.
基金Project(KZCX2-YW-T12)supported by the Chinese Academy of Science,China
文摘In the traditional strength reduction method,the cohesion and the friction angle adopt the same reduction parameter,resulting in equivalent proportional reduction.This method does not consider the different effects of the cohesion and friction angle on the stability of the same slope and is defective to some extent.Regarding this defect,a strength reduction method based on double reduction parameters,which adopts different reduction parameters,is proposed.The core of the double-parameter reduction method is the matching reduction principle of the slope with different angles.This principle is represented by the ratio of the reduction parameter of the cohesion to that of the friction angle,described as η.With the increase in the slopeangle,ηincreases; in particular,when the slope angle is 45°,tηis 1.0.Through the matching reduction principle,different safety margin factors can be calculated for the cohesion and friction angle.In combination with these two safety margin factors,a formula for calculating the overall safety factor of the slope is proposed,reflecting the different contributions of the cohesion and friction angle to the slope stability.Finally,it is shown that the strength reduction method based on double reduction parameters acquires a larger safety factor than the classic limit equilibrium method,but the calculation results are very close to those obtained by the limit equilibrium method.
基金supported by the National Natural Science Foundation of China(No.51379039)the Excellent Young Scientists Fund(No.51222904)
文摘This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.
基金supported by the National Special Fund for Major Research Instrument Development(2011YQ140145)111 Project(B07009)+1 种基金the National Natural Science Foundation of China(11002013)Defense Industrial Technology Development Program(A2120110001 and B2120110011)
文摘A new numerical technique named as fuzzy finite difference method is proposed to solve the heat conduction problems with fuzzy uncertainties in both the phys- ical parameters and initial/boundary conditions. In virtue of the level-cut method, the difference discrete equations with fuzzy parameters are equivalently transformed into groups of interval equations. New stability analysis theory suited to fuzzy difference schemes is developed. Based on the parameter perturbation method, the interval ranges of the uncertain temperature field can be approximately predicted. Subsequently, fuzzy solutions to the original difference equations are obtained by the fuzzy resolution theorem. Two numerical examples are given to demonstrate the feasibility and efficiency of the presented method for solving both steady-state and transient heat conduction problems.
文摘The inverse heat conduction method is one of methods to identify the casting simulation parameters. A new inverse method was presented according to the Tikhonov regularization theory. One appropriate regularized functional was established, and the functional was solved by the sensitivity coefficient and Newtonaphson iteration method. Moreover, the orthogonal experimental design was used to estimate the appropriate initial value and variation domain of each variable to decrease the number of iteration and improve the identification accuracy and efficiency. It illustrated a detailed case of AlSiTMg sand mold casting and the temperature measurement experiment was done. The physical properties of sand mold and the interracial heat transfer coefficient were identified at the meantime. The results indicated that the new regularization method was efficient in overcoming the ill-posedness of the inverse heat conduction problem and improving the stability and accuracy of the solutions.
基金supported by the Program for New Century Excellent Talents in University(NCET11-0086)the National Natural Science Foundation of China(10902024)+1 种基金the Doctoral Program of Higher Education of China(20130092120039)the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD-1105007001)
文摘A hybrid numerical-experimental approach to identify elastic modulus of a textile composite panel using vibration test data is proposed and investi- gated. Homogenization method is adopted to predict the initial values of elastic parameters of the composite, and parameter identification is transformed to an optimization problem in which the objective function is the minimization of the discrepancies between the experimental and numerical modal data. Case study is conducted employing a woven fabric reinforced composite panel. Three parameters （Ell, E22, G12） with higher sensitivities are selected to be identified. It is shown that the elastic parameters can be accurately identified from experimental modal data.
文摘The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time.These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation.Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters calculation without any need of adjoint flux. The most straightforward scheme is known as the ‘‘prompt method'' and has been used widely in literature. The main objective of this article is dedicated to point kinetic parameters calculation in Tehran research reactor(TRR) using deterministic as well as probabilistic techniques. WIMS-D5B and CITATION codes have been used in deterministic calculation of forward and adjoint fluxes in the TRR core. On the other hand, the MCNP Monte Carlo code has been employed in the ‘‘prompt method''scheme for effective delayed neutron fraction evaluation.Deterministic results have been cross-checked with probabilistic ones and validated with SAR and experimental data. In comparison with experimental results, the relativedifferences of deterministic as well as probabilistic methods are 7.6 and 3.2%, respectively. These quantities are10.7 and 6.4%, respectively, in comparison with SAR report.
基金supported by the National Natural Science Foundation of China (10972151)
文摘This paper focuses on studying the integration method of a generalized Birkhoffian system.The method of variation on parameters for the dynamical equations of a generalized Birkhoffian system is presented.The procedure for solving the problem can be divided into two steps：the first step,a system of auxiliary equations is constructed and its general solution is given;the second step,the parameters are varied,and the solution of the problem is obtained by using the properties of generalized canonical transformation.The method of variation on parameters for the generalized Birkhoffian system is of universal significance,and we take a nonholonomic system and a nonconservative system as examples to illustrate the application of the results of this paper.
文摘The origin and movement of groundwater are the fundamental questions that address both the temporal and spatial aspects of ground water run and water supply related issues in hydrological systems.As groundwater flows through an aquifer,its composition and temperature may variation dependent on the aquifer condition through which it flows.Thus,hydrologic investigations can also provide useful information about the subsurface geology of a region.But because such studies investigate processes that follow under the Earth's shallow,obtaining the information necessary to answer these questions is not continuously easy.Springs,which discharge groundwater table directly,afford to study subsurface hydrogeological processes.The present study of estimation of aquifer factors such as transmissivity(T)and storativity(S)are vital for the evaluation of groundwater resources.There are several methods to estimate the accurate aquifer parameters(i.e.hydrograph analysis,pumping test,etc.).In initial days,these parameters are projected either by means of in-situ test or execution test on aquifer well samples carried in the laboratory.The simultaneous information on the hydraulic behavior of the well(borehole)that provides on this method,the reservoir and the reservoir boundaries,are important for efficient aquifer and well data management and analysis.The most common in-situ test is pumping test performed on wells,which involves the measurement of the fall and increase of groundwater level with respect to time.The alteration in groundwater level(drawdown/recovery)is caused due to pumping of water from the well.Theis(1935)was first to propose method to evaluate aquifer parameters from the pumping test on a bore well in a confined aquifer.It is essential to know the transmissivity(T=Kb,where b is the aquifer thickness;pumping flow rate,Q=TW(dh/dl)flow through an aquifer)and storativity(confined aquifer:S=bS_s,unconfined:S=S_y),for the characterization of the aquifer parameters in an unknown area so as to predict the rate of drawdown of the groundwater table/potentiometric surface throughout the pumping test of an aquifer.The determination of aquifer's parameters is an important basis for groundwater resources evaluation,numerical simulation,development and protection as well as scientific management.For determining aquifer's parameters,pumping test is a main method.A case study shows that these techniques have been fast speed and high correctness.The results of parameter's determination are optimized so that it has important applied value for scientific research and geology engineering preparation.
基金Financial support from the Sino-Danish Center for Education and Research(SDC)the Hempel Foundation to CoaST(The Hempel Foundation Coatings Science and Technology Centre)Hempel A/S。
文摘The Hansen solubility parameters(HSP)are frequently used for solvent selection and characterization of polymers,and are directly related to the suspension behavior of pigments in solvent mixtures.The performance of currently available group contribution(GC)methods for HSP were evaluated and found to be insufficient for computer-aided product design(CAPD)of paints and coatings.A revised and,for this purpose,improved GC method is presented for estimating HSP of organic compounds,intended for organic pigments.Due to the significant limitations of GC methods,an uncertainty analysis and parameter confidence intervals are provided in order to better quantify the estimation accuracy of the proposed approach.Compared to other applicable GC methods,the prediction error is reduced significantly with average absolute errors of 0.45 MPa^(1/2),1.35 MPa^(1/2),and 1.09 MPa^(1/2) for the partial dispersion(δD),polar(δP)and hydrogen-bonding(δH)solubility parameters respectively for a database of 1106 compounds.The performance for organic pigments is comparable to the overall method performance,with higher average errors forδD and lower average errors forδP andδH.
基金the Natural Sciences and Engineering Research Council of Canada(Grant No.498020-16)Hydro-Quebec(NC525700)the Mitacs Accelerate Program(Grant Ref.IT10008)
文摘Among the methods used for evaluating the potential hydraulic erodibility of rock,the most common methods are those based on the correlation between the force of flowing water and the capacity of a rock to resist erosion,such as Annandale’s and Pells’methods.The capacity of a rock to resist erosion is evaluated based on erodibility indices that are determined from specific geomechanical parameters of a rock mass.These indices include unconfined compressive strength(UCS)of rock,rock block size,joint shear strength,a block’s shape and orientation relative to the direction of flow,joint openings,and the nature of the surface to be potentially eroded.However,it is difficult to determine the relevant geomechanical parameters for evaluating the hydraulic erodibility of rock.The assessment of eroded unlined spillways of dams has shown that the capacity of a rock to resist erosion is not accurately evaluated.Using more than 100 case studies,we develop a method to determine the relevant geomechanical parameters for evaluating the hydraulic erodibility of rock in unlined spillways.The UCS of rock is found not to be a relevant parameter for evaluating the hydraulic erodibility of rock.On the other hand,we find that the use of three-dimensional(3D)block volume measurements,instead of the block size factor used in Annandale’s method,improves the rock block size estimation.Furthermore,the parameter representing the effect of a rock block’s shape and orientation relative to the direction of flow,as considered in Pells’method,is more accurate than the parameter adopted by Annandale’s method.
文摘Particle swarm optimization （PSO）, like other evolutionary algorithms is a population-based stochastic algorithm inspired from the metaphor of social interaction in birds, insects, wasps, etc. It has been used for finding promising solutions in complex search space through the interaction of particles in a swarm. It is a well recognized fact that the performance of evolutionary algorithms to a great extent depends on the choice of appropriate strategy/operating parameters like population size, crossover rate, mutation rate, crossover operator, etc. Generally, these parameters are selected through hit and trial process, which is very unsystematic and requires rigorous experimentation. This paper proposes a systematic based on Taguchi method reasoning scheme for rapidly identifying the strategy parameters for the PSO algorithm. The Taguchi method is a robust design approach using fractional factorial design to study a large number of parameters with small number of experiments. Computer simulations have been performed on two benchmark functionsiRosenbrock function and Griewank functionito validate the approach.
基金supported by the National Special Fund for Major Research Instrument Development(2011YQ140145)111 Project (B07009)+1 种基金the National Natural Science Foundation of China(11002013)Defense Industrial Technology Development Program(A2120110001 and B2120110011)
文摘A new numerical technique named interval finite difference method is proposed for the steady-state temperature field prediction with uncertainties in both physical parameters and boundary conditions. Interval variables are used to quantitatively describe the uncertain parameters with limited information. Based on different Taylor and Neumann series, two kinds of parameter perturbation methods are presented to approximately yield the ranges of the uncertain temperature field. By comparing the results with traditional Monte Carlo simulation, a numerical example is given to demonstrate the feasibility and effectiveness of the proposed method for solving steady-state heat conduction problem with uncertain-but-bounded parameters.
基金This work was supported by the National Natural Science Foundation of China
文摘This paper discusses a kind of implicit iterative methods with some variable parameters, which are called control parameters, for solving ill-posed operator equations. The theoretical results show that the new methods always lead to optimal convergence rates and have some other important features, especially the methods can be implemented parallelly.
基金supported by the National Natural Science Foundation of China(No.52008395).
文摘During the service lives of ballasted tracks,the ballast experiences degradation,such as breakage and fragmentation,which reduces the stability of the tracks.Ballast breakage is directly related to the shapes,sizes,and stress states of particles.The key to solving this problem is to determine the breakage parameters of railway ballast.In this study,through uniaxial compression and uniaxial shear tests,the factors associated with ballast breakage for a simple stress state were obtained.Then,a refined discrete element model of railway ballast was established.The Box-Behnken method was used to design a simulation test,and a response surface method was used to obtain the optimal ballast breakage parameters.Lastly,the results of a direct shear test were compared with the simulation results to verify the correctness of the parameters and to study the shear resistance of the ballast aggregates.The results showed that for ballast particles with sizes of 22.4–63.0 mm,the characteristic stresses of compression and shear were not significantly affected by the size of the particles,but were greatly affected by their shape.In particular,long particles were more likely to break.During the direct shear test,the cumulative ballast breakage ratio gradually increased to 13.97%.The optimal breakage parameters of the ballast determined by this approach have high application value in the management of ballasted tracks.
