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.展开更多
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 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.展开更多
In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introduc...In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introducing only one additional evaluation of the function in each step.Furthermore,some new efficient methods with a higher-order of convergence are obtained by using only a single matrix inversion in each iteration.Analyses of convergence properties and computational efficiency of these new methods are made and testified by several numerical problems.By comparison,the new schemes are more efficient than the corresponding existing ones,particularly for large problem sizes.展开更多
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 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,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.展开更多
The current theory in NF EN 1995-1-1/NA of Eurocode 5, which is based on maximum deflection, has been investigated on softwoods. Therefore, this theory is not adapted for slender glulam beam columns made of tropical h...The current theory in NF EN 1995-1-1/NA of Eurocode 5, which is based on maximum deflection, has been investigated on softwoods. Therefore, this theory is not adapted for slender glulam beam columns made of tropical hardwood species from the Congo Basin. This maximum deflection is caused by a set of loads applied to the structure. However, Eurocode 5 doesn’t provide how to predict this deflection in case of long-term load for such structures. This can be done by studying load-displacement (P-Δ) behaviour of these structures while taking into account second order effects. To reach this goal, a nonlinear analysis has been performed on a three-dimensional beam column embedded on both ends. Since conducting experimental investigations on large span structural products is time-consuming and expensive especially in developing countries, a numerical model has been implemented using the Newton-Raphson method to predict load-displacement (P-Δ) curve on a slender glulam beam column made of tropical hardwood species. On one hand, the beam has been analyzed without wood connection. On the other hand, the beam has been analyzed with a bolted wood connection and a slotted-in steel plate. The load cases considered include self-weight and a uniformly applied long-term load. Combinations of serviceability limit states (SLS) and ultimate limit states (ULS) have also been considered, among other factors. A finite-element software RFEM 5 has been used to implement the model. The results showed that the use of steel can reduce displacement by 20.96%. Additionally, compared to the maximum deflection provided by Eurocode 5 for softwoods, hardwoods can exhibit an increasing rate of 85.63%. By harnessing the plastic resistance of steel, the bending resistance of wood can be increased by 32.94%.展开更多
The Stirling engine,as a closed-cycle power machine,exhibits excellent emission characteristics and broad energy adaptability.Second-order analysis methods are extensively used during the foundational design and therm...The Stirling engine,as a closed-cycle power machine,exhibits excellent emission characteristics and broad energy adaptability.Second-order analysis methods are extensively used during the foundational design and thermodynamic examination of Stirling engines,owing to their commendable model precision and remarkable efficiency.To scrutinize the effect of Stirling engine design parameters on the cyclical work output and efficiency,this study formulates a series of differential equations for the Stirling cycle by employing second-order analysis methods,subsequently augmenting the predictive accuracy by integrating considerations of loss mechanisms.In addition,an iterative method for the convergence of the average pressure was introduced.The predictive capability of the established model was validated using GPU-3 and RE-1000 experimental data.According to the model,parameters such as the operational fluid,porosity of the regenerator,and diameter of the wire mesh and their influence on the resulting work output and cyclic efficiency of the Stirling engine were analyzed,thereby facilitating a broader understanding of the engine's functional characteristics.These findings suggest that hydrogen,owing to its lower dynamic viscosity coefficient,can provide superior output power.The loss due to flow resistance tends to increase with the rotational speed.Additionally,under conditions of elevated rotational speed,the loss from flow resistance declines in cases of increased porosity,and the enhancement of the porosity to diminish flow resistance losses can boost both the output work and the cyclic efficiency of the engine.As the porosity increased further,the hydraulic diameter and dead volume in the regenerator continued to expand,causing the pressure drop within the engine to become the dominant factor in the gradual reduction of output power.Furthermore,extending the length of the regenerator results in a decrease in the output work,although the thermal cycle efficiency initially increases before eventually decreasing.Based on these insights,this study pursues the optimal designs for Stirling engines.展开更多
Current damage detection methods based on model updating and sensitivity Jacobian matrixes show a low convergence ratio and computational efficiency for online calculations.The aim of this paper is to construct a real...Current damage detection methods based on model updating and sensitivity Jacobian matrixes show a low convergence ratio and computational efficiency for online calculations.The aim of this paper is to construct a real-time automated damage detection method by developing a theory-assisted adaptive mutiagent twin delayed deep deterministic(TA2-MATD3)policy gradient algorithm.First,the theoretical framework of reinforcement-learning-driven damage detection is established.To address the disadvantages of traditional mutiagent twin delayed deep deterministic(MATD3)method,the theory-assisted mechanism and the adaptive experience playback mechanism are introduced.Moreover,a historical residential house built in 1889 was taken as an example,using its 12-month structural health monitoring data.TA2-MATD3 was compared with existing damage detection methods in terms of the convergence ratio,online computing efficiency,and damage detection accuracy.The results show that the computational efficiency of TA2-MATD3 is approximately 117–160 times that of the traditional methods.The convergence ratio of damage detection on the training set is approximately 97%,and that on the test set is in the range of 86.2%–91.9%.In addition,the main apparent damages found in the field survey were identified by TA2-MATD3.The results indicate that the proposed method can significantly improve the online computing efficiency and damage detection accuracy.This research can provide novel perspectives for the use of reinforcement learning methods to conduct damage detection in online structural health monitoring.展开更多
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.
Iterative methods based on finite element simulation are effective approaches to design mold shape to compensate springback in sheet metal forming. However, convergence rate of iterative methods is difficult to improv...Iterative methods based on finite element simulation are effective approaches to design mold shape to compensate springback in sheet metal forming. However, convergence rate of iterative methods is difficult to improve greatly. To increase the springback compensate speed of designing age forming mold, process of calculating springback for a certain mold with finite element method is analyzed. Springback compensation is abstracted as finding a solution for a set of nonlinear functions and a springback compensation algorithm is presented on the basis of quasi Newton method. The accuracy of algorithm is verified by developing an ABAQUS secondary development program with MATLAB. Three rectangular integrated panels of dimensions 710 mmx750 mm integrated panels with intersected ribs of 10 mm are selected to perform case studies. The algorithm is used to compute mold contours for the panels with cylinder, sphere and saddle contours respectively and it takes 57%, 22% and 33% iterations as compared to that of displacement adjustment （DA） method. At the end of iterations, maximum deviations on the three panels are 0.618 4 mm, 0.624 1 mm and 0.342 0 mm that are smaller than the deviations determined by DA method （0.740 8 mm, 0.740 8 mm and 0.713 7 mm respectively）. In following experimental verification, mold contour for another integrated panel with 400 ram^380 mm size is designed by the algorithm. Then the panel is age formed in an autoclave and measured by a three dimensional digital measurement devise. Deviation between measuring results and the panel＇s design contour is less than 1 mm. Finally, the iterations with different mesh sizes （40 mm, 35 mm, 30 mm, 25 mm, 20 mm） in finite element models are compared and found no considerable difference. Another possible compensation method, Broyden-Fletcher-Shanmo method, is also presented based on the solving nonlinear fimctions idea. The Broyden-Fletcher-Shanmo method is employed to compute mold contour for the second panel. It only takes 50% iterations compared to that of DA. The proposed method can serve a faster mold contour compensation method for sheet metal forming.展开更多
A fourth-order convergence method of solving roots for nonlinear equation, which is a variant of Newton＇s method given. Its convergence properties is proved. It is at least fourth-order convergence near simple roots ...A fourth-order convergence method of solving roots for nonlinear equation, which is a variant of Newton＇s method given. Its convergence properties is proved. It is at least fourth-order convergence near simple roots and one order convergence near multiple roots. In the end, numerical tests are given and compared with other known Newton and Newton-type methods. The results show that the proposed method has some more advantages than others. It enriches the methods to find the roots of non-linear equations and it is important in both theory and application.展开更多
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.展开更多
The box constrained variational inequality problem can be reformulated as a nonsmooth equation by using median operator.In this paper,we present a smoothing Newton method for solving the box constrained variational in...The box constrained variational inequality problem can be reformulated as a nonsmooth equation by using median operator.In this paper,we present a smoothing Newton method for solving the box constrained variational inequality problem based on a new smoothing approximation function.The proposed algorithm is proved to be well defined and convergent globally under weaker conditions.展开更多
基金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.
基金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.
文摘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.
基金Supported by the National Natural Science Foundation of China(12061048)NSF of Jiangxi Province(20232BAB201026,20232BAB201018)。
文摘In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introducing only one additional evaluation of the function in each step.Furthermore,some new efficient methods with a higher-order of convergence are obtained by using only a single matrix inversion in each iteration.Analyses of convergence properties and computational efficiency of these new methods are made and testified by several numerical problems.By comparison,the new schemes are more efficient than the corresponding existing ones,particularly for large problem sizes.
基金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 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.
基金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.
文摘The current theory in NF EN 1995-1-1/NA of Eurocode 5, which is based on maximum deflection, has been investigated on softwoods. Therefore, this theory is not adapted for slender glulam beam columns made of tropical hardwood species from the Congo Basin. This maximum deflection is caused by a set of loads applied to the structure. However, Eurocode 5 doesn’t provide how to predict this deflection in case of long-term load for such structures. This can be done by studying load-displacement (P-Δ) behaviour of these structures while taking into account second order effects. To reach this goal, a nonlinear analysis has been performed on a three-dimensional beam column embedded on both ends. Since conducting experimental investigations on large span structural products is time-consuming and expensive especially in developing countries, a numerical model has been implemented using the Newton-Raphson method to predict load-displacement (P-Δ) curve on a slender glulam beam column made of tropical hardwood species. On one hand, the beam has been analyzed without wood connection. On the other hand, the beam has been analyzed with a bolted wood connection and a slotted-in steel plate. The load cases considered include self-weight and a uniformly applied long-term load. Combinations of serviceability limit states (SLS) and ultimate limit states (ULS) have also been considered, among other factors. A finite-element software RFEM 5 has been used to implement the model. The results showed that the use of steel can reduce displacement by 20.96%. Additionally, compared to the maximum deflection provided by Eurocode 5 for softwoods, hardwoods can exhibit an increasing rate of 85.63%. By harnessing the plastic resistance of steel, the bending resistance of wood can be increased by 32.94%.
基金supported by Sichuan Science and Technology Program(No.24NSFSC4579)National Natural Science Foundation of China(No.12305193)+2 种基金Sichuan Science and Technology Program(No.23NSFSC6149)National Natural Science Foundation of China(No.12305194)Technology on Reactor System Design Technology Laboratory Stable support Funding(No.2023_JCJQ_LB_003).
文摘The Stirling engine,as a closed-cycle power machine,exhibits excellent emission characteristics and broad energy adaptability.Second-order analysis methods are extensively used during the foundational design and thermodynamic examination of Stirling engines,owing to their commendable model precision and remarkable efficiency.To scrutinize the effect of Stirling engine design parameters on the cyclical work output and efficiency,this study formulates a series of differential equations for the Stirling cycle by employing second-order analysis methods,subsequently augmenting the predictive accuracy by integrating considerations of loss mechanisms.In addition,an iterative method for the convergence of the average pressure was introduced.The predictive capability of the established model was validated using GPU-3 and RE-1000 experimental data.According to the model,parameters such as the operational fluid,porosity of the regenerator,and diameter of the wire mesh and their influence on the resulting work output and cyclic efficiency of the Stirling engine were analyzed,thereby facilitating a broader understanding of the engine's functional characteristics.These findings suggest that hydrogen,owing to its lower dynamic viscosity coefficient,can provide superior output power.The loss due to flow resistance tends to increase with the rotational speed.Additionally,under conditions of elevated rotational speed,the loss from flow resistance declines in cases of increased porosity,and the enhancement of the porosity to diminish flow resistance losses can boost both the output work and the cyclic efficiency of the engine.As the porosity increased further,the hydraulic diameter and dead volume in the regenerator continued to expand,causing the pressure drop within the engine to become the dominant factor in the gradual reduction of output power.Furthermore,extending the length of the regenerator results in a decrease in the output work,although the thermal cycle efficiency initially increases before eventually decreasing.Based on these insights,this study pursues the optimal designs for Stirling engines.
基金supported by National Key Research and Development Program of China(2023YFF0906100)National Natural Science Foundation of China(52408008)Key Research and Development Program of Jiangsu Province(BE2022833).
文摘Current damage detection methods based on model updating and sensitivity Jacobian matrixes show a low convergence ratio and computational efficiency for online calculations.The aim of this paper is to construct a real-time automated damage detection method by developing a theory-assisted adaptive mutiagent twin delayed deep deterministic(TA2-MATD3)policy gradient algorithm.First,the theoretical framework of reinforcement-learning-driven damage detection is established.To address the disadvantages of traditional mutiagent twin delayed deep deterministic(MATD3)method,the theory-assisted mechanism and the adaptive experience playback mechanism are introduced.Moreover,a historical residential house built in 1889 was taken as an example,using its 12-month structural health monitoring data.TA2-MATD3 was compared with existing damage detection methods in terms of the convergence ratio,online computing efficiency,and damage detection accuracy.The results show that the computational efficiency of TA2-MATD3 is approximately 117–160 times that of the traditional methods.The convergence ratio of damage detection on the training set is approximately 97%,and that on the test set is in the range of 86.2%–91.9%.In addition,the main apparent damages found in the field survey were identified by TA2-MATD3.The results indicate that the proposed method can significantly improve the online computing efficiency and damage detection accuracy.This research can provide novel perspectives for the use of reinforcement learning methods to conduct damage detection in online structural health monitoring.
文摘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.
文摘Iterative methods based on finite element simulation are effective approaches to design mold shape to compensate springback in sheet metal forming. However, convergence rate of iterative methods is difficult to improve greatly. To increase the springback compensate speed of designing age forming mold, process of calculating springback for a certain mold with finite element method is analyzed. Springback compensation is abstracted as finding a solution for a set of nonlinear functions and a springback compensation algorithm is presented on the basis of quasi Newton method. The accuracy of algorithm is verified by developing an ABAQUS secondary development program with MATLAB. Three rectangular integrated panels of dimensions 710 mmx750 mm integrated panels with intersected ribs of 10 mm are selected to perform case studies. The algorithm is used to compute mold contours for the panels with cylinder, sphere and saddle contours respectively and it takes 57%, 22% and 33% iterations as compared to that of displacement adjustment （DA） method. At the end of iterations, maximum deviations on the three panels are 0.618 4 mm, 0.624 1 mm and 0.342 0 mm that are smaller than the deviations determined by DA method （0.740 8 mm, 0.740 8 mm and 0.713 7 mm respectively）. In following experimental verification, mold contour for another integrated panel with 400 ram^380 mm size is designed by the algorithm. Then the panel is age formed in an autoclave and measured by a three dimensional digital measurement devise. Deviation between measuring results and the panel＇s design contour is less than 1 mm. Finally, the iterations with different mesh sizes （40 mm, 35 mm, 30 mm, 25 mm, 20 mm） in finite element models are compared and found no considerable difference. Another possible compensation method, Broyden-Fletcher-Shanmo method, is also presented based on the solving nonlinear fimctions idea. The Broyden-Fletcher-Shanmo method is employed to compute mold contour for the second panel. It only takes 50% iterations compared to that of DA. The proposed method can serve a faster mold contour compensation method for sheet metal forming.
基金Foundation item： Supported by the National Science Foundation of China（10701066） Supported by the National Foundation of the Education Department of Henan Province（2008A110022）
文摘A fourth-order convergence method of solving roots for nonlinear equation, which is a variant of Newton＇s method given. Its convergence properties is proved. It is at least fourth-order convergence near simple roots and one order convergence near multiple roots. In the end, numerical tests are given and compared with other known Newton and Newton-type methods. The results show that the proposed method has some more advantages than others. It enriches the methods to find the roots of non-linear equations and it is important in both theory and application.
基金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 NNSF of China（11071041）Supported by the Fujian Natural Science Foundation（2009J01002）Supported by the Fujian Department of Education Foundation（JA11270）
文摘The box constrained variational inequality problem can be reformulated as a nonsmooth equation by using median operator.In this paper,we present a smoothing Newton method for solving the box constrained variational inequality problem based on a new smoothing approximation function.The proposed algorithm is proved to be well defined and convergent globally under weaker conditions.
