A METHOD FOR SOLVING THE DYNAMICS OF MULTIBODY SYSTEMS WITH RHEONOMIC AND NONHOLONOMIC CONSTRAINTS

ＡＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＴＨＥＤＹＮＡＭＩＣＳＯＦＭＵＬＴＩＢＯＤＹＳＹＳＴＥＭＳＷＩＴＨＲＨＥＯＮＯＭＩＣＡＮＤＮＯＮＨＯＬＯＮＯＭＩＣＣＯＮＳＴＲＡＩＮＴＳ￥ＳｈｕｉＸｉａｏｐｉｎｇ（水小平）ＺｈａｎｇＹｏｎｇｆａ（张永发）（Ｄｅｐａｒ...     

PRECONDITIONED GAUSS-SEIDEL TYPE ITERATIVE METHOD FOR SOLVING LINEAR SYSTEMS

The preconditioned Gauss-Seidel type iterative method for solving linear systems, with the proper choice of the preconditioner, is presented. Convergence of the preconditioned method applied to Z-matrices is discussed. Also the optimal parameter is presented. Numerical results show that the proper choice of the preconditioner can lead to effective by the preconditioned Gauss-Seidel type iterative methods for solving linear systems.     

ENERGY THEOREMS FOR SOLVING EQUATIONS OF DEFLECTIONS

In this paper the unit-dummy-load method is ge-neralized on the basis of Castigliano’s Theorem.Onthese grounds the general equations of deflection sur-faces of the structures.such as a kind of beams,platesand shells.are directly derived by the force method.We derived the equations of the deflection sur-faces of the rectangular thin plates and thick platesconsidering the effect of transverse shearing defor-mations with the inhomogeneous displacement boundaryconditions.At the same time we give the equations ofdeflection axes of the corresponding straight beams.The applications of the reciprocal theorem arealso generalized.Three simple calculated examples are given.     

THE RELATIONSHIP BETWEEN THE STABILITY AND THE OPTIMALITY OF LINEAR SYSTEMS(Ⅱ)

Different from the inverse problem put forward by R.E.Kalman, another kind ofinverse problem of linear optimal control is proposed and discussed in[1] as follows:Givenan asymptotically stable linear constant system and a nonnegative quadratic performanceindex, when can a state-feedback be separated from the stable system so that this state-feedback control law is optimal for the given index? In this paper this problem is extended.Similar conclusions are obtained for linear discrete systems and linear time-variablesystems. According to these conclusions we can say that the correspondence between theasymptotically stable system and the optimal feedback system is the inherent character ofall kinds of linear systems.     

THE PERTURBATION FINITE ELEMENT METHOD FOR SOLVING THE PLANE PROBLEM IN CONSIDERATION OF LINEAR CREEP

In this paper, amethod (PFMC) for solving plane problem of linear creep is presented by using perturbation finite element. It can be used in plane problem in consideration of creep, such as reinforced concrete beam, prestressed concrete beam, reinforced concrete cylinder and reinforced concrete tunnel in elastic or visco-elastic medium, as well as underground building and so on. In the presented method, the assumption made in the general increment method that variables remain constant in a divtded time interval is not taken. The accuracy is improved and the length of time step becomes larger. The computer storage can be reduced and the calculating efficiency can be increased. Perturbation finite element formulae for four-node quadrilateral isoparametric element including reinforcement are established and five numerical examples are given. As contrasted with the analytical solution, the accuracy is satisfactory.     

FUZZY ARITHMETIC AND SOLVING OF THE STATIC GOVERNING EQUATIONS OF FUZZY FINITE ELEMENT METHOD

The key component of finite element analysis of structures with fuzzy parameters, which is associated with handling of some fuzzy information and arithmetic relation of fuzzy variables , was the solving of the governing equations of fuzzy finite element method. Based on a given interval representation of fuzzy numbers, some arithmetic rules of fuzzy numbers and fuzzy variables were developed in terms of the properties of interval arithmetic. According to the rules and by the theory of interval finite element method, procedures for solving the static governing equations of fuzzy finite element method of structures were presented. By the proposed procedure, the possibility distributions of responses of fuzzy structures can be generated in terms of the membership functions of the input fuzzy numbers. It is shown by a numerical example that the computational burden of the presented procedures is low and easy to implement. The effectiveness and usefulness of the presented procedures are also illustrated.     

IMPULSIVE CONTROL FOR THE STABILIZATION AND SYNCHRONIZATION OF LUR'E SYSTEMS

IntroductionInpractice ,thereexistmanyexamplesofimpulsivecontrolsystems (seeRefs.[1 ,2 ,3 ] ) .Recently ,impulsivecontrolhasbeenwidelyusedtostabilizeandsynchronizechaoticsystems(seeRefs.[4～ 1 5] ) .Forexample,thepredictivePoincarecontrol[4 ]andtheoccasionalproportionalfeedbackcontrol[5 ]aretwoimpulsivecontrolschemeswithvaryingimpulsiveintervals.Yangetal.inRefs.[7,1 1 ]studied ,respectively ,thestabilizationandsynchronizationofaclassofchaoticsystemscalledLorenzsystems .YangandChuainRef.[1 …     

THE ANALYTIC RESOLUTIONS AND APPLICATIONS OF THE NON－LINEAR SEEPAGE FLOW EQUATIONS OF COAL INFUSION

ＴＨＥＡＮＡＬＹＴＩＣＲＥＳＯＬＵＴＩＯＮＳＡＮＤＡＰＰＬＩＣＡＴＩＯＮＳＯＦＴＨＥＮＯＮ－ＬＩＮＥＡＲＳＥＥＰＡＧＥＦＬＯＷＥＱＵＡＴＩＯＮＳＯＦＣＯＡＬＩＮＦＵＳＩＯＮＺｈａｎｇＹａｎｓｏｎｇ（张延松）（ＲｅｃｅｉｖｅｄＪａｎ．１,１９９４．Ｃｏ...     

APPLICATION 0F THE THEORY OF BRANCHING OF SOLUTIONS OF NONLINEAR EQUATIONS IN THE ESTIMATION OF POLE PERTURBANCE OF LINEAR CLOSED-LOOP SYSTEMS

When the state and input matrices of a multivariable lineartime-invariant system are perturbed.the problem of the es-timation of pole perturbance of closed-loop system is con-sidered by making use of the theory of branching of solutionsfor nonlinear equations.     

THE THEOREM OF THE STABILITY OF LINEAR NONAUTONOMOUS SYSTEMS UNDER THE FREQUENTLY-ACTING PERTURBATION

In the paper a theorem on the stability of linear nonautonomous system under thefrequently-acting perturbation has been given and proved on the basis of Malkin＇sTheorem.     

Asynchronous relaxed iterative methods for solving linear systems of equations

I.AsynchronousRelaxedIterativeMethodsMoreandmorelargeorverylargescaleproblemsofscientificcomputationhavebeenproposedandarebeingproposedinmanyimportantfieldsofscienceandengineering.Manyoftheseproblemsresultinsolvinglargeorverylargelinearalgebraicsystemsofe…     

Local and parallel finite element algorithms for time-dependent convection-diffusion equations

Local and parallel finite element algorithms based on two-grid discretization for the time-dependent convection-diffusion equations are presented. These algorithms are motivated by the observation that, for a solution to the convection-diffusion problem, low frequency components can be approximated well by a relatively coarse grid and high frequency components can be computed on a fine grid by some local and parallel proce- dures. Hence, these local and parallel algorithms only involve one small original problem on the coarse mesh and some correction problems on the local fine grid. One technical tool for the analysis is the local a priori estimates that are also obtained. Some numerical examples are given to support our theoretical analvsis.     

求解非线性方程组的混合遗传算法

非线性方程组的求解是数值计算领域中最困难的问题。大多数的数值求解算法例如牛顿法的收敛性和性能特征在很大程度上依赖于初始点。但是对于很多非线性方程组,选择好的初始点是一件非常困难的事情。本文结合遗传算法和经典算法的优点,提出了一种用于求解非线性方程组的混合遗传算法。该混合算法充分发挥了遗传算法的群体搜索和全局收敛性,有效地克服了经典算法的初始点敏感问题;同时在遗传算法中引入经典算法(Powell法、拟牛顿迭代法)作局部搜索,克服了遗传算法收敛速度慢和精度差的缺点。选择了几个典型非线性方程组,从收敛可靠性、计算成本和适用性等指标分析对比了不同算法。计算结果表明所设计的混合遗传算法有着可靠的收敛性和较高的收敛速度和精度,是求解非线性方程组的一种成功算法。     

The correlative potential function and a new method for solving Maxwell equations

This paper is a continuation of Paper [1]. 1. A new potential ψ which is defined as the correlative potential is developed in this paper. The potential ψ is different from the classical scalar potential φ and the vector potential developed by Helmholtz. The new formulae of the solution of eqs.are given in terms of ψ· 2. In the time-varying electromagnetic field, two new retarded potentials, the electric type retarded correlative potential ψ_e and the magnetic-type retarded correlative potential ψ_m, which are distinct from the classical retarded potentials and φ,are used to solve Maxwell equations. The new formulae of solution of Maxwell equation are given in terms of ψ_e and ψ_m 3. The methods for constructing a rotational field with given curl function (vorticity) is proposed.     

求解非线性方程组的混沌粒子群算法及应用

针对非线性方程组的求解在工程上具有广泛的实际意义,经典的数值算法如牛顿法存在其收敛性依赖于初值而实际计算中初值难确定的问题,提出以混沌粒子群算法求解非线性方程。它通过将混沌搜索机制有机地引入粒子群算法,使每个粒子从混沌搜索机制与粒子群算法搜索机制中获得适当的搜索方向,以混沌变量的遍历性增强粒子的搜索性能与更全面地应用目标函数的信息,并反映到逐代更新的个体极值和群体极值中,可更有效地调整粒子的移向并最终获得最优解。测试结果表明这一尝试的有效性。最后将所提的方法用于建立复合材料结构的疲劳寿命与应力、温度、湿度的关系模型。     

基于物理信息的深度学习求解矩形薄板力学正反问题

建立基于物理信息的神经网络框架,利用深度学习求解矩形薄板力学正反问题.力学正问题为已知矩形薄板的基本参数、边界条件和受力情况,求薄板各点挠度;反问题为已知薄板部分点的挠度、基本参数和受力情况等,识别边界条件.基于物理信息的神经网络模型中,损失函数除基于数据驱动模型的挠度数据拟合部分以外,还引入薄板弯曲基本方程和应力应变...     

Weak formulation and first order form of the equations of motion for a class of constrained mechanical systems

Some new theoretical results are presented on modeling the dynamic response of a class of discrete mechanical systems subject to equality motion constraints. Both the development and presentation are facilitated by employing some fundamental concepts of differential geometry. At the beginning, the equations of motion of the corresponding unconstrained system are presented on a configuration manifold with general properties, first in strong and then in a primal weak form, using Newton׳s law of motion as a foundation. Next, the final weak form is obtained by performing a crucial integration by parts step, involving a covariant derivative. This step required the clarification and enhancement of some concepts related to the variations employed in generating the weak form. The second part of this work is devoted to systems involving holonomic and non-holonomic scleronomic constraints. The equations of motion derived in a recent study of the authors are utilized as a basis. The novel characteristic of these equations is that they form a set of second order ordinary differential equations (ODEs) in both the coordinates and the Lagrange multipliers associated to the constraint action. Based on these equations, the corresponding weak form is first obtained, leading eventually to a consistent first order ODE form of the equations of motion. These equations are found to appear in a form resembling the form obtained after application of the classical Hamilton׳s canonical equations. Finally, the new theoretical findings are illustrated by three representative examples.     

Some existence theorems of common and coincidence solutions for a class of systems of functional equations arising in dynamic programming

Some existence theorems of common and coincidence solutions for a class of moregeneral systems of functional equations arising in dynamic programming are shown.Theresults presented in this paper not only contain the corresponding results of[6,7]as specialcases,but also give an existence theorem of solutions for a class of functional equationssuggested by Wang recently.     

The method of mixed boundary condition for a kion of linear and nonlinear composite structure

This paper deals with the axisymmetrical deformation of shallow shells in large deflection which are in conjunction with linear elastic structures at the boundary: A method of mixed boundary condition for this problem is introduced, then the problem of a composite structure is transformed into a problem of a single structure and the integral equations are given. The perturbation method is used to obtain the solutions and an example of composite structure consisting of a shallow spherical and a cylindrical shell is presented.Communicated by Yeh Kai-yuan