地震作用下水-轴对称柱体相互作用的子结构分析方法

针对水-轴对称柱体动力相互作用问题,提出了一种地震作用下水-结构相互作用的时域子结构分析方法.基于三维不可压缩水体的波动方程和边界条件,利用分离变量法将其转换为环向解析、竖向和径向数值的二维模型;基于比例边界有限元推导了截断边界处无限域水体的动力刚度方程,并将水体内域有限元方程和人工边界处的动水压力进行耦合,从而得到结...     

恒压下树脂在双尺度多孔介质中的流动特性

基于复合材料液态模塑(LCM)工艺过程中存在半饱和区域的实验现象以及对预制体双尺度效应的逐步认识, 一些学者提出用沉浸模型来研究双尺度多孔介质的不饱和流动。通过体积均匀化方法描述了双尺度多孔介质复合材料液态模塑工艺模型的特征, 得到含有沉浸项的双尺度多孔介质的质量守恒方程, 并采用有限元法对方程进行数值求解, 通过具体算例计算了考虑双尺度效应时恒压树脂注射下不同时段的压力分布状态, 得到树脂在填充过程中流动前沿半饱和区域从出现到消失的过程, 采用不同注射压力进行模拟并比较。结果表明, 与单尺度多孔介质模型不同, 双尺度多孔介质模型更能反映实际树脂填充过程中出现的半饱和区域现象。     

A new finite element technique for the analysis of steady viscous flow problems

A new finite element technique for two-dimensional viscous incompressible fluid flow problems is presented in this paper. The vorticity transport equation is integrated in a small control volume, which results in the conservation law of vorticity. The finite element technique is applied to this equation together with the continuity equation, where simple linear triangular elements with three nodes are used for the formulation. Resulting sets of algebraic equations are solved by the use of a kind of relaxation method. Numerical results for viscous flow past a cavity show good agreement with experimental results.     

考虑构件弹性的电动Stewart并联机构刚度建模与仿真

基于集总解析建模方法和构件有限元分析建立包含驱动副、被动万向铰链和运动杆件弹性变形以及预载作用下的Stewart机构刚度矩阵模型。采用添加虚拟铰链等效构件弹性的方式,将分支等效为一系列刚性构件经由主、被动副以及虚拟铰链连接的形式,给出了运动关节和虚拟铰链变量对机构末端位姿的运动学Jacobian矩阵的数值计算方法,应用虚功原理得到静平衡方程,最终建立了机构无预载以及预载下的刚度矩阵模型。该模型不仅考虑了控制环路刚度,还将构件柔性的有限元分析结果与解析建模相结合,在降低计算成本的同时保证了精度。通过一机构分析实例,考察了两种模型下刚度分布的差异。     

树脂传递模塑工艺中的非饱和流动过程模拟与实验研究

针对编织类纤维增强体的纤维束之间与纤维束内孔隙的双尺度特点,建立了平纹织物的细观结构模型,并推导了汇函数的数学表达式。建立了局部细观流动特征的非饱和流动控制方程,利用有限元/控制体积方法求解,得到了局部饱和度分布。与实验进行比较,吻合较好。      

液体模塑成型工艺二维径向非饱和流动数值模拟

通过引入沉浸函数建立了双尺度多孔介质非饱和流动模型,并采用有限元/控制体积法实现了恒压及恒流注射条件下液体模塑成型（LCM）工艺二维径向非饱和流动的数值模拟,得到了不同注射条件下纤维织物内的压力场分布及半饱和区域长度随时间的变化规律,并将双尺度非饱和理论结果与单尺度饱和理论结果进行对比。结果表明：非饱和流动过程中,半饱和区域内的压力和压力梯度明显下降;半饱和区域长度随时间逐渐增加随后保持稳定,当流动前沿到达出口后半饱和区域长度开始逐渐减小;当两个主方向渗透率不同时,沿主方向半饱和区域长度也不同,渗透率越大该方向的半饱和区域长度也越大,纤维织物完全浸润时间取决于较小的渗透率。研究结果对合理预测树脂填充过程中压力分布及纤维预制件的浸润具有指导意义。     

INFLUENCE OF GEOMETRIC NON-LINEARITIES ON THE FREE VIBRATIONS OF ORTHOTROPIC OPEN CYLINDRICAL SHELLS

This paper presents a general approach to predict the influence of geometric non-linearities on the free vibration of elastic, thin, orthotropic and non-uniform open cylindrical shells. The open shells are assumed to be freely simply supported along their curved edges and to have arbitrary straight edge boundary conditions. The method is a hybrid of finite element and classical thin shell theories. The solution is divided into two parts. In part one, the displacement functions are obtained from Sanders' linear shell theory and the mass and linear stiffness matrices are obtained by the finite element procedure. In part two, the modal coefficients derived from the Sanders–Koiter non-linear theory of thin shells are obtained for these displacement functions. Expressions for the second- and third-order non-linear stiffness matrices are then determined through the finite element method. The non-linear equation of motion is solved by the fourth-order Runge–Kutta numerical method. The linear and non-linear natural frequency variations are determined as a function of shell amplitudes for different cases. The results obtained reveal that the frequencies calculated by this method are in good agreement with those obtained by other authors. © 1997 by John Wiley & Sons, Ltd.     

基于能量有限元法的功能梯度梁振动分析

功能梯度材料(Functionally Graded Material,FGM)由于其优良的结构性能和重要的应用价值,近些年来得到了广泛的研究和关注。采用能量有限元法对功能梯度梁和耦合梁的弯曲振动特性进行研究,推导了功能梯度材料梁的能量密度控制方程、能量有限元矩阵方程以及耦合梁的能量有限元方程,从而得到梁中的能量密度和能量流。以一简支功能梯度梁为例,分别采用该方法和传统有限元法计算了梁弯曲振动时的能量密度,通过对比验证了能量有限元法求解的准确性。在此基础上进一步对耦合功能梯度梁结构的能量密度和能量流进行了求解,得到其能量分布特征。该研究为基于能量有限元法分析复杂功能梯度材料结构的振动特性提供了理论基础。     

Bruce Irons: A non-conforming engineering scientist to be remembered and rediscovered

A method is described which calculates the fluid flow pattern for any generalized Newtonian fluid forced through the narrow annular space between two concentric cylindrical surfaces. The inner surface is engraved with a number of identical equally spaced spiral channels of small constant pitch angle but of varying cross-sectional area. For equal flow into the spiral passages, the pressure distribution along them is described by a quasi-linear ordinary difference-differential equation. A finite element approximation is employed to solve this equation for the pressure and resultant flow field. The method is illustrated by calculating and comparing flow patterns through two different spiral geometries.     

Topology optimization for submerged buoyant structures

This paper presents an evolutionary structural topology optimization method for the design of completely submerged buoyant modules with design-dependent fluid pressure loading. This type of structure is used to support offshore rig installation and pipeline transportation at all water depths. The proposed optimization method seeks to identify the buoy design that has the highest stiffness, allowing it to withstand deepwater pressure, uses the least material and has a minimum prescribed buoyancy. Laplace's equation is used to simulate underwater fluid pressure, and a polymer buoyancy module is considered to be linearly elastic. Both domains are solved with the finite element method. Using an extended bi-directional evolutionary structural optimization (BESO) method, the design-dependent pressure loads are modelled in a straightforward manner without any need for pressure surface parametrization. A new buoyancy inequality constraint sets a minimum required buoyancy effect, measured by the joint volume of the structure and its interior voids. Solid elements with low strain energy are iteratively removed from the initial design domain until a certain prescribed volume fraction. A test case is described to validate the optimization problem, and a buoy design problem is used to explore the features of the proposed method.     

树脂膜熔渗工艺充模过程的模拟与分析

树脂膜熔渗工艺(RFI)是一种新型的复合材料成型工艺.为了更深入了解树脂膜熔渗工艺过程中充模阶段的控制参数对制品质量的影响,避免制品出现空斑、充模不完全等问题,针对该工艺过程中树脂在复杂形状预制件中的流动行为进行了分析,在达西定律基础上建立了二维等温流动控制方程,利用有限元/控制体方法建立了数值分析模型,编制了FORTRAN程序进行模拟运算,并讨论了流动过程中施加的压力对充模时间的影响.由计算实例可见,所编制程序能够很好地预测树脂膜熔渗工艺过程中充模时间、各个时刻树脂的流动前沿位置及模腔中的压力分布.     

Stabilized finite element method for viscoplastic flow: formulation with state variable evolution

A stabilized, mixed finite element formulation for modelling viscoplastic flow, which can be used to model approximately steady‐state metal‐forming processes, is presented. The mixed formulation is expressed in terms of the velocity, pressure and state variable fields, where the state variable is used to describe the evolution of the material's resistance to plastic flow. The resulting system of equations has two sources of well‐known instabilities, one due to the incompressibility constraint and one due to the convection‐type state variable equation. Both of these instabilities are handled by adding mesh‐dependent stabilization terms, which are functions of the Euler–Lagrange equations, to the usual Galerkin method. Linearization of the weak form is derived to enable a Newton–Raphson implementation into an object‐oriented finite element framework. A progressive solution strategy is used for improving convergence for highly non‐linear material behaviour, typical for metals. Numerical experiments using the stabilization method with hierarchic shape functions for the velocity, pressure and state variable fields in viscoplastic flow and metal‐forming problems show that the stabilized finite element method is effective and efficient for non‐linear steady forming problems. Finally, the results are discussed and conclusions are inferred. Copyright © 2002 John Wiley & Sons, Ltd.     

Hybrid finite element/volume method for shallow water equations

A hybrid numerical scheme based on finite element and finite volume methods is developed to solve shallow water equations. In the recent past, we introduced a series of hybrid methods to solve incompressible low and high Reynolds number flows for single and two‐fluid flow problems. The present work extends the application of hybrid method to shallow water equations. In our hybrid shallow water flow solver, we write the governing equations in non‐conservation form and solve the non‐linear wave equation using finite element method with linear interpolation functions in space. On the other hand, the momentum equation is solved with highly accurate cell‐center finite volume method. Our hybrid numerical scheme is truly a segregated method with primitive variables stored and solved at both node and element centers. To enhance the stability of the hybrid method around discontinuities, we introduce a new shock capturing which will act only around sharp interfaces without sacrificing the accuracy elsewhere. Matrix‐free GMRES iterative solvers are used to solve both the wave and momentum equations in finite element and finite volume schemes. Several test problems are presented to demonstrate the robustness and applicability of the numerical method. Copyright © 2010 John Wiley & Sons, Ltd.     

A control volume mixed finite element formulation for tetrahedral elements

The control‐volume mixed finite element method is formulated for and applied to a computational domain consisting of a tetrahedral partitioning to solve the steady groundwater flow equations. Test functions consistent with piecewise constant and piecewise linear pressure distributions are used in the formulation. Comparisons are made with a standard mixed finite element formulation using lowest‐order Raviart Thomas basis functions. Results suggest that the control‐volume based formulation is a viable alternative to the standard formulation. Copyright © 2015 John Wiley & Sons, Ltd.     

Application of the control volume mixed finite element method to a triangular discretization

A two‐dimensional control volume mixed finite element method is applied to the elliptic equation. Discretization of the computational domain is based in triangular elements. Shape functions and test functions are formulated on the basis of an equilateral reference triangle with unit edges. A pressure support based on the linear interpolation of elemental edge pressures is used in this formulation. Comparisons are made between results from the standard mixed finite element method and this control volume mixed finite element method. Published 2011. This article is a US Government work and is in the public domain in the USA.     

考虑剪滞变形时箱形梁广义力矩的数值分析

为了简化变截面箱梁等复杂结构的剪滞效应分析,在明确定义相应于剪滞位移的广义力矩和有关几何特性的基础上,提出一种梁段有限元数值分析方法。选取控制微分方程的齐次解作为单元位移函数,以各积分常数为中间转换变量,推导梁段单元刚度矩阵和等效节点力向量的具体表达式,并给出用单元节点力直接计算应力的一般公式。编制了箱梁梁段有限元程序,对简支、悬臂、连续箱梁3个有机玻璃模型进行计算并与实测结果对比,验证了该文方法及公式的正确性。用所编程序对箱梁的剪滞广义力矩进行数值分析,并揭示了其变化规律。研究表明,在竖向荷载作用下,剪滞力矩与弯矩具有相似的分布规律,而且数值大小也接近。     

Thermal analysis of induction motors using 3D finite element method

Abstract

To design a reliable and economical induction motor, it is necessary to be able to predict accurately the temperature distribution within the motor. In this paper, a 3D thermal model of an induction motor is presented. Except for providing a more accurate representation of the problem, the proposed model can also reduce computer memory and time. The finite element method (FEM) is used to analyze the three dimensional (3D) heat flow equation which describes the thermal model. Galerkin's procedure is used to derive the element equations and first order tetrahedral elements are used to discretize the field region. Galerkin's time‐stepping scheme is employed to treat time differential terms. Values of surface heat transfer coefficients are obtained from the empirical formula and heat losses are revised by the factory test. Application of the proposed method to the analysis of a 9,000 HP induction motor yields temperature distribution very close to the experimental data.     

Elastic buckling analysis of 3-D trusses and frames with thin-walled members

 A finite element method is presented for the determination of the elastic buckling load of three-dimensional trusses and frames with rigid joints. The beam element stiffness matrix is constructed on the basis of the exact solution of the governing equations describing the coupled flexural-torsional buckling behaviour of a three-dimensional beam with an open thin-walled section in the framework of a small deformation theory. Large deformation effects are taken into account approximately through consideration of P−Δ effects. The structural stiffness matrix is obtained by an appropriate superposition of the various element stiffness matrices. The axial force distribution in the members is obtained iteratively for every value of the externally applied loading and the vanishing of the determinant of the structural stiffness matrix is the criterion used to numerically determine the elastic buckling load of the structure. The effect of initial member imperfections is also included in the formulation. Comparisons of accuracy and efficiency of the present exact finite element method against the conventional approximate finite element method are made. Cases where the axial force distribution determination can be done without iterations are also identified. The effect of neglecting the warping stiffness of some mono-symmetric sections is also investigated. Numerical examples involving simple and complex three-dimensional trusses and frames are presented to illustrate the method and demonstrate its merits. Received: 2 May 2000 / Accepted: 15 July 2002     

Explicit finite element artificial boundary scheme for transient scalar waves in two‐dimensional unbounded waveguide

To simulate the transient scalar wave propagation in a two‐dimensional unbounded waveguide, an explicit finite element artificial boundary scheme is proposed, which couples the standard dynamic finite element method for complex near field and a high‐order accurate artificial boundary condition (ABC) for simple far field. An exact dynamic‐stiffness ABC that is global in space and time is constructed. A temporal localization method is developed, which consists of the rational function approximation in the frequency domain and the auxiliary variable realization into time domain. This method is applied to the dynamic‐stiffness ABC to result in a high‐order accurate ABC that is local in time but global in space. By discretizing the high‐order accurate ABC along artificial boundary and coupling the result with the standard lumped‐mass finite element equation of near field, a coupled dynamic equation is obtained, which is a symmetric system of purely second‐order ordinary differential equations in time with the diagonal mass and non‐diagonal damping matrices. A new explicit time integration algorithm in structural dynamics is used to solve this equation. Numerical examples are given to demonstrate the effectiveness of the proposed scheme. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite element solution to an elastica problem of beams

In this study, finite element solution procedures are developed for an elastica problem of inextensible beams. The element stiffness matrices are obtained by using Galerkin's method. Results of a numerical example compare reasonably well with those obtained by using the elliptical integral. Extensions of this research are currently being made to include membrane effects of beams and to develop plate element stiffness matrices for elastica problems of plate structures.