SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection (EEP) method for computation of super-convergent results in one-dimensional finite element method (FEM), the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.     

微重力环境下Marangoni对流的有限元数值模拟

本文以悬浮区为背景研究液桥中气液交界面上由表面张力所驱动的流体对流。我们采用有限元方法对轴对称、定常运动方程,能量方程,自由面切向、法向应力平衡条件迭代求解,首次考虑了边界形状的影响,得到了自洽的流区位形和流场、温度场、表面压强分布。一般说来,流区自由面呈弯月形。结果表明,只要表面张力数 C_a<1,静态平衡的流区位形就是动态情况的良好近似。本文还分析了C_a数及G_r 数对流区位形的影响,得到了不同 M_a 数及散热条件下的温度场和流函数分布。     

用有限元自适网格控制结构应力分析的精度

本文涉及有限元结构应务分析精度控制的三个基本方面：１）自适应的有限元网格自动生成技术；２）有限元应力分析的误差估计：３）用单元尺寸场将误差分布用于形成新网格。     

A METHOD OF DETERMINING BUCKLED STATES OF THIN PLATES AT A DOUBLE EIGENVALUE

A method of determining bifurcation directions at a double eigenvalue is presented by combining the finite element method with the perturbation method. By using the present method, the buckled states of rectangular plates at a double eigenvalue are numerically analyzed. The results show that this method is effective.     

压电材料椭圆切口的力学分析

在线性压电本构方程框架下，用复势函数方法对椭圆切口模型进行了精确的数值计算。完整地考虑了各向异性力电耦合效应以及切口内不同电介质的介电性质。给出了切口内部不同介电性质对压电材料内部应力的影响。指出了Sosa文章里的一些计算错误。由于现在文献中很少有关于电导通边界条件下理论解的数值结果，所以本文同时提供了不同电边界条件下理论解的数值结果，所以本文同时提供了不同电边界条件下的理论解的数值结果。最后通过最小势能原理建立了8结点有限元模型，对椭圆切口问题进行了计算并与理论解进行了仔细比较。     

点式压电智能结构的模型修正与振动主动控制

提出了在点式压电智能结构中应用摄动有限元方法对结构的有限元模型进行修正 ,从而达到提高建模精度 ,改善实际结构振动主动控制效果的目的。通过对一悬臂梁在模型修正前后进行振动主动控制的不同的控制效果验证了该方法的有效性     

基于拉格朗日乘子法和群论的具有任意位移边界条件的旋转周期对称结构有限元分析

对任意位移边界条件下的旋转周期对称结构，由拉格朗日乘子法建立有限元方程。在对称适应的坐标系下，由结构刚度矩阵的块循环性质，利用群变换给出一种新的求解方法。数值验证给出令人满意的结果。     

非连续边界元——有限元耦合方法分析

对边界元－有限元耦合方法进行了分析，采用非连续元离散边界积分方程，解决了耦合分析中自由度约束问题，给出了非连续边界元同有限元耦合的具体实施步骤，通过对二维弹性力学和流＝固耦合问题分析，表明了该文方法的有效性。     

非线性动力学常微分方程组高精度数值积分方法

建立了一种求解非线性动力学常微分方程组初值问题的新方法．若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解．引入一般自治或非自治非线性系统的首次近似Jacobi矩阵，不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel卷积分解析表达式，在时间步长内进行数值积分选代求解，在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解．积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法．数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率．算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题．基于以上理论和数值方法，计算了线性非线性算例并进行了分析，验证了本方法简捷而有广泛的适应性，可以有足够的精确性．     

Structural dynamics and stability analysis of 2D-FG microbeams with two-directional porosity distribution and variable material length scale parameter

Abstract

Since the two-directional functionally graded (2D-FG) materials can satisfy the new requirements raised based on the elimination of the stress concentration, delamination and cracking problems accompanying with the low cost and lightweight on the structures without sacrificing the stiffness and strength, the structural analyses of these structures become more important than ever. Moreover, the usage of the micro-electromechanical systems composed of 2D-FG materials has been increasing in automotive, military, space, biomedical, and nuclear energy industries. Within this study, the free vibration and buckling behaviors of 2D-FG porous microbeams are investigated based on the modified couple stress theory by employing a transverse shear-normal deformation beam theory and using finite element method. The effects of the thickness to material length scale parameter (MLSP) accompanying with the micro-porosity volume fraction ratio, boundary condition, aspect ratio, and gradient index on the dimensionless fundamental frequencies and dimensionless critical buckling loads of the 2D-FG porous microbeams are investigated. Moreover, with assumption of the variable material length scale parameters (VMLSP), the computed results are compared with ones obtained by employing constant MLSP. It is found that VMLSP increases the stiffness of the 2D-FG porous microbeams and effects the free vibration and buckling responses of these structures.