This paper explores model order reduction(MOR)methods for discrete linear and discrete bilinear systems via discrete pulse orthogonal functions(DPOFs).Firstly,the discrete linear systems and the discrete bilinear syst...This paper explores model order reduction(MOR)methods for discrete linear and discrete bilinear systems via discrete pulse orthogonal functions(DPOFs).Firstly,the discrete linear systems and the discrete bilinear systems are expanded in the space spanned by DPOFs,and two recurrence formulas for the expansion coefficients of the system’s state variables are obtained.Then,a modified Arnoldi process is applied to both recurrence formulas to construct the orthogonal projection matrices,by which the reduced-order systems are obtained.Theoretical analysis shows that the output variables of the reducedorder systems can match a certain number of the expansion coefficients of the original system’s output variables.Finally,two numerical examples illustrate the feasibility and effectiveness of the proposed methods.展开更多
Routh order reduction method of the relativistic Birkhoffian equations is studied. For a relativistic Birkhoffian system, the cyclic integrals can be found by using the perfect differential method. Through these cycli...Routh order reduction method of the relativistic Birkhoffian equations is studied. For a relativistic Birkhoffian system, the cyclic integrals can be found by using the perfect differential method. Through these cyclic integrals, the order of the system can be reduced. If the relativistic Birkhoffian system has a cyclic integral, then the Birkhoffian equations can be reduced at least by two degrees and the Birkhoffian form can be kept. The relations among the relativistic Birkhoffian mechanics, the relativistic Hamiltonian mechanics, and the relativistic Lagrangian mechanics are discussed, and the Routh order reduction method of the relativistic Lagrangian system is obtained. And an example is given to illustrate the application of the result.展开更多
Telemedicine plays an important role in Corona Virus Disease 2019(COVID-19).The virtual surgery simulation system,as a key component in telemedicine,requires to compute in real-time.Therefore,this paper proposes a rea...Telemedicine plays an important role in Corona Virus Disease 2019(COVID-19).The virtual surgery simulation system,as a key component in telemedicine,requires to compute in real-time.Therefore,this paper proposes a realtime cutting model based on finite element and order reduction method,which improves the computational speed and ensure the real-time performance.The proposed model uses the finite element model to construct a deformation model of the virtual lung.Meanwhile,a model order reduction method combining proper orthogonal decomposition and Galerkin projection is employed to reduce the amount of deformation computation.In addition,the cutting path is formed according to the collision intersection position of the surgical instrument and the lesion area of the virtual lung.Then,the Bezier curve is adopted to draw the incision outline after the virtual lung has been cut.Finally,the simulation system is set up on the PHANTOM OMNI force haptic feedback device to realize the cutting simulation of the virtual lung.Experimental results show that the proposed model can enhance the real-time performance of telemedicine,reduce the complexity of the cutting simulation and make the incision smoother and more natural.展开更多
Structural components may enter an initial-elastic state,a plastic-hardening state and a residual-elastic state during strong seismic excitations.In the residual-elastic state,structural components keep in an unloadin...Structural components may enter an initial-elastic state,a plastic-hardening state and a residual-elastic state during strong seismic excitations.In the residual-elastic state,structural components keep in an unloading/reloading stage that is dominated by a tangent stiffness,thus structural components remain residual deformations but behave in an elastic manner.It has a great potential to make model order reduction for such structural components using the tangent-stiffness-based vibration modes as a reduced order basis.In this paper,an adaptive substructure-based model order reduction method is developed to perform nonlinear seismic analysis for structures that have a priori unknown damage distribution.This method is able to generate time-varying substructures and make nonlinear model order reduction for substructures in the residual-elastic phase.The finite element program OpenSees has been extended to provide the adaptive substructure-based nonlinear seismic analysis.At the low level of OpenSees framework,a new abstract layer is created to represent the time-varying substructures and implement the modeling process of substructures.At the high level of OpenSees framework,a new transient analysis class is created to implement the solving process of substructure-based governing equations.Compared with the conventional time step integration method,the adaptive substructure-based model order reduction method can yield comparative results with a higher computational efficiency.展开更多
Model Order Reduction (MOR) plays more and more imp or tant role in complex system simulation, design and control recently. For example , for the large-size space structures, VLSI and MEMS (Micro-ElectroMechanical Sys...Model Order Reduction (MOR) plays more and more imp or tant role in complex system simulation, design and control recently. For example , for the large-size space structures, VLSI and MEMS (Micro-ElectroMechanical Systems) etc., in order to shorten the development cost, increase the system co ntrolling accuracy and reduce the complexity of controllers, the reduced order model must be constructed. Even in Virtual Reality (VR), the simulation and d isplay must be in real-time, the model order must be reduced too. The recent advances of MOR research are overviewed in the article. The MOR theor y and methods may be classified as Singular Value decomposition (SVD) based, the Krylov subspace based and others. The merits and demerits of the different meth ods are analyzed, and the existed problems are pointed out. Moreover, the applic ation’s fields are overviewed, and the potential applications are forecaste d. After the existed problems analyzed, the future work is described. There are som e problems in the traditional methods such as SVD and Krylov subspace, they are that it’s difficult to (1)guarantee the stability of the original system, (2) b e adaptive to nonlinear system, and (3) control the modeling accuracy. The f uture works may be solving the above problems on the foundation of the tradition al methods, and applying other methods such as wavelet or signal compression.展开更多
Numerical solutions could not perform rapid system-level simulation of the behavior of micro-electro-mechanical systems（MEMS） and analytic solutions for the describing partial differential equations are only availab...Numerical solutions could not perform rapid system-level simulation of the behavior of micro-electro-mechanical systems（MEMS） and analytic solutions for the describing partial differential equations are only available for simple geometries.Model order reduction（MOR） can extract approximate low-order model from the original large scale system.Conventional model order reduction algorithm is based on first-order system model,however,most structure mechanical MEMS systems are naturally second-order in time.For the purpose of solving the above problem,a direct second-order system model order reduction approach based on Krylov subspace projection for the coupled dynamic study of electrostatic torsional micromirrors is presented.The block Arnoldi process is applied to create the orthonormal vectors to construct the projection matrix,which enables the extraction of the low order model from the discretized system assembled through finite element analysis.The transfer functions of the reduced order model and the original model are expanded to demonstrate the moment-matching property of the second-order model reduction algorithm.The torsion and bending effect are included in the finite element model,and the squeeze film damping effect is considered as well.An empirical method considering relative error convergence is adopted to obtain the optimal choice of the order for the reduced model.A comparison research between the full model and the reduced model is carried out.The modeling accuracy and computation efficiency of the presented second-order model reduction method are confirmed by the comparison research results.The research provides references for MOR of MEMS.展开更多
Starting with the governing equations in terms of displacements of 3D elastic media, the solutions to displacement components and their first derivatives are obtained by the application of a double Fourier transform a...Starting with the governing equations in terms of displacements of 3D elastic media, the solutions to displacement components and their first derivatives are obtained by the application of a double Fourier transform and an order reduction method based on the Cayley-Hamilton theorem. Combining the solutions and the constitutive equations which connect the displacements and stresses, the transfer matrix of a single soil layer is acquired. Then, the state space solution to multilayered elastic soils is further obtained by introducing the boundary conditions and continuity conditions between adjacent soil layers. The numerical analysis based on the present theory is carried out, and the vertical displacements of multilayered foundation with a weak and a hard underlying stratums are compared and discussed.展开更多
Airframe structural optimization at different design stages results in new mass and stiffness distributions which modify the critical design loads envelop. Determination of aircraft critical loads is an extensive anal...Airframe structural optimization at different design stages results in new mass and stiffness distributions which modify the critical design loads envelop. Determination of aircraft critical loads is an extensive analysis procedure which involves simulating the aircraft at thousands of load cases as defmed in the certification requirements. It is computationally prohibitive to use a GFEM （Global Finite Element Model） for the load analysis, hence reduced order structural models are required which closely represent the dynamic characteristics of the GFEM. This paper presents the implementation of CMS （Component Mode Synthesis） method for the generation of high fidelity ROM （Reduced Order Model） of complex airframes. Here, sub-structuring technique is used to divide the complex higher order airframe dynamical system into a set of subsystems. Each subsystem is reduced to fewer degrees of freedom using matrix projection onto a carefully chosen reduced order basis subspace. The reduced structural matrices are assembled for all the subsystems through interface coupling and the dynamic response of the total system is solved. The CMS method is employed to develop the ROM of a Bombardier Aerospace business jet which is coupled with aerodynamic model for dynamic aeroelasticity loads analysis under gust turbulence. Another set of dynamic aeroelastic loads is also generated employing a stick model of same aircraft. Stick model is the reduced order modelling methodology commonly used in the aerospace industry based on stiffness generation by unitary loading application. The extracted aeroelastic loads from both models are compared against those generated employing the GFEM. Critical loads modal participation factors and modal characteristics of the different ROMs are investigated and compared against those of the GFEM. Results obtained show that the ROM generated using Craig Bampton CMS reduction process has a superior dynamic characteristics compared to the stick model.展开更多
The oscillatory stability analysis of multi-converter-fed systems(MCFSs)with modest computational resources needs a precise parametric reduced-order impedance model(PROIM).However,the traditional Krylov subspace based...The oscillatory stability analysis of multi-converter-fed systems(MCFSs)with modest computational resources needs a precise parametric reduced-order impedance model(PROIM).However,the traditional Krylov subspace based parametric model order reduction(KS-PMOR)method has difficulty in building precise PROIM for MCFSs.This is because the factors related to the errors of PROIM are complicated and coupled.To fill this gap,the factors associated with the accuracy of the KS-PMOR method are estimated by defining three indicators:the convergence error,cumulative error,and rank of projection matrix.Using the three indicators,a frequency-domain adaptive parametric model order reduction(FDA-PMOR)method is developed to form the precise PROIM of MCFSs for the accurate and fast oscillatory stability analysis.The accuracy of the obtained PROIM using the proposed FDA-PMOR method and its efficiency in actual oscillatory stability analysis are validated by three MCFSs with different scales,i.e.,a small-scale MCFS with four paralleled converter-based renewable energy generators(CREGs),a real-time simulation-based MCFS with eighteen paralleled CREGs,and a larger MCFS with ninety paralleled CREGs.展开更多
In this paper, we present an accelerated simulation approach on waveform relaxation using Krylov subspace for a large time-dependent system composed of some subsystems. This approach first allows these subsystems to b...In this paper, we present an accelerated simulation approach on waveform relaxation using Krylov subspace for a large time-dependent system composed of some subsystems. This approach first allows these subsystems to be decoupled by waveform relaxation. Then the Arnoldi procedure based on Krylov subspace is provided to accelerate the simulation of the decoupled subsystems independently. For the new approach, the convergent conditions on waveform relaxation are derived. The robust behavior is also successfully illustrated via numerical examples.展开更多
Model order reduction of interconnect circuits is an important technique to reduce the circuit complexity and improve the efficiency of post-layout verification process in the nanometer VLSI design. Existing works usi...Model order reduction of interconnect circuits is an important technique to reduce the circuit complexity and improve the efficiency of post-layout verification process in the nanometer VLSI design. Existing works using the Krylov subspace method are very efficient, but the resulting models are less compact and lack global accuracy. Also, existing methods cannot handle interconnect circuits with large input and output ports. Recent advances in reduction techniques using non-Krylov subspace techniques such as truncated balanced realization （TBR） hold some promise to solve these problems. In this paper, we first review the classic TBR-based reduction methods and then present the recent developments in fast TBR-based reduction and techniques such as PMTBR, SBPOR, and ETBR methods. These newly proposed methods try to avoid the expensive computing steps in traditional TBR methods at some cost to accuracy to boost efficiency and scalability, which is critical to reduce large interconnect parasitics modeled as RLCK circuits. The ETBR method can also reduce circuits with massive ports by considering the input signals. We show the pros and cons of each method and compare them on a set of large interconnect circuits, and finally point to some new research directions for this area.展开更多
Adaptive mesh refinement (AMR) is fairly practiced in the context of high-dimensional, mesh-based computational models. However, it is in its infancy in that of low-dimensional, generalized-coordinate-based computatio...Adaptive mesh refinement (AMR) is fairly practiced in the context of high-dimensional, mesh-based computational models. However, it is in its infancy in that of low-dimensional, generalized-coordinate-based computational models such as projection-based reduced-order models. This paper presents a complete framework for projection-based model order reduction (PMOR) of nonlinear problems in the presence of AMR that builds on elements from existing methods and augments them with critical new contributions. In particular, it proposes an analytical algorithm for computing a pseudo-meshless inner product between adapted solution snapshots for the purpose of clustering and PMOR. It exploits hyperreduction—specifically, the energy-conserving sampling and weighting hyperreduction method—to deliver for nonlinear and/or parametric problems the desired computational gains. Most importantly, the proposed framework for PMOR in the presence of AMR capitalizes on the concept of state-local reduced-order bases to make the most of the notion of a supermesh, while achieving computational tractability. Its features are illustrated with CFD applications grounded in AMR and its significance is demonstrated by the reported wall-clock speedup factors.展开更多
This paper presents a design of continuous-time sliding mode control for the higher order systems via reduced order model. It is shown that a continuous-time sliding mode control designed for the reduced order model g...This paper presents a design of continuous-time sliding mode control for the higher order systems via reduced order model. It is shown that a continuous-time sliding mode control designed for the reduced order model gives similar performance for thc higher order system. The method is illustrated by numerical examples. The paper also introduces a technique for design of a sliding surface such that the system satisfies a cost-optimality condition when on the sliding surface.展开更多
Reduction of complex protein networks models is of great importance.The accuracy of a passivity preserving algorithm (PRIMA) for model order reduction (MOR) is here tested on protein networks,introducing innovative va...Reduction of complex protein networks models is of great importance.The accuracy of a passivity preserving algorithm (PRIMA) for model order reduction (MOR) is here tested on protein networks,introducing innovative variations of the standard PRIMA method to fit the problem at hand.The reduction method does not require to solve the complete system,resulting in a promising tool for studying very large-scale models for which the full solution cannot be computed.The mathematical structure of the considered kinetic equations is preserved.Keeping constant the reduction factor,the approximation error is lower for larger systems.展开更多
In a previous paper,a technique was suggested to avoid order reduction with any explicit exponential Runge-Kutta method when integrating initial boundary value nonlinear problems with time-dependent boundary condition...In a previous paper,a technique was suggested to avoid order reduction with any explicit exponential Runge-Kutta method when integrating initial boundary value nonlinear problems with time-dependent boundary conditions.In this paper,we significantly simplify the full discretization formulas to be applied under conditions which are nearly always satisfied in practice.Not only a simpler linear combination of'j-functions is given for both the stages and the solution,but also the information required on the boundary is so much simplified that,in order to get local order three,it is no longer necessary to resort to numerical differentiation in space.In many cases,even to get local order 4.The technique is then shown to be computationally competitive against other widely used methods with high enough stiff order through the standard method of lines.展开更多
A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordin...A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.展开更多
In recent years,the Active Flutter Suppression(AFS)employing Linear ParameterVarying(LPV)framework has become a hot spot in the research field.Nevertheless,the flutter suppression technique is facing two severe challe...In recent years,the Active Flutter Suppression(AFS)employing Linear ParameterVarying(LPV)framework has become a hot spot in the research field.Nevertheless,the flutter suppression technique is facing two severe challenges.On the one hand,due to the fatal risk of flight test near critical airspeed,it is hard to obtain the accurate mathematical model of the aeroelastic system from the testing data.On the other hand,saturation of the actuator may degrade the closed-loop performance,which was often neglected in the past work.To tackle these two problems,a new active controller design procedure is proposed to suppress flutter in this paper.Firstly,with the aid of LPV model order reduction method and State-space Model Interpolation of Local Estimates(SMILE)technique,a set of high-fidelity Linear Time-Invariant(LTI)models which are usually derived from flight tests at different subcritical airspeeds are reduced and interpolated to construct an LPV model of an aeroelastic system.And then,the unstable aeroelastic dynamics beyond critical airspeed are‘predicted’by extrapolating the resulting LPV model.Secondly,based on the control-oriented LPV model,an AFS controller in LPV framework which is composed of a nominal LPV controller and an LPV anti-windup compensator is designed to suppress the aeroelastic vibration and overcome the performance degradation caused by actuator saturation.Although the nominal LPV controller may have superior performance in linear simulation in which the saturation effect is ignored,the results of the numerical simulations show that the nominal LPV controller fails to suppress the Body Freedom Flutter(BFF)when encountering the actuator saturation.However,the LPV anti-windup compensator not only enhances the nominal controller’s performance but also helps the nominal controller to stabilize the unstable aeroelastic system whenencountering serious actuator saturation.展开更多
Model order reduction(MOR)is considered as a good alternative to reduce the computational scale for electro-magnetic problems.The aim of this work is to introduce the use of dynamic mode decomposition(DMD)as a promisi...Model order reduction(MOR)is considered as a good alternative to reduce the computational scale for electro-magnetic problems.The aim of this work is to introduce the use of dynamic mode decomposition(DMD)as a promising tool for MOR to analyze its effectiveness in creating a fast model-based design platform for the permanent magnet motor design for ur-ban aerial vehicles(UAVs).Using a singular value decomposition(SVD)based DMD,the design process is constructed and verified against different scenarios.展开更多
This article is devoted to the problem of composite control design for continuous nonlinear singularly perturbed(SP)system using approximate feedback linearization(AFL)method.The essence of AFL method lies in the feed...This article is devoted to the problem of composite control design for continuous nonlinear singularly perturbed(SP)system using approximate feedback linearization(AFL)method.The essence of AFL method lies in the feedback linearization only of a certain part of the original nonlinear system.According to AFL approach,we suggest to solve feedback linearization problems for continuous nonlinear SP system by reducing it to two feedback linearization problems for slow and fast subsystems separately.The resulting AFL control is constructed in the form of asymptotic composition(composite control).Standard procedure for the composite control design consists of the following steps:1)system decomposition,2)solution of control problem for fast subsystem,3)solution of control problem for slow subsystem,4)construction of the resulting control in the form of the composition of slow and fast controls.The main difficulty during system decomposition is associated with dynamics separation condition for nonlinear SP system.To overcome this,we propose to change the sequence of the design procedure:1)solving the control problem for fast state variables part,2)system decomposition,3)solving the control problem for slow state variables part,4)construction of the resulting composite control.By this way,fast feedback linearizing control is chosen so that the dynamics separation condition would be met and the fast subsystem would be stabilizable.The application of the proposed approach is illustrated through several examples.展开更多
基金supported by Natural Science Foundation of Xinjiang Uygur Autonomous Region of China“Research on model order reduction methods based on the discrete orthogonal polynomials”(2023D01C163)The Tianchi Talent Introduction Plan Project of Xinjiang Uygur Autonomous Region of China“Research on orthogonal decomposition model order reduction methods for discrete control systems”.
文摘This paper explores model order reduction(MOR)methods for discrete linear and discrete bilinear systems via discrete pulse orthogonal functions(DPOFs).Firstly,the discrete linear systems and the discrete bilinear systems are expanded in the space spanned by DPOFs,and two recurrence formulas for the expansion coefficients of the system’s state variables are obtained.Then,a modified Arnoldi process is applied to both recurrence formulas to construct the orthogonal projection matrices,by which the reduced-order systems are obtained.Theoretical analysis shows that the output variables of the reducedorder systems can match a certain number of the expansion coefficients of the original system’s output variables.Finally,two numerical examples illustrate the feasibility and effectiveness of the proposed methods.
基金The project supported by National Natural Science Foundation of China under Grant Nos, 10372053 and 10472040, the Natural Science Foundation of Hunan Province under Grant No. 03JJY3005, the Scientific Research Foundation of Eduction Burean of Hunan Province under Grant No. 02C033 and the 0utstanding Young Talents Training Fund of Liaoning Province under Grant No. 3040005
文摘Routh order reduction method of the relativistic Birkhoffian equations is studied. For a relativistic Birkhoffian system, the cyclic integrals can be found by using the perfect differential method. Through these cyclic integrals, the order of the system can be reduced. If the relativistic Birkhoffian system has a cyclic integral, then the Birkhoffian equations can be reduced at least by two degrees and the Birkhoffian form can be kept. The relations among the relativistic Birkhoffian mechanics, the relativistic Hamiltonian mechanics, and the relativistic Lagrangian mechanics are discussed, and the Routh order reduction method of the relativistic Lagrangian system is obtained. And an example is given to illustrate the application of the result.
基金supported,in part,by the Natural Science Foundation of Jiangsu Province under Grant Numbers BK20201136,BK20191401in part,by the National Nature Science Foundation of China under Grant Numbers 61502240,61502096,61304205,61773219in part,by the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)fund.
文摘Telemedicine plays an important role in Corona Virus Disease 2019(COVID-19).The virtual surgery simulation system,as a key component in telemedicine,requires to compute in real-time.Therefore,this paper proposes a realtime cutting model based on finite element and order reduction method,which improves the computational speed and ensure the real-time performance.The proposed model uses the finite element model to construct a deformation model of the virtual lung.Meanwhile,a model order reduction method combining proper orthogonal decomposition and Galerkin projection is employed to reduce the amount of deformation computation.In addition,the cutting path is formed according to the collision intersection position of the surgical instrument and the lesion area of the virtual lung.Then,the Bezier curve is adopted to draw the incision outline after the virtual lung has been cut.Finally,the simulation system is set up on the PHANTOM OMNI force haptic feedback device to realize the cutting simulation of the virtual lung.Experimental results show that the proposed model can enhance the real-time performance of telemedicine,reduce the complexity of the cutting simulation and make the incision smoother and more natural.
基金supported by the National Nature Science Foundation of China(No.51678210)National Key Research and Development Program of China(No.2016YFC0701400).
文摘Structural components may enter an initial-elastic state,a plastic-hardening state and a residual-elastic state during strong seismic excitations.In the residual-elastic state,structural components keep in an unloading/reloading stage that is dominated by a tangent stiffness,thus structural components remain residual deformations but behave in an elastic manner.It has a great potential to make model order reduction for such structural components using the tangent-stiffness-based vibration modes as a reduced order basis.In this paper,an adaptive substructure-based model order reduction method is developed to perform nonlinear seismic analysis for structures that have a priori unknown damage distribution.This method is able to generate time-varying substructures and make nonlinear model order reduction for substructures in the residual-elastic phase.The finite element program OpenSees has been extended to provide the adaptive substructure-based nonlinear seismic analysis.At the low level of OpenSees framework,a new abstract layer is created to represent the time-varying substructures and implement the modeling process of substructures.At the high level of OpenSees framework,a new transient analysis class is created to implement the solving process of substructure-based governing equations.Compared with the conventional time step integration method,the adaptive substructure-based model order reduction method can yield comparative results with a higher computational efficiency.
文摘Model Order Reduction (MOR) plays more and more imp or tant role in complex system simulation, design and control recently. For example , for the large-size space structures, VLSI and MEMS (Micro-ElectroMechanical Systems) etc., in order to shorten the development cost, increase the system co ntrolling accuracy and reduce the complexity of controllers, the reduced order model must be constructed. Even in Virtual Reality (VR), the simulation and d isplay must be in real-time, the model order must be reduced too. The recent advances of MOR research are overviewed in the article. The MOR theor y and methods may be classified as Singular Value decomposition (SVD) based, the Krylov subspace based and others. The merits and demerits of the different meth ods are analyzed, and the existed problems are pointed out. Moreover, the applic ation’s fields are overviewed, and the potential applications are forecaste d. After the existed problems analyzed, the future work is described. There are som e problems in the traditional methods such as SVD and Krylov subspace, they are that it’s difficult to (1)guarantee the stability of the original system, (2) b e adaptive to nonlinear system, and (3) control the modeling accuracy. The f uture works may be solving the above problems on the foundation of the tradition al methods, and applying other methods such as wavelet or signal compression.
基金supported by National Natural Science Foundation of China (Grant No. 50775201)National Science & Technology Major Project of China (Grant No. 2009ZX04014-031)PhD Programs Foundation of Ministry of Education of China (Grant No. 200803350031)
文摘Numerical solutions could not perform rapid system-level simulation of the behavior of micro-electro-mechanical systems（MEMS） and analytic solutions for the describing partial differential equations are only available for simple geometries.Model order reduction（MOR） can extract approximate low-order model from the original large scale system.Conventional model order reduction algorithm is based on first-order system model,however,most structure mechanical MEMS systems are naturally second-order in time.For the purpose of solving the above problem,a direct second-order system model order reduction approach based on Krylov subspace projection for the coupled dynamic study of electrostatic torsional micromirrors is presented.The block Arnoldi process is applied to create the orthonormal vectors to construct the projection matrix,which enables the extraction of the low order model from the discretized system assembled through finite element analysis.The transfer functions of the reduced order model and the original model are expanded to demonstrate the moment-matching property of the second-order model reduction algorithm.The torsion and bending effect are included in the finite element model,and the squeeze film damping effect is considered as well.An empirical method considering relative error convergence is adopted to obtain the optimal choice of the order for the reduced model.A comparison research between the full model and the reduced model is carried out.The modeling accuracy and computation efficiency of the presented second-order model reduction method are confirmed by the comparison research results.The research provides references for MOR of MEMS.
文摘Starting with the governing equations in terms of displacements of 3D elastic media, the solutions to displacement components and their first derivatives are obtained by the application of a double Fourier transform and an order reduction method based on the Cayley-Hamilton theorem. Combining the solutions and the constitutive equations which connect the displacements and stresses, the transfer matrix of a single soil layer is acquired. Then, the state space solution to multilayered elastic soils is further obtained by introducing the boundary conditions and continuity conditions between adjacent soil layers. The numerical analysis based on the present theory is carried out, and the vertical displacements of multilayered foundation with a weak and a hard underlying stratums are compared and discussed.
文摘Airframe structural optimization at different design stages results in new mass and stiffness distributions which modify the critical design loads envelop. Determination of aircraft critical loads is an extensive analysis procedure which involves simulating the aircraft at thousands of load cases as defmed in the certification requirements. It is computationally prohibitive to use a GFEM （Global Finite Element Model） for the load analysis, hence reduced order structural models are required which closely represent the dynamic characteristics of the GFEM. This paper presents the implementation of CMS （Component Mode Synthesis） method for the generation of high fidelity ROM （Reduced Order Model） of complex airframes. Here, sub-structuring technique is used to divide the complex higher order airframe dynamical system into a set of subsystems. Each subsystem is reduced to fewer degrees of freedom using matrix projection onto a carefully chosen reduced order basis subspace. The reduced structural matrices are assembled for all the subsystems through interface coupling and the dynamic response of the total system is solved. The CMS method is employed to develop the ROM of a Bombardier Aerospace business jet which is coupled with aerodynamic model for dynamic aeroelasticity loads analysis under gust turbulence. Another set of dynamic aeroelastic loads is also generated employing a stick model of same aircraft. Stick model is the reduced order modelling methodology commonly used in the aerospace industry based on stiffness generation by unitary loading application. The extracted aeroelastic loads from both models are compared against those generated employing the GFEM. Critical loads modal participation factors and modal characteristics of the different ROMs are investigated and compared against those of the GFEM. Results obtained show that the ROM generated using Craig Bampton CMS reduction process has a superior dynamic characteristics compared to the stick model.
基金supported by Open Fund of State Key Laboratory of Operation and Control of Renewable Energy&Storage Systems(China Electric Power Research Institute)(No.NYB51202201695)National Natural Science Foundation of China(No.51677050)111 Project(No.BP0719039).
文摘The oscillatory stability analysis of multi-converter-fed systems(MCFSs)with modest computational resources needs a precise parametric reduced-order impedance model(PROIM).However,the traditional Krylov subspace based parametric model order reduction(KS-PMOR)method has difficulty in building precise PROIM for MCFSs.This is because the factors related to the errors of PROIM are complicated and coupled.To fill this gap,the factors associated with the accuracy of the KS-PMOR method are estimated by defining three indicators:the convergence error,cumulative error,and rank of projection matrix.Using the three indicators,a frequency-domain adaptive parametric model order reduction(FDA-PMOR)method is developed to form the precise PROIM of MCFSs for the accurate and fast oscillatory stability analysis.The accuracy of the obtained PROIM using the proposed FDA-PMOR method and its efficiency in actual oscillatory stability analysis are validated by three MCFSs with different scales,i.e.,a small-scale MCFS with four paralleled converter-based renewable energy generators(CREGs),a real-time simulation-based MCFS with eighteen paralleled CREGs,and a larger MCFS with ninety paralleled CREGs.
基金This work was supported by the Natural Science Foundation of China（NSFC） under grant 11071192 and the International Science and Technology Cooperation Program of China under grant 2010DFA14700.
文摘In this paper, we present an accelerated simulation approach on waveform relaxation using Krylov subspace for a large time-dependent system composed of some subsystems. This approach first allows these subsystems to be decoupled by waveform relaxation. Then the Arnoldi procedure based on Krylov subspace is provided to accelerate the simulation of the decoupled subsystems independently. For the new approach, the convergent conditions on waveform relaxation are derived. The robust behavior is also successfully illustrated via numerical examples.
基金Supported in part by National Science Foundation (NSF) (Nos.CCF-0448534 and OISE-0929699)in part by the National Natural Science Foundation of China (No. 60828008)
文摘Model order reduction of interconnect circuits is an important technique to reduce the circuit complexity and improve the efficiency of post-layout verification process in the nanometer VLSI design. Existing works using the Krylov subspace method are very efficient, but the resulting models are less compact and lack global accuracy. Also, existing methods cannot handle interconnect circuits with large input and output ports. Recent advances in reduction techniques using non-Krylov subspace techniques such as truncated balanced realization （TBR） hold some promise to solve these problems. In this paper, we first review the classic TBR-based reduction methods and then present the recent developments in fast TBR-based reduction and techniques such as PMTBR, SBPOR, and ETBR methods. These newly proposed methods try to avoid the expensive computing steps in traditional TBR methods at some cost to accuracy to boost efficiency and scalability, which is critical to reduce large interconnect parasitics modeled as RLCK circuits. The ETBR method can also reduce circuits with massive ports by considering the input signals. We show the pros and cons of each method and compare them on a set of large interconnect circuits, and finally point to some new research directions for this area.
基金support by the Air Force Office of Scientific Research under Grant No.FA9550-20-1-0358 and Grant No.FA9550-22-1-0004.
文摘Adaptive mesh refinement (AMR) is fairly practiced in the context of high-dimensional, mesh-based computational models. However, it is in its infancy in that of low-dimensional, generalized-coordinate-based computational models such as projection-based reduced-order models. This paper presents a complete framework for projection-based model order reduction (PMOR) of nonlinear problems in the presence of AMR that builds on elements from existing methods and augments them with critical new contributions. In particular, it proposes an analytical algorithm for computing a pseudo-meshless inner product between adapted solution snapshots for the purpose of clustering and PMOR. It exploits hyperreduction—specifically, the energy-conserving sampling and weighting hyperreduction method—to deliver for nonlinear and/or parametric problems the desired computational gains. Most importantly, the proposed framework for PMOR in the presence of AMR capitalizes on the concept of state-local reduced-order bases to make the most of the notion of a supermesh, while achieving computational tractability. Its features are illustrated with CFD applications grounded in AMR and its significance is demonstrated by the reported wall-clock speedup factors.
文摘This paper presents a design of continuous-time sliding mode control for the higher order systems via reduced order model. It is shown that a continuous-time sliding mode control designed for the reduced order model gives similar performance for thc higher order system. The method is illustrated by numerical examples. The paper also introduces a technique for design of a sliding surface such that the system satisfies a cost-optimality condition when on the sliding surface.
文摘Reduction of complex protein networks models is of great importance.The accuracy of a passivity preserving algorithm (PRIMA) for model order reduction (MOR) is here tested on protein networks,introducing innovative variations of the standard PRIMA method to fit the problem at hand.The reduction method does not require to solve the complete system,resulting in a promising tool for studying very large-scale models for which the full solution cannot be computed.The mathematical structure of the considered kinetic equations is preserved.Keeping constant the reduction factor,the approximation error is lower for larger systems.
基金supported by Junta de Castilla y León and Feder(Project VA169P20).
文摘In a previous paper,a technique was suggested to avoid order reduction with any explicit exponential Runge-Kutta method when integrating initial boundary value nonlinear problems with time-dependent boundary conditions.In this paper,we significantly simplify the full discretization formulas to be applied under conditions which are nearly always satisfied in practice.Not only a simpler linear combination of'j-functions is given for both the stages and the solution,but also the information required on the boundary is so much simplified that,in order to get local order three,it is no longer necessary to resort to numerical differentiation in space.In many cases,even to get local order 4.The technique is then shown to be computationally competitive against other widely used methods with high enough stiff order through the standard method of lines.
基金The National Natural Science Foundation of China (No10471023)
文摘A numerical simulation for a model of wood drying process is considered. The model is given by a couple of nonlinear differential equations. One is a nonlinear parabolic equation and the other one is a nonlinear ordinary equation. A difference scheme is derived by the method of reduction of order. First, a new variable is introduced and the original problem is rewritten into a system of the first-order differential equations. Secondly, a difference scheme is constructed for the later problem. The solvability, stability and convergence of the difference scheme are proved by the energy method. The convergence order of the difference scheme is secondorder both in time and in space. A prior error estimate is put forward. The new variable is put aside to reduce the computational cost. A numerical example testifies the theoretical result.
基金the National Natural Science Foundation of China(No.61573289)Space Science and Technology Fund,and Natural Science Basic Research Plan in Shaanxi Province of China(No.2019JM042)Fundamental Research Funds for the Central Universities of China(No.3102019ZDHKY11)。
文摘In recent years,the Active Flutter Suppression(AFS)employing Linear ParameterVarying(LPV)framework has become a hot spot in the research field.Nevertheless,the flutter suppression technique is facing two severe challenges.On the one hand,due to the fatal risk of flight test near critical airspeed,it is hard to obtain the accurate mathematical model of the aeroelastic system from the testing data.On the other hand,saturation of the actuator may degrade the closed-loop performance,which was often neglected in the past work.To tackle these two problems,a new active controller design procedure is proposed to suppress flutter in this paper.Firstly,with the aid of LPV model order reduction method and State-space Model Interpolation of Local Estimates(SMILE)technique,a set of high-fidelity Linear Time-Invariant(LTI)models which are usually derived from flight tests at different subcritical airspeeds are reduced and interpolated to construct an LPV model of an aeroelastic system.And then,the unstable aeroelastic dynamics beyond critical airspeed are‘predicted’by extrapolating the resulting LPV model.Secondly,based on the control-oriented LPV model,an AFS controller in LPV framework which is composed of a nominal LPV controller and an LPV anti-windup compensator is designed to suppress the aeroelastic vibration and overcome the performance degradation caused by actuator saturation.Although the nominal LPV controller may have superior performance in linear simulation in which the saturation effect is ignored,the results of the numerical simulations show that the nominal LPV controller fails to suppress the Body Freedom Flutter(BFF)when encountering the actuator saturation.However,the LPV anti-windup compensator not only enhances the nominal controller’s performance but also helps the nominal controller to stabilize the unstable aeroelastic system whenencountering serious actuator saturation.
基金This work was supported by Dong-A University research fund.(Corresponding author:J.Chang)
文摘Model order reduction(MOR)is considered as a good alternative to reduce the computational scale for electro-magnetic problems.The aim of this work is to introduce the use of dynamic mode decomposition(DMD)as a promising tool for MOR to analyze its effectiveness in creating a fast model-based design platform for the permanent magnet motor design for ur-ban aerial vehicles(UAVs).Using a singular value decomposition(SVD)based DMD,the design process is constructed and verified against different scenarios.
基金supported by Russian Foundation for Basic Research(No.15-08-06859a)and by the Ministry of Education and Science of the Russian Federation in the framework of the basic part of the state order(No.2.8629.2017).
文摘This article is devoted to the problem of composite control design for continuous nonlinear singularly perturbed(SP)system using approximate feedback linearization(AFL)method.The essence of AFL method lies in the feedback linearization only of a certain part of the original nonlinear system.According to AFL approach,we suggest to solve feedback linearization problems for continuous nonlinear SP system by reducing it to two feedback linearization problems for slow and fast subsystems separately.The resulting AFL control is constructed in the form of asymptotic composition(composite control).Standard procedure for the composite control design consists of the following steps:1)system decomposition,2)solution of control problem for fast subsystem,3)solution of control problem for slow subsystem,4)construction of the resulting control in the form of the composition of slow and fast controls.The main difficulty during system decomposition is associated with dynamics separation condition for nonlinear SP system.To overcome this,we propose to change the sequence of the design procedure:1)solving the control problem for fast state variables part,2)system decomposition,3)solving the control problem for slow state variables part,4)construction of the resulting composite control.By this way,fast feedback linearizing control is chosen so that the dynamics separation condition would be met and the fast subsystem would be stabilizable.The application of the proposed approach is illustrated through several examples.
