A step-by-step load was utilized to mimic the load history of the backfill column in the in-situ curing process.The inner damage of the specimen during curing and uniaxial compressive testing was monitored by electric...A step-by-step load was utilized to mimic the load history of the backfill column in the in-situ curing process.The inner damage of the specimen during curing and uniaxial compressive testing was monitored by electrical resistivity and ultrasonic equipment.Results show that:1)Uniaxial compressive strength(UCS)and elastic modulus(EM)of the samples curing under pressure are higher than those of the control samples without pressure,ranging in ratio from 0.5%to 20.2%and 7.1%to 52.3%,respectively,and are influenced by the initial loading age(ILA)and stress strength ratio(SSR).The SSR during curing should not exceed 80%.2)The earlier the ILA is,the higher the total strain becomes.The higher the SSR applies,the larger the total strain gets.The creep strain increases with the increase of SSR and can be described by Burger’s viscoelastic creep model.When SSR is less than 80%,the earlier the ILA is,the smaller the creep strain becomes after the last step-loading.3)The stability of the early age backfill column under pressure can be monitored based on the change of ultrasonic pulse velocity(UPV)and electrical resistivity.展开更多
In recent years,oil spill accidents occur frequently in the marine area of China.Finding out the spilled oil source is a key step in the relevant investigation.In this paper,a step-by-step fingerprinting identificatio...In recent years,oil spill accidents occur frequently in the marine area of China.Finding out the spilled oil source is a key step in the relevant investigation.In this paper,a step-by-step fingerprinting identification method was used in a spilled oil accident in the Bohai Sea in 2002.Advanced chemical fingerprinting and data interpretation techniques were used to characterize the chemical composition and determine the possible sources of two spilled oil samples.The original gas chromatography -flame ionization detec-tion (GC-FID) chromatogram of saturated hydrocarbons was compared.The gas chromatography-mass spectrometry (GC/MS) chromatograms of aromatic hydrocarbons terpane and sterane,n-alkane and poly-aromatic hydrocarbons (PAHs) were analyzed.The correlation analysis on diagnostic ratios was performed with Student’s t-test.It is found that the oil fingerprinting of the spilled oil (designated as sz1) from the polluted sand beach was identical with the suspected oil (designated as ky1) from a nearby crude oil refinery factory.They both showed the fingerprinting character of mixed oil.The oil fingerprinting of the spilled oil (designated as ms1) collected from the port was significantly different from oil ky1 and oil sz1 and was with a lubricating oil fingerprint character.The identification result not only gave support for the spilled oil investigation,but also served as an example for studying spilled oil accidents.展开更多
A three-dimensional(3D)step-by-step inversion strategy for audio magnetotellurics(AMT)is investigated in this study.The objective function is minimized by iteratively solving the Gauss-Newton normal equation,and the i...A three-dimensional(3D)step-by-step inversion strategy for audio magnetotellurics(AMT)is investigated in this study.The objective function is minimized by iteratively solving the Gauss-Newton normal equation,and the inversion region is discretized with unstructured tetrahedral elements.The inversion proceeds step-by-step from a coarse mesh to a fine mesh.In the inversion iteration process,a mesh is adaptively optimized according to the spatial gradient information about the model resistivity to fine delineate the boundaries of abnormal bodies.In the early stage of inversion execution,a coarse mesh is used for inversion,and the inversion stability is improved by reducing the number of inversion elements.In addition,mesh refinement is performed in the iterative inversion process.The inversion results obtained from the previous mesh are used as the reference and initial models for the next mesh iterative inversion.The step-by-step inversion strategy can ensure that the inversion is performed in the correct direction,improving the inversion stability and results gradually.Synthetic results show that the step-by-step inversion strategy with a Gauss-Newton method for 3D AMT inversion is stable and reliable,which lays a foundation for further practical 3D AMT data inversion.展开更多
We applied near-infrared(NIR)spectroscopy with chemometrics for the rapid and reagent-fee analysis of serum urea nitrogen(SUN).The modeling is based on the average effect of multiple sample partitions to achieve param...We applied near-infrared(NIR)spectroscopy with chemometrics for the rapid and reagent-fee analysis of serum urea nitrogen(SUN).The modeling is based on the average effect of multiple sample partitions to achieve parameter selection with stability.A multiparameter optimization platform with Norris derivative filter-partial least squares(Norris-PLS)was developed to select the most suitable mode(d=2,s=33,g=15).Using equidistant combination PLS(EC-PLS)with four parameters(initial wavelength I,number of wavelengths N,number of wavelength gaps G and latent variables LV),we performed wavelength screening after eliminating high-absorption wavebands.The optimal EC-PLS parameters were I=1228 nm,N=26,G=16 and LV=12.The root-mean square error(SEP),correlation coefficient(R_(p))for prediction and ratio of performance-to-deviation(RPD)for validation were 1.03 mmol L^(-1),0.992 and 7.6,respectively.We proposed the wavelength step-by-step phase-out PLS(WSP-PLS)to remove redun-dant wavelengths in the top 100 EC-PLS models with improved prediction performance.The combination of 19 wavelengths was identifed as the optimal model for SUN.The SEP,Rp and RPD in validation were 1.01 mmol L^(-1),0.992 and 7.7,respectively.The prediction effect and wavelength complexity were better than those of EC-PIS.Our results showed that NIR spectroscopy combined with the EC-PLS and WSP-PLS methods enabled the high-precision analysis ofSUN.WSP-PLS is a secondary optimization method that can further optimize any wavelength moc odel obtained through other continuous or discrete strategies to establish a simple and better model.展开更多
Objective: To evaluate the treatment outcomes of patients with pelvic ring injury by applying step-by-step external pelvic fixation and circular external fixation device.Methods: A total of 28 patients suffering from ...Objective: To evaluate the treatment outcomes of patients with pelvic ring injury by applying step-by-step external pelvic fixation and circular external fixation device.Methods: A total of 28 patients suffering from disintegrated pelvic ring injury are involved in the study. Fourteen patients(the treatment group) underwent step-by-step external pelvic fixation by applying anterior(anti-shock) and posterior modules. For the rest 14 patients(the experimental group), the osteosynthesis were conducted by means of a circular external fixation device. The long-term outcomes were evaluated in a year after the injury.Results: The residual deformity of 5(4–7) mm was observed in 10 patients(71.4%) from the experimental group. In the treatment group, the residual deformity was evident only in 4(28.6%) cases being 2.5(2–3) mm(P = 0.000 319) on the average. The functional result(according to the Majeed scale) was statistically better in the treatment group(P = 0.000 319). Nine(64.3%) and five(35.7%) patients in treatment group showed excellent and positive results, respectively. The excellent result was demonstrated by 3patients(21.4%) of the experimental group, the positive outcomes were observed in 6 cases(42.9%) and the unsatisfactory one was displayed by 1 patient(7.1%) of the same group.Conclusions: The modular approach applied is the advantage of the transosseous osteosynthesis allowing for a separate anterior(anti-shock) fixation and final posterior reposition of the pelvic ring preceded by the stabilization of vital functions. The above mentioned method gives an opportunity to increase the amount of techniques applied for the pelvic external fixation in polytrauma cases.展开更多
A method is presented that coordinates the calculation of the displacement, velocity and acceleration of structures within the time-steps of different types of step-by-step integration. The dynamic equation is solved ...A method is presented that coordinates the calculation of the displacement, velocity and acceleration of structures within the time-steps of different types of step-by-step integration. The dynamic equation is solved using an energy equation and the calculating data of the original method. The method presented is better than the original method in terms of calculating postulations and is in better conformity with the system's movement. Take the Wilson-θ method as an example. By using the coordination process, the calculation precision has been greatly im proved (reducing the errors by approximately 90% ), and the greater part of overshooting of the calculation result has been eliminated. The study suggests that the mal-coordination of the motion parameters within the time-step is the major factor that contributes to the result errors of step-by-step integration for the dynamic equation.展开更多
A new active control method was proposed, in which the analytical control law was deduced by using a step by step integral method to differential equation of motion under the condition of static error being zero. This...A new active control method was proposed, in which the analytical control law was deduced by using a step by step integral method to differential equation of motion under the condition of static error being zero. This control law is terse in mathematical expression and convenient for practical use. The simulation results demonstrate that the proposed method can provide much more remarkable peak response reduction of seismically excited structures than the classical LQR method.展开更多
In order to improve the design and research and development (R & D) efficiency of the pressure- compensating drip irrigation emitter,a step-by-step computational fluid dynamics (CFD) design method was proposed bas...In order to improve the design and research and development (R & D) efficiency of the pressure- compensating drip irrigation emitter,a step-by-step computational fluid dynamics (CFD) design method was proposed based on CFD theory combined with the finite element method. By analyzing its hydraulic performance through the step-by-step CFD method,the prediction pressure-flow curve(p-Q curve) of the pressure-compensating emitter was obtained. Then the test samples were fabricated using rapid prototype and manufacturing(RP & M) technology. The emitters' hydraulic performance experiment was carried out and the experimental p-Q curve was obtained. The step-by-step CFD design method was verified by comparing the experimental p-Q curve with the prediction values,which showed that the prediction values met the experimental results well within the normal range of the emitter's working pressure. On this basis,the effect of the emitter structure on its pressure-compensating performance was studied,which showed that the height of the pressure-compensating region had significant effects on the emitter's pressure-compensating performance. Series products of the pressure-compensating emitter could be designed by changing the region's height.展开更多
A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentration...A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentrations less than 10-6 mol·L-1 do not activate any chemical reaction. So, within each step, the time steplength Δt of the algorithm is determined from the fastest reaction rate maxR by the formula Δt = 10-6mol·L-1/max R. All the reversible elementary reactions occur simultaneously;however, by a simple book-keeping technique, the updating of species concentrations, within each step of the algorithm, is performed within each elementary reaction separately. The above proposed simple algorithm for Chemical Kinetics is applied to a simple model for hydrogen combustion with only five reversible elementary reactions (Initiation, Propagation, First and Second Branching, Termination by wall destruction) with six species (H2, O2, H, O, HO, H2O). These five reversible reactions are recommended in the literature as the most significant elementary reactions of hydrogen combustion [1] [2]. Based on the proposed here simple algorithm for Chemical Kinetics, applied to the global mechanism of proposed five reversible elementary reactions for hydrogen combustion, a simple and short computer program has been developed with only about 120 Fortran instructions. By this proposed program, the following are obtained: 1) The total species concentration of hydrogen combustion, starting from the sum of initial reactants concentrations [H2] + [O2], gradually diminishes, due to termination reaction by wall destruction, and tends to the final concentration of the product [H2O], that is to the 2/3 of its initial value, in accordance to the established overall stoichiometric reaction of hydrogen combustion 2H2 + O2 → 2H2O. 2) Time-histories for concentrations of main species H2, O2, H, H2O of hydrogen combustion, in explosion and equilibrium regions, obtained by the proposed program, are compared to corresponding ones obtained by accurate computational studies of [3]. 3) In the first step of the algorithm, the only nonzero species concentrations are those of reactants [H2], [O2]. So, the maximum reaction rate is that of the forward initiation reaction max R = Rif = kif[H2] [O2], where the rate constant kif is very slow. Thus, the first time steplength Δt1 = 10-6mol·L-1/max R results long in sec. After the first step, the sequences of all the following Δt’s are very short, in μsec. So, the first time steplength Δt1 can be considered as ignition delay time. 4) It is assumed that explosion corresponds to ignition delay time Δt1 t1 = 10 sec., can be considered as explosion limit curve. This curve is compared to the corresponding one obtained by the accurate computational studies of [2].展开更多
对码长n=2m-1纠两个错误的二进制BCH码,本文介绍了一种新的step by step译码算法。译码算法能直接确定接收的比特是否正确而不需要知道接收矢量的错误数目,也不需要知道相应的错误位置多项式。基于这种译码算法的译码器的优点在于占用...对码长n=2m-1纠两个错误的二进制BCH码,本文介绍了一种新的step by step译码算法。译码算法能直接确定接收的比特是否正确而不需要知道接收矢量的错误数目,也不需要知道相应的错误位置多项式。基于这种译码算法的译码器的优点在于占用硬件资源最少,实现简单快速。此外,这种译码方法特别适合于长码。展开更多
Multi-way principal component analysis (MPCA) had been successfully applied to monitoring the batch and semi-batch process in most chemical industry. An improved MPCA approach, step-by-step adaptive MPCA (SAMPCA), usi...Multi-way principal component analysis (MPCA) had been successfully applied to monitoring the batch and semi-batch process in most chemical industry. An improved MPCA approach, step-by-step adaptive MPCA (SAMPCA), using the process variable trajectories to monitoring the batch process is presented in this paper. It does not need to estimate or fill in the unknown part of the process variable trajectory deviation from the current time until the end. The approach is based on a MPCA method that processes the data in a sequential and adaptive manner. The adaptive rate is easily controlled through a forgetting factor that controls the weight of past data in a summation. This algorithm is used to evaluate the industrial streptomycin fermentation process data and is compared with the traditional MPCA. The results show that the method is more advantageous than MPCA, especially when monitoring multi-stage batch process where the latent vector structure can change at several points during the batch.展开更多
The vehicle routing problem(VRP)is a typical discrete combinatorial optimization problem,and many models and algorithms have been proposed to solve the VRP and its variants.Although existing approaches have contribute...The vehicle routing problem(VRP)is a typical discrete combinatorial optimization problem,and many models and algorithms have been proposed to solve the VRP and its variants.Although existing approaches have contributed significantly to the development of this field,these approaches either are limited in problem size or need manual intervention in choosing parameters.To solve these difficulties,many studies have considered learning-based optimization(LBO)algorithms to solve the VRP.This paper reviews recent advances in this field and divides relevant approaches into end-to-end approaches and step-by-step approaches.We performed a statistical analysis of the reviewed articles from various aspects and designed three experiments to evaluate the performance of four representative LBO algorithms.Finally,we conclude the applicable types of problems for different LBO algorithms and suggest directions in which researchers can improve LBO algorithms.展开更多
As all generators are distributed in different areas among large scale power systems,the cooperative manipulation of the multi-generator system cannot operate well without consideration of the distance information of ...As all generators are distributed in different areas among large scale power systems,the cooperative manipulation of the multi-generator system cannot operate well without consideration of the distance information of the generators.A distributed step-by-step finite-time consensus scheme for the heterogeneous battery energy storage system(BESS)is proposed in this paper,where the coordinated consensus can come into reality within a limited time,which is appealing for the electrical engineering community.To be concrete,at first,all BESSs are classified into several clusters according to their locations,and in each cluster,there is an active leader in charge of information receiving from outside.Then,in order to coordinate the multi BESSs,five inputs,which are function oriented,were used to achieve energy level balancing,active/reactive power sharing,and the consensus of voltage and frequency of the multi BESSs.Furthermore,the frequency and voltage restoration to the nominal values of the main grid were made possible by the introduction of a virtual leader,which is actually an external leader.Compared with the centralized methods,this control scheme is entirely distributed,and each BESS only utilizes the information of its own and its neighbors.In addition,this control is robust to the load perturbation and the plug-and-play of the communication topology.Finally,some simulation experiments are executed on the modified IEEE 57-bus system to verify the suggested scheme.展开更多
A new approach which is a direct integration method with integral model (DIM-IM) to solve dynamic governing equations is presented. The governing equations ore integrated into the integral equations. An algorithm with...A new approach which is a direct integration method with integral model (DIM-IM) to solve dynamic governing equations is presented. The governing equations ore integrated into the integral equations. An algorithm with explicit and predict-correct and self-starting and fourth-order accuracy to integrate the integral equations is given. Theoretical analysis and numerical examples shaw that DIM-IM described in this paper suitable for strong nonlinear and non-conservative system have higher accuracy than central difference, Houbolt, Newmark and Wilson-Theta methods.展开更多
The time integration method with four-order accuracy, self-starting and implicit for the diffuse chemical reaction kinetics equation or the transient instantaneous temperature filed equation was presented. The example...The time integration method with four-order accuracy, self-starting and implicit for the diffuse chemical reaction kinetics equation or the transient instantaneous temperature filed equation was presented. The examples show that both accuracy and stability are better than Runge-Kutta method with four-order. The coefficients of the equation are stored with sparse matrix pattern, so an algorithm is presented which combines a compact storage scheme with reduced computation cost. The computation of the competitive and consecutive reaction in the rotating packed bed, taken as examples, shows that the method is effective.展开更多
In the case of increasing fragmentation of wetlands, the study of the relationship between wetland landscape characteristics and total nitrogen(TN) in water is of great significance to reveal the mechanism of wetland ...In the case of increasing fragmentation of wetlands, the study of the relationship between wetland landscape characteristics and total nitrogen(TN) in water is of great significance to reveal the mechanism of wetland water purification. Taking the Naoli River(NR) wetlands in Northeast China as the research object, 10 uniformly distributed sampling sites in the study area were sampled in August 2015 to test the TN concentration and interpret the images of NR wetlands in the same period. Taking the sampling site as the control point, the whole wetlands were divided into 10 regions, and the landscape index of each region was extracted. In order to reveal whether the landscape characteristics are related to the TN concentration in the wetlands water body, the landscape index and the TN concentration in the control point water body were analyzed by correlation analysis, step-by-step elimination analysis and path analysis to reveal whether the landscape characteristics are related to the TN concentration under wetlands receiving agricultural drainages. The results showed that the correlation coefficients between four area indexes or eight shape indexes and TN concentration did not reach a significant correlation level(P > 0.05), indicating that TN removal was not only determined by a single landscape index. The path coefficient of edge density(ED) index is –0.41, indicating that wetland patch connectivity is the primary factor of TN removal, and there is no relationship between the larger patch area and the higher TN removal. The removal of TN in wetlands is restricted by the synergistic effect of landscape area and shape characteristics.展开更多
Two important extensions of a technique to perform a nonlinear error propagation analysis for an explicit pseudodynamic algorithm (Chang, 2003) are presented. One extends the stability study from a given time step t...Two important extensions of a technique to perform a nonlinear error propagation analysis for an explicit pseudodynamic algorithm (Chang, 2003) are presented. One extends the stability study from a given time step to a complete step-by-step integration procedure. It is analytically proven that ensuring stability conditions in each time step leads to a stable computation of the entire step-by-step integration procedure. The other extension shows that the nonlinear error propagation results, which are derived for a nonlinear single degree of freedom (SDOF) system, can be applied to a nonlinear multiple degree of freedom (MDOF) system. This application is dependent upon the determination of the natural frequencies of the system in each time step, since all the numerical properties and error propagation properties in the time step are closely related to these frequencies. The results are derived from the step degree of nonlinearity. An instantaneous degree of nonlinearity is introduced to replace the step degree of nonlinearity and is shown to be easier to use in practice. The extensions can be also applied to the results derived from a SDOF system based on the instantaneous degree of nonlinearity, and hence a time step might be appropriately chosen to perform a pseudodynamic test prior to testing.展开更多
In this paper, to begin with. the nonlinear differential equations of a truncaled shallow spherical shell with variable thickness under uniformal distributed load are linearized by step-by-step loading method. The lin...In this paper, to begin with. the nonlinear differential equations of a truncaled shallow spherical shell with variable thickness under uniformal distributed load are linearized by step-by-step loading method. The linear differential equations can be solved by spline collocanon method. Critical loads have been obtained accordingly.展开更多
A simplified model is proposed for an easy understanding of the coarse-grained technique and for achieving a first approximation to the behavior of gases. A mole of a gas substance, within a cubic container, is repres...A simplified model is proposed for an easy understanding of the coarse-grained technique and for achieving a first approximation to the behavior of gases. A mole of a gas substance, within a cubic container, is represented by six particles symmetrically moving. The impacts of particles on container walls, the inter-particle collisions, as well as the volume of particles and the inter-particle attractive forces, obeying a Lennard-Jones curve, are taken into account. Thanks to the symmetry, the problem is reduced to the nonlinear dynamic analysis of a SDOF oscillator, which is numerically solved by a step-by-step time integration algorithm. Five applications of proposed model, on Carbon Dioxide, are presented: 1) Ideal gas in STP conditions. 2) Real gas in STP conditions. 3) Condensation for small molar volume. 4) Critical point. 5) Iso-kinetic energy curves and iso-therms in the critical point region. Results of the proposed model are compared with test data and results of the Van der Waals model for real gases.展开更多
基金Project(51974192)supported by the National Natural Science Foundation of ChinaProject(201803D31044)supported by the Program for Key Research Project of Shanxi Province in the Field of Social Development,ChinaProject(201801D121092)supported by the Applied Basic Research Project of Shanxi Province,China。
文摘A step-by-step load was utilized to mimic the load history of the backfill column in the in-situ curing process.The inner damage of the specimen during curing and uniaxial compressive testing was monitored by electrical resistivity and ultrasonic equipment.Results show that:1)Uniaxial compressive strength(UCS)and elastic modulus(EM)of the samples curing under pressure are higher than those of the control samples without pressure,ranging in ratio from 0.5%to 20.2%and 7.1%to 52.3%,respectively,and are influenced by the initial loading age(ILA)and stress strength ratio(SSR).The SSR during curing should not exceed 80%.2)The earlier the ILA is,the higher the total strain becomes.The higher the SSR applies,the larger the total strain gets.The creep strain increases with the increase of SSR and can be described by Burger’s viscoelastic creep model.When SSR is less than 80%,the earlier the ILA is,the smaller the creep strain becomes after the last step-loading.3)The stability of the early age backfill column under pressure can be monitored based on the change of ultrasonic pulse velocity(UPV)and electrical resistivity.
文摘In recent years,oil spill accidents occur frequently in the marine area of China.Finding out the spilled oil source is a key step in the relevant investigation.In this paper,a step-by-step fingerprinting identification method was used in a spilled oil accident in the Bohai Sea in 2002.Advanced chemical fingerprinting and data interpretation techniques were used to characterize the chemical composition and determine the possible sources of two spilled oil samples.The original gas chromatography -flame ionization detec-tion (GC-FID) chromatogram of saturated hydrocarbons was compared.The gas chromatography-mass spectrometry (GC/MS) chromatograms of aromatic hydrocarbons terpane and sterane,n-alkane and poly-aromatic hydrocarbons (PAHs) were analyzed.The correlation analysis on diagnostic ratios was performed with Student’s t-test.It is found that the oil fingerprinting of the spilled oil (designated as sz1) from the polluted sand beach was identical with the suspected oil (designated as ky1) from a nearby crude oil refinery factory.They both showed the fingerprinting character of mixed oil.The oil fingerprinting of the spilled oil (designated as ms1) collected from the port was significantly different from oil ky1 and oil sz1 and was with a lubricating oil fingerprint character.The identification result not only gave support for the spilled oil investigation,but also served as an example for studying spilled oil accidents.
文摘A three-dimensional(3D)step-by-step inversion strategy for audio magnetotellurics(AMT)is investigated in this study.The objective function is minimized by iteratively solving the Gauss-Newton normal equation,and the inversion region is discretized with unstructured tetrahedral elements.The inversion proceeds step-by-step from a coarse mesh to a fine mesh.In the inversion iteration process,a mesh is adaptively optimized according to the spatial gradient information about the model resistivity to fine delineate the boundaries of abnormal bodies.In the early stage of inversion execution,a coarse mesh is used for inversion,and the inversion stability is improved by reducing the number of inversion elements.In addition,mesh refinement is performed in the iterative inversion process.The inversion results obtained from the previous mesh are used as the reference and initial models for the next mesh iterative inversion.The step-by-step inversion strategy can ensure that the inversion is performed in the correct direction,improving the inversion stability and results gradually.Synthetic results show that the step-by-step inversion strategy with a Gauss-Newton method for 3D AMT inversion is stable and reliable,which lays a foundation for further practical 3D AMT data inversion.
基金supported by the Science and Technology Project of Guangdong Province of China(Nos.2014A020213016,2014A020212445)the University-enterprise Joint Research Project"Intelligent detection network technology joint research centre"(No.40115031).
文摘We applied near-infrared(NIR)spectroscopy with chemometrics for the rapid and reagent-fee analysis of serum urea nitrogen(SUN).The modeling is based on the average effect of multiple sample partitions to achieve parameter selection with stability.A multiparameter optimization platform with Norris derivative filter-partial least squares(Norris-PLS)was developed to select the most suitable mode(d=2,s=33,g=15).Using equidistant combination PLS(EC-PLS)with four parameters(initial wavelength I,number of wavelengths N,number of wavelength gaps G and latent variables LV),we performed wavelength screening after eliminating high-absorption wavebands.The optimal EC-PLS parameters were I=1228 nm,N=26,G=16 and LV=12.The root-mean square error(SEP),correlation coefficient(R_(p))for prediction and ratio of performance-to-deviation(RPD)for validation were 1.03 mmol L^(-1),0.992 and 7.6,respectively.We proposed the wavelength step-by-step phase-out PLS(WSP-PLS)to remove redun-dant wavelengths in the top 100 EC-PLS models with improved prediction performance.The combination of 19 wavelengths was identifed as the optimal model for SUN.The SEP,Rp and RPD in validation were 1.01 mmol L^(-1),0.992 and 7.7,respectively.The prediction effect and wavelength complexity were better than those of EC-PIS.Our results showed that NIR spectroscopy combined with the EC-PLS and WSP-PLS methods enabled the high-precision analysis ofSUN.WSP-PLS is a secondary optimization method that can further optimize any wavelength moc odel obtained through other continuous or discrete strategies to establish a simple and better model.
文摘Objective: To evaluate the treatment outcomes of patients with pelvic ring injury by applying step-by-step external pelvic fixation and circular external fixation device.Methods: A total of 28 patients suffering from disintegrated pelvic ring injury are involved in the study. Fourteen patients(the treatment group) underwent step-by-step external pelvic fixation by applying anterior(anti-shock) and posterior modules. For the rest 14 patients(the experimental group), the osteosynthesis were conducted by means of a circular external fixation device. The long-term outcomes were evaluated in a year after the injury.Results: The residual deformity of 5(4–7) mm was observed in 10 patients(71.4%) from the experimental group. In the treatment group, the residual deformity was evident only in 4(28.6%) cases being 2.5(2–3) mm(P = 0.000 319) on the average. The functional result(according to the Majeed scale) was statistically better in the treatment group(P = 0.000 319). Nine(64.3%) and five(35.7%) patients in treatment group showed excellent and positive results, respectively. The excellent result was demonstrated by 3patients(21.4%) of the experimental group, the positive outcomes were observed in 6 cases(42.9%) and the unsatisfactory one was displayed by 1 patient(7.1%) of the same group.Conclusions: The modular approach applied is the advantage of the transosseous osteosynthesis allowing for a separate anterior(anti-shock) fixation and final posterior reposition of the pelvic ring preceded by the stabilization of vital functions. The above mentioned method gives an opportunity to increase the amount of techniques applied for the pelvic external fixation in polytrauma cases.
文摘A method is presented that coordinates the calculation of the displacement, velocity and acceleration of structures within the time-steps of different types of step-by-step integration. The dynamic equation is solved using an energy equation and the calculating data of the original method. The method presented is better than the original method in terms of calculating postulations and is in better conformity with the system's movement. Take the Wilson-θ method as an example. By using the coordination process, the calculation precision has been greatly im proved (reducing the errors by approximately 90% ), and the greater part of overshooting of the calculation result has been eliminated. The study suggests that the mal-coordination of the motion parameters within the time-step is the major factor that contributes to the result errors of step-by-step integration for the dynamic equation.
文摘A new active control method was proposed, in which the analytical control law was deduced by using a step by step integral method to differential equation of motion under the condition of static error being zero. This control law is terse in mathematical expression and convenient for practical use. The simulation results demonstrate that the proposed method can provide much more remarkable peak response reduction of seismically excited structures than the classical LQR method.
基金The National Natural Science Fund(No.50975227)The National High-tech R & D Program("863"Program)(No.2011AA100507-04)
文摘In order to improve the design and research and development (R & D) efficiency of the pressure- compensating drip irrigation emitter,a step-by-step computational fluid dynamics (CFD) design method was proposed based on CFD theory combined with the finite element method. By analyzing its hydraulic performance through the step-by-step CFD method,the prediction pressure-flow curve(p-Q curve) of the pressure-compensating emitter was obtained. Then the test samples were fabricated using rapid prototype and manufacturing(RP & M) technology. The emitters' hydraulic performance experiment was carried out and the experimental p-Q curve was obtained. The step-by-step CFD design method was verified by comparing the experimental p-Q curve with the prediction values,which showed that the prediction values met the experimental results well within the normal range of the emitter's working pressure. On this basis,the effect of the emitter structure on its pressure-compensating performance was studied,which showed that the height of the pressure-compensating region had significant effects on the emitter's pressure-compensating performance. Series products of the pressure-compensating emitter could be designed by changing the region's height.
文摘A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentrations less than 10-6 mol·L-1 do not activate any chemical reaction. So, within each step, the time steplength Δt of the algorithm is determined from the fastest reaction rate maxR by the formula Δt = 10-6mol·L-1/max R. All the reversible elementary reactions occur simultaneously;however, by a simple book-keeping technique, the updating of species concentrations, within each step of the algorithm, is performed within each elementary reaction separately. The above proposed simple algorithm for Chemical Kinetics is applied to a simple model for hydrogen combustion with only five reversible elementary reactions (Initiation, Propagation, First and Second Branching, Termination by wall destruction) with six species (H2, O2, H, O, HO, H2O). These five reversible reactions are recommended in the literature as the most significant elementary reactions of hydrogen combustion [1] [2]. Based on the proposed here simple algorithm for Chemical Kinetics, applied to the global mechanism of proposed five reversible elementary reactions for hydrogen combustion, a simple and short computer program has been developed with only about 120 Fortran instructions. By this proposed program, the following are obtained: 1) The total species concentration of hydrogen combustion, starting from the sum of initial reactants concentrations [H2] + [O2], gradually diminishes, due to termination reaction by wall destruction, and tends to the final concentration of the product [H2O], that is to the 2/3 of its initial value, in accordance to the established overall stoichiometric reaction of hydrogen combustion 2H2 + O2 → 2H2O. 2) Time-histories for concentrations of main species H2, O2, H, H2O of hydrogen combustion, in explosion and equilibrium regions, obtained by the proposed program, are compared to corresponding ones obtained by accurate computational studies of [3]. 3) In the first step of the algorithm, the only nonzero species concentrations are those of reactants [H2], [O2]. So, the maximum reaction rate is that of the forward initiation reaction max R = Rif = kif[H2] [O2], where the rate constant kif is very slow. Thus, the first time steplength Δt1 = 10-6mol·L-1/max R results long in sec. After the first step, the sequences of all the following Δt’s are very short, in μsec. So, the first time steplength Δt1 can be considered as ignition delay time. 4) It is assumed that explosion corresponds to ignition delay time Δt1 t1 = 10 sec., can be considered as explosion limit curve. This curve is compared to the corresponding one obtained by the accurate computational studies of [2].
文摘对码长n=2m-1纠两个错误的二进制BCH码,本文介绍了一种新的step by step译码算法。译码算法能直接确定接收的比特是否正确而不需要知道接收矢量的错误数目,也不需要知道相应的错误位置多项式。基于这种译码算法的译码器的优点在于占用硬件资源最少,实现简单快速。此外,这种译码方法特别适合于长码。
基金Supported by the National High-tech Program of China (No. 2001 AA413110).
文摘Multi-way principal component analysis (MPCA) had been successfully applied to monitoring the batch and semi-batch process in most chemical industry. An improved MPCA approach, step-by-step adaptive MPCA (SAMPCA), using the process variable trajectories to monitoring the batch process is presented in this paper. It does not need to estimate or fill in the unknown part of the process variable trajectory deviation from the current time until the end. The approach is based on a MPCA method that processes the data in a sequential and adaptive manner. The adaptive rate is easily controlled through a forgetting factor that controls the weight of past data in a summation. This algorithm is used to evaluate the industrial streptomycin fermentation process data and is compared with the traditional MPCA. The results show that the method is more advantageous than MPCA, especially when monitoring multi-stage batch process where the latent vector structure can change at several points during the batch.
文摘The vehicle routing problem(VRP)is a typical discrete combinatorial optimization problem,and many models and algorithms have been proposed to solve the VRP and its variants.Although existing approaches have contributed significantly to the development of this field,these approaches either are limited in problem size or need manual intervention in choosing parameters.To solve these difficulties,many studies have considered learning-based optimization(LBO)algorithms to solve the VRP.This paper reviews recent advances in this field and divides relevant approaches into end-to-end approaches and step-by-step approaches.We performed a statistical analysis of the reviewed articles from various aspects and designed three experiments to evaluate the performance of four representative LBO algorithms.Finally,we conclude the applicable types of problems for different LBO algorithms and suggest directions in which researchers can improve LBO algorithms.
基金supported in part by the National Natural Science Foundation of China(61340041,61374079 and 61903126)the Natural Science Foundation of Henan Province(182300410112).
文摘As all generators are distributed in different areas among large scale power systems,the cooperative manipulation of the multi-generator system cannot operate well without consideration of the distance information of the generators.A distributed step-by-step finite-time consensus scheme for the heterogeneous battery energy storage system(BESS)is proposed in this paper,where the coordinated consensus can come into reality within a limited time,which is appealing for the electrical engineering community.To be concrete,at first,all BESSs are classified into several clusters according to their locations,and in each cluster,there is an active leader in charge of information receiving from outside.Then,in order to coordinate the multi BESSs,five inputs,which are function oriented,were used to achieve energy level balancing,active/reactive power sharing,and the consensus of voltage and frequency of the multi BESSs.Furthermore,the frequency and voltage restoration to the nominal values of the main grid were made possible by the introduction of a virtual leader,which is actually an external leader.Compared with the centralized methods,this control scheme is entirely distributed,and each BESS only utilizes the information of its own and its neighbors.In addition,this control is robust to the load perturbation and the plug-and-play of the communication topology.Finally,some simulation experiments are executed on the modified IEEE 57-bus system to verify the suggested scheme.
文摘A new approach which is a direct integration method with integral model (DIM-IM) to solve dynamic governing equations is presented. The governing equations ore integrated into the integral equations. An algorithm with explicit and predict-correct and self-starting and fourth-order accuracy to integrate the integral equations is given. Theoretical analysis and numerical examples shaw that DIM-IM described in this paper suitable for strong nonlinear and non-conservative system have higher accuracy than central difference, Houbolt, Newmark and Wilson-Theta methods.
文摘The time integration method with four-order accuracy, self-starting and implicit for the diffuse chemical reaction kinetics equation or the transient instantaneous temperature filed equation was presented. The examples show that both accuracy and stability are better than Runge-Kutta method with four-order. The coefficients of the equation are stored with sparse matrix pattern, so an algorithm is presented which combines a compact storage scheme with reduced computation cost. The computation of the competitive and consecutive reaction in the rotating packed bed, taken as examples, shows that the method is effective.
基金Under the auspices of National Natural Science Foundations of China(No.41620104005,31500307,41601263,41771120)Technology Development Program of Jilin Province(No.20180101082JC,20180520085JC,20190201256JC,20190201018JC)+1 种基金Natural Science Foundation of Changchun Normal University(No.2016-009)Northeast Institute of Geography and Agroecology,Chinese Academy of Sciences(No.IGA-135-05)
文摘In the case of increasing fragmentation of wetlands, the study of the relationship between wetland landscape characteristics and total nitrogen(TN) in water is of great significance to reveal the mechanism of wetland water purification. Taking the Naoli River(NR) wetlands in Northeast China as the research object, 10 uniformly distributed sampling sites in the study area were sampled in August 2015 to test the TN concentration and interpret the images of NR wetlands in the same period. Taking the sampling site as the control point, the whole wetlands were divided into 10 regions, and the landscape index of each region was extracted. In order to reveal whether the landscape characteristics are related to the TN concentration in the wetlands water body, the landscape index and the TN concentration in the control point water body were analyzed by correlation analysis, step-by-step elimination analysis and path analysis to reveal whether the landscape characteristics are related to the TN concentration under wetlands receiving agricultural drainages. The results showed that the correlation coefficients between four area indexes or eight shape indexes and TN concentration did not reach a significant correlation level(P > 0.05), indicating that TN removal was not only determined by a single landscape index. The path coefficient of edge density(ED) index is –0.41, indicating that wetland patch connectivity is the primary factor of TN removal, and there is no relationship between the larger patch area and the higher TN removal. The removal of TN in wetlands is restricted by the synergistic effect of landscape area and shape characteristics.
基金NSC, Chinese Taipei Under Grant No. NSC-97-2221-E-027-036-MY2
文摘Two important extensions of a technique to perform a nonlinear error propagation analysis for an explicit pseudodynamic algorithm (Chang, 2003) are presented. One extends the stability study from a given time step to a complete step-by-step integration procedure. It is analytically proven that ensuring stability conditions in each time step leads to a stable computation of the entire step-by-step integration procedure. The other extension shows that the nonlinear error propagation results, which are derived for a nonlinear single degree of freedom (SDOF) system, can be applied to a nonlinear multiple degree of freedom (MDOF) system. This application is dependent upon the determination of the natural frequencies of the system in each time step, since all the numerical properties and error propagation properties in the time step are closely related to these frequencies. The results are derived from the step degree of nonlinearity. An instantaneous degree of nonlinearity is introduced to replace the step degree of nonlinearity and is shown to be easier to use in practice. The extensions can be also applied to the results derived from a SDOF system based on the instantaneous degree of nonlinearity, and hence a time step might be appropriately chosen to perform a pseudodynamic test prior to testing.
文摘In this paper, to begin with. the nonlinear differential equations of a truncaled shallow spherical shell with variable thickness under uniformal distributed load are linearized by step-by-step loading method. The linear differential equations can be solved by spline collocanon method. Critical loads have been obtained accordingly.
文摘A simplified model is proposed for an easy understanding of the coarse-grained technique and for achieving a first approximation to the behavior of gases. A mole of a gas substance, within a cubic container, is represented by six particles symmetrically moving. The impacts of particles on container walls, the inter-particle collisions, as well as the volume of particles and the inter-particle attractive forces, obeying a Lennard-Jones curve, are taken into account. Thanks to the symmetry, the problem is reduced to the nonlinear dynamic analysis of a SDOF oscillator, which is numerically solved by a step-by-step time integration algorithm. Five applications of proposed model, on Carbon Dioxide, are presented: 1) Ideal gas in STP conditions. 2) Real gas in STP conditions. 3) Condensation for small molar volume. 4) Critical point. 5) Iso-kinetic energy curves and iso-therms in the critical point region. Results of the proposed model are compared with test data and results of the Van der Waals model for real gases.