求解具有弹性接头桩基动力学大变形的微分—代数方法

本文采用弧坐标首先建立了求解具有弹性接头的桩基大变形分析的非线性动力学微分方程,其中, 广义Winkler模型用来模拟土对桩基的抗力.其次,在空间域内应用微分求积单元法来离散非线性微分方程组,并给出了处理弹性接头处连接条件的微分求积单元公式,得到了时间域内的一组微分-代数方程,采用二阶向后差分来代替二阶时间导数离散微分-代数方程组,得到一组离散化的非线性代数方程,应用Newton-Raphson方法求解了该非线性代数方程组.最后给出了数值算例,得到了桩基在顶部处受到组合动载荷作用时的响应,考察了弹性接头的刚度、位置对桩基动力学行为的影响.     

基础-饱和土地基耦合系统分析的微分求积单元法

研究了基础-饱和土地基耦合系统的动力学特性.首先根据多孔介质理论,在小变形的假设下,分别建立了耦合系统中不可压流体饱和土地基和弹性基础的运动微分方程以及相应的边界条件,连接条件和初始条件;然后在空间域采用微分求积单元法对基础-饱和土地基的控制方程进行离散,并提供了正确处理耦合系统界面之间连接条件和间断性条件的方法,从而得到时间域内的一组代数-微分方程;接着运用隐式二阶向后差分格式处理了代数-微分方程组;最后利用牛顿迭代法数值求解了该方程组,得到了耦合系统的数值解,考察了所布节点数和参数对数值结果的影响.     

弹性或弹塑性土体中桩基的大变形分析

采用弧坐标,首先建立了位于弹性地基或弹塑性地基上并具有初始位移的桩基大变形行为的非线性微分方程组,并采用Winkeler模型来模拟地基对桩基的抗力;其次,应用微分求积方法离散非线性微分方程组,得到一组离散化的非线性代数方程,并给出了利用Newn-Raphson方法求解非线性代数方程的步骤;作为应用给出了数值算例,得到了桩顶受组合载荷作用时,变形后桩基的构形、弯矩和剪力,考察了土的弹性和弹塑性性质、桩基初始位移、荷载等参数对桩基力学行为的影响.最后将模型进行简化,得到了小变形理论的解析解,并比较了由大变形理论与小变形理论所得结果的差别.     

非线性转子系统周期响应的分析

基于时域的时间有限元法,将描述转子系统动力学特征的非线性微分方程组离散成一组非线性代数方程,然后应用吴消去法的特征列思维对所得到的非线性代数方程组进行降维求解,进而得到待求节点位移响应的解形式,并据此对一具有非线性支撑的柔性Jeffcott转子模型响应的性质进行了分析。     

微分求积方法对于桩基大变形分析的应用

本文基于大变形的理论,采用弧坐标首先建立了具有初始位移的桩基的非线性数学模型,一组强非线性的微分-积分方程,其中,地基的抗力采用了Winkeler模型;其次,引入变数变换将微分-积分方程转化为一组非线性微分方程,并用微分求积方法离散了方程组,得到一组离散化的非线性代数方程;最后用Newton-Raphson迭代方法对离散化方程进行了求解,得到了桩基变形前后的构形、弯矩和剪力.计算中选取了两种不同类型的初始位移,并考察了它们对桩基大变形力学行为的影响.     

单桩横向非线性惯性响应简化分析计算方法

根据已建立的地基土-单桩系统横向非线性动力相互作用的简化分析模型,在动力Winkler模型的基础上,研究单桩横向非线性惯性响应。简化模型中．地基土对桩轴的横向非线性作用力．由非线性迟滞退化弹簧模拟,地基土在运动过程中的能量耗散,由与弹簧并联的粘壶模拟。利用有限元空间离散技术和Newmark-β法,讨论了横向惯性荷载作用下．单桩运动方程在时域中的求解方法。计算结果表明,该模型能够反映单桩非线性惯性响应的主要特征。能够揭示各种非线性因素,如：桩周土的迟滞效应、退化、屈服,以及土-桩界面的分离等,对单桩惯性响应的影响。计算结果与试验结果相当吻合。     

一种求解双曲方程新的有限差分算法

本文提出一种适于求解一阶双向系统的新的差分格式。它的建立方法是：将所要求解的方程与解的空间导数所满足的微分方程同时离散化，然后再通过插值函数构成封闭的离散变量代数方程。在线性情况下的误差分析表明：该格式的幅值与位相误差均小于常用的一、二阶差分格式；当其应用于非线性气动方程求解时，基本上可以消除数值扩散与振荡这两种非正常现象。     

一种求解双曲方程新有有限差分算法

本文提出一种适于求解一阶双曲系统的新的差分格式。它的建立方法是：将所要求解的方程与解的空间导数所满足的微分方程同时离散化，然后再通过插值函数构成封闭的离散变量代数方程。在线性情况下的误差分析表明：该格式的幅值与位相误差均小于常用的一，二阶差分格式，当其应于非线性气动方程求解时，基本上可以消除数值扩散与振荡这两种非正常现象。     

任意圈层径向非均质土中桩的纵向振动特性

研究三维轴对称条件下径向非均质土中桩的纵向振动. 首先将桩周土体沿径向分为任意圈层来考虑土体的径向非均质性,每个圈层土体均质; 然后结合边界条件和相邻圈层土体之间接触面上位移和应力连续条件,对任意圈层土体动力平衡方程由外而内逐层求解,进而利用桩-土完全耦合条件求解桩动力平衡方程,得到桩顶的频域响应解析解和时域响应半解析解; 最后通过对土体主要控制参数的研究,得出了土体径向非均匀性对桩-土动力响应的影响.      

多体系统动力学约束Hamilton方程多步块方法

针对多体系统动力学Hamilton体系正则微分-代数方程,基于等距节点、Chebyshev节点、Legendre节点等非等距节点,在时间区间上建立r-级多步块求解格式,得到单步区间各节点的非线性代数方程,求解过程利用Newton迭代。以双连杆机械臂系统为例,使用r-级多步块格式进行数值仿真,通过改变步长、仿真时间以及节点数分别对指标-1、-2、-3的Hamilton系统正则微分-代数方程进行数值验证。数值结果表明,本文方法稳定性好、精度高、计算效率高,能很好保持Hamilton能量误差以及各级约束误差,特别是对于长时间仿真,Hamilton能量误差不会产生漂移,适合长时间多体系统动力学仿真。     

Differential-algebraic approach to large deformation analysis of frame structures subjected to dynamic loads

A nonlinear mathematical model for the analysis of large deformation of frame structures with discontinuity conditions and initial displacements, subject to dynamic loads is formulated with arc-coordinates. The differential quadrature element method （DQEM） is then applied to discretize the nonlinear mathematical model in the spatial domain, An effective method is presented to deal with discontinuity conditions of multivariables in the application of DQEM. A set of DQEM discretization equations are obtained, which are a set of nonlinear differential-algebraic equations with singularity in the time domain. This paper also presents a method to solve nonlinear differential-algebra equations. As application, static and dynamical analyses of large deformation of frames and combined frame structures, subjected to concentrated and distributed forces, are presented. The obtained results are compared with those in the literatures. Numerical results show that the proposed method is general, and effective in dealing with disconti- nuity conditions of multi-variables and solving differential-algebraic equations. It requires only a small number of nodes and has low computation complexity with high precision and a good convergence property.     

Vibration of fluid-conveying pipe with nonlinear supports at both ends

The axial fluid-induced vibration of pipes is very widespread in engineering applications. The nonlinear forced vibration of a viscoelastic fluid-conveying pipe with nonlinear supports at both ends is investigated. The multi-scale method combined with the modal revision method is formulated for the fluid-conveying pipe system with nonlinear boundary conditions. The governing equations and the nonlinear boundary conditions are rescaled simultaneously as linear inhomogeneous equations and linear inhomogeneous boundary conditions on different time-scales. The modal revision method is used to transform the linear inhomogeneous boundary problem into a linear homogeneous boundary problem. The differential quadrature element method (DQEM) is used to verify the approximate analytical results. The results show good agreement between these two methods. A detailed analysis of the boundary nonlinearity is also presented. The obtained results demonstrate that the boundary nonlinearities have a significant effect on the dynamic characteristics of the fluid-conveying pipe, and can lead to significant differences in the dynamic responses of the pipe system.

     

Differential quadrature element method for static analysis of Reissner–Mindlin polar plates

In this paper, a new numerical technique, the differential quadrature element method (DQEM) , has been developed for static analysis of the two-dimensional polar Reissner–Mindlin plate in the polar coordinate system by integrating the domain decomposition method (DDM) with the differential quadrature method (DQM) . The detailed formulations for the sectorial DQEM plate bending element and the compatibility conditions between each element are presented. The convergence properties and the accuracy of the DQEM for bending of thick polar plates are investigated through a number of numerical computations. Consequently, the DQEM has been successfully applied to analyze several annular sector plates with discontinuous loading and boundary conditions and cutouts to illustrate the simplicity and flexibility of this method for solving Reissner–Mindlin plates in polar coordinate system which are not solvable directly using the differential quadrature method. The numerical results are verified by the existing exact solutions or the FEM solutions obtained using the software package ANSYS (Version 5.3) .     

Nonlinear dynamical characteristics of piles under horizontal vibration

The pile-soil system is regarded as a visco-elastic half-space embedded pile. Based on the method of continuum mechanics, a nonlinear mathematical model of pile-soil interaction was established-a coupling nonlinear boundary value problem. Under the case of horizontal vibration, the nonlinearly dynamical characteristics of pile applying the axis force were studied in horizontal direction in frequency domain. The effects of parameters, especially the axis force on the stiffness were studied in detail. The numerical results suggest that it is possible that the pile applying an axis force will lose its stability. So,the effect of the axis force on the pile is considered.     

A unified approach to direct and inverse boundary layer solutions

A unified approach is presented for solving the two-dimensional incompressible boundary layer equations. Solutions are obtained for direct and inverse options using the same equation formulation by a simple interchange of boundary conditions. A modified form of the mechul function scheme obtains inverse solutions with specification of transformed wall shear, skin friction coefficient or displacement thickness distributions. Direct solutions may be obtained without altering the block tridiagonal structure of the system by simply requiring no corrections on the streamwise pressure gradient parameter. Fourth-order spline discretization approximates normal derivatives with two- and three-point backward differences approximating streamwise derivatives, yielding a fully implicit solution method. The resulting spline/finite difference equations are solved by Newton-Raphson iteration together with partial pivoting. The results of the study demonstrate the importance of proper linearization of all equations. The successful use of spline discretization is also tied to the use of strong two-point boundary conditions at the wall for cases involving reversed flow. Numerical solutions are presented for several non-similar flows and compared with published results.     

同震断层三维裂纹扩展波动时域超奇异积分法

应用波动时域超奇异积分法将P波、S波和磁电热弹多场耦合作用下同震断层任意形状三维裂纹扩展问题转化为求解以广义位移间断率为未知函数的超奇异积分方程组问题;定义了广义应力强度因子,得到裂纹前沿广义奇异应力增量解析表达式;应用波动时域有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以三维矩形裂纹扩展问题为例,通过典型算例,研究了广义应力强度因子随裂纹位置变化规律;分析了同震断层裂纹扩展中力、磁、电场辐射规律.      

基于弹塑性模型的群桩-土-结构动力非线性数值分析

提出一种群桩-土弹塑性模型,结合动力文克尔理论,推导出了与桩(筏)-土属性及SSI体系频率相关的各项弹簧-阻尼单元动力阻抗,建立了三维框架土-结构相互作用有限元简化分析模型。针对不同地震激励,在不同桩-土条件下对模型进行了动力非线性时程分析,结果表明,在某些地震动和土-基础条件下,上部结构非线性作用效应结果可能大于固基假定情形,且桩-土弹塑性模型对上部结构柔弱层位置产生影响。应用本文简化方法可以快速、较准确有效地进行复杂的上下部结构动力时程分析及抗震评估。     

Hybrid finite element formulations for elastodynamic analysis in the frequency domain

Three alternative sets of hybrid formulations to solve linear elastodynamic problems by the finite element method are presented. They are termed hybrid–mixed, hybrid and hybrid–Trefftz and differ essentially on the field conditions that the approximation functions are constrained to satisfy locally. Two models, namely the displacement and the stress models, are obtained for each formulation depending on whether the tractions or the boundary displacements are the field chosen to implement interelement continuity. A Fourier time discretization is used to uncouple the solving system in the frequency domain. The basic space discretization criterion is implemented directly on the fundamental relations of elastodynamics and used to derive the stress and displacement models of the hybrid–mixed formulation. The hybrid and hybrid–Trefftz formulations are presented in sequence as the variants of the hybrid–mixed formulation obtained by progressively increasing the constraints on the approximation bases. Numerical implementation aspects are briefly discussed and the performance of the finite element models is illustrated with numerical applications.     

Stability and vibration of empty and fluid-filled circular cylindrical shells under static and periodic axial loads

In the present study, the dynamic stability of simply supported, circular cylindrical shells subjected to dynamic axial loads is analysed. Geometric nonlinearities due to finite-amplitude shell motion are considered by using the Donnell’s nonlinear shallow-shell theory. The effect of structural damping is taken into account. A discretization method based on a series expansion involving a relatively large number of linear modes, including axisymmetric and asymmetric modes, and on the Galerkin procedure is developed. Axisymmetric modes are included; indeed, they are essential in simulating the inward deflection of the mean oscillation with respect to the equilibrium position and in describing the axisymmetric deflection due to axial loads. A finite length, simply supported shell is considered; the boundary conditions are satisfied, including the contribution of external axial loads acting at the shell edges. The effect of a contained liquid is investigated. The linear dynamic stability and nonlinear response are analysed by using continuation techniques and direct simulations.     

Nonlinear bending analysis of ring-stiffened circular and annular general angle-ply laminated plates with various boundary conditions

In this study, elastic large deflection analysis of axisymmetric ring-stiffened circular and annular general angle-ply laminated plates subjected to transverse uniform load is studied. Based on first order shear deformation theory (FSDT) and large deflection von-Karman relations, the governing equations are derived. The dynamic relaxation (DR) method in conjunction with the central finite difference discretization technique is used to solve the nonlinear equilibrium equations. A detailed parametric study is carried out to investigate the influences of plate thicknesses, stiffener width, stiffener depth, fiber orientation, stacking sequence and different types of boundary conditions. Also, some linear and nonlinear analysis is provided to consider the effect of nonlinearity on the results.