斜拉桥分析的索膜单元理论及其应用

本文提出分析斜拉桥的索膜单元理论,建立起索膜单元概念与计算模型,推导出索膜单元刚度矩阵,说明了用索膜单元分析斜拉桥的方法细节,应用索膜单元对144+360+144米公路钢斜拉桥方案进行了分析计算,得出了主梁弯矩影响面与扭矩图,讨论了索梁锚固点与截面水平主惯性轴偏心对主梁弯矩的影响。     

基于能量平衡原理的有限单元耦合方法

结构非线性数值计算分析应真实反映局部损伤破坏细节,以作为损伤演化全过程分析的依据。对同类构件,有限单元耦合方法可以解决破坏细节与整体模拟的空间尺度差异问题。基于能量平衡原理,建立了梁与实体单元、梁与壳单元以及壳与实体单元的耦合方程,适用于结构的损伤数值计算。对某RC框架结构原位推覆试验的损伤数值分析表明,有限单元耦合模型能正确反映整体结构的承载力和变形性能,并且能准确反映局部损伤破坏细节。     

单元矩阵的自动生成

本文对计算单元矩阵的过程进行了仔细地分析,根据运算中的特点,提出了一种自动生成单元矩阵的方法。根据该方法编制的计算程序,经过例题检验表明准确可靠。由不同的插值函数和单元几何形状,即可导出各种单元矩阵,将其与总体分析程序接口,可完成不同问题的有限元法数值分析。     

曲厚壳单元在斜板拱桥分析中的应用

介绍了利用曲厚壳单元对斜板拱桥进行空间受力分析,并将曲厚壳单元所得结果与采用实体单元分析的结果和模型实验所得结果进行了比较,证明利用曲厚壳单元来分析斜板拱桥是行之有效的。     

基于Abaqus的剪切型阻尼器单元开发与应用

基于Abaqus软件平台,利用FORTRAN语言二次开发了一种剪切型金属阻尼器用户子单元(UEL)。通过与Abaqus软件库中用于近似模拟阻尼器的单元进行比较,验证了该用户子单元的正确性,为Abaqus软件二次开发提供了参考。与软件库自带的单元相比,该用户子单元参数输入简单,无需试算,动力分析下的滞回性能更明显。同时,本文对软件中的不同单元进行了探讨,为剪切型金属阻尼器的模拟提供了新的思路。     

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

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

16自由度的中厚板弯曲矩形单元

基于有限条带思想,引入结点扭率自由度,利用深梁单元的位移模式建立了一个4结点16自由度中厚板弯曲高阶单元,此单元是薄板单元BFS-16的推广形式,其特点是单元的横向位移、转角位移、剪应变位移模式直接构造,在边界上位移模式与深梁单元一致,方便与梁单元叠加,适应于带加劲肋的板弯曲问题分析,用于薄壁结构时可考虑翘曲。实例计算显示,此单元精度高,计算稳定,收敛快,无剪切闭锁现象,能较好地反映中厚板的边界效应。     

基于应力分析的桁架结构双向拓扑优化方法

结构拓扑优化问题的研究多是采用基结构的思路,通过删除在设计区域内的不必要单元来得到结构的最优拓扑构型。本文探索了一种增加单元与删除单元相结合的双向拓扑优化方法,采用了网格与杆件两类单元对桁架进行分析:在高应力杆件单元周围生成新网格单元,并且删除低应力的杆单元,结构逐渐进化,从而得到优化的拓扑构型。文章最后通过数值算例,表明该方法是可行的。     

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

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

含界面裂纹Reissner板弯曲问题分析的奇异单元

首先,采用特征函数渐近展开法,推导了Reissner板弯曲界面裂纹尖端附近位移场渐近展开的前两阶显式表达式,并利用所获得的位移场渐近表达式构造了一种可用于Reissner板弯曲界面裂纹分析的奇异单元。然后,将该奇异单元与外部的常规有限单元相结合,开展了含界面裂纹Reissner板弯曲断裂问题的数值分析。奇异单元可以较好地描述裂纹尖端附近的内力场与位移场,其优势是它与常规单元进行连接时不需要使用过渡单元,并且可以直接给出应力强度因子等断裂参数的高精度数值结果。最后,通过两个数值算例验证了本文方法的有效性。     

From Single-Phase to Compositional Flow: Applicability of Mixed Finite Elements

In this paper we discuss the formulation of the governing equations that describe flow of fluids in porous media. Various types of fluid flow, ranging from single-phase flow to compositional flow, are considered. It is shown that all the differential equations governing these types of flow can be effectively rewritten in a fractional flow formulation; i.e., in terms of a global pressure and saturation (or saturations), and that mixed finite element methods can be accurately exploited to solve the pressure equation. Numerical results are presented to see the performance of the mixed methods for the flow equations in three space dimensions.     

A new decoupled finite element algorithm for viscoelastic flow. Part 1: Numerical algorithm and sample results

This paper presents an algorithm for two-dimensional Steady viscoelastic flow Simulation in which the Solution of the momentum and continuity equations is decoupled from that of the constitutive equations. The governing equations are discretized by the finite element method, with 3 × 3 element subdivision for the stress field approximation. Non-consistent Streamline upwinding is also used. Results are given for flow through a converging channel and through an abrupt planar 4:1 contraction.     

计算离心泵叶轮流场的扫描流束有限元方法

本文提出一种求解离心式叶轮流场的数值方法，将流动求解区域离散为有限个由流线构成其边界的单元，采用伽辽金法建立的单元方程在一条流束上集合为方程组，流线上的节点坐标亦作为未知量包含在有限元方程中，通过扫描计算，逐步解得流线位置及流动参数。本文应用叶轮的通流理论流动模型，采用扫描流速有限元方法对离心泵叶轮流场进行了计算，并与有关文献作了比较。     

Simulation of Coupled Viscous and Porous Flow Problems

Theoretical and numerical formulations are presented for the conjugate problem involving incompressible flow and flow in a saturated porous medium. The major focus of the work is the development of a generally applicable finite element method for the simulation of both fixed interface and evolving porous interface problems. The available alternatives for coupling Darcy and non-Darcy models to the Navier-Stokes equations have been studied and evaluated in a mixed finite element framework. Questions regarding convergence of the finite element method for porous flow models have been addressed. Numerical experiments on simple flow geometries have revealed the shortcomings of both the Darcy and Brinkman models. Application of the more realistic models to practical, multidimensional, flow studies has also been demonstrated.     

On the use of characteristic‐based split meshfree method for solving flow problems

This study presents characteristic‐based split (CBS) algorithm in the meshfree context. This algorithm is the extension of general CBS method which was initially introduced in finite element framework. In this work, the general equations of flow have been represented in the meshfree context. A new finite element and MFree code is developed for solving flow problems. This computational code is capable of solving both time‐dependent and steady‐state flow problems. Numerical simulation of some known benchmark flow problems has been studied. Computational results of MFree method have been compared to those of finite element method. The results obtained have been verified by known numerical, analytical and experimental data in the literature. A number of shape functions are used for field variable interpolation. The performance of each interpolation method is discussed. It is concluded that the MFree method is more accurate than FEM if the same numbers of nodes are used for each solver. Meshfree CBS algorithm is completely stable even at high Reynolds numbers. Copyright © 2007 John Wiley & Sons, Ltd.     

A Quasilinear Based Procedure for Saturated/Unsaturated Water Movement in Soils

The quasilinear form of Richards equation for one-dimensional unsaturated flow in soils can be readily solved for a wide variety of conditions. However, it cannot explain saturated/unsaturated flow and the constant diffusivity assumption, used to linearise the transient quasilinear equation, can introduce significant error. This paper presents a quasi-analytical solution to transient saturated/unsaturated flow based on the quasilinear equation, with saturated flow explained by a transformed Darcy's equation. The procedure presented is based on the modified finite analytic method. With this approach, the problem domain is divided into elements, with the element equations being solutions to a constant coefficient form of the governing partial differential equation. While the element equations are based on a constant diffusivity assumption, transient diffusivity behaviour is incorporated by time stepping. Profile heterogeneity can be incorporated into the procedure by allowing flow properties to vary from element to element. Two procedures are presented for the temporal solution; a Laplace transform procedure and a finite difference scheme. An advantage of the Laplace transform procedure is the ability to incorporate transient boundary condition behaviour directly into the analytical solutions. The scheme is shown to work well for two different flow problems, for three soil types. The technique presented can yield results of high accuracy if the spatial discretisation is sufficient, or alternatively can produce approximate solutions with low computational overheads by using large sized elements. Error was shown to be stable, linearly related to element size.     

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)     

A unified mesh refinement method with applications to porous media flow

A unified algorithm is presented for the refinement of finite element meshes consisting of tensor product Lagrange elements in any number of space dimensions. The method leads to repeatedly refined n-irregular grids with associated constraint equations. Through an object-oriented implementation existing solvers can be extended to handle mesh refinements without modifying the implementation of the finite element equations. Various versions of the refinement procedure are investigated in a porous media flow problem involving singularities around wells. A domain decomposition-type finite element method is also proposed based on the refinement technique. This method is applied to flow in heterogeneous porous media. © 1998 John Wiley & Sons, Ltd.     

On the implementation of symmetric and antisymmetric periodic boundary conditions for incompressible flow

In this paper we consider symmetric and antisymmetric periodic boundary conditions for flows governed by the incompressible Navier-Stokes equations. Classical periodic boundary conditions are studied as well as symmetric and antisymmetric periodic boundary conditions in which there is a pressure difference between inlet and outlet. The implementation of this type of boundary conditions in a finite element code using the penalty function formulation is treated and also the implementation in a finite volume code based on pressure correction. The methods are demonstrated by computation of a flow through a staggered tube bundle.     

Numerical analysis of density flows in adiabatic two-phase fluids through characteristic finite element method

The purpose of this study is to analyze the density flow in adiabatic two-phase fluids through the characteristic finite element method. The fluids are assumed to be liquids. The equations of conservations of mass and momentum for the adiabatic flows and the Birch–Murnaghan equation of state are employed as the governing equations. The employed finite element method is a combination of the characteristic method and the implicit method. The governing equations are divided into two parts: the advection part and the non-advection part. The characteristic method is applied to the advection part. The Hermite interpolation function, which is based on the complete third-order polynomial interpolation using triangular finite element is employed for the interpolation of both velocity and density. Using the discontinuity conditions, an interface translocation method can be derived. The interface of the two flow densities are interpolated through the third-order spline function, using which the curvature of the interface can be directly computed. For the numerical study, the development of density flow over the Tokyo bay is presented. It is detected out that high density area is abruptly diffused over the whole area. According to the differences in the two densities, various flow patterns are computed.