Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin (SUPG) method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid are performed by the Galerkin method. The second-order semiimplicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are linearized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effiectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.     

基于等几何分析的比例边界有限元方法

提出了一种具有比例边界有限元的半解析特性和等几何分析的几何特性的新方法。该新方法是在比例边界有限元框架中用NURBS曲线或曲面精确描述域边界几何形状,同时域边界位移场采用描述几何形状的NURBS形函数等参构造。这种新方法具有比例边界有限元固有的径向解析特性和NURBS的高阶连续性的优点。数值算例显示,与传统的比例边界有限元相比,基于等几何分析的比例边界有限元方法提高了域边界单元和域内应力场的连续性,减少了计算自由度。应用此方法可以用较少的计算自由度获得更高连续阶和更高精度的位移、应力和应变场。     

Control structure interaction in the nonlinear analysis of flexible mechanical systems

The effect of the control structure interaction on the feedforward control law as well as the dynamics of flexible mechanical systems is examined in this investigation. An inverse dynamics procedure is developed for the analysis of the dynamic motion of interconnected rigid and flexible bodies. This method is used to examine the effect of the elastic deformation on the driving forces in flexible mechanical systems. The driving forces are expressed in terms of the specified motion trajectories and the deformations of the elastic members. The system equations of motion are formulated using Lagrange's equation. A finite element discretization of the flexible bodies is used to define the deformation degrees of freedom. The algebraic constraint equations that describe the motion trajectories and joint constraints between adjacent bodies are adjoined to the system differential equations of motion using the vector of Lagrange multipliers. A unique displacement field is then identified by imposing an appropriate set of reference conditions. The effect of the nonlinear centrifugal and Coriolis forces that depend on the body displacements and velocities are taken into consideration. A direct numerical integration method coupled with a Newton-Raphson algorithm is used to solve the resulting nonlinear differential and algebraic equations of motion. The formulation obtained for the flexible mechanical system is compared with the rigid body dynamic formulation. The effect of the sampling time, number of vibration modes, the viscous damping, and the selection of the constrained modes are examined. The results presented in this numerical study demonstrate that the use of the driving forees obtained using the rigid body analysis can lead to a significant error when these forces are used as the feedforward control law for the flexible mechanical system. The analysis presented in this investigation differs significantly from previously published work in many ways. It includes the effect of the structural flexibility on the centrifugal and Coriolis forces, it accounts for all inertia nonlinearities resulting from the coupling between the rigid body and elastic displacements, it uses a precise definition of the equipollent systems of forces in flexible body dynamics, it demonstrates the use of general purpose multibody computer codes in the feedforward control of flexible mechanical systems, and it demonstrates numerically the effect of the selected set of constrained modes on the feedforward control law.     

Finite element analysis for consolidation in interaction between structure and saturated soil foundation

The consolidation analysis of interaction between structure and saturated soil foundation is discussed. With the use of substructure technique, the structure is condensed onto the interface of the soil, and then the consolidation governing equations to describe the interaction between soil and structure are derived, The solution with non-iterative algorithm is proposed in this paper. The pressure Master-Slave relation method is used to deal with the non-permeability conditions on soil boundaries. A numerical example is illustrated. Based on this paper, the interactive consolidation analysis between large structure and soil has been more practical.     

Numerical analysis of one-dimensional nonlinear large-strain consolidation by the finite element method

The nonlinear partial differential equation model of Gibson et al. which governs one-dimensional large-strain consolidation is solved numerically using a semi-discrete formulation involving a Galerkin weighted residual approach. The use of quadratic Lagrange basis functions usually complicates the task of solving the system of time-dependent ordinary differential equations that are obtained with the semi-discrete Galerkin procedure. However, an efficient algorithm has been discovered yielding the advantages of quadratic interpolation without undue computational burden.Although considerable effort has already been made to solve the PDE of large-strain consolidation by numerical methods, a satisfactory set of benchmarks is still needed to assess accuracy. To fill this need, three procedures are reported which allow numerical solutions of the large-strain model to be reliably evaluated. One involves the use of perturbation methodology to provide a solution when only self-weight effects are present. A second utilizes an analytical solution developed by Philip when self-weight effects are absent and the third involves the exact calculation of the discharge flux through the upper boundary of a deposit consolidating through self-weight effects alone. All three are restricted to early-time consolidation and are illustrated in the context of the finite element method.     

Numerical modeling of coupled seabed scour and pipe interaction

Scouring of the seabed by ice masses poses great threat to structural integrity and safety of buried structures such as oil and gas pipes. Large ridge motions and large seabed deformations complicate study of the seabed scouring and pipe interaction using classical Lagrangian methods. We present a new numerical approach for modeling the coupled seabed scour problem and its interaction with embedded marine pipes. In our work, we overcome the common issues associated with finite deformation inherently present within the seabed scour problem using a rheological approach for soil flow. The seabed is modeled as a viscous non-Newtonian fluid, and its interface described implicitly through a level-set function. The history-independent model allows for large displacements and deformations of the ridge and seabed, respectively. Numerical examples are presented to demonstrate the capabilities of the method with regard to seabed scouring and the ensuing large soil deformation.     

波导本征问题分析的比例边界有限元方法

引入了一种求解波导本征值问题的高效而精确算法-比例边界有限元方法SBFEM (Scaled Boundary Finite Element Method).该方法的一个特点是只需在边界上进行离散,问题降低一维,使计算工作量大大减少;另一特点是所建立的控制方程为二阶常微分方程,可以解析地求解,使计算精度得到了保证.论文利用变分原理并通过比例边界坐标变换,推导了TE波和TM波波导的比例边界有限元频域方程以及波导动剐度方程,同时给出了波导动刚度矩阵的连分式解形式,通过引入辅助变量进一步得出波导特征值方程并求出波导本征值.以矩形、L形波导和叶型加载矩形波导的本征问题分析为例,通过与解析解及其他数值方法比较,结果表明,此方法具有精度高、计算工作量小的优点,而且随着连分式阶数增加收敛速度快.进一步分析了一类角切四脊正方形波导的传输特性.     

On a new extended finite element method for dislocations: Core enrichment and nonlinear formulation

A recently developed finite element method for the modeling of dislocations is improved by adding enrichments in the neighborhood of the dislocation core. In this method, the dislocation is modeled by a line or surface of discontinuity in two or three dimensions. The method is applicable to nonlinear and anisotropic materials, large deformations, and complicated geometries. Two separate enrichments are considered: a discontinuous jump enrichment and a singular enrichment based on the closed-form, infinite-domain solutions for the dislocation core. Several examples are presented for dislocations constrained in layered materials in 2D and 3D to illustrate the applicability of the method to interface problems.     

几何精确NURBS有限元中边界条件施加方式对精度影响的三维计算分析

非均匀有理B样条(NURBS)有限元法把计算机辅助几何设计(CAGD)中的NURBS几何构形方法与有限元方法有机结合起来,有效消除了有限元离散模型的几何误差,提高了计算精度。但是由于NURBS基函数不是插值函数,直接在控制节点上施加位移边界条件会引起较大误差。本文详细讨论了NURBS基函数的插值特性,在NURBS有限元分析中采用罚函数法施加位移边界条件,提高了收敛率和计算精度。结合典型三维弹性力学问题,对两种施加位移边界条件的方法进行了对比和分析。计算结果表明,直接施加位移边界条件会导致收敛率和精度的明显降低,而基于罚函数法的NURBS有限元分析则能达到最优收敛率,并具有更高的精度。     

流体-结构耦合问题的有限元并行计算研究

流体-结构耦合问题广泛存在于各种工程领域,本文采用ALE显式有限元法求解该类问题,并对该方法的并行性进行讨论。同时根据流体-结构耦合问题与ALE显式有限元的计算特点,在坐标递归分区方法的基础上设计并程序实现了基于流体-结构耦合均衡的分区算法。通过与坐标递归分区方法的计算结果相比较,对于流体-结构耦合问题的求解,耦合均衡并行分区方法具有更好的加速比和并行效率。     

裂纹扩展过程中线性内聚力模型计算的半解析有限元法

提出了求解基于线性内聚力模型的平面裂纹扩展问题的半解析有限元法,利用弹性平面扇形域哈密顿体系的方程,通过分离变量法及共轭辛本征函数向量展开法,推导了一个环形和一个圆形奇异超级解析单元列式,组装这两个超级单元能准确地描述裂纹表面作用有双线性内聚力的平面裂纹尖端场。将该解析元与有限元相结合,构成半解析的有限元法,可求解任意几何形状和载荷的基于线性内聚力模型的平面裂纹扩展问题。典型算例的计算结果表明本文方法简单有效,具有令人满意的精度。     

Infinite spline boundary element method and its application to structure-foundation interaction

In view of the infinity behaviors of3-D Kelvin solution,we constructed an infinite spline boundary element which has fine precision in the analysis of the half space foundation subjected to uniform pressure on the circular domain.We also analysed a square plate resting on elastic half space foundation.The results indicate that this model not only fits for the coupled analysis of foundation and structures but also has the advantage of fewer degrees of freedom and fine precision.     

ALE FINITE ELEMENT ANALYSIS OF THE OPENING AND CLOSING PROCESS OF THE ARTIFICIAL MECHANICAL VALVE

ＡＬＥＦＩＮＩＴＥＥＬＥＭＥＮＴＡＮＡＬＹＳＩＳＯＦＴＨＥＯＰＥＮＩＮＧＡＮＤＣＬＯＳＩＮＧＰＲＯＣＥＳＳＯＦＴＨＥＡＲＴＩＦＩＣＩＡＬＭＥＣＨＡＮＩＣＡＬＶＡＬＶＥ（张建海）（陈大鹏）（邹盛铨）ＡＬＥＦＩＮＩＴＥＥＬＥＭＥＮＴＡＮＡＬＹＳＩＳＯＦＴ...     

A comparison of boundary element and finite element methods for modeling axisymmetric polymeric drop deformation

A modified boundary element method (BEM) and the DEVSS‐G finite element method (FEM) are applied to model the deformation of a polymeric drop suspended in another fluid subjected to start‐up uniaxial extensional flow. The effects of viscoelasticity, via the Oldroyd‐B differential model, are considered for the drop phase using both FEM and BEM and for both the drop and matrix phases using FEM. Where possible, results are compared with the linear deformation theory. Consistent predictions are obtained among the BEM, FEM, and linear theory for purely Newtonian systems and between FEM and linear theory for fully viscoelastic systems. FEM and BEM predictions for viscoelastic drops in a Newtonian matrix agree very well at short times but differ at longer times, with worst agreement occurring as critical flow strength is approached. This suggests that the dominant computational advantages held by the BEM over the FEM for this and similar problems may diminish or even disappear when the issue of accuracy is appropriately considered. Fully viscoelastic problems, which are only feasible using the FEM formulation, shed new insight on the role of viscoelasticity of the matrix fluid in drop deformation. Copyright © 2001 John Wiley & Sons, Ltd.     

Treatment of the small time instability in the finite element analysis of fluid structure interaction problems

In this paper, the fluid–structure interaction problem in mechanical systems in which a high frequency vibrating solid structure interacts with the surrounding fluid flow is considered. Such a situation normally appears in many microelectromechanical systems like a wide variety of microfluidic devices. A different implementation of the residual‐based variational multiscale flow method is employed within the arbitrary Lagrangian–Eulerian formulation. The combination of the variational multiscale method with appropriate stabilization parameters is used to handle the so‐called small time step instability in the finite element analysis of the fluid part in the coupled fluid–structure interaction problem. The capability of the employed approach has been demonstrated through finite element study of a benchmark example and FEM simulation of two different mechanical micropumping devices. High frequency vibrations of the solid membrane are used to derive the fluid flow in these micropumps. Results of FEM simulations are shown to be in good agreement with available experimental data.Copyright © 2012 John Wiley & Sons, Ltd.     

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

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

悬索桥结构分析中的索-鞍座单元

为在悬索桥结构的有限元分析中真实、简洁、高效地模拟索鞍,本文建立了一类新的单元。新单元包括索段的一端固定在与其接触部分为单一半径圆弧的索鞍上,另一端分别位于索鞍两侧的两节点“左索-鞍座单元”和“右索-鞍座单元”,以及索段两端点分别位于索鞍两侧,中间一点固定于鞍座上的三节点“索-鞍座单元”,后者的鞍槽可为两不同半径圆弧的组合。根据要求的成桥状态几何参数确定结构的无应力状态时,可利用前二者进行悬索桥的单跨分析。新单元通过自动调整索与鞍座的脱离点而处于平衡状态,从而简化了计算。单元算法的推导基于有限元分析的基本原理和弹性悬链线的精确解,并利用了处于平衡状态时索与鞍座之间的内力关系。新单元可象常规单元一样直接用于成桥状态或施工过程中悬索桥结构的有限元分析。设计的算例验证了新单元的正确性,并举例说明了新单元在悬索桥结构分析中的应用。     

Plotting isoclinics for hybrid photoelasticity and finite element analysis

Displacement-based finite element method formulations are coupled with stress-based photoelasticity analysis. As the stress field is discontinuous at the interelement boundaries, the introduced smoothing procedure enables the generation of high-quality digital images acceptable for hybird experimental-numerical techniques. The proposed methods are applicable for the analysis of static and dynamic results of experimental photoelasticity.     

Plastic limit analysis of incompatible finite element method

This paper describes an incompatible finite element model satisfying the consistencycondition of energy to solve the numerical precision problem of finite element solution inperfectly plastic analysis.In this paper the reason and criterion of the application of themodel to plastic limit analysis are discussed,and an algorithm of computing plastic limitload is given.