Motion of a vertical wall fixed on springs under the action of surface waves

A two-dimensional unsteady hydroelastic problem of interaction between surface waves and a moving vertical wall fixed on springs is studied. An analytical solution of the problem is constructed using a linear approximation, and a numerical solution within the framework of a nonlinear model of a potential fluid flow is found by a complex boundary element method. By means of analysis of the linear and nonlinear solutions, it is found that the linear solution can be used to predict some important characteristics of the wall motion and the fluid flow in the case of moderate wave amplitudes.     

A finite element approximation of steady flow through a rotating non-aligned straight tube

A finite element solution is developed for a penalty function formulation of the equations which govern the steady motion of a Newtonian fluid through a pipe that rotates about an axis not parallel to its own. The motion in this system is driven by the Coriolis acceleration, which has components in the axial direction as well as in the transverse plane of the pipe. The relative magnitudes of these components significantly affect the qualitative and quantitative nature of the primary and secondary flow field. The present results compare favourably with those of previously reported experimental and theoretical studies over a wide range of flow regimes.     

A layerwise boundary integral equation model for layers and layered media

A hybrid method is presented for the analysis of layers, plates, and multilayered systems consisting of isotropic and linear elastic materials. The problem is formulated for the general case of a multilayered system using a total potential energy formulation. The layerwise laminate theory of Reddy is employed to develop a layerwise, two-dimensional, displacement-based, hybrid boundary element model that assumes piecewise continuous distribution of the displacement components through the system's thickness. A one-dimensional finite element model is used for the analysis of the multilayered system through its thickness, and integral Fourier transforms are used to obtain the exact solution for the in-plane problem. Explicit expressions are obtained for the fundamental solution of a typical infinite layer (element) assuming linear displacement distribution through its thickness. This fundamental solution is given in a closed form in the cartesian space, and it can be applied in the two-dimensional boundary integral equation model to analyze layered structures with finite dimensions. The proposed method provides a simple, efficient, and versatile model for a three-dimensional analysis of thick plates or multilayered systems.     

A finite element analysis of a free surface drainage problem of two immiscible fluids

A sharp interface problem arising in the flow of two immiscible fluids, slag and molten metal in a blast furnace, is formulated using a two-dimensional model and solved numerically. This problem is a transient two-phase free or moving boundary problem, the slag surface and the slag–metal interface being the free boundaries. At each time step the hydraulic potential of each fluid satisfies the Laplace equation which is solved by the finite element method. The ordinary differential equations determining the motion of the free boundaries are treated using an implicit time-stepping scheme. The systems of linear equations obtained by discretization of the Laplace equations and the equations of motion of the free boundaries are incorporated into a large system of linear equations. At each time step the hydraulic potential in the interior domain and its derivatives on the free boundaries are obtained simultaneously by solving this linear system of equations. In addition, this solution directly gives the shape of the free boundaries at the next time step. The implicit scheme mentioned above enables us to get the solution without handling normal derivatives, which results in a good numerical solution of the present problem. A numerical example that simulates the flow in a blast furnace is given.     

Hydroelastic behavior of a complex structure floating on waves

A numerical method is developed to solve the plane problem of the hydroelastic behavior of a complex structure floating on the surface of an ideal incompressible fluid of finite depth. The motion of the structure described by a deflection function is considered steady-state under the action of incident waves. The hydrodynamic part of the problem is solved using the proposed approach based on the normal-mode method for homogeneous plates. The problem is reduced to a system of linear algebraic equations by means of a transition matrix between representations of the required deflection in the form of expansion in the vibration eigenfunctions of the structure and the plate. It is shown that the results of the calculation performed are in good agreement with available calculation results for a two-part hinged structure at wavelengths comparable to the length of the structure.     

The Motion of a Fluid–Rigid Disc System at the Zero Limit of the Rigid Disc Radius

We consider the two-dimensional motion of the coupled system of a viscous incompressible fluid and a rigid disc moving with the fluid, in the whole plane. The fluid motion is described by the Navier–Stokes equations and the motion of the rigid body by conservation laws of linear and angular momentum. We show that, assuming that the rigid disc is not allowed to rotate, as the radius of the disc goes to zero, the solution of this system converges, in an appropriate sense, to the solution of the Navier–Stokes equations describing the motion of only fluid in the whole plane. We also prove that the trajectory of the centre of the disc, at the zero limit of its radius, coincides with a fluid particle trajectory.     

Modeling micro-inertia in heterogeneous materials under dynamic loading

A continuum model including micro-inertia for heterogeneous materials under dynamic loading is proposed using a micro-mechanics method. The macro strain and stress are defined as the volume averages of the strain and stress fields in the representative volume element (RVE). The macro equations of motion are derived by using Hamilton’s principle together with the strain energy density and kinetic energy density involving the micro-inertia terms. The new macro equations of motion are used to study harmonic and transient wave propagation in layered media. Using a simple linear displacement field for the RVE, the dispersion curves obtained from the present model agree with the exact solutions very well for a range of wavelengths. The present model is also applied to analyze the transient response of layered media subjected to a triangular pulse loading. Comparison is made between the results of the present model and a finite element analysis.     

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

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

Finite element technique to solve the elastic strain for Leonov fluid flow

The finite element method is used to find the elastic strain (and thus the stress) for given velocity fields of the Leonov model fluid. With a simple linearization technique and the Galerkin formulation, the quasi-linear coupled first-order hyperbolic differential equations together with a non-linear equality constraint are solved over the entire domain based on a weighted residual scheme. The proposed numerical scheme has yielded efficient and accurate convective integrations for both the planar channel and the diverging radial flows for the Leonov model fluid. Only the strain in the inflow plane is required to be prescribed as the boundary conditions. In application, it can be conveniently incorporated in an existing finite element algorithm to simulate the Leonov viscoelastic fluid flow with more complex geometry in which the velocity field is not known a priori and an iterative procedure is needed.     

Highly elastic solutions for Oldroyd-B and Phan-Thien/Tanner fluids with a finite volume/element method: planar contraction flows

A hybrid finite volume/element method is analysed through the computation of creeping flows of viscoelastic fluids in plane 4:1 sharp and rounded-corner contraction geometries. Simulations are presented for three models: a constant viscosity Oldroyd-B fluid, and Phan-Thien/Tanner (PTT) shear thinning fluids of exponential and linear approximation form. A Taylor–Galerkin/pressure-correction scheme is implemented as the base time-stepping framework. The momentum equations are solved by a finite element method, whilst the constitutive equations are solved by a finite volume approach. Mesh convergence is analysed via refinement around the contraction to capture boundary layers and flow structure. Pressure drop is shown to increase with flow rate for a fixed fluid. For the Oldroyd-B model, singular behaviour is reported in the main stress component as one approaches the corner in the rounded, as with the sharp geometry. Velocity components display an asymptotic trend with a positive slope. Higher values of Weissenberg numbers (We) are reached with these finite volume schemes compared to their finite element counterparts, attributing this to superior accuracy properties.     

流固耦合地震波动问题的显式谱元模拟方法

流固耦合地震波动问题主要研究由流体和固体构成的复杂系统中地震波传播特性及其规律. 传统模拟方法中一般以声波方程、弹性波方程的数值解分别描述理想流体和弹性固体中的波动, 并实时地处理两种不同性质介质之间的相互耦合作用, 数值格式复杂且限制数值模拟精度与计算效率. 本文采用谱元法结合多次透射公式人工边界条件实现了一种流固耦合地震波动问题的高阶显式数值计算方法. 该方法利用了流固耦合问题统一计算框架,可将饱和多孔介质的Biot波动方程分别退化为理想流体的声波方程和弹性固体的弹性波方程. 通过P波垂直入射的水平成层理想流体-饱和多孔介质-弹性固体场地模型、P波斜入射的不规则层状界面以及任意形状界面的理想流体-饱和多孔介质-弹性固体场地模型等三个算例, 与传递函数法解析解以及集中质量有限元法计算结果进行对比分析, 证明了本文方法的正确性与有效性. 数值模拟结果表明, 本文方法相较传统有限元法可以少得多的节点数量获得更高的数值精度, 并且在较宽的频率范围内都能可靠地模拟出流固耦合系统的动力响应, 充分体现出本文方法兼顾高精度、计算效率和复杂场地建模灵活的特点.      

Finite elasto-plastic deformation—I: Theory and numerical examples

It is demonstrated that the problem of elasto-plastic finite deformation is governed by a quasi-linear model irrespective of deformation magnitude. This feature follows from the adoption of a rate viewpoint toward finite deformation analysis in an Eulerian reference frame. Objectivity of the formulation is preserved by introduction of a frame-invariant stress rate.Equations for piece-wise linear incremental finite element analysis are developed by application of the Galerkin method to the instantaneously linear governing differential equations of the quasi-linear model. Numerical solution capability has been established for problems of plane strain and plane stress. The accuracy of the numerical analysis is demonstrated by consideration of a number of problems of homogeneous finite deformation admiting comparative analytic solution. It is shown that accurate, objective numerical solutions can be obtained for problems involving dimensional changes of an order of magnitude and rotations of a full radian.     

双层旋转壳流固耦合振动有限元分析

本文采用旋转壳体元和流体元分析了双层旋转壳体流固耦合振动特性,基于Novozhilov壳体理论、水弹性理论以及Hamilton变分原理推导出耦合系统运动方程,计算结果与实测值比较表明本分析方法有较好的精度。     

The application of 2D base force element method (BFEM) to geometrically non-linear analysis

Based on the concept of the base forces by Gao, a new finite element method – the base force element method (BFEM) on complementary energy principle for two-dimensional geometrically non-linear problems is presented. A 4-mid-node plane element model of the BFEM for geometrically non-linear problem is derived by assuming that the stress is uniformly distributed on each sides of a plane element. The explicit formulations of the control equations for the BFEM are derived using the modified complementary energy principle. The BFEM is naturally universal for small displacement and large displacement problems. A number of example problems are solved using the BFEM and the results are compared with corresponding analytical solutions and those obtained from the standard displacement finite element method. A good agreement of the results, and better performance of the BFEM, compared to the displacement model, in the large displacement and large rotation calculations, is observed.     

Numerical analysis of finite amplitude motion of waves and a moored floating body under severe storm conditions

The motion of a moored floating body under the action of wave forces, which is influenced by fluid forces, shape of the floating body and mooring forces, should be analysed as a complex coupled motion system. Especially under severe storm conditions or resonant motion of the floating body it is necessary to consider finite amplitude motions of the waves, the floating body and the mooring lines as well as non-linear interactions of these finite amplitude motions. The problem of a floating body has been studied on the basis of linear wave theory by many researchers. However, the finite amplitude motion under a correlated motion system has rarely been taken into account. This paper presents a numerical method for calculating the finite amplitude motion when a floating body is moored by non-linear mooring lines such as chains and cables under severe storm conditions.     

Three‐dimensional sloshing: A consistent finite element approach

This paper presents a new computational methodology based on Legendre's polynomials to predict the slosh and acoustic motion in nearly incompressible fluids in both rigid and flexible structures with free surface. Here, we have used a finite element formulation based on Lagrangian frame of reference to model the fluid motion derived using Hamiltonian equation of the fluid system. We formulated three hexahedral finite elements based on strain fields expressed in terms of extended Legendre's polynomials. Sloshing and acoustic motion of liquid is investigated using these newly formulated elements and inf–sup test is performed on these new elements to check the performance of these elements in modeling sloshing under two severe constraints, namely incompressibility and irrotationality. Comparisons of slosh and acoustic frequencies, and mode shapes with exact solutions are given. Dynamic analysis with earthquake and harmonic kind of forcing function is carried out to validate the formulated hexahedral elements to analyze the sloshing response. Numerical results obtained with these new finite elements, and with the present finite element formulation of the mathematical model agree well with the exact solution and as well as with published experimental literature. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

An implicit three‐dimensional numerical model to simulate transport processes in coastal water bodies

A three‐dimensional baroclinic numerical model has been developed to compute water levels and water particle velocity distributions in coastal waters. The numerical model consists of hydrodynamic, transport and turbulence model components. In the hydrodynamic model component, the Navier–Stokes equations are solved with the hydrostatic pressure distribution assumption and the Boussinesq approximation. The transport model component consists of the pollutant transport model and the water temperature and salinity transport models. In this component, the three‐dimensional convective diffusion equations are solved for each of the three quantities. In the turbulence model, a two‐equation k–ϵ formulation is solved to calculate the kinetic energy of the turbulence and its rate of dissipation, which provides the variable vertical turbulent eddy viscosity. Horizontal eddy viscosities can be simulated by the Smagorinsky algebraic sub grid scale turbulence model. The solution method is a composite finite difference–finite element method. In the horizontal plane, finite difference approximations, and in the vertical plane, finite element shape functions are used. The governing equations are solved implicitly in the Cartesian co‐ordinate system. The horizontal mesh sizes can be variable. To increase the vertical resolution, grid clustering can be applied. In the treatment of coastal land boundaries, the flooding and drying processes can be considered. The developed numerical model predictions are compared with the analytical solutions of the steady wind driven circulatory flow in a closed basin and of the uni‐nodal standing oscillation. Furthermore, model predictions are verified by the experiments performed on the wind driven turbulent flow of an homogeneous fluid and by the hydraulic model studies conducted on the forced flushing of marinas in enclosed seas. Copyright © 2000 John Wiley & Sons, Ltd.     

The perturbation finite element method for solving the plane problem in consideration of linear creep

In this paper, a method (PFMC) for solving plane problem of linear creep is presented by using perturbation finite element. It can be used in plane problem in consideration of creep, such as reinforced concrete beam, presiressed concrete beam, reinforced concrete cylinder and reinforced concrete tunnel in elastic or visco-elastic medium, as well as underground building and so on.In the presented method, the assumption made in the general increment method that variables remain constant in a divided time interval is not taken. The accuracy is improved and the length of time step becomes larger. The computer storage can be reduced and the calculating efficiency can be increased.Perturbation finite element formulae for four-node quadrilateral isoparametric element including reinforcement are established and five numerical examples are given. As contrasted with the analytical solution, the accuracy is satisfactory.