非线性问题和分叉问题及其数值方法

本文给出了一个一般性的分叉定义，说明了伪弧长算法在分叉计算中的应用，概述了静分叉点定位、用单纯形算法准确确定静分叉后各分叉解枝初始方向的算法，以及Ｈｏｐｆ分叉点定位和大范围连续追踪周期解轨道的数值方法。     

爆轰波强间断问题的伪弧长算法及其人为解验证

针对该问题开展了伪弧长数值算法研究,通过引入弧长参数,使网格按照一定的形式自适应移动,达到在强间断区域自动加密的效果,从而提高网格分辨率。基于伪弧长算法编写了二维程序,并对程序进行人为解方法验证。将伪弧长算法和直接有限体积法的数值结果进行对比,通过误差分析,显示出伪弧长算法能有效提高计算精度。最后将伪弧长算法应用于气相爆轰波在二维管道中的传播问题,研究了波阵面的捕捉效果和爆轰波胞格结构的形成过程。     

计算动力学中的伪弧长方法研究

动力学问题通常采用微分方程来描绘,但由于工程实际问题的复杂性,微分方程模型常伴随着解的不连续性、刚性或激波间断奇异性特点,传统方法很难求解,奇异性问题是计算动力学难点,同时也是国内外学者研究的热点.伪弧长数值算法是针对计算动力学中的奇异性问题所提出的,其基本思想为通过在解曲线上引入伪弧长参数,并增加一个约束方程,在伪弧长参数作用下,使得原始离散单元发生扭曲形变,从而达到消除或减弱奇异性的目的.本文首先介绍伪弧长方法求解定常对流-扩散方程的奇异性问题,并提出针对双曲守恒定律的局部伪弧长算法,其思想在于首先通过间断解的梯度变换来确定强间断所处位置,进而通过局部网格点重构以及数值修正来达到强间断处奇异性消除与降低的目的.针对高维问题,提出全局伪弧长方法,通过对整个计算区域内的网格点进行重构,使得所有网格点向奇异间断点处移动,从而降低间断点的影响域,达到降低奇异性的目的.重点讨论了三维全局伪弧长算法问题的计算难点,即三维空间网格扭曲大变形导致的数值算法不收敛,并提出在算法设计过程中采用分块重构与整体计算相结合的策略,实现了三维空间中的伪弧长数值算法,最后通过数值实验来验证伪弧长算法对于奇异性问题的有效性.     

双曲偏微分方程的局部伪弧长方法研究

重点研究了局部伪弧长方法在处理偏微分方程,尤其是双曲型偏微分方程出现激波间断的奇异性问题,对比分析了全局伪弧长方法空间转化的形式及其网格自适应的性质。为提高求解效率,提出了局部伪弧长方法,利用激波间断的性质,给出了判断奇异点位置以及模板选择的方法,涉及如何处理激波振荡,如何引入弧长参数,以及怎样求解间断等问题。通过数值算例验证了局部伪弧长在激波捕捉和追踪方面的可行性,通过比较局部伪弧长方法与Godunov方法处理不同初值条件的双曲问题,显示出局部伪弧长方法处理双曲偏微分方程的优越性,为伪弧长方法应用到物理问题奠定基础。     

含参数的非线性方程组的数值解法

本文对含参数λ的非线性方程组f(α,λ)=O提出了在位移-载荷空间R~n+1中以弧长为参数追踪解曲线c(s)的弧长增量法,并统一地处理了搜寻极值点(extre-mum point)与分叉点(bifurcation point)的问题,我们将此法成功地用于弹性旋转薄壳的非线性分析。     

爆炸冲击波强间断问题的高阶伪弧长算法研究

为了提高对冲击波强间断处的分辨率,通过引入弧长参数,使网格自适应地朝着间断处移动,并结合高精度WENO数值格式,进而达到了对大梯度物理量的高分辨率捕捉。针对网格移动造成的非均匀和非正交现象,通过坐标变换,使得计算过程在均匀正交的计算空间中进行。通过和有限体积下的数值结果对比,结合数值误差分析,可以看到高阶伪弧长数值算法不仅保证了高精度而且对间断的捕捉更加明显,在间断附近解的整体光滑性较好,网格的自适应移动使得解的奇异性得到了削弱,因此可以削弱高阶格式容易引起数值振荡这个缺点。最后采用高阶伪弧长算法计算了化学反应流问题,结果表明高阶伪弧长算法有着较快的收敛率,对于解决爆炸与冲击强间断问题有着较为明显的优势。     

受轴向冲击有限长弹性直杆中应力波引起的分叉问题

研究了有限长理想直杆受阶跃载荷作用的弹性动力屈曲问题。将直杆的屈问题归结为由轴向应力波传播而导致的杆的分叉问题，并考虑了应力波反射的影响。给出了分叉发生的临界条件并对具体实例进行了计算。最后对阶跃载荷及脉冲载荷对杆的动力屈曲的影响进行了讨论。     

结构后屈曲分析中的混合Newton-Lanczos算法

本文提出一种适于结构非线性后屈曲分析的混合Newton-Lanczos算法。与当前流行的弧长法不同,本文提出的算法采用传统的载荷增量法进行逐步求解,可求出给定载荷下的结构变形且适于任意外加载荷。对于临界载荷附近的迭代应用了Lanczos法求解方程及相应变载技巧。文中给出的若干数值计算结果表明了该算法在结构非线性后屈曲分析中的适用性。     

弹性地基上HDAJ接头桩的非线性稳定性分析

本文采用弧坐标首先建立了弹性地基中受轴向载荷作用的高柔性抗震拼接头桩(High Ductility Aseis-matic Joint Spliced Pile)的非线性数学模型，并假定土(基础)对桩基的反作用力服从Winkler模型；在此基础上对该模型进行了线性化，并得到HDAJ接头桩的临界载荷。最后根据分叉理论的观点和方法，讨论了HDAJ接头桩在临界载荷处的稳定性问题。研究结果表明HDAJ接头桩在临界载荷附近必发生分叉，且分叉解是唯一的，稳定的，并且给出了分叉解的渐近表达式。物理上，这表示HDAJ接头桩的平衡构形在临界载荷处必然发生改变，并且从一个稳定的平衡构形变化到另一个稳定的平衡构形。同时考察了土的液化对临界载荷的影响，说明液化的影响是非常明显的。当考虑土的液化时，桩基的临界载荷低于不考虑土的液化时桩基的临界载荷。     

有限变形的弧长算法

近些年来,人们提出了很多方法来解决结构静力非线性跟踪分析问题,其中,弧长算法应用最为广泛,但是,其中仍存在很多问题,本文针对梁板壳结构计算中的有限变形弧长算法,首先引入了将位移自由度与转动自由度分离提法,在此前提下对前人已有的处法加以改造,建立一个N＋1维的增量弧长方程组进行跟踪求解,本文同时引入了无量纲化,增量板长自动调节系数,奇点的判定准则,最终提出一个实用的弧长算法,本文在结尾将给出几个算例以显示该算法良好的跟踪性能。     

Continuation of limit cycles near saddle homoclinic points using splines in phase space

We present a new algorithm for continuation of limit cycles of autonomous systems as a system parameter is varied. The algorithm works in phase space with an ordered set of points on the limit cycle, along with spline interpolation. Currently popular algorithms in bifurcation analysis packages compute time-domain approximations of limit cycles using either shooting or collocation. The present approach seems useful for continuation near saddle homoclinic points, where it encounters a corner while time-domain methods essentially encounter a discontinuity (a relatively short period of rapid variation). Other phase space-based algorithms use rescaled arclength in place of time, but subsequently resemble the time-domain methods. Compared to these, we introduce additional freedom through a variable stretching of arclength based on local curvature, through the use of an auxiliary index-based variable. Several numerical examples are presented. Comparisons with results from the popular package, MATCONT, are favorable close to saddle homoclinic points.     

Bifurcations and chaotic dynamics in suspended cables under simultaneous parametric and external excitations

The bifurcations and chaotic dynamics of parametrically and externally excited suspended cables are investigated in this paper. The equations of motion governing such systems contain quadratic and cubic nonlinearities, which may result in two-to-one and one-to-one internal resonances. The Galerkin procedure is introduced to simplify the governing equations of motion to ordinary differential equations with two-degree-of-freedom. The case of one-to-one internal resonance between the modes of suspended cables, primary resonant excitation, and principal parametric excitation of suspended cables is considered. Using the method of multiple scales, a parametrically and externally excited system is transformed to the averaged equations. A pseudo arclength scheme is used to trace the branches of the equilibrium solutions and an investigation of the eigenvalues of the Jacobian matrix is used to assess their stability. The equilibrium solutions experience pitchfork, saddle-node, and Hopf bifurcations. A detailed bifurcation analysis of the dynamic (periodic and chaotic) solutions of the averaged equations is presented. Five branches of dynamic solutions are found. Three of these branches that emerge from two Hopf bifurcations and the other two are isolated. The two Hopf bifurcation points, one is supercritical Hopf bifurcation point and another is primary Hopf bifurcation point. The limit cycles undergo symmetry-breaking, cyclic-fold, and period-doubling bifurcations, whereas the chaotic attractors undergo attractor-merging, boundary crises. Simultaneous occurrence of the limit cycle and chaotic attractors, homoclinic orbits, homoclinic explosions and hyperchaos are also observed.     

A local pseudo arc-length method for hyperbolic conservation laws

A local pseudo arc-length method（LPALM）for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.     

静电激励MEMS微结构吸合电压尺寸效应研究

基于应变梯度弹性理论,研究了静电激励MEMS微结构吸合电压的尺寸效应。利用最小势能原理分别推导出含尺寸效应的一维梁模型和二维板模型的高阶控制方程。采用广义微分求积法和拟弧长算法对控制方程进行了数值求解。结果表明,随着结构尺寸的降低,新模型所预测的归一化的吸合电压呈非线性增长,表现出尺寸效应(特别是当结构尺寸与内禀常数在同一数量级时尺寸效应更加强烈);而相应的经典理论模型并不能预测此尺寸效应。两种新模型可视为相应经典理论的推广。本文有助于研究MEMS微结构的特性并对微结构的设计有潜在的应用价值。     

Non-linear vibrations of doubly curved shallow shells

Large amplitude (geometrically non-linear) vibrations of doubly curved shallow shells with rectangular base, simply supported at the four edges and subjected to harmonic excitation normal to the surface in the spectral neighbourhood of the fundamental mode are investigated. Two different non-linear strain-displacement relationships, from the Donnell's and Novozhilov's shell theories, are used to calculate the elastic strain energy. In-plane inertia and geometric imperfections are taken into account. The solution is obtained by Lagrangian approach. The non-linear equations of motion are studied by using (i) a code based on arclength continuation method that allows bifurcation analysis and (ii) direct time integration. Numerical results are compared to those available in the literature and convergence of the solution is shown. Interaction of modes having integer ratio among their natural frequencies, giving rise to internal resonances, is discussed. Shell stability under static and dynamic load is also investigated by using continuation method, bifurcation diagram from direct time integration and calculation of the Lyapunov exponents and Lyapunov dimension. Interesting phenomena such as (i) snap-through instability, (ii) subharmonic response, (iii) period doubling bifurcations and (iv) chaotic behaviour have been observed.     

Secondary buckling of an elastic column with spring-supports at clamped ends

Summary The postbuckling behavior of an elastic column with spring supports of equal stiffness of extensional type at both clamped ends is studied. Attention is focused on those of spring stiffnesses near the critical value at which, under axial load, the column becomes critical with respect to two buckling modes simultaneously. By using the Liapunov-Schmidt-Koiter approach, we show that there are precisely two secondary bifurcation points on each primary postbuckling state for the spring stiffness greater than the critical value. The bifurcation takes place at one of the two least buckling loads. The corresponding secondary postbuckling states connect all the secondary bifurcation points in a loop. For the spring stiffness less than the critical value, no secondary bifurcation occurs. Asymptotic expansions of the primary and secondary postbuckling states are constructed. The stability analysis indicates that the primary postbuckling state for the spring stiffness greater than the critical value is bifurcating from the first buckling load and becomes unstable from a stable state via the secondary bifurcation, i.e., secondary buckling occurs. Received 22 April 1997; accepted for publication 22 December 1997     

Nonlinear Nonplanar Dynamics of Parametrically Excited Cantilever Beams

The nonlinear nonplanar response of cantilever inextensional metallic beams to a principal parametric excitation of two of its flexural modes, one in each plane, is investigated. The lowest torsional frequencies of the beams considered are much larger than the frequencies of the excited modes so that the torsional inertia can be neglected. Using this condition as well as the inextensionality condition, we develop a Lagrangian whose variation leads to two integro-partial-differential equations governing the motions of the beams. The method of time-averaged Lagrangian is used to derive four first-order nonlinear ordinary-differential equations governing the modulation of the amplitudes and phases of the two interacting modes. These modulation equations exhibit symmetry properties. A pseudo arclength scheme is used to trace the branches of the equilibrium solutions and an investigation of the eigenvalues of the Jacobian matrix is used to assess their stability. The equilibrium solutions experience pitchfork, saddle-node, Hopf, and codimension-2 bifurcations. A detailed bifurcation analysis of the dynamic solutions of the modulation equations is presented. Five branches of dynamic (periodic and chaotic) solutions were found. Two of these branches emerge from two Hopf bifurcations and the other three are isolated. The limit cycles undergo symmetry-breaking, cyclic-fold, and period-doubling bifurcations, whereas the chaotic attractors undergo attractor-merging and boundary crises.     

复合材料层合板的二次屈曲和二次分枝点分析

为了研究复合材料层合板的二次分叉特性 ,利用能量变分原理和非线性几何方程建立了具有弹性约束的复合材料层合板在面内载荷作用下的非线性稳定性控制方程组。控制方程组用广义傅立叶级数法进行求解 ,并得到载荷 -挠度曲线。基于分叉理论中的 Lerray-Schaulder定理 ,采用小挠动法 ,直接导出了复合材料层合板的二次失稳方程。研究结果表明 ,非对称层板也可能存在分叉 ,弹性转动支持系数和铺层等因素对二次分叉有很重要的影响。随着弹性系数的增大 ,二次失稳载荷值与初次失稳载荷值之比下降     

Sensitivity analysis of an optimized bar-spring model with hill-top branching

Summary Characteristics of optimal solutions under nonlinear buckling constraints are investigated by using a bar-spring model. It is demonstrated that optimization under buckling constraints of a symmetric system often leads to a structure with hill-top branching, where a limit point and bifurcation points coincide. A general formulation is derived for imperfection sensitivity of the critical load factor corresponding to a hill-top branching point. It is shown that the critical load is not imperfection-sensitive even for the case where an asymmetric bifurcation point exists at the limit point.