On vibration of hemispherical shell by using finite element method

This paper establishes the finite element equation for the spherical shell.The resonantfrequencies of the above shell under different boundary conditions are also discussed andcalculated.     

A finite element method study of Eulerian droplets impingement models

To compute droplet impingement on airfoils, an Eulerian model for air flows containing water droplets is proposed as an alternative to the traditional Lagrangian particle tracking approach. Appropriate boundary conditions are presented for the droplets equations, with a stability analysis of the solution near the airfoil surface. Several finite element formulations are proposed to solve the droplets equations, based on conservative and non‐conservative forms and using different stabilization terms. Numerical results on single and multi‐elements airfoils for droplets of mean volume diameter, as well as for a Langmuir distribution of diameters, are presented and validated against measured values. Copyright © 1999 John Wiley & Sons, Ltd.     

A direct finite element method study of generalized thermoelastic problems

By using the principle of virtual work, the finite element equations corresponding to the generalized thermoelasticity with two relaxation times (i.e., the G–L theory) are derived. In order to improve the accuracy of integral transform method, especially in two/three-dimension, the equations are solved directly in time-domain. As numerical examples, the developed method is used to investigate the generalized thermoelastic behavior of a slim strip and a semi-infinite plate subject to thermal shock loading. The results demonstrate that the developed method can faithfully predict the deformation of the structure and most importantly the delicate second sound effect of heat conduction in both one and two-dimensional solids whilst it is usually difficult to model, especially in two-dimensional, by using the transform method. This method can also be used to study generalized piezothermoelastic and magnetothermoelastic problems.     

基于小波有限元的车辆-轨道-桥梁系统竖向振动分析

基于小波有限元建立了车辆-轨道-桥梁系统竖向运动方程。将车辆、轨道和桥梁作为一个整体系统,钢轨和桥梁采用区间B样条小波单元离散,钢轨与桥梁之间的钢轨基础采用均布的弹簧和阻尼模拟,采用虚功原理建立了基于小波有限元的四轴车辆-轨道-桥梁竖向振动分析模型。结果表明,采用区间B样条小波单元可较大程度上减小系统的自由度数,缩减计算量,节省计算时间。     

Large amplitude free flexural vibration of rings using finite element approach

Here, the large amplitude free flexural vibrations of isotropic/laminated orthotropic rings are investigated, using a shear flexible curved beam element based on field consistency principle. A laminated refined beam theory is introduced for developing the element, which satisfies the interface transverse shear stress and displacement continuity, and has a vanishing shear stress on the inner and outer surfaces of the beam. The formulation includes in-plane and rotary inertia effects, and the non-linearity due to the finite deformation of the ring. The governing equations obtained using Lagrange's equations of motion are solved through the direct integration technique. Amplitude-frequency relationships evaluated from the dynamic response history are examined. Detailed numerical results are presented considering various parameters such as radius-to-thickness ratio, circumferential wave number and ovality for isotropic and laminated orthotropic rings. The nature and degree of the participation of various modes in non-linear asymmetric vibration of oval ring brought out through the present study are useful for accurate modelling of the closed non-circular structures.     

Reliability of elasto-plastic structure using finite element method

A solution of probabilistic FEM for elastic-plastic materials is presented based on the incremental theory of plasticity and a modified initial stress method. The formulations are deduced through a direct differentiation scheme. Partial differentiation of displacement, stress and the performance function can be iteratively performed with the computation of the mean values of displacement and stress. The presented method enjoys the efficiency of both the perturbation method and the finite difference method, but avoids the approximation during the partial differentiation calculation. In order to improve the efficiency, the adjoint vector method is introduced to calculate the differentiation of stress and displacement with respect to random variables. In addition, a time-saving computational method for reliability index of elastic-plastic materials is suggested based upon the advanced First Order Second Moment (FOSM) and by the usage of Taylor expansion for displacement. The suggested method is also applicable to 3-D cases. The project supported by the Research Grant Council of Hong Kong (HKUST 722196E, 6039197E), the National Natural Science Foundation of China(59809003) and the Foundation of University Key Teacher by the Chinese Ministry of Education     

Nonlinear vibration of symmetrically laminated composite skew plates by finite element method

Here, the large amplitude free flexural vibration behavior of symmetrically laminated composite skew plates is investigated using the finite element method. The formulation includes the effects of shear deformation, in-plane and rotary inertia. The geometric non-linearity based on von Kármán's assumptions is introduced. The nonlinear matrix amplitude equation obtained by employing Galerkin's method is solved by direct iteration technique. Time history for the nonlinear free vibration of composite skew plate is also obtained using Newmark's time integration technique to examine the accuracy of matrix amplitude equation. The variation of nonlinear frequency ratios with amplitudes is brought out considering different parameters such as skew angle, fiber orientation and boundary condition.     

Parametric study on supersonic flutter of angle-ply laminated plates using shear deformable finite element method

The influence of fiber orientation,flow yaw angle and length-to-thickness ratio on flutter characteristics of angle-ply laminated plates in supersonic flow is studied by finite element approach.The structural model is established using the Reissner-Mindlin theory in which the transverse shear deformation is considered.The aerodynamic pressure is evaluated by the quasi-steady first-order piston theory.The equations of motion are formulated based on the principle of virtual work.With the harmonic motion assumption,the flutter boundary is determined by solving a series of complex eigenvalue problems.Numerical study shows that (1) The flutter dynamic pressure and the coalescence of flutter modes depend on fiber orientation,flow yaw angle and length-to-thickness ratio;(2) The laminated plate with all fibers aligned with the flow direction gives the highest flutter dynamic pressure,but a slight yawing of the flow from the fiber orientation results in a sharp decrease of the flutter dynamic pressure;(3) The angle-ply laminated plate with fiber orientation angle equal to flow yaw angle gives high flutter dynamic pressure,but not the maximum flutter dynamic pressure;(4) With the decrease of length-to-thickness ratio,an adverse effect due to mode transition on the flutter dynamic pressure is found.     

A finite element method to solve the Navier-Stokes equations using the method of characteristics

We present a simple and efficient finite element method to solve the Navier-Stokes equations in primitive variables V, p. It uses (a) an explicit advection step, by upwind differencing. Improvement with regard to the classical upwind differencing scheme of the first order is realized by accurate calculation of the characteristic curve across several elements, and higher order interpolation; (b) an implicit diffusion step, avoiding any theoretical limitation on the time increment, and (c) determination of the pressure field by solving the Poisson equation. Two laminar flow calculations are presented and compared to available numerical and experimental results.     

A 3D finite element ALE method using an approximate Riemann solution

Arbitrary Lagrangian–Eulerian finite volume methods that solve a multidimensional Riemann‐like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangian–Eulerian SGH method that incorporates a multidimensional Riemann‐like problem. Two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well‐known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problem results are presented. Copyright © 2016 John Wiley & Sons, Ltd.     

A finite element study on the large amplitude flexural vibration characteristics of FGM plates under aerodynamic load

The large amplitude flexural vibration characteristics of functionally graded material (FGM) plates are investigated here using a shear flexible finite element approach. Material properties of the plate are assumed to be graded in the thickness direction according to a simple power-law distribution in terms of volume fractions of the constituents. The effective material properties are then evaluated based on the rule of mixture. The FGM plate is modeled using the first-order shear deformation theory based on exact neutral surface position and von Kármán’s assumptions for large displacement. The third-order piston theory is employed to evaluate the aerodynamic pressure. The governing equations of motion are solved by harmonic balance method to study the vibration amplitude of FGM plates under supersonic air flow. Thereafter, the non-linear equations of motion are solved using Newmark’s time integration technique to understand the flexural vibration behavior of FGM plates in time domain (simple harmonic or periodic or quasi-periodic). This work is new in the sense that it deals with the non-linear flutter characteristics of FGM plates under high supersonic airflow accounting for both the geometric and aerodynamic non-linearities. Some parametric study is conducted to understand the influence of these non-linearities on the flutter characteristics of FGM plates.     

Exact finite element method

In this paper,a new method,exact element method for constructing finite element,ispresented.It can be applied to solve nonpositive definite or positive definite partialdifferential equation with arbitrary variable coefficient under arbitrary boundarycondition.Its convergence is proved and its united formula for solving partial differentialequation is given.By the present method,a noncompatible element can be obtained and thecompatibility conditions between elements can be treated very easily.Comparing the exactelement method with the general finite element method with the same degrees of freedom,the high convergence rate of the high order derivatives of solution can be obtained.Threenumerical examples are given at the end of this paper,which indicate all results canconverge to exact solution and have higher numerical precision.     

Buckling and post-buckling of piezoelectric plates using finite element method

The finite element equations considering the geometrical nonlinearity of piezoelectric smart structures are derived based on the total Lagrange method under the assumption of weak coupling between electricity and mechanics. Buckling and post-buckling of piezoelectric-plate with various boundary conditions are investigated. The calculated results show that piezoelectric effects and external voltage can hardly affect the buckling and post-buckling characteristics of piezoelectric-plate under uniaxial pressure while the buckling caused by displacement in-plane has much to do with the electric field. The present work is supported by the National Natural Science Foundation of China (No. 59635140).     

A finite element method on the plane elastic material analysis

Based on the finite element displacement method, a finite element method on the analysis of mechanical behaviour of plane elastic materials is proposed in this paper. By using this method and the corresponding computational program, the material behaviour of any unknown plane elastic material can be determined and all the elastic constants can be calculated.     

A finite element method for liquid with memory

Steady, two-dimensional flows of viscoelastic liquid with memory are analyzed by means of a nonlinear integral constitutive equation, novel streamlined finite elements, Galerkin's method of weighted residuals, and Newton iteration. By their relation to streamlines the new elements allow, for the first time, full Newton iteration of the algebraic equation set to which the governing integrodifferential system reduces. Streamlines are computed simultaneously with velocity components and pressure, the primary unknowns; to track the history of liquid particles the deformation equations are solved analytically in a Protean system of coordinates that conforms to the streamlines. In the cases studied the Newton iteration converges quadratically to a creeping flow state of Weissenberg number (product of upstream wall shear rate and relaxation time of the liquid) up to about 20 in channel flow, 10 in film flow, and 2 in die-swell flow. Shear-thinning shrinks the domain of convergence; it also reduces die-swell. A slip boundary condition in the vicinity of the contact line extends the domain of convergence.     

Self-adaptive one-dimensional nonlinear finite element method based on element energy projection method

The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.     

拟纯概率有限元法的研究

本文以提高复杂结构可靠性分析精度为主要目的，通过拟纯概率思想建立的三参数可靠性模型，对拟纯概率有限元方法的基本理论及计算进行了研究。文中以材料性能参数、几何参数和外载荷等为随机变量，导出了三维概率有限元法基本计算表达式。应用本文方法计算了纯弯曲下标准缺口试件的可靠度，并与实验结果和常规的二参数可靠性分析计算结果进行了对比。结果表明，本文提出的方法具有更高的计算精度。