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

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

Application of the modified method of multiple scales to solving the problem of a thin clamped circular plate of a very large deflection

In this paper,asymptotic behaviour of the solution to the problem of a thin clamped circular plate under uniform normal pressure at very large deflection is restudied by means of the modified method of multiple scales given in[1?].The result presented herein is in good agreement with the one obtained by professor Chien Wei-zang who first proposed the method of composite expansions to solve this problem in[3].However,by contrast,the advantage of the modified method of multiple scales it seems to be relatively simpler than the method used in[3].It is also shown that the restriction of the method of paper[1-2]pointed out in paper[4]is not essential,and several computation errors in[3]are corrected as well.     

Limit cycle analysis of a class of strongly nonlinear oscillation equations

The limit cycle of a class of strongly nonlinear oscillation equations of the form % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiGc9yrFr0xXdbba91rFfpec8Eeeu0x% Xdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9pgea0dXdar-Jb9hs% 0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaacaqabeaadaqaaqGaaO% qaaiqadwhagaWaaiabgUcaRmXvP5wqonvsaeHbbjxAHXgiofMCY92D% aGqbciab-DgaNjab-HcaOiaadwhacqWFPaqkcqWF9aqpcqaH1oqzca% WGMbGaaiikaiaadwhacaGGSaGabmyDayaacaGaaiykaaaa!50B8!\[\ddot u + g(u) = \varepsilon f(u,\dot u)\] is investigated by means of a modified version of the KBM method, where is a positive small parameter. The advantage of our method is its straightforwardness and effectiveness, which is suitable for the above equation, where g(u) need not be restricted to an odd function of u, provided that the reduced equation, corresponding to =0, has a periodic solution. A specific example is presented to demonstrate the validity and accuracy of our 09 method by comparing our results with numerical ones, which are in good agreement with each other even for relatively large .     

The auto-adjustable damping method for solving nonlinear equations

I.IntroductionTherearealotofnonlinearproblemsinengineering.Forexamples,theheattrallsfel'with'nonlinearradiation-convectionboundaryalldlimitedheatcapacity,thenonlinearhydrodynamicsandaerodynamics,thenonlillearstrtlcturaldynamicresponsesandtheIihe.Itisalmostimpossibletosolvethe11onlinearequationsbydirectmethodexceptfi)rsomespecialnonlinearequations.Thegeneralapproachislinearizingtheequationsandformingvariousiterationprocedures.Forthestronglynonlillearproblems,thesolutionobtainedintheitcrativep…     

Computer computation of the method of multiple scales—dirichlet problem for a class of system of nonlinear differential equations

The method of boundary layer with multiple scales and computer algebra were applied to study the asymptotic behavior of solution of boundary value problems for a class of system of nonlinear differential equations. The asymptotic expansions of solution were constructed. The remainders were estimated. And an example was analysed. It provides a new foreground for the application of the method of boundary layer with multiple scales. Contributed by Jiang Fu-ru, Original Member of Editorial Committe, AMM Biography: Xie La-bing (1976∼); Jiang Fu-ru(1927∼)     

Comparison of multiple scales and KBM methods for strongly nonlinear oscillators with slowly varying parameters

Two versions of multiple scales and Krylov–Bogoliubov–Mitropolsky (KBM) methods are extended to obtain the asymptotic solutions of a class of strongly nonlinear oscillators with slowly varying parameters. A comparison of the two methods is made to show that the multiple scales method is equivalent to the KBM method for the first order approximation. An example of pendulum with slowly varying length is given to illustrate the comparison of the two methods.     

Initial value problem for a class of fourth-order nonlinear wave equations

In this paper, existence and uniqueness of the generalized global solution and the classical global solution to the initial value problem for a class of fourth-order nonlinear wave equations are studied in the fractional order Sobolev space using the contraction mapping principle and the extension theorem. The sufficient conditions for the blow up of the solution to the initial value problem are given.     

CHAOTIC BEHAVIOUR OF FORCED OSCILLATOR CONTAINING A SQUARE NONLINEAR TERM ON PRINCIPAL RESONANCE CURVES

ＣＨＡＯＴＩＣＢＥＨＡＶＩＯＵＲＯＦＦＯＲＣＥＤＯＳＣＩＬＬＡＴＯＲＣＯＮＴＡＩＮＩＮＧＡＳＱＵＡＲＥＮＯＮＬＩＮＥＡＲＴＥＲＭＯＮＰＲＩＮＣＩＰＡＬＲＥＳＯＮＡＮＣＥＣＵＲＶＥＳＰｅｉＱｉｎ－ｙｕａｎ（裴钦元）（ＣｈａｎｇｓｈａＲａｉｌｗａｙＵｎｉ...     

求解变分型积分方程的一种新型数值方法——有限变分法

将作者提出的多虚拟裂纹扩展法(MVCE法)拓展为求解变分型积分方程问题的一种新型数值方法——有限变分法(FVM)。它的基本思想是,给定有限个(N个)局部变分模式,将所求解的未知量用适当的方法离散化,针对这N个局部变分模式列出N个方程,求解N个未知系数,从而求得未知量。单一未知变量FVM的最终方程组的系数矩阵通常是一个对称的窄带矩阵,对角元是大数,有很好的数值计算性能。用FVM求解了三维I型裂纹前缘的应力强度因子(SIF)分布。利用基于FVM的通用权函数法计算程序,可以高精度和高效率地求解表面力、体积力和温度载荷共同作用情况下三维裂纹前缘SIF的分布及其时间历程。FVM可以被推广到更广泛的领域,是一个求解变分型积分方程问题的普遍适用的新型数值方法。     

A field method for solving the equations of motion of nonholonomic systems

In this paper, the field method^[1] for solving the equations of motion of holonomic nonconservative systems is extended to nonholonomic systems with constant mass and with variable mass. Two examples are given to illustrate its application. The project supported by the National Natural Science Foundation of China.     

The first integral method for solving some important nonlinear partial differential equations

Exact solutions of some important nonlinear partial differential equations are obtained by using the first integral method. The efficiency of the method is demonstrated by applying it for two selected equations.     

A deterministic vortex method for solving the Navier-Stokes equations

In this paper, a new 2-D vortex method is developed, which treats the vorticity diffusion in a deterministical way. The Laplacian operator, which describes vorticity diffusion, is approximated by a contour integral. The numerical results of two model problems show that this method has a good accuracy. A primary error estimation is given, and the self-adaptive vortex blob and the boundary conditions are discussed. The project supported by the National Natural Science Foundation of China     

NONOSCILLATION RESULTS FOR A CLASS OF THIRD ORDER NONLINEAR DIFFERENTIAL EQUATIONS

ＩｎｔｒｏｄｕｃｔｉｏｎＦｉｎｄｉｎｇｓｕｆｆｉｃｉｅｎｔｃｏｎｄｉｔｉｏｎｓｆｏｒｎｏｎｏｓｃｉｌｌａｔｏｒｙｏｆｓｏｌｕｔｉｏｎｓｉｓａｐｒｏｂｌｅｍｏｆｇｅｎｅｒａｌｉｎｔｅｒｅｓｔｉｎｔｈｅｔｈｅｏｒｙｏｆｏｒｄｉｎａｒｙａｎｄｄｅｌａｙｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓ.Ｉｎｔｈｉｓｗｏｒｋｗｅｃｏｎｓｉｄｅｒｄｄｓ σ(ｓ) ｄｄｓ ｒ(ｓ) ｄｘｄｓ +ｑ(ｓ) ｄｘｄｓβ +ｐ(ｓ)ｘα =ｆ(ｓ) ,(1)ｗｈｅｒｅσ ,ｒ,ｑ,ｐａｎｄｆａｒｅｒｅａｌ＿ｖａｌｕｅｄｃｏｎｔｉｎｕｏｕｓｆｕｎｃｔｉ…     

Initial value problems for a class of nonlinear integro-partial differential equations

１ＰｒｏｂｌｅｍｓａｎｄＭａｉｎＲｅｓｕｌｔｓＩｎｔｈｉｓｐａｐｅｒ,ｗｅｓｔｕｄｙｔｈｅｎｏｎｌｉｎｅａｒｖｉｂｒａｔｉｏｎｓｏｆｉｎｆｉｎｉｔｅｒｏｄｓｗｉｔｈｖｉｓｃｏｅｌａｓｔｉｃｉｔｙ．Ｔｈｅｃｏｎｓｔｉｔｕｔｉｏｎｌａｗｏｆｔｈｅｒｏｄｓ...     

一类非线性发展方程的交替分组的并行数值算法

给出了一类非线性发展方程的交替分组4点格式,并得到该方法的无条件稳定性及具有并行本性兼顾的结果。能够适合在并行计算系统上使用。文中还进行了并行计算的数值实验。     

EXISTENCE THEOREMS OF EXTREMAL SOLUTIONS FOR A CLASS OF NONLINEAR INTEGRO－DIFFERENTIAL EQUATIONS

ＥＸＩＳＴＥＮＣＥＴＨＥＯＲＥＭＳＯＦＥＸＴＲＥＭＡＬＳＯＬＵＴＩＯＮＳＦＯＲＡＣＬＡＳＳＯＦＮＯＮＬＩＮＥＡＲＩＮＴＥＧＲＯ－ＤＩＦＦＥＲＥＮＴＩＡＬＥＱＵＡＴＩＯＮＳＳｏｎｇＧｕａｎｇｘｉｎｇ（宋光兴）（ＲｅｃｅｉｖｅｄＭａｙ４，１９９５；Ｃｏｍ...     

Duhamel项的精细积分方法在非线性微分方程数值求解中的应用

基于Duhamel项的精细积分方法,构造了几种求解非线性微分方程的数值算法。首先将非线性微分方程在形式上划分为线性部分和非线性部分,对非线性部分进行多项式近似,利用Duhamel积分矩阵,导出了非线性方程求解的一般格式。然后结合传统的数值积分技术,例如Adams线性多步法等,构造了基于精细积分方法的相应算法。本文算法利用了精细积分方法对线性部分求解高度精确的优点,大大提高了传统算法的数值精度和稳定性,尤其是对于刚性问题。本文构造的算法不需要对线性系统矩阵求逆,可以方便的考察不同的线性系统矩阵对算法性能的影响。数值算例验证了本文算法的有效性,并表明非线性系统的线性化矩阵作为线性部分是比较合理的选择。     

APPLICATION OF THE MODIFIED METHOD OF MULTIPLE SCALES TO THE BENDING PROBLEMS FOR CIRCULAR THIN PLATE AT VERY LARGE DEFLECTION AND THE ASYMPTOTICS OF SOLUTIONS(Ⅱ)

This paper is a continuation of part (Ⅰ), on the asymptotics behaviors of the series solutions investigated in (Ⅰ). The remainder terms of the series solutions are estimated by the maximum norm.     

Analysis of guided waves with a nonlinear boundary condition caused by internal resonance using the method of multiple scales

We theoretically investigated the cumulative nonlinear guided waves caused by internal resonance, using the method of multiple scales (MMS), which can construct better approximations to the solutions of perturbation problems. In this study, we consider nonlinearity only on the boundary instead of material nonlinearity or geometric nonlinearity. We showed nonlinear effects on the amplitudes of a lower mode and a higher mode depending on the propagation length. Also, we examined effects of wavenumber detuning from a phase matching condition of the two modes. If the wavenumber detuning is exactly equal to zero, the mechanical energy of the lower mode is transferred through nonlinear coupling to the energy of the higher mode, unilaterally. However, if a wavenumber detuning is not equal to zero, amplitude of the two modes change in a cyclic fashion during wave propagation. The amount of this amplitude variation and its cycle length are determined by the eigenfunctions of the two modes, the nonlinear parameter and the wavenumber detuning.     

Effect of gravity on nonlinear oscillations of a horizontal,immovable-end beam

A critical problem in designing large structures for space applications, such as space stations and parabolic antennas, is the limitation of testing these structures and their substructures on earth. These structures will exhibit very high flexibilities due to the small loads expected to be encountered in orbit. It has been reported in the literature that the gravitational sag effect under dead weight is of extreme importance during ground tests of space-station structural components [1–4]. An investigation of a horizontal, pinned-pinned beam with complete axial restraint and undergoing large-amplitude oscillations about the statically deflected position is presented here. This paper presents a solution for the frequency-amplitude relationship of the nonlinear free oscillations of a horizontal, immovable-end beam under the influence of gravity.The governing equation of motion used for the analysis is the Bernoulli-Euler type modified to include the effects of mid-plane stretching and gravity. Boundary conditions are simply supported such that at both ends there is no bending moment and no transverse and axial displacements. These boundary conditions give rise to an initial tension in the statically deflected shape. The displacement function consists of an assumed space mode using a simple sine function and unknown amplitude which is a function of time. This assumption provides for satisfaction of the boundary conditions and leads to an ordinary differential equation which is nonlinear, containing both quadratic and cubic functions of the amplitude. The perturbation method of multiple scales is used to provide an approximate solution for the fundamental frequency-amplitude relationship.Since the beam is initially deflected the small-amplitude fundamental natural frequency always increases relative to the free vibration situation provided in zero gravity. The nonlinear equation provides for interactions between frequency and amplitude in that both hardening and softening effects arise. The coefficient of the quadratic term in the nonlinear equation arises from the static (dead load) portion of the deflection. This quadratic term, depending upon its magnitude, introduces a softening effect that overcomes the hardening term (due to initial axial tension developed by deflection) for large slenderness ratios.For very large slender, immovable-end beams, the fundamental natural frequency is greater than that of beams without axial constraints undergoing small amplitude oscillations. This phenomenon is attributed to the stiffening effect of the statically-induced axial tension. However, the stiffening effect of axial tension in beams with slenderness ratios greater than approximately 392 undergoing large-amplitude symmetric-mode oscillations is overpowered by the presence of gravitational loading.Nomenclature A amplitude of the first harmonic - A ₁ cross-sectional area of beam - a(t) vibratory amplitude - E Young's modulus - g acceleration due to gravity - g ₁,g ₂,g ₃ constants defined in equations (8) - I area moment of inertia of cross-section - L length of beam - N axial tension force induced by gravitational loading