时滞对于参数激励系统周期运动的影响

本文以结构物承受周期激励地震波和具有时滞的弹性地基为背景，研究时滞对于结构物振动响应的影响规律。问胚的数学模型是一个非线性参数激励时滞系统，采用中心流形和平均法得到Hopf分岔方程，研究结果表明时滞量对结构物的地震响应有重要影响，它可使响应振幅增大，意味着结构物遭破坏的危险性增大，同时，结果也表明可以通过阻尼控制器克服时滞引起的危险性。最后，对于不同的时滞量，比较理论分析与数值模拟结果基本吻合，说明了本文结果的可靠性。     

非线性Winkler地基上矩形薄板在移动荷载作用下的非线性动力分析

分析了非线性Winkler地基上矩形薄板在车辆移动荷载作用下的非线性动力特性。考虑地基反力的存在,基于Hamilton能量变分原理,建立了车辆、板、地基耦合系统非线性振动的控制微分方程;并将方程进行了量纲归一化处理,构造了满足周边自由矩形薄板全部边界条件的试探函数;运用伽辽金法和谐波平衡法对耦合系统控制方程进行了求解,讨论了板参数、地基参数、车辆系统参数等变化对耦合系统板振动幅频曲线的影响。结果表明:该耦合系统振动的频率都随板振幅的增大而增大;当板振动的幅值一定时,系统振动频率随着板厚、地基反应模量、车辆运行速度、车体刚度的增大而增大,但随着车体质量的增大而减小。因此,适当增加地基的反应模量可优化地基板的振动,并且从行车舒适性角度考虑,适当控制车速和车体刚度是有益的。     

横向荷载作用下弹性地基上梁的主共振分析

基于Winkler地基模型及Euler-Bernoulli梁理论,建立了弹性地基上有限长梁的非线性运动方程.运用Galerkin方法对运动方程进行一阶模态截断,并利用多尺度法求得该系统主共振的一阶近似解.分析了长细比、地基刚度、外激励幅值和阻尼系数等参数对系统主共振幅频响应的影响,然后通过与非共振硬激励情况对比分析主共振对其动力响应的影响.结果表明:主共振幅频响应存在跳跃和滞后现象;阻尼对主共振响应有抑制作用;主共振显著增大系统稳态动力响应位移.     

具有光滑中心的波纹圆板的非线性振动

基于各向同性及各向异性圆板的大挠度理论,研究了具有光滑中心的波纹圆板的非线性自由振动。以波纹板的中心最大振幅为摄动参数,采用摄动变分法,得出了波纹板的二次近似非线性固有频率。     

地基土—多自由度非线性结构相互作用体系的非平稳随机动力响应分析

本文应用一个改进的等效线性化方法探讨了地基土——多自由度非线性结构相互作用体系在随机地震荷载作用下的非平稳动力响应问题。对非线性恢复力模型采用一个具有三线性滞回特性的非线性系统。最后将在平稳Gauss过滤白噪声激励下的非平稳随机响应与Monte-carlo法的统计结果进行了比较,得出了较为满意的结果。     

矩形移动荷载作用下饱和-非饱和土双层地基的动力响应分析

考虑地基为饱和-非饱和土双层半空间,利用连续介质力学和多相孔隙介质理论,构建双层地基的统一动力控制方程并进行耦合求解。利用Dirac-delta函数和Heaviside阶跃函数将矩形移动荷载作用描述为时间和空间坐标的解析函数,将荷载函数代入地基动力控制方程,采用三重Fourier变换以及降阶法进行求解,并对推导结果进行退化验证及对饱和-非饱和土双层地基的动力响应进行分析。研究表明,当荷载移动速度小于瑞利波速时,竖向振动峰值很小,振幅随速度的增大发生小幅增涨,但当荷载速度达到瑞利波速时,竖向振动发生激增;随着速度进一步增大,竖向位移多次出现峰点。非饱和土的饱和度及土层厚度也对地基振幅存在显著影响。     

一种新的覆冰导线舞动非线性有限元分析方法

本文基于Hamilton原理,建立了耦合三个平动自由度和一个扭转自由度以及偏心覆冰作用下导线的非线性舞动方程,提出了一种新的覆冰导线舞动分析非线性有限元模型,将相邻跨导线和绝缘子等效为线性弹簧单元,同时考虑覆冰导线非线性气动力和大幅舞动的几何非线性,采用Newmark-β时间积分结合修正的Newton-Raphson非线性迭代求解舞动有限元方程。本文数值解与D形覆冰导线舞动实测值无论在振幅还是频率方面均吻合得非常好,证明了本文新提出的导线舞动有限元模型的准确性。本文研究还表明：舞动是一种以上下运动为主的低频振动,通常发生在一阶上下自振频率附近,其振幅和频率由输电线路的物理参数和风载荷唯一确定,与初始运动状态无关。     

机床非线性颤振的描述函数分析

本文采用描述函数方法分析机床的非线性颤振.以描述函数刻划切削过程复杂的非线性特性,并将非线性振动分析所依赖的复平面推广到三维空间进行稳定性分析。从而求出颤振振幅和颤振频率。该法不仅可揭示切削过程等效刚度系数和等效阻尼系数的变化,而且还可考察各切削参数对颤振频率和振幅的影响.结果表明,描述函数方法是分析机床颤振这个极其复杂的非线性系统的一个有效的手段.     

样条边界元法分析圆形基础板非线性非光滑接触

本文用样条边界元法分析轴对称圆形基础板的非线性非光滑接触。除事先给出外力与地基表面位移的对应关系外,地基的模型、板的厚薄(Reissner/Kirchhoff型)和外荷载沿径向的分布均可是任意的。文中将地基反力看作板的外载,对变量沿径向采用不等距祥条插值,并用无地基效应的板基本解。算例表明本文方法只要用少量自由度,便可得到令人满意的结果。     

轴对称弹性圆锥扁壳非线性自由振动直接摄动解法

本文采用弹性圆锥扁壳中心无量纲振幅和壳体母线的倾角为参数，将挠度、应力函数的导数以及自由振动频率展开为双参量的幂级数形式．用直接摄动法获得各级递推线性偏微分方程．应用变分法求得各级递推方程的近似解答．从而给出弹性圆锥扁壳非线性自由振动频率的基本公式。     

Velocity–pressure integrated versus penalty finite element methods for high-Reynolds-number flows

Velocity–pressure integrated and consistent penalty finite element computations of high-Reynolds-number laminar flows are presented. In both methods the pressure has been interpolated using linear shape functions for a triangular element which is contained inside the biquadratic flow element. It has been shown previously that the pressure interpolation method, when used in conjunction with the velocity-pressure integrated method, yields accurate computational results for high-Reynolds-number flows. It is shown in this paper that use of the same pressure interpolation method in the consistent penalty finite element method yields computational results which are comparable to those of the velocity–pressure integrated method for both the velocity and the pressure fields. Accuracy of the two finite element methods has been demonstrated by comparing the computational results with available experimental data and/or fine grid finite difference computational results. Advantages and disadvantages of the two finite element methods are discussed on the basis of accuracy and convergence nature. Example problems considered include a lid-driven cavity flow of Reynolds number 10 000, a laminar backward-facing step flow and a laminar flow through a nest of cylinders.     

A coupled immersed boundary‐lattice Boltzmann method with smoothed point interpolation method for fluid‐structure interaction problems

The immersed boundary‐lattice Boltzmann method has been verified to be an effective tool for fluid‐structure interaction simulation associated with thin and flexible bodies. The newly developed smoothed point interpolation method (S‐PIM) can handle the largely deformable solids owing to its softened model stiffness and insensitivity to mesh distortion. In this work, a novel coupled method has been proposed by combining the immersed boundary‐lattice Boltzmann method with the S‐PIM for fluid‐structure interaction problems with large‐displacement solids. The proposed method preserves the simplicity of the lattice Boltzmann method for fluid solvers, utilizes the S‐PIM to establish the realistic constitutive laws for nonlinear solids, and avoids mesh regeneration based on the frame of the immersed boundary method. Both two‐ and three‐dimensional numerical examples have been carried out to validate the accuracy, convergence, and stability of the proposed method in consideration of comparative results with referenced solutions.     

Overlapping grids and multigrid methods for three-dimensional unsteady flow calculations in IC engines

A new computational methodology with emphasis on using an overlapping grid technique and a multigrid method has been developed. The main feature of the present overlapping-grid system is of extended flexibility to deal with three-dimensional complex multicomponent geometries. The multigrid method is incorporated into this technique to accelerate the convergence of the numerical solution. The current scheme has been applied for computations of the laminar flows in the multicomponent configuration of internal combusion engines. The flow is governed by three-dimensional, time-dependent, incompressible Navier-Stokes equations with the continuity equation. A time-independent grid system is constructed for the moving boundary, i.e. the moving piston in the engine. This grid system is entirely different from others for the same problem in previous works. The performance of the present method has been validated by comparing the results with those from an equivalent, single-grid method and those from experiments. In addition, the flexibility and potential of the method has been demonstrated by calculating several cases which would be very difficult to be handled by other schemes.     

Thermal solution of cylindrical composite systems using meshless method

This paper presents the thermal solution of cylindrical composite systems using meshless element free Galerkin (EFG) method. The EFG method utilizes the moving least square approximants, which are constructed by using a weight function, a basis function and a set of non-constant coefficients to approximate the unknown function of temperature. Dirichlet (essential) boundary conditions have been enforced using Lagrange multiplier and penalty methods. Existing rational weight function has been modified and used in the present analysis. MATLAB codes have been developed to obtain the numerical solution. The EFG results have been obtained using cubicspline, quarticspline, Gaussian, quadratic, hyperbolic, exponential, rational and cosine weight functions for a model problem. The results obtained using different EFG weight functions are also compared with those obtained by finite element method. The effect of scaling and penalty parameters has also been studied in detail.     

超音速粘性流动的SUPG有限元数值解法

本文构造了准简化 N-S 方程组的 SUPG(Streamline Upwind/Petrov-Galerk-in)加权剩余式,并利用该方法对 Burgers 方程、无粘性激波反射问题、以及超音速平板和压缩拐角的层流流动作了数值求解。计算结果表明,本文方法是精确、收敛和稳定的。     

金属基纳米复合材料等效弹性模量的均匀化方法数值模拟

均匀化理论利用位移场双尺度渐近展开建立有限元列式，本文将其与有限元通用程序相结合，应用于金属基复合材料的弹性本构数值模拟。通过对不同尺度增强相金属基复合材料等效模量的数值模拟，考察了均匀化方法的适用情况。数值计算结果表明，对常规尺度增强相金属基复合材料，均匀化方法可以较准确地预测其等效弹性模量；对纳米增强相金属基复合材料，该方法仍可给出较好的预测，但存在某种程度的系统偏差。通过对纳米尺度增强机理的分析讨论，认为纳米增强相与基体材料问的界面效应可能有别于连续介质假设，指出可以考虑采用离散原子-连续介质耦合模型改进数值模拟结果。     

The determination of an equivalent elastic modulus for bodies with cracks by boundary element method

The equivalent elastic modulus of cracked bodies with orderly distributed cracks was computed with the boundary element method. A practical self-consistent scheme has been proposed in consideration of the mutual interaction effects of the cracks. The influence of friction coefficients and orientation of cracks has been investigated. Some computational examples have been given, and the results show that the proposed method is adequate and the scheme is efficient.This project is supported by the National Natural Science Foundation of China.     

Nonlinear bending of shallow spherical shells on the Pasternak foundations

Based on Kármán’s nonlinear fundamental differential equations,the new approach,which combines modified iteration method with Galerkin’s one,has been put forward tosolve nonlinear bending of shallow spherical shell with concave base and clamped edges onthe Pasternak foundation under uniform loads in this paper.Mathematical expression ofload-deflection has been given;furthermore,results obtained are in good agreement withexistent ones.     

空间刚架结构的广义逆力法

在结构计算中，根据算法中所采用的基本未知量的不同．结构分析方法可以分为力法、位移法和混合法。其中位移法由于适宜计算机处理而在结构计算领域得到了广泛的应用．经典力法相比之下应用就远不如位移法普遍，虽然力法本身在力学上有其独特的优势。广义逆矩阵做为一种较新的数学工具，自二十世纪五十年代诞生以来正日益表现出越来越旺盛的生命力。广义逆力法就是一种基于力法和广义逆矩阵理论的新的迭代算法。这种算法是一种完全适合计算机处理的力法方法。该算法的思路以及对于求解线弹性空间刚架结构问题的具体公式均在文中给出并给出了算例。从算例计算结果可以看到广义逆力法有着较好的计算效率和计算精度。该算法的提出为力法在计算机计算领域的应用开拓了新的发展空间。该算法在材料非线性问题和结构并行计算方面也有着较好的应用前景。     

Numerical methods for theoretical analysis of foil bearing dynamics

The bearing under consideration is a self-acting, gas bearing with a flexibly supported foil. In order to conduct the theoretical analysis, three systems were identified: a rotor, a gas film and a flexible structure. Each of these elements has been mathematically analysed and analytical equations have been formulated. It was found, that the gas flow in the bearing can be described by the Reynolds equation, whereas a spring-damper model was selected for the structural analysis. The Reynolds equation is a differential equation the exact solution to which is unknown. The work describes the finite difference method in detail, where the partial derivatives in the Reynolds equation are replaced by a system of algebraic equations. In order to solve the resulting system, the Alternating Direction Implicit method (ADI) was used. Based on that and on the analysis of interactions between the three systems the equations have been computed. The accuracy of methods has been verified by means of series of numerical tests. Chosen results are described in this paper.