P型有限元线法分析板弯曲问题

Ｐ型有限元线法分析板弯曲问题袁驷（北京清华大学，１０００８４）关键词有限元线法，板弯曲问题，中厚板理论互方法简介有限元线法（ＦＥＭＯＬ）是一种新型的以常微分方程（ＯＤＥ）求解器为支撑软件的半解析数值方法，其基本理论可参见文【１］【２］及有关文献，这里...     

薄板弯曲单元的精化矩阵

本文提出一个改进薄板弯曲单元计算精度的一般性方法：用精化矩阵对原有的单元刚度阵进行修正。实践表明，这种方法可以有效地改变单元的刚柔特性，从而达到提高单元精度的目的。文中证明了，单刚中加进精化矩阵不会破坏原来单元的收敛性质。     

细胞自动机方法解小挠度薄板弯曲问题

固体力学中的细胞自动机方法是一种分析固体力学问题的新方法,它将结构的静平衡认为是在力和位移边界条件下,物体元胞问的自组织过程,它不需要象现有有限元那样形成整体刚度矩阵和求解整体平衡方程,对计算机容量要求低,可望成为分析大型结构的有效方法,且特别适宜于并行计算.为拓展其应用范围,本文将其应用于小挠度弹性薄板弯曲的静力计算以检验其可行性.文中给出了详细的计算步骤,以实际算例为基础,探讨了该方法的有效性,并对其优点和目前存在的问题进行了讨论.     

一种新型的薄板弯曲单元

本文将薄板弯曲问题的4阶微分方程,化为两个2阶微分方程,并提出了一个与之相应的三变量广义变分原理,然后从这个原理出发,进行有限元求解。该方法不但成功地将C~1连续性条件降为C~0问题,而且降低了板单元的自由度数目,同时单元的精度也是令人满意的。     

矩形薄板弹性弯曲的精确解析解法

本文的解析解法和常用的迭加法原则上是相同的,但本文的方法简单,易掌握,能求解任意荷载作用,任意边界条件及板结构问题。     

薄板弯曲问题的神经网络方法

为发展神经网络方法在求解薄板弯曲问题中的应用,基于Kirchhoff板理论,提出一种采用全连接层求解薄板弯曲四阶偏微分控制方程的神经网络方法.首先在求解域、边界中随机生成数据点作为神经网络输入层的参数,由前向传播系统求出预测解;其次计算预测解在域内及边界处的误差,利用反向传播系统优化神经网络系统的计算参数;最后,不断训...     

薄板弹塑性弯曲的X样条有限条方法

薄板弹塑性弯曲的Ｘ样条有限条方法叶金铎，杨海元（天津大学冶金分校，３００４００）（天津大学，３０００７２）关键词Ｘ样条有限条，弹塑性，分层方案，不分展方案，不规则区域１前言采用有限元法解弹塑性问题，计算量大、计算费用高，采用传统有限条法解弹塑性问题又...     

P-version hybrid analytical finite element method for plate bending problems

Two sets of trial functions with different variables are constructed for the admissible space of the finite element analysis. The trial functions satisfy the equilibrium differential equation inside elements, while the deflections and rotations on the edges of the elements are approximated by the Peano hierarchical interpolation functions. Then, a generalized variational principle is applied to set up the p-version hybrid analytical finite element method for plate bending problems. The accuracy of finite element computation can be improved by increasing the order of the interpolation polynomials with fixed mesh. In the finite element formulation, to obtain the stiffness matrices and the load vectors, it is only necessary to perform quadrature over the edges of the elements. These matrices and vectors possess an embedding structure. The conformability between the elements can be controlled automatically.This work is supported by the Natural Science Foundation of China and the Aeronautical Science Foundation of China.     

斜交板梁桥的有限单元算法

借助斜坐标与直角坐标的关系,构造了斜交板梁单元的位移函数,由此提出了斜交板梁桥的实用有限单元算法。     

薄板动态热耦合弯曲问题的拉普拉斯变换有限元法

以薄板动力热耦合弯曲问题的变分原理为基础，提出了一种解决该问题的拉普拉斯变换有限元法，该方法在拉氏域内建立和求解有限元方程，再将结果数值逆变换至时域，得到各给定时刻的解答，周边绝热简支方板的热冲击问题算例说明方法是成功的，该方法为研究薄板此类问题的可能性，其计算结果可供工程设计参考。     

A NEW IHIERARCHICAL BOUNDARY ELEMENT METHOD AND ITS ADAPTIVE PROCESSES FOR PLATE BENDING PROBLEM

In this paper a new hierarchical boundary element method is introducedfor solving the problem of plate bending.Exact solutions of the governing equationsare used inside the domain together with independent deflections and rotations on theboundary.A generalized variational principle is employed to achieve the boundaryelement formulation.Further,the adaptive processes of the method are also discussed.By virtue of the error estimate technique proposed by Zienkiewicz et al.areasonable error indicator and adaptive scheme are suggested.Numerical examplesillustrate the high accuracy of the new elements and show the excellent efficiency ofthe adaptive computation described in this paper.     

多变量样条有限元法

本文提出了一种基于胡海昌－鹫津久一郎三类变量广义变分原理［１，２］的多变量样条有限元法，文中应用乘积型二元三次Ｂ样条插值函数来构造板壳的广义变量场函数，由三类变量广义变分原理来建立多变理样有限元模型。在计算种种场变量时，既不必求导，又无需用应力应变关系式，可直接算得其结果，因而对各种场变量均有足够的精度，文中，还解算了振动与稳定特征值问题，均得到了精度较好的数值结果。     

HIERARCHICAL STOCHASTIC FINITE ELEMENT METHOD FOR STRUCTURAL ANALYSIS

In this paper, the hierarchical approach is adopted for series representation of the stochastic nodal displacement vector using the hierarchical basis vectors, while the Karhunen-Lòeve series expansion technique is employed to discretize the random field into a set of random variables. A set of hierarchical basis vectors are defined to approximate the stochastic response quantities. The stochastic variational principle instead of the projection scheme is adopted to develop a hierarchical stochastic finite element method (HSFEM) for stochastic structures under stochastic loads. Simplified expressions of coefficients of governing equations and the first two statistical moments of the response quantities in the schemes of the HSFEM are developed, so that the time consumed for computation can be greatly reduced. Investigation in this paper suggests that the HSFEM yields a series of stiffness equations with similar dimensionality as the perturbation stochastic finite element method (PSFEM). Two examples are presented for numerical study on the performance of the HSFEM in elastic structural problems with stochastic Young’s Modulus and external loads. Results show that the proposed method can achieve higher accuracy than the PSFEM for cases with large coefficients of variation, and yield results agreeing well with those obtained by the Monte Carlo simulation (MCS).     

牛顿迭代一致性算法及其在板弹塑性有限元分析中的应用

本文简略讨论了有限载荷增量弹塑性有限元分析中传统切线刚度法丧失精度和牛顿迭代平方收敛速度的原因,并提出保持牛顿迭代平方收敛速度、保证一阶精度和无条件稳定性的一致性算法.一致性算法具备以下两个特征:1)采用路径无关计算格式;2)采用一致弹塑性切线模量。根据一致性算法构造出以弯矩和曲率为基本变量的弹塑性板弯曲有限元NIDKQ元。数值结果表明NIDKQ元具有令人满意的精度,同时验证了有限载荷增量下牛顿迭代一致性算法的平方收敛率特性,而传统切线刚度法随着塑性区的扩展将大大降低收敛速度。     

GURTIN VARIATIONAL PRINCIPLE AND FINITE ELEMENT SIMULATION FOR DYNAMICAL PROBLEMS OF FLUID-SATURATED POROUS MEDIA

Based on the theory of porous media,a general Gurtin variational principle for theinitial boundary value problem of dynamical response of fluid-saturated elastic porous media isdeveloped by assuming infinitesimal deformation and incompressible constituents of the solid andfluid phase.The finite element formulation based on this variational principle is also derived.Asthe functional of the variational principle is a spatial integral of the convolution formulation,thegeneral finite element discretization in space results in symmetrical differential-integral equationsin the time domain.In some situations,the differential-integral equations can be reduced to sym-metrical differential equations and,as a numerical example,it is employed to analyze the reflectionof one-dimensional longitudinal wave in a fluid-saturated porous solid.The numerical results canprovide further understanding of the wave propagation in porous media.     

AN EULERIAN FINITE ELEMENT METHOD FOR TIME-DEPENDENT FREE SURFACE PROBLEMS IN HYDRODYNAMICS

A numerical method based on the finite element method is presented for simulating the two-dimensional transient motion of a viscous liquid with free surfaces. For ease of numerical treatment of the free surface expressed by a multiple-valued function, the marker particle method is employed. Numerous virtual particles are spread over all regions occupied by liquid. They move about on a fixed finite element mesh with the liquid velocity at their positions. These particles contribute nothing to the dynamics of the liquid and only serve as markers of liquid regions. The velocity field within liquid regions is calculated by solving the Navier– Stokes equations and the equation of continuity by the finite element method based on quadrilateral elements. A detailed discussion is given of the methodological problems arising in the implementation of the marker particle method on an unstructured finite element mesh and of the solutions to these problems. The proposed method is demonstrated on three sample problems: the broken dam problem, the impact of a falling liquid drop on a still liquid and the entry of a rigid block into water. Good agreement has been obtained in the comparison of the present numerical results with available experimental data.     

PARAMETRIC VARIATIONAL PRINCIPLE BASED ELASTIC-PLASTIC ANALYSIS OF HETEROGENEOUS MATERIALS WITH VORONOI FINITE ELEMENT METHOD

The Voronoi cell finite element method (VCFEM) is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.