THE SOLUTION FOR THE GENERALIZED RICCATI ALGEBRAIC EQUATIONS OF LINEAR EQUALITY CONSTRAINT SYSTEM

I.IntroductionTheestablishmentofanalogytheorybetweencomputationalstructuralmechanicsalldoptimalcontrolisbasedonthesubstructuralchaintheoryinstructuralmechanicsandthelillearquadratic(LQ)controlprobleminoptimalcontrolllJ.Whendealillgwiththecontrolproblenlbyapplyingtheabovethe.ry,wecanrealizetheadvantagesofmethodsinstrLlcturalInechanicsl2'3].Howevermostworkhasbeenlimitedtothenon-constraintsystemandlittleworkhasbeendonefortheconstraintsystemsofLlr.InRef.I4],basedonthegeneralizedvariationalprin…     

多体系统动力学方程的无违约数值计算方法

多体系统动力学方程为3阶微分代数方程,已有的约束违约稳定法存在位移违约问题,数值仿真准确性和稳定性不足。本文将求解高阶微分代数方程的降阶理论、ε嵌入处理方式与隐式龙格库塔法相结合,提出了直接满足位移约束条件的多体系统动力学方程的无违约算法,避免了约束违约问题。该方法先将多体动力学方程转化为2阶微分代数方程,并与位移约束方程联立;再应用ε嵌入隐式龙格库塔法进行数值求解。应用两种方法分别对单摆机构进行数值仿真,结果表明本文的方法不仅能适应较大步长,且准确性和稳定性均优于约束违约稳定法。     

Numerical solutions of linear quadratic control for time-varying systems via symplectic conservative perturbation

Optimal control system of state space is a conservative system, whose approximate method should be symplectic conservation. Based on the precise integration method, an algorithm of symplectic conservative perturbation is presented. It gives a uniform way to solve the linear quadratic control (LQ control) problems for linear time-varying systems accurately and efficiently, whose key points are solutions of differential Riccati equation (DRE) with variable coefficients and the state feedback equation. The method is symplectic conservative and has a good numerical stability and high precision. Numerical examples demonstrate the effectiveness of the proposed method.     

Numerical solutions of linear quadratic control for time-varying systems via symplectic par conservative perturbation

Optimal control system of state space is a conservative system, whose approximate method should be symplectic conservation. Based on the precise integration method, an algorithm of symplectic conservative perturbation is presented.It gives a uniform way to solve the linear quadratic control (LQ control) problems for linear time-varying systems accurately and efficiently, whose key points are solutions of differential Riccati equation (DRE) with variable coefficients and the state feedback equation.The method is symplectic conservative and has a good numerical stability and high precision. Numerical examples demonstrate the effectiveness of the proposed method.     

An automatic constraint violation stabilization method for differential/ algebraic equations of motion in multibody system dynamics

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｄｙｎａｍｉｃｅｑｕａｔｉｏｎｓｏｆｍｏｔｉｏｎｏｆｍｕｌｔｉｂｏｄｙｓｙｓｔｅｍｓｗｉｔｈｃｏｎｓｔｒａｉｎｔｓａｒｅｔｈｅｆｏｌｌｏｗｉｎｇｄｉｆｆｅｒｅｎｔｉａｌ／ａｌｇｅｂｒａｉｃｅｑｕａｔｉｏｎｓ,ｉ．ｅ．,Ｅ...     

Equivalent method for accurate solution to linear interval equations

Based on linear interval equations, an accurate interval finite element method for solving structural static problems with uncertain parameters in terms of optimization is discussed.On the premise of ensuring the consistency of solution sets, the original interval equations are equivalently transformed into some deterministic inequations.On this basis, calculating the structural displacement response with interval parameters is predigested to a number of deterministic linear optimization problems.The results are proved to be accurate to the interval governing equations.Finally, a numerical example is given to demonstrate the feasibility and efficiency of the proposed method.     

A linear solver based on algebraic multigrid and defect correction for the solution of adjoint RANS equations

A new solution approach for discrete adjoint Reynolds‐averaged Navier–Stokes equations is suggested. The approach is based on a combination of algebraic multigrid and defect correction. The efficient interplay of these two methods results in a fast and robust solution technique. Algebraic multigrid deals very well with unstructured grids, and defect correction completes it in such a way that a required second‐order accuracy of the solution is achieved. The performance of the suggested solution approach is demonstrated on a number of representative benchmarks.Copyright © 2014 John Wiley & Sons, Ltd.     

边界元法中线性方程组的拟波阵分块求解

以边界单元法形成的线性方程组为例提出了一种分块求解具有非线称满系数矩阵方程组的拟波阵方法。这种在于高斯消去法基本原理的分块求解法,采用了拟波阵求解技术,使常驻计算机内存的仅为二行分块的系数子矩阵。     

Nonconforming local projection stabilization for generalized Oseen equations

A new method of nonconforming local projection stabilization for the gen- eralized Oseen equations is proposed by a nonconforming inf-sup stable element pair for approximating the velocity and the pressure. The method has several attractive features. It adds a local projection term only on the sub-scale （H ≥ h）. The stabilized term is simple compared with the residual-free bubble element method. The method can handle the influence of strong convection. The numerical results agree with the theoretical expectations very well.     

基于等效系统假定的变厚度板弯曲广义积分变换解

针对等厚度薄板的弯曲问题,研究人员已给出了基于不同数值算法的经典数值解。针对变厚度薄板弯曲问题的解答较少,且以有限元数值模拟计算为主,计算耗时较大。本文基于广义积分变换原理建立了求解变厚度等效系统的广义积分变换算法,分析了线性和二次变化的变厚度板在多种边界条件下的弯曲问题,利用文献已发表结果同本文建立的广义积分变换解进行验证。计算结果表明,本文建立的基于广义积分变换的变厚度板弯曲求解方法具有较高准确性。同时,通过参数化分析手段,分别利用广义积分变换方法和有限元数值模拟方法讨论了不同边界约束和长宽比等条件对中心点处挠度的影响,计算结果具有较好的一致性,证明本文建立的广义积分变换方法可用于求解变厚度板弯曲问题,且具有较高的准确性。     

Periodicity and strict oscillation for generalized lyness equations

ＩｎｔｒｏｄｕｃｔｉｏｎＣｏｎｓｉｄｅｒｔｈｅｆｏｌｌｏｗｉｎｇｄｅｌａｙｄｉｆｆｅｒｅｎｃｅｅｑｕａｔｉｏｎｘｎ 1=ｘｎ(ａ ｂｘｎ)ｘｎ- 1,　　ｎ=0 ,1 ,2 ,… ,(1 )ｗｈｅｒｅａ ,ｂ∈ [0 ,∞ )ｗｉｔｈａ ｂ>0 (2 )ａｎｄｗｈｅｒｅｔｈｅｉｎｉｔｉａｌｖａｌｕｅｓｘ- 1ａｎｄｘ0 ａｒｅａｒｂｉｔｒａｒｙｐｏｓｉｔｉｖｅｎｕｍｂｅｒｓ.Ｅｑ .(1 )ｉｓｒｅｇａｒｄｅｄａｓａｇｅｎｅｒａｌｉｚｅｄＬｙｎｅｓｓｅｑｕａｔｉｏｎｂｙＧ .Ｌａｄａｓｉｎ [1 ] .Ｏｂｖｉｏｕｓｌｙ ,ｕｎｄｅｒｃｏｎｄｉｔｉｏｎ (2 ) ,…     

A new algorithm for solving differential/algebraic equations of multibody system dynamics

I.Intr0ductionTheequationsofmoti0nofconstrainedmultibodysystemdynamicscanbewrittenasfollows11lwherelistimeparameter*q6R'isgeneralizedcoordinatesvectorofsystem,M(q,t):RnXR-R"xnisgeneralizedmassn1atrixofsystem,k6RmisLagrangemultipliervector,fo(q,t):RnXR-R"(…     

多体系统动力学方程违约修正的数值计算方法

多体系统动力学方程为微分代数方程,一般将其转化成常微分方程组进行数值计算,在数值积分的过程中约束方程的违约会逐渐增大.本文对具有完整、定常约束的多体系统,在修改的带乘子Lagrange正则形式的方程的基础上,根据Baumgarte提出的违约修正的方法,给出了一种多体系统微分代数方程违约修正法和系统的动力学方程的矩阵表达式.通过对曲柄-滑块机构的数值仿真,计算结果表明本文给出的方法在计算精度和计算效率上好于Baumgarte提出的两种违约修正的方法.     

BIFURCATIONS OF TRAVELLING WAVE SOLUTIONS FOR GENERALIZED DRINFELD-SOKOLOV EQUATIONS

Ansatz method and the theory of dynamical systems are used to study the traveling wave solutions for the generalized Drinfeld-Sokolov equations. Under two groups of the parametric conditions, more solitary wave solutions, kink and anti-kink wave solutions and periodic wave solutions are obtained. Exact explicit parametric representations of these travelling wave solutions are given.     

The addendum for the generalized d'Alembert equations of motion

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｏｓｔｂａｓｉｃｃｏｎｄｉｔｉｏｎｓｆｏｒｆｅｅｄｂａｃｋｃｌｏｓｅｄ＿ｌｏｏｐｃｏｎｔｒｏｌｏｆｍａｎｉｐｕｌａｔｏｒｓａｒｅｔｈｅｆｉｎｅｓｔｒｕｃｔｕｒｅａｎｄｔｈｅｅｆｆｉｃｉｅｎｔｒｅａｌ＿ｔｉｍｅｃｏｍｐｕ...     

The small disturbance problems associated with the algebraic riccati equation of continuous linear time-invariant systems

The stability problem of the disturbed algebraic Riccati equation of continuous linear time-invariant systems is discussed in this paper. Through matrix norm analysis the estimation (expressed in terms of the disturbance range of the system parameters) of the disturbance range in the solution of the disturbed algebraic Riccati equation is established. Apparently this method is quite convenient for the practical computational purposes.     

Parallel Multisplitting AOR method for solving a class of system of nonlinear algebraic equations

ＰＡＲＡＬＬＥＬＭＵＬＴＩＳＰＬＩＴＴＩＮＧＡＯＲＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＡＣＬＡＳＳＯＦＳＹＳＴＥＭＯＦＮＯＮＬＩＮＥＡＲＡＬＧＥＢＲＡＩＣＥＱＵＡＴＩＯＮＳＢａｉＺｈｏｎｇｚｈｉ（白中治）（ＩｎｓｔｉｔｕｉｅｏｆＭａｔｈｅｍａｔｉｃｓＦｕ...     

Potential symmetries and conservation laws for generalized quasilinear hyperbolic equations

Based on the Lie group method, the potential symmetries and invariant solutions for generalized quasilinear hyperbolic equations are studied. To obtain the invariant solutions in an explicit form, the physically interesting situations with potential symmetries are focused on, and the conservation laws for these equations in three physically interesting cases are found by using the partial Lagrangian approach.     

Exponential stability of stochastic generalized porous media equations with jump

Stochastic generalized porous media equation with jump is considered. The aim is to show the moment exponential stability and the almost certain exponential stability of the stochastic equation.     

Instability of solution for the fourth order linear differential equation with varied coefficient

In this paper,we give some sufficient conditions of the instability for the fourth order linear differential equation with varied coefficient,at least one of the characteristic roots of which has positive real part,by means of Liapunov’s second method.