节点梯度光滑有限元配点法

配点法构造简单、计算高效,但需要用到数值离散形函数的高阶梯度,而传统有限元形函数的梯度在单元边界处通常仅具有C0连续性,因此无法直接用于配点法分析.本文通过引入有限元形函数的光滑梯度,提出了节点梯度光滑有限元配点法.首先基于广义梯度光滑方法,定义了有限元形函数在节点处的一阶光滑梯度值,然后以有限元形函数为核函数构造了有...     

Application of Kalman Filter Finite Element Method and AIC

This paper presents a numerical simulation for application of the Kalman filter finite element method. The Kalman filter is employed frequently for the solution of time series analysis including observation and system noises. Applying the Kalman filter to the finite element method, the present method is capable of the estimation in time and space directions. In this method, the matrix generated by the finite element method is applied to the state transition matrix. Using the Kalman filter finite element method, the characteristics of both the Kalman filter and the finite element method can be strengthened. In this paper, the state transition matrix is based on the shallow water equations which are approximated by the finite element method. This method can estimate the tidal current not only in time but also in space directions.     

Finite element analysis of non-Newtonian fluid flow in 2-d branching channel

This paper presents finite element analysis of non-Newtonian fluid flow in 2-d branching channel. The Galerkin method and mixed finite element method are used. Here the fluid is considered as incompressible, non-Newtonian fluid with Oldyord differential-type constitutive equation. The non-linear algebraic equation system which is formulated with finite element method is solved by means of continuous differential method. The results show that finite element method is suitable for the analysis of non-Newtonian fluid flow with complex geometry.     

基于扩展有限元的应力强度因子的位移外推法

针对平面裂纹问题,阐述了扩展有限元法的单元位移模式、推导了扩展有限元法的控制方程、介绍了特殊单元的数值积分技术.基于最小二乘法,建立了应力强度因子位移外推法的计算公式.利用MATLAB编写计算程序,对平面裂纹问题用扩展有限元法进行了计算.基于扩展有限元法的计算结果,分别利用位移外推法和相互作用积分法,对平面裂纹的应力强度因子进行了计算.计算结果表明,位移外推法比相互作用积分法能更方便和准确地计算平面裂纹的应力强度因子.     

RTM充模过程数值模拟的隐式有限元算法

建立了基于欧拉方法描述树脂传递模塑(RTM)工艺充模过程的基本数学方程,并采用有限元隐式时间积分方法对基本方程进行了数值求解．编制了基于隐式有限元算法及传统有限元控制体算法的程序,通过具体算例比较了这两种算法的优缺点．与传统的有限元控制体法相比,该文提出的隐式有限元算法能节省计算时间,特别适合于单元、节点数目多的情况．隐式有限元算法是一种纯有限元方法,不需要使用控制体积技术,采用该算法计算出的流动前沿与时间步长无关。     

应用分形有限元方法于外域声场计算

 应用二级分形有限元方法计算了外域声场. 用一人工边界把外域声场分为两 部分，人工边界以内使用常规有限元方法，人工边界以外的无限大区域使用分形有限元方法. 使用分形有限元方法的优点是：一方面形成几何自相似网格使得相邻层之间的单元刚度矩阵 和质量矩阵具有非常简单的关系；另一方面引用自动满足无限远辐射条件的全域插值函数把 节点自由度变换为一组广义坐标，因而计算量可以大大减少. 数值算例表明：该方法对于计 算无限大外域声场是有效的.     

Adaptive finite element method for high-speed flow-structure interaction

An adaptive finite element method for high-speed flow-structure interaction is presented. The cell-centered finite element method is combined with an adaptive meshing technique to solve the Navier-Stokes equations for high-speed compressible flow behavior. The energy equation and the quasi-static structural equations for aerodynamically heated structures are solved by applying the Galerkin finite element method. The finite element formulation and computational procedure are described. Interactions between the high-speed flow, structural heat transfer, and deformation are studied by two applications of Mach 10 flow over an inclined plate, and Mach 4 flow in a channel. The project supported by the Thailand Research Fund (TRF)     

局部有限元模型边界条件的直接传递方法

提出一种将整体分析得到的节点力或节点位移直接传递到精细化局部有限元模型的方法,即部分混合单元法。沿精细化局部有限元模型周边建立一组过渡单元,该组过渡单元采用与整体模型一致的单元类型和模拟方式,其外侧边界上的节点与整体模型节点的相对坐标对应,内侧边界与精细化局部有限元模型采用基于面约束的方式连接。在外侧边界上根据节点坐标对应施加整体分析获得的节点力或节点位移,过渡单元就可直接将边界条件传递到精细化局部有限元模型。通过贵州红水河特大桥钢-混结合段的精细化有限元分析,验证了本文方法的实用性和有效性。     

Coupling of assumed stress finite element and boundary element methods with stress-traction equilibrium

In this study, the stress based finite element method is coupled with the boundary element method in two different ways. In the first one, the ordinary distribution matrix is used for coupling. In the second one, the stress traction equilibrium is used at the interface line of both regions as a new coupling process. This new coupling procedure is presented without a distribution matrix. Several case studies are solved for the validation of the developed coupling procedure. The results of case studies are compared with the distribution matrix coupling, displacement based finite element method, assumed stress finite element method, boundary element method, ANSYS and analytical results whenever possible. It is shown that the coupling of the stress traction equilibrium with assumed stress finite elements gives as accurate results as those by the distribution matrix coupling.     

Mixed time discontinuous space-time finite element method for convection diffusion equations

A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.     

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

薄板小波有限元理论及其应用

利用样条小波尺度函数构造了常用的三角形和矩形薄板单元的位移函数,得到了利用小波函数表示的形函数。采用合理的局部坐标,对单元进行压缩,使单元在局部坐标区间上有其值,成功地推导出了分域的三角形和矩形薄板小波有限元列式。在此基础上,提出了弹性地基薄板的小波有限元求解方法。通过两个算例对薄板的挠度和弯矩进行了计算,数值结果表明,求解结果具有收敛快、精度高的特点。     

激光辐照下圆柱壳结构的稳定性分析

针对航天工程中常用的承受轴压作用的薄壁圆柱壳,分别采用解析方法与特征值屈曲有限元方法分析了圆柱壳结构在均匀轴压作用下的稳定性能,得到了屈曲承载力.并进行了对比,验证了有限元模型的合理性.采用线性屈曲特征值有限元方法分析了强激光辐照作用下圆柱壳的稳定性能,分析表明激光辐照导致壳体局部温度上升并由此带来材料参数的改变,是因为激光辐照在壳体中引起了热应力与热应变,使轴压圆柱壳的屈曲承载力明显降低.论文还着重对加筋圆柱壳结构的稳定性进行了研究,数值分析结果表明,加筋能有效提高圆柱壳结构的抗压承载力,激光辐照作用下,加筋对轴压作用下圆柱壳的屈曲承载力的提高作用更为明显.     

地下水地表水联合模拟的解析单元法理论及应用

本文介绍了解析单元法的基本理论及适用条件 ,并运用该方法 (GFL OW1软件 )计算了某水库蓄水过程中库首结晶岩区风化裂隙水渗流场变化 ,取得良好效果.     

Galerkin-Petrov least squares mixed element method for stationary incompressible magnetohydrodynamics

The Galerkin-Petrov least squares method is combined with the mixed finite element method to deal with the stationary, incompressible magnetohydrodynamics system of equations with viscosity. A Galerkin-Petrov least squares mixed finite element format for the stationary incompressible magnetohydrodynamics equations is presented. And the existence and error estimates of its solution are derived. Through this method, the combination among the mixed finite element spaces does not demand the discrete Babuska-Brezzi stability conditions so that the mixed finite element spaces could be chosen arbitrartily and the error estimates with optimal order could be obtained.     

Comparison of methods for vibration analysis of electrostatic precipitators

The paper presents two methods for the formulation of free vibration analysis of collecting electrodes of precipitators.The first,called the hybrid finite element method, combines the finit element method used for calculations of spring deformations with the rigid finite element method used to reflect mass and geometrical features,which is called the hybrid finite element method.As a result,a model with a diagonal mass matrix is obtained.Due to a specific geometry of the electrodes,which are long plates of complicated shapes,the second method proposed is the strip method which is a semi-analytical method.The strip method allows us to formulate the equations of motion with a considerably smaller number of generalized coordinates.Results of numerical calculations obtained by both methods are compared with those obtained using commercial software like ANSYS and ABAQUS.Good compatibility of results is achieved.     

加肋开孔柱壳混合有限条——有限元屈曲分析

本文提出了一种分析环加肋开孔柱壳屈曲问题的混合有限条──有限元法。环加肋柱壳作为一个构造上的正交各向异性壳处理，柱壳非开孔区采用有限条元离散，开孔区采用有限单元离散。在有限条元与有限单元交界面上，根据位移协调条件建立条元和单元的耦合关系，并据此构造一种特殊的过渡单元、联接条元和单元，进行整体分析。算例表明，这一方法对开孔柱壳屈曲问题的分折十分有效。     

结构动力分析自适应有限元方法综述

结构动力分析自适应有限元方法主要研究有限元动力分析的误差估计理论,建立适用于复杂结构动力分析的有限元网格自适应过程．介绍了结构动力问题自适应有限元方法的重要发展,包括固有振动和动响应分析的误差估计及相应的自适应策略;且简要介绍了几种现有的网格生成技术及其特点．最后指出这种方法存在的问题和今后的研究方向．     

基于互补理论和扩展有限元法的摩擦接触裂纹应力强度因子计算

从虚功方程出发,结合扩展有限元离散技术与接触条件的非线性互补表述,建立了摩擦接触裂纹问题的扩展有限元非线性互补模型,将不等式接触条件转化为非线性互补类的非光滑方程组,并采用基于广义导数的非光滑阻尼牛顿法求解方程组,无需引入任何额外人工变量以及迭代求解。以含中心倾斜裂纹平板和边裂纹平板为例,运用相互作用积分法计算摩擦接触裂纹的应力强度因子,将其结果与理论解进行对比分析,该方法都能给出精确的计算结果;基于扩展有限元方法对单轴压缩作用下倾斜裂纹扩展过程进行了数值模拟,计算结果表明,受压裂纹数值结果与实验结果比较吻合,从而验证了本文方法的有效性与正确性。     

拟协调元研究综述

拟协调元是有限元中十分重要的、具有特色的一种列式体系. 拟协调元列式简单、灵活, 统一了协 调元、非协调元等列式方法. 在列式中, 拟协调元将几何方程和平衡方程同时弱化, 并强调基函数在有限元空 间中的重要作用; 借助对位移和应变离散精度的控制, 拟协调元保障了单元的收敛性, 并可以利用泰勒展开校 核进行简便直接的收敛性分析. 研究者们利用拟协调元已经构造了大量的优秀的单元, 并广泛地应用到结构问 题、流体流动问题、非线性分析、稳定性和破坏分析等方面. 这些工作集中体现了拟协调元的理论价值和工程 应用价值. 对拟协调列式方法、列式理论和已发表文献中的主要拟协调单元进行了总结. 最后对拟协调的研究 工作进行了展望.