In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is deriv...In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.展开更多
The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction di...The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.展开更多
In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transp...In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.展开更多
The paper studies the parametric stochastic roll motion in the random waves.The differential equation of the ship parametric roll under random wave is established with considering the nonlinear damping and ship speed....The paper studies the parametric stochastic roll motion in the random waves.The differential equation of the ship parametric roll under random wave is established with considering the nonlinear damping and ship speed.Random sea surface is treated as a narrow-band stochastic process,and the stochastic parametric excitation is studied based on the effective wave theory.The nonlinear restored arm function obtained from the numerical simulation is expressed as the approximate analytic function.By using the stochastic averaging method,the differential equation of motion is transformed into Ito’s stochastic differential equation.The steady-state probability density function of roll motion is obtained,and the results are validated with the numerical simulation and model test.展开更多
This paper puts forward a complex inner product averaging method for calculating normal form of ODE. Compared with conventional averaging method, the theoretic analytical process has such simple forms as to realize co...This paper puts forward a complex inner product averaging method for calculating normal form of ODE. Compared with conventional averaging method, the theoretic analytical process has such simple forms as to realize computer program easily. Results can be applied in both autonomous and non-autonomous systems. At last, an example is resolved to verify the method.展开更多
The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equation...The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equations were expanded with respect to a snmll parameter to get the zeroth- and first-order approximate equations. By using the differenl2al quadrature method with only a few grid points, the high-accurate numerical results were obtained.展开更多
Considering the rainfall’s importance in hydrological modeling, the objective of this study was the performance comparison, in convergence terms, of techniques often used to estimate the average rainfall over an area...Considering the rainfall’s importance in hydrological modeling, the objective of this study was the performance comparison, in convergence terms, of techniques often used to estimate the average rainfall over an area: Thiessen Polygon (TP) Method;Reciprocal Distance Squared (RDS) Method;Kriging Method (KM) and Multiquadric Equations (ME) Method. The comparison was done indirectly, using GORE and BALANCE index to assess the convergence results from each method by increasing the rain gauges density in a region, through six scenarios. The Coremas/Mae D’água Watershed employed as study area, with an area of 8385 km2, is situated on Brazilian semi-arid. The results showed the TP, as RDS and ME techniques to be employed successfully to obtain the average rainfall over an area, highlighting the MEM. On the other hand, KM, using two variograms models, had an unstable behavior, pointing the prior study of data and variogram’s choice as a need to practical applying.展开更多
We consider the persistence of affine periodic solutions for perturbed affine periodic systems.Such(Q,T)-affine periodic solutions have the form x(t+T)=Qx(t)for all t∈R,where T>0 is fixed and Q is a nonsingular ma...We consider the persistence of affine periodic solutions for perturbed affine periodic systems.Such(Q,T)-affine periodic solutions have the form x(t+T)=Qx(t)for all t∈R,where T>0 is fixed and Q is a nonsingular matrix.These are a kind of spatiotemporal symmetric solutions,e.g.spiral waves.We give the averaging method for the existence of affine periodic solutions in two situations:one in which the initial values of the affine periodic solutions of the unperturbed system form a manifold,and another that does not rely on the structure of the initial values of the unperturbed system's affine periodic solutions.The transversal condition is determined using the Brouwer degree.We also present a higher order averaging method for general degenerate systems by means of the Brouwer degree and a Lyapunov-Schmidt reduction.展开更多
The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elasti...The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspon- dence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale （SOTS） computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calcu- lating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement ap- proximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.展开更多
In this paper,a stochastic second-order two-scale(SSOTS)method is proposed for predicting the non-deterministic mechanical properties of composites with random interpenetrating phase.Firstly,based on random morphology...In this paper,a stochastic second-order two-scale(SSOTS)method is proposed for predicting the non-deterministic mechanical properties of composites with random interpenetrating phase.Firstly,based on random morphology description functions(RMDF),the randomness of the material properties of the constituents as well as the correlation among these random properties are fully characterized through the topologies of the constituents.Then,by virtue of multiscale asymptotic analysis,the random effective quantities such as stiffness parameters and strength parameters along with their numerical computation formulae are derived by a SSOTS strategy combined with the Monte-Carlo method.Finally,the SSOTS method developed in this paper shows an excellent computational accuracy,and therefore present an important advance towards computationally efficient multiscale modeling frameworks considering microstructure uncertainties.展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.
A multi-objective linear programming problem is made from fuzzy linear programming problem. It is due the fact that it is used fuzzy programming method during the solution. The Multi objective linear programming probl...A multi-objective linear programming problem is made from fuzzy linear programming problem. It is due the fact that it is used fuzzy programming method during the solution. The Multi objective linear programming problem can be converted into the single objective function by various methods as Chandra Sen’s method, weighted sum method, ranking function method, statistical averaging method. In this paper, Chandra Sen’s method and statistical averaging method both are used here for making single objective function from multi-objective function. Two multi-objective programming problems are solved to verify the result. One is numerical example and the other is real life example. Then the problems are solved by ordinary simplex method and fuzzy programming method. It can be seen that fuzzy programming method gives better optimal values than the ordinary simplex method.展开更多
The energy approach is used to theoretically verify that the average acceleration method （AAM）, which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical ...The energy approach is used to theoretically verify that the average acceleration method （AAM）, which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical nonlinear damping, including the special case of velocity power type damping with a bilinear restoring force model. Based on the energy approach, the stability of the AAM is proven for SDOF structures using the mathematical features of the velocity power function and for MDOF structures by applying the virtual displacement theorem. Finally, numerical examples are given to demonstrate the accuracy of the theoretical analysis.展开更多
Numerical properties of the time integration method proposed by the first author of this paper in 2007 are the same as those of the constant average acceleration method （AAM） for linear elastic systems, except that ...Numerical properties of the time integration method proposed by the first author of this paper in 2007 are the same as those of the constant average acceleration method （AAM） for linear elastic systems, except that the capability to capture dynamic loading was not explored. It was found that there were different quadrature equations to predict the next step displacement increment. A modified quadrature equation of this method was derived so that the equation to determine the next step displacement was numerically equivalent to the equation used in the constant AAM. It was verified that the original form of this method, in general, had a better capability to capture dynamic loadings than the constant AAM. This excellent property, in addition to computational efficiency, will help to make this method competitive with general secondorder accurate integration methods.展开更多
The heat transfer of a magnetohydrodynamics nanofluid inside an annulus considering the second-order slip condition and nanoparticle migration is theoret-ically investigated. A second-order slip condition, which appro...The heat transfer of a magnetohydrodynamics nanofluid inside an annulus considering the second-order slip condition and nanoparticle migration is theoret-ically investigated. A second-order slip condition, which appropriately represents the non-equilibrium region near the interface, is prescribed rather than the no-slip condition and the linear Navier slip condition. To impose different temperature gradients, the outer wall is subjected to q2, the inner wall is subjected to q1, and q1 〉 q2. A modified two-component four-equation non-homogeneous equilibrium model is employed for the nanofiuid, which have been reduced to two-point ordinary boundary value differential equations in the consideration of the thermally and hydrodynamically fully developed flow. The homotopy analysis method （HAM） is employed to solve the equations, and the h-curves are plotted to verify the accuracy and efficiency of the solutions. Moreover, the effects of the physical factors on the flow and heat transfer are discussed in detail, and the semi-analytical relation between NUB and NBT is obtained.展开更多
Neutrons have been extensively used in many fields,such as nuclear physics,biology,geology,medical science,and national defense,owing to their unique penetration characteristics.Gamma rays are usually accompanied by t...Neutrons have been extensively used in many fields,such as nuclear physics,biology,geology,medical science,and national defense,owing to their unique penetration characteristics.Gamma rays are usually accompanied by the detection of neutrons.The capability to discriminate neutrons from gamma rays is important for evaluating plastic scintillator neutron detectors because similar pulse shapes are generated from both forms of radiation in the detection system.The pulse signals measured by plastic scintillators contain noise,which decreases the accuracy of n-y discrimination.To improve the performance of n-y discrimination,the noise of the pulse signals should be filtered before the n-y discrimination process.In this study,the influences of the Fourier transform,wavelet transform,moving-average filter,and Kalman algorithm on the charge comparison method,fractal spectrum method,and back-propagation neural network methods were studied.It was found that the Fourier transform filtering algorithm exhibits better adaptability to the charge comparison method than others,with an increasing accuracy of 6.87%compared to that without the filtering process.Meanwhile,the Kalman filter offers an improvement of 3.04%over the fractal spectrum method,and the adaptability of the moving-average filter in backpropagation neural network discrimination is better than that in other methods,with an increase in 8.48%.The Kalman filtering algorithm has a significant impact on the peak value of the pulse,reaching 4.49%,and it has an insignificant impact on the energy resolution of the spectrum measurement after discrimination.展开更多
A stochastic averaging method for predicting the response of quasi partially integrable and non-resonant Hamiltoniansystems to fractional Gaussian noise （fGla） with the Hurst index 1/2〈H〈l is proposed. The average...A stochastic averaging method for predicting the response of quasi partially integrable and non-resonant Hamiltoniansystems to fractional Gaussian noise （fGla） with the Hurst index 1/2〈H〈l is proposed. The averaged stochastic differential equa-tions （SDEs） for the first integrals of the associated Hamiltonian system are derived. The dimension of averaged SDEs is less thanthat of the original system. The stationary probability density and statistics of the original system are obtained approximately fromsolving the averaged SDEs numerically. Two systems are worked out to illustrate the proposed stochastic averaging method. It isshown that the results obtained by using the proposed stochastic averaging method and those from digital simulation of originalsystem agree well, and the computational time for the former results is less than that for the latter ones.展开更多
A stochastic averaging method of quasi integrable and resonant Hamiltonian systems under excitation of fractional Gaussian noise （fGn） with the Hurst index 1/2 〈 H 〈 1 is proposed. First, the definition and the ba...A stochastic averaging method of quasi integrable and resonant Hamiltonian systems under excitation of fractional Gaussian noise （fGn） with the Hurst index 1/2 〈 H 〈 1 is proposed. First, the definition and the basic property of fGn and related fractional Brownian motion （iBm） are briefly introduced. Then, the averaged fractional stochastic differential equations （SDEs） for the first integrals and combinations of angle variables of the associated Hamiltonian systems are derived. The stationary probability density and statistics of the original systems are then obtained approximately by simulating the averaged SDEs numerically. An example is worked out to illustrate the proposed stochastic averaging method. It is shown that the results obtained by using the proposed stochastic averaging method and those from digital simulation of original system agree well.展开更多
The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial second-order fluid was studied for modeling wet granular materials using the discrete element method. Based on the Reyn...The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial second-order fluid was studied for modeling wet granular materials using the discrete element method. Based on the Reynolds' lubrication theory, the small parameter method was introduced to approximately analyze velocity field and stress distribution between the two disks. Then a similar procedure was carried out for analyzing the normal interaction between two nearly touching, arbitrary rigid spheres to obtain the pressure distribution and the resulting squeeze force. It has been proved that the solutions can be reduced to the case of a Newtonian fluid when the non-Newtonian terms are neglected.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.11971303 and 11871330)。
文摘In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.
基金supported by the National Natural Science Foundation of China(Grant No.91130013)the Open Foundation of State Key Laboratory of HighPerformance Computing of China
文摘The energy preserving average vector field （AVF） method is applied to the coupled Schr6dinger-KdV equations. Two energy preserving schemes are constructed by using Fourier pseudospectral method in space direction discretization. In order to accelerate our simulation, the split-step technique is used. The numerical experiments show that the non-splitting scheme and splitting scheme are both effective, and have excellent long time numerical behavior. The comparisons show that the splitting scheme is faster than the non-splitting scheme, but it is not as good as the non-splitting scheme in preserving the invariants.
基金supported in part by a grant from the Research Grants Council of the Hong Kong Special Administrative Region,China (Project No.CityU 102204).
文摘In this paper,we study the electromagnetic scattering from a two dimen- sional large rectangular open cavity embedded in an infinite ground plane,which is modelled by Helmholtz equations.By introducing nonlocal transparent boundary con- ditions,the problem in the open cavity is reduced to a bounded domain problem.A hypersingular integral operator and a weakly singular integral operator are involved in the TM and TE cases,respectively.A new second-order Toeplitz type approximation and a second-order finite difference scheme are proposed for approximating the hyper- singular integral operator on the aperture and the Helmholtz in the cavity,respectively. The existence and uniqueness of the numerical solution in the TE case are established for arbitrary wavenumbers.A fast algorithm for the second-order approximation is pro- posed for solving the cavity model with layered media.Numerical results show the second-order accuracy and efficiency of the fast algorithm.More important is that the algorithm is easy to implement as a preconditioner for cavity models with more general media.
基金the State Administration of Science,Technology and Industry for National Defense of China(Grant No.B2420132001).
文摘The paper studies the parametric stochastic roll motion in the random waves.The differential equation of the ship parametric roll under random wave is established with considering the nonlinear damping and ship speed.Random sea surface is treated as a narrow-band stochastic process,and the stochastic parametric excitation is studied based on the effective wave theory.The nonlinear restored arm function obtained from the numerical simulation is expressed as the approximate analytic function.By using the stochastic averaging method,the differential equation of motion is transformed into Ito’s stochastic differential equation.The steady-state probability density function of roll motion is obtained,and the results are validated with the numerical simulation and model test.
文摘This paper puts forward a complex inner product averaging method for calculating normal form of ODE. Compared with conventional averaging method, the theoretic analytical process has such simple forms as to realize computer program easily. Results can be applied in both autonomous and non-autonomous systems. At last, an example is resolved to verify the method.
文摘The two-dimensional steady flow of an incompressible second-order viscoelastic fluid between two parallel plates was studied in terms of vorticity, the stream function and temperature equations. The governing equations were expanded with respect to a snmll parameter to get the zeroth- and first-order approximate equations. By using the differenl2al quadrature method with only a few grid points, the high-accurate numerical results were obtained.
文摘Considering the rainfall’s importance in hydrological modeling, the objective of this study was the performance comparison, in convergence terms, of techniques often used to estimate the average rainfall over an area: Thiessen Polygon (TP) Method;Reciprocal Distance Squared (RDS) Method;Kriging Method (KM) and Multiquadric Equations (ME) Method. The comparison was done indirectly, using GORE and BALANCE index to assess the convergence results from each method by increasing the rain gauges density in a region, through six scenarios. The Coremas/Mae D’água Watershed employed as study area, with an area of 8385 km2, is situated on Brazilian semi-arid. The results showed the TP, as RDS and ME techniques to be employed successfully to obtain the average rainfall over an area, highlighting the MEM. On the other hand, KM, using two variograms models, had an unstable behavior, pointing the prior study of data and variogram’s choice as a need to practical applying.
基金supported by the National Natural Science Foundation of China(1237119112071175)+4 种基金supported by the NSFC(1207117511901080)supported by the NSFC(12071175)the Fundamental Research Funds For the Central Universities(2412023YQ003)the Natural Science Foundation of Jilin Province(20200201253JC)。
文摘We consider the persistence of affine periodic solutions for perturbed affine periodic systems.Such(Q,T)-affine periodic solutions have the form x(t+T)=Qx(t)for all t∈R,where T>0 is fixed and Q is a nonsingular matrix.These are a kind of spatiotemporal symmetric solutions,e.g.spiral waves.We give the averaging method for the existence of affine periodic solutions in two situations:one in which the initial values of the affine periodic solutions of the unperturbed system form a manifold,and another that does not rely on the structure of the initial values of the unperturbed system's affine periodic solutions.The transversal condition is determined using the Brouwer degree.We also present a higher order averaging method for general degenerate systems by means of the Brouwer degree and a Lyapunov-Schmidt reduction.
基金Project supported by the National Natural Science Foundation of China(No.11471262)
文摘The correspondence principle is an important mathematical technique to compute the non-ageing linear viscoelastic problem as it allows to take advantage of the computational methods originally developed for the elastic case. However, the correspon- dence principle becomes invalid when the materials exhibit ageing. To deal with this problem, a second-order two-scale （SOTS） computational method in the time domain is presented to predict the ageing linear viscoelastic performance of composite materials with a periodic structure. First, in the time domain, the SOTS formulation for calcu- lating the effective relaxation modulus and displacement approximate solutions of the ageing viscoelastic problem is formally derived. Error estimates of the displacement ap- proximate solutions for SOTS method are then given. Numerical results obtained by the SOTS method are shown and compared with those by the finite element method in a very fine mesh. Both the analytical and numerical results show that the SOTS computational method is feasible and efficient to predict the ageing linear viscoelastic performance of composite materials with a periodic structure.
基金partially supported by China Postdoctoral Science Foundation(2018M643573)National Natural Science Foundation of Shaanxi Province(2019JQ-048)+2 种基金National Natural Science Foundation of China(51739007,61971328,11301392 and 11961009)of ChinaShanghai Peak Discipline Program for Higher Education Institutions(ClassⅠ)–Civil EngineeringFundamental Research Funds for the Central Universities(No.22120180529)。
文摘In this paper,a stochastic second-order two-scale(SSOTS)method is proposed for predicting the non-deterministic mechanical properties of composites with random interpenetrating phase.Firstly,based on random morphology description functions(RMDF),the randomness of the material properties of the constituents as well as the correlation among these random properties are fully characterized through the topologies of the constituents.Then,by virtue of multiscale asymptotic analysis,the random effective quantities such as stiffness parameters and strength parameters along with their numerical computation formulae are derived by a SSOTS strategy combined with the Monte-Carlo method.Finally,the SSOTS method developed in this paper shows an excellent computational accuracy,and therefore present an important advance towards computationally efficient multiscale modeling frameworks considering microstructure uncertainties.
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
文摘In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
文摘In this paper, we construct a uniform second-order difference scheme for a class of boundary value problems of fourth-order ordinary differential equations. Finally, a numerical example is given.
文摘A multi-objective linear programming problem is made from fuzzy linear programming problem. It is due the fact that it is used fuzzy programming method during the solution. The Multi objective linear programming problem can be converted into the single objective function by various methods as Chandra Sen’s method, weighted sum method, ranking function method, statistical averaging method. In this paper, Chandra Sen’s method and statistical averaging method both are used here for making single objective function from multi-objective function. Two multi-objective programming problems are solved to verify the result. One is numerical example and the other is real life example. Then the problems are solved by ordinary simplex method and fuzzy programming method. It can be seen that fuzzy programming method gives better optimal values than the ordinary simplex method.
基金National Natural Science Foundation of ChinaUnder Grant No. 50578047, 50338020 China Ministry ofEducation (Program for New Century Excellent Talents inUniversity) China Ministry of Science and Technology UnderGrant No.2003AA602150
文摘The energy approach is used to theoretically verify that the average acceleration method （AAM）, which is unconditionally stable for linear dynamic systems, is also unconditionally stable for structures with typical nonlinear damping, including the special case of velocity power type damping with a bilinear restoring force model. Based on the energy approach, the stability of the AAM is proven for SDOF structures using the mathematical features of the velocity power function and for MDOF structures by applying the virtual displacement theorem. Finally, numerical examples are given to demonstrate the accuracy of the theoretical analysis.
基金Science Council (NSC),Chinese Taipei Under Grant No.NSC-96-2221-E-027-030
文摘Numerical properties of the time integration method proposed by the first author of this paper in 2007 are the same as those of the constant average acceleration method （AAM） for linear elastic systems, except that the capability to capture dynamic loading was not explored. It was found that there were different quadrature equations to predict the next step displacement increment. A modified quadrature equation of this method was derived so that the equation to determine the next step displacement was numerically equivalent to the equation used in the constant AAM. It was verified that the original form of this method, in general, had a better capability to capture dynamic loadings than the constant AAM. This excellent property, in addition to computational efficiency, will help to make this method competitive with general secondorder accurate integration methods.
基金Project supported by the National Natural Science Foundation of China(Nos.51476191 and51406008)
文摘The heat transfer of a magnetohydrodynamics nanofluid inside an annulus considering the second-order slip condition and nanoparticle migration is theoret-ically investigated. A second-order slip condition, which appropriately represents the non-equilibrium region near the interface, is prescribed rather than the no-slip condition and the linear Navier slip condition. To impose different temperature gradients, the outer wall is subjected to q2, the inner wall is subjected to q1, and q1 〉 q2. A modified two-component four-equation non-homogeneous equilibrium model is employed for the nanofiuid, which have been reduced to two-point ordinary boundary value differential equations in the consideration of the thermally and hydrodynamically fully developed flow. The homotopy analysis method （HAM） is employed to solve the equations, and the h-curves are plotted to verify the accuracy and efficiency of the solutions. Moreover, the effects of the physical factors on the flow and heat transfer are discussed in detail, and the semi-analytical relation between NUB and NBT is obtained.
基金supported by the Key Natural Science Projects of the Sichuan Education Department(No.18ZA0067)the Key Science and Technology Projects of Leshan(No.19SZD117)。
文摘Neutrons have been extensively used in many fields,such as nuclear physics,biology,geology,medical science,and national defense,owing to their unique penetration characteristics.Gamma rays are usually accompanied by the detection of neutrons.The capability to discriminate neutrons from gamma rays is important for evaluating plastic scintillator neutron detectors because similar pulse shapes are generated from both forms of radiation in the detection system.The pulse signals measured by plastic scintillators contain noise,which decreases the accuracy of n-y discrimination.To improve the performance of n-y discrimination,the noise of the pulse signals should be filtered before the n-y discrimination process.In this study,the influences of the Fourier transform,wavelet transform,moving-average filter,and Kalman algorithm on the charge comparison method,fractal spectrum method,and back-propagation neural network methods were studied.It was found that the Fourier transform filtering algorithm exhibits better adaptability to the charge comparison method than others,with an increasing accuracy of 6.87%compared to that without the filtering process.Meanwhile,the Kalman filter offers an improvement of 3.04%over the fractal spectrum method,and the adaptability of the moving-average filter in backpropagation neural network discrimination is better than that in other methods,with an increase in 8.48%.The Kalman filtering algorithm has a significant impact on the peak value of the pulse,reaching 4.49%,and it has an insignificant impact on the energy resolution of the spectrum measurement after discrimination.
基金supported by the National Natural Science Foundation of China(Nos.11172259,11272279,11321202,and 11432012)
文摘A stochastic averaging method for predicting the response of quasi partially integrable and non-resonant Hamiltoniansystems to fractional Gaussian noise （fGla） with the Hurst index 1/2〈H〈l is proposed. The averaged stochastic differential equa-tions （SDEs） for the first integrals of the associated Hamiltonian system are derived. The dimension of averaged SDEs is less thanthat of the original system. The stationary probability density and statistics of the original system are obtained approximately fromsolving the averaged SDEs numerically. Two systems are worked out to illustrate the proposed stochastic averaging method. It isshown that the results obtained by using the proposed stochastic averaging method and those from digital simulation of originalsystem agree well, and the computational time for the former results is less than that for the latter ones.
基金supported by the National Natural Science Foundation of China under grants nos.:11272279,11321202 and 11432012
文摘A stochastic averaging method of quasi integrable and resonant Hamiltonian systems under excitation of fractional Gaussian noise （fGn） with the Hurst index 1/2 〈 H 〈 1 is proposed. First, the definition and the basic property of fGn and related fractional Brownian motion （iBm） are briefly introduced. Then, the averaged fractional stochastic differential equations （SDEs） for the first integrals and combinations of angle variables of the associated Hamiltonian systems are derived. The stationary probability density and statistics of the original systems are then obtained approximately by simulating the averaged SDEs numerically. An example is worked out to illustrate the proposed stochastic averaging method. It is shown that the results obtained by using the proposed stochastic averaging method and those from digital simulation of original system agree well.
文摘The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial second-order fluid was studied for modeling wet granular materials using the discrete element method. Based on the Reynolds' lubrication theory, the small parameter method was introduced to approximately analyze velocity field and stress distribution between the two disks. Then a similar procedure was carried out for analyzing the normal interaction between two nearly touching, arbitrary rigid spheres to obtain the pressure distribution and the resulting squeeze force. It has been proved that the solutions can be reduced to the case of a Newtonian fluid when the non-Newtonian terms are neglected.
