NUMERICAL SIMULATION OF STRESS WAVE PROPAGATION WITH UNSTRUCTURED FINITE VOLUME TIME DOMAIN METHOD

An unstructured finite volume time domain method (UFVTDM) is proposed to simulate stress wave propagation. The original variables of displacement and stress are solved based on the dynamic equilibrium equations. An Euler explicit and unstructured finite volume method is used for time and spacial terms respectively. The displacements are stored on the cell vertex and a vertex based finite volume is formed with the integral surface and the stress is assumed as uniform in the cell. This is some similar with the stager grid method in computational fluid dynamics. Several cases are used to show the capability of the algorithm.     

H 1 space-time discontinuous finite element method for convection-diffusion equations

An H~1 space-time discontinuous Galerkin (STDG) scheme for convectiondiffusion equations in one spatial dimension is constructed and analyzed. This method is formulated by combining the H~1 Galerkin method and the space-time discontinuous finite element method that is discontinuous in time and continuous in space. The existence and the uniqueness of the approximate solution are proved. The convergence of the scheme is analyzed by using the techniques in the finite difference and finite element methods. An optimal a-priori error estimate in the L~∞ (H~1 ) norm is derived. The numerical exper- iments are presented to verify the theoretical results.     

Three-dimensional steady-state closed form solution for multilayered fluid-saturated anisotropic finite media due to surface/internal point source

A three-dimensional(3 D)steady-state solution of fluid saturated anisotropic finite media is presented.The eigenequation method and the pseudo-Stroh formalism are used to obtain the exact solution for homogeneous saturated finite media.The propagator matrix method is introduced to deal with the corresponding multilayered poroelastic media.The poroelastic solutions due to surface or internal point fluid source are obtained.The comparison of the results of the saturated isotropic media in a half space and those obtained by the finite element method is given to illustrate the accuracy of the solution in a finite domain.Numerical solutions of a sandwich poroelastic medium are presented to analyze its hydromechanical behaviors.Two ratios of the horizontal permeability to vertical permeability and different source positions are investigated.The results show that the fluid parameters and source positions have great influence on the hydromechanical behaviors of the layered media.     

ANALYSIS AND EXPERIMENT ON TRANSIENT DYNAMIC RESPONSE IN FINITE MINDLIN PLATE

The transient wave propagation in the finite rectangular Mindlin plate is investigated by the analytical and experimental methods. The generalized ray method (GRM) which has been successfully applied to study the transient responses of beams, planar trusses, space frames and infinite layered media is extended to investigate the transient wave propagation and early short time transient response in finite Mindlin plate. Combining the wave solution, the shock source and the boundary conditions, the ray groups transmitted in the finite rectangular plate can be determined. Numerical simulations and experiments are performed and compared with each other. The results show that the transient wave propagation and early short time transient responses in the finite plate can be studied using the GRM. The early short time transient accelerations are very large for the finite plate subjected to the unit impulse, while the early short time transient displacements are very small. The early short time transient accelerations under the unit impulse are much larger than those under the unit step impulse. The thickness and material characteristics have remarkable effects on the early short time transient responses.     

Local and parallel finite element algorithms based on two-grid discretization for steady Navier-Stokes equations

Local and parallel finite element algorithms based on two-grid discretization for Navier-Stokes equations in two dimension are presented. Its basis is a coarse finite element space on the global domain and a fine finite element space on the subdomain. The local algorithm consists of finding a solution for a given nonlinear problem in the coarse finite element space and a solution for a linear problem in the fine finite element space, then droping the coarse solution of the region near the boundary. By overlapping domain decomposition, the parallel algorithms are obtained. This paper analyzes the error of these algorithms and gets some error estimates which are better than those of the standard finite element method. The numerical experiments are given too. By analyzing and comparing these results, it is shown that these algorithms are correct and high efficient.     

Interval finite element method and its application on anti-slide stability analysis

The problem of interval correlation results in interval extension is discussed by the relationship of interval-valued functions and real-valued functions. The methods of reducing interval extension are given. Based on the ideas of the paper, the formulas of sub-interval perturbed finite element method based on the elements are given. The sub-interval amount is discussed and the approximate computation formula is given. At the same time, the computational precision is discussed and some measures of improving computational efficiency are given. Finally, based on sub-interval perturbed finite element method and anti-slide stability analysis method, the formula for computing the bounds of stability factor is given. It provides a basis for estimating and evaluating reasonably anti-slide stability of structures.     

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.     

Finite volume element method for analysis of unsteady reaction–diffusion problems

A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element method together. The finite volume method is used to discretize the unsteady reaction-diffusion equation, while the finite element method is applied to estimate the gradient quantities at cell faces. Robustness and efficiency of the combined method have been evaluated on uniform rectangular grids by using available numerical solutions of the two-dimensional reaction-diffusion problems. The numerical solutions demonstrate that the combined method is stable and can provide accurate solution without spurious oscillation along the high-gradient boundary layers.     

Analysis of bonded anisotropic wedges with interface crack under anti-plane shear loading

The antiplane stress analysis of two anisotropic finite wedges with arbitrary radii and apex angles that are bonded together along a common edge is investigated. The wedge radial boundaries can be subjected to displacement-displacement boundary condi- tions, and the circular boundary of the wedge is free from any traction. The new finite complex transforms are employed to solve the problem. These finite complex transforms have complex analogies to both kinds of standard finite Mellin transforms. The traction free condition on the crack faces is expressed as a singular integral equation by using the exact analytical method. The explicit terms for the strength of singularity are extracted, showing the dependence of the order of the stress singularity on the wedge angle, material constants, and boundary conditions. A numerical method is used for solving the resul- tant singular integral equations. The displacement boundary condition may be a general term of the Taylor series expansion for the displacement prescribed on the radial edge of the wedge. Thus, the analysis of every kind of displacement boundary conditions can be obtained by the achieved results from the foregoing general displacement boundary condition. The obtained stress intensity factors （SIFs） at the crack tips are plotted and compared with those obtained by the finite element analysis （FEA）.     

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.     

On Taylor weak statement finite element methods for computational fluid dynamics

A Taylor series augmentation of a weak statement (a ‘Taylor weak statement’ or ‘Taylor-Galerkin’ method) is used to systematically reduce the dispersion error in a finite element approximation of the one-dimensional transient advection equation. A frequency analysis is applied to determine the phase velocity of semi-implicit linear, quadratic and cubic basis one-dimensional finite element methods and of several comparative finite difference/finite volume algorithms. The finite element methods analysed include both Galerkin and Taylor weak statements. The frequency analysis is used to obtain an improved linear basis Taylor weak statement finite element algorithm. Solutions are reported for verification problems in one and two dimensions and are compared with finite volume solutions. The improved finite element algorithms have sufficient phase accuracy to achieve highly accurate linear transient solutions with little or no artificial diffusion.     

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

配点法构造简单、计算高效, 但需要用到数值离散形函数的高阶梯度,而传统有限元形函数的梯度在单元边界处通常仅具有C$^{0}$连续性,因此无法直接用于配点法分析. 本文通过引入有限元形函数的光滑梯度,提出了节点梯度光滑有限元配点法. 首先基于广义梯度光滑方法,定义了有限元形函数在节点处的一阶光滑梯度值,然后以有限元形函数为核函数构造了有限元形函数的一阶光滑梯度,进而对一阶光滑梯度直接求导并用一阶光滑梯度替换有限元形函数的标准梯度,即完成了有限元形函数二阶光滑梯度的构造.文中以线性有限元形函数为基础的理论分析表明,其光滑梯度不仅满足传统线性有限元形函数梯度对应的一阶一致性条件,而且在均布网格假定下满足更高一阶的二阶一致性条件.因此与传统线性有限元法相比,基于线性形函数的节点梯度光滑有限元法的$L_{2}$和$H_{1}$误差均具有二次精度,即其$H_{1}$误差收敛阶次比传统有限元法高一阶, 呈现超收敛特性.文中通过典型算例验证了节点梯度光滑有限元配点法的精度和收敛性,特别是其$H_{1}$或能量误差的精度和收敛率都明显高于传统有限元法.      

关于无振荡、无自由参数有限元格式的研究

利用双曲守恒律方程的Ｔａｙｌｏｒ弱解表达式，建立了有限元法修正方程，选择合适的展开式系数能得到一系列数值格式．通过稳定性分析研究了格式的稳定性、色散误差与有限元修正方程导数项系数之间的关系，该关系与差分法的ＮＮＤ格式一致．在选定格式下，通过ＣＦＬ数可控制有限元离散解的振荡而使格式不含自由参数．最后，用数值算例验证了这一关系，并在二、三维欧拉方程作了推广应用．     

带源参数的二维热传导反问题的无网格方法

利用无网格有限点法求解带源参数的二维热传导反问题，推导了相应的离散方程. 与 其它基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不 需要划分网格，用配点法离散控制方程，可以直接施加边界条件，不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.     

Finite difference and finite element methods for mhd channel flows

A Galerkin finite element method and two finite difference techniques of the control volume variety have been used to study magnetohydrodynamic channel flows as a function of the Reynolds number, interaction parameter, electrode length and wall conductivity. The finite element and finite difference formulations use unequally spaced grids to accurately resolve the flow field near the channel wall and electrode edges where steep flow gradients are expected. It is shown that the axial velocity profiles are distorted into M-shapes by the applied electromagnetic field and that the distortion increases as the Reynolds number, interaction parameter and electrode length are increased. It is also shown that the finite element method predicts larger electromagnetic pinch effects at the electrode entrance and exit and larger pressure rises along the electrodes than the primitive-variable and streamfunction–vorticity finite difference formulations. However, the primitive-variable formulation predicts steeper axial velocity gradients at the channel walls and lower axial velocities at the channel centreline than the streamfunction–vorticity finite difference and the finite element methods. The differences between the results of the finite difference and finite element methods are attributed to the different grids used in the calculations and to the methods used to evaluate the pressure field. In particular, the computation of the velocity field from the streamfunction–vorticity formulation introduces computational noise, which is somewhat smoothed out when the pressure field is calculated by integrating the Navier–Stokes equations. It is also shown that the wall electric potential increases as the wall conductivity increases and that, at sufficiently high interaction parameters, recirculation zones may be created at the channel centreline, whereas the flow near the wall may show jet-like characteristics.     

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

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

无穷域势流问题的有限元／差分线法混合求解

提出了一种将有限元和差分线法相结合求解无穷域势流问题的算法。用两同心圆将求解域划分为存在重叠的有限和无限两个区域,在有限和无限域上分别用有限元和差分线法求解Laplace方程边值问题。用差分线法推导出的关系式修正有限元方程,求解该方程组从而得到原问题的解。本算法将求解无穷域问题转化为代数特征值问题和有限域内线性方程组的...     

OBJECT-ORIENTED FINITE ELEMENT ANALYSIS AND PROGRAMMING IN VC++

ＩｎｔｒｏｄｕｃｔｉｏｎＳｉｎｃｅｆｉｎｉｔｅｅｌｅｍｅｎｔａｎａｌｙｓｉｓｐｒｏｇｒａｍｓｄｅａｌｗｉｔｈｔｈｅｔｈｅｏｒｉｅｓａｎｄｍｅｔｈｏｄｓｏｆｔｈｒｅｅｆｉｅｌｄｓｏｆｍｅｃｈａｎｉｃｓ,ａｐｐｌｉｅｄｍａｔｈｅｍａｔｉｃｓａｎｄｃｏｍｐｕｔｅｒｓｃｉｅｎｃｅ,ｔｈｅｙａｒｅｃｏｍｐｌｉｃａｔｅｄａｎｄｅｒｒｏｒ＿ｐｒｏｎｅ .Ｄｕｒｉｎｇｄｅｖｅｌｏｐｍｅｎｔｐｅｏｐｌｅａｌｌｔｈｒｏｕｇｈｓｅｅｋｔｈｅｍｅｔｈｏｄｔｏｉｍｐｒｏｖｅｐｒｏｇｒａｍｓ.Ｔｈｅｏｂｊｅｃｔ＿ｏｒｉｅｎｔｅ…     

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)