用光力学和有限元法确定钢壳的动态响应

本文提出了壳体在冲击载荷作用下三维位移分布的光力学方法的实验解和有限元数值解。采用一种新型的电磁脉冲加载装置,实现了对试件施加重复性很好的冲击载荷,利用测得的力一时间关系曲线作为外力边界条件,用逐步积分法对壳体在冲击载荷下的动态响应作了分析。光力学方法的实验解和有限元数值解吻合较好。     

可变网格的有限单元分析在边坡稳定计算中的应用

本文在非线性有限单元法的基础上，引进八节点的狭长单元来描述边坡中的滑动面，同时对于有限元单元网格可调节点坐标进行变更形成新的网格，于是利用原有基本网格架，通过多次有限元分析求出边坡中可能滑动面的稳定系数曲线Ｋ＝ｆ（α），并进而确定最危险的滑动面，在此基础上开展一系列稳定分析工作。     

基于静力响应面的结构有限元模型修正方法

提出了基于静力响应面的结构有限元模型修正方法.运用响应面方法,将结构静力响应和结构参数之间复杂的隐式关系用显式函数近似表达出来;在此响应面模型(函数)基础上,通过优化计算对结构有限元模型参数进行修正.阐述了基于静力响应面的结构有限元模型修正方法的基本理论和一般实现过程.对两跨连续梁结构的静力模型修正数值算例分析结果表明:基于静力响应面的有限元模型修正方法可以减少结构有限元计算的次数、提高模型修正的优化效率,结构有限元模型修正结果具有可接受的精度.     

Finite element-boundary integral formulation for electromagnetic scattering

A new formulation is described which combines the most robust attributes of the volume finite element and surface integral equation approaches to electromagnetic boundary value solutions. The result is a numerical technique which may be applied to scattering problems involving configurations having metallic surfaces and inhomogeneous penetrable material situated in open spatial regions. This is accomplished by way of coupling internal region finite element modal field solutions to equivalent currents on the surrounding boundary surface through an appropriate surface integral equation. The method is demonstrated for the special case of scattering by axisymmetric inhomogeneous penetrable objects. Example numerical calculations are presented for validation of the procedure and potential problem areas are discussed.     

SEMI-ELLIPTIC SURFACE CRACK IN AN ELASTIC SOLID WITH FINITE SIZE UNDER IMPACT LOADING

In this paper a semi-elliptic surface crack problem in an elastic solid of finite size under impact loading is investigated. An analysis is performed by means of fracture dynamics and the finite element method, and a three-dimensional finite element program is developed to compute the dynamic stress intensity factor. The results reveal that the effects of the solid's boundary surface, crack surface, material inertia and stress wave interactions play significant roles in dynamic fracture.     

Shape optimization of a body located in incompressible flow using adjoint method and partial control algorithm

The purpose of this study is to obtain an optimal shape of a body located in an incompressible viscous flow. The optimal shape of the body is defined so as to minimize the fluid forces acting on it by determining the surface coordinates based on the finite element method and the optimal control theory. The performance function, which is used to judge the optimality of a shape, is defined as the square sum of the drag and lift forces. The minimization problem is solved using an adjoint equation method. The gradient in the adjoint equation is affected by the finite element configuration. The use of a finite element mesh whose shape is appropriate for the procedure is important in shape optimization. If the finite element mesh used is not suitable for computations, the exact gradient is not calculated. Therefore, a structured mesh is used for the adjacent area of the body and all finite element meshes are refined using the Delaunay triangulation at each iteration computation. The weighted gradient method is applied as the minimization technique. Using an algorithm in which all nodal coordinates on the surface of the body are employed and starting from a circle as an initial shape, a front‐edged and rear‐round shape is obtained because of the vortices at the back of the body. To overcome this difficulty, we introduced the partial control algorithm, in which some of the nodal coordinates on the surface of the body are updated. From four cases of computational studies, we reveal that the optimal shape has both sharp front and sharp rear edges. All computations are conducted at Reynolds number Re=250. The minimum value of the performance function is obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamic modeling and control of a rotating piezo-bonded beam-mass system

A variational finite element method is used to derive the dynamic equations for a high speed rotating beam-mass system with embedded piezoelectric materials. Two piezoelectric layers are attached on the beam surfaces. The layer on the top surface is used for sensing and monitoring. The layer on the bottom surface is used for feedback control. The objective is to develop a simple finite element based dynamic model which can be used for structural characterization, monitoring and control. The structural dynamics introduced by the gyroscopic effect is included in the model. The nodal displacements of the structure are monitored by the output voltage level of the sensing layer. The voltage induced by the sensing layer is negatively fed back to the control layer for the adjustment of nodal displacements. Several simulation examples are presented to validate the effectiveness of the proposed finite element based control model.     

Surface cracks in a plate of finite width under extension or bending

In this paper the problem of a finite plate containing collinear surface cracks is considered. The problem is solved by using the line spring model with plane elasticity and Reissner's plate theory. The main purpose of the study is to investigate the effect of interaction between two cracks or between cracks and stress-free plate boundaries on the stress intensity factors and to provide extensive numerical results which may be useful in applications. First, some sample results are obtained and are compared with the existing finite element results. Then the problem is solved for a single (internal) crack, two collinear cracks and two corner cracks for wide range of relative dimensions. Particularly in corner cracks the agreement with the finite element solution is surprisingly very good. The results are obtained for semielliptic and rectangular crack profiles which may, in practice, correspond to two limiting cases of the actual profile of a subcritically growing surface crack.     

Dual variational formulation for Trefftz finite element method of elastic materials

A dual variational principle is presented for Trefftz finite element analysis. The proof of the stationary conditions of the variational functional and the theorem on the existence of extremum are provided in this paper. They are boundary displacement condition, surface traction condition and interelement continuity condition. Based on the assumed intraelement and frame fields, element stiffness matrix equation is obtained which can easily be implemented into computer programs for numerical analysis with Trefftz finite element method. Two numerical examples are considered to illustrate the effectiveness and applicability of the proposed element model.     

Topology optimization of piezoelectric nanostructures

We present an extended finite element formulation for piezoelectric nanobeams and nanoplates that is coupled with topology optimization to study the energy harvesting potential of piezoelectric nanostructures. The finite element model for the nanoplates is based on the Kirchoff plate model, with a linear through the thickness distribution of electric potential. Based on the topology optimization, the largest enhancements in energy harvesting are found for closed circuit boundary conditions, though significant gains are also found for open circuit boundary conditions. Most interestingly, our results demonstrate the competition between surface elasticity, which reduces the energy conversion efficiency, and surface piezoelectricity, which enhances the energy conversion efficiency, in governing the energy harvesting potential of piezoelectric nanostructures.     

A Newton method for three-dimensional fretting problems

The present paper concerns the numerical treatment of fretting problems using a finite element analysis. The governing equations resulting from a formal finite element discretization of an elastic body with a potential contact surface are considered in a quasi-static setting. The constitutive equations of the potential contact surface are Signorinis contact conditions, Coulombs law of friction and Archards law of wear. Using a backward Euler time discretization and an approach based on projections, the governing equations are written as an augmented Lagrangian formulation which is implemented and solved using a Newton algorithm for three-dimensional fretting problems of didactic nature. Details concerning the implementation are provided.     

爆炸冲击波作用下高层建筑基础与地基间的滑移与翘离

本文采用能反映基础与地基间接触面力学特性的弹塑性薄层有限连接单元与二维有限元相结合的方法,探讨了高层建筑在爆炸空气波作用下的动力反应,着重研究了基础与地基间产生的滑移与翘离效应。还探讨了有限元边界条件,结沟、基础与地基材料等因素对结构的动力反应、基础与地基间的滑移与翘离效应的影响。     

针对有砟道床动力特性分析的嵌入式离散元-有限元耦合方法

在离散元-有限元耦合方法中,离散元和有限元交界面处的耦合方式对整体有砟道床的力学行为影响显著.采用基于球形单元的镶嵌单元或粘结单元模拟有砟道床时,由于球形单元和有限单元表面的自锁能力较差,使道砟层在列车载荷作用下容易产生侧向滑移,导致数值模型不稳定.此外,在实际铁路道床中,底部道砟均不同程度地嵌入路堤.为此,发展了一种...     

铸件凝固温度场有限元分析中界面热阻的处理

提出一种处理铸件与铸型界面热阻问题的虚拟界面单元法，并给出了有限元计算公式。由于该公式不显含单元厚度(△l)，故该单元厚度△l可取任意值。当△l取为零时，使问题处理变得极为方便。针对某一具体金属型铝合金活塞的铸造凝固过程，按考虑和不考虑铸件与铸型间热阻影响两种方法作了有限元计算，通过与实测值相比较，本文提出的算法其计算精度远高于不考虑铸件与铸型间热阻影响的计算结果。另外．该方法使有限元建模方便、通用性强。     

Linear and Nonlinear Algorithms for Stress Separation in Photoelasticity

An experimental-numerical hybrid method for the stress separation in photoelasticity is proposed in this study. In the proposed method, boundary conditions for a local finite element model, that is, tractions along boundaries are inversely determined from photoelastic fringes. Two algorithms are proposed for determining the boundary condition. One is a linear algorithm in which the tractions are obtained by the method of linear least-squares from both principal stress difference and principal direction. Another is the nonlinear algorithm in which the tractions are determined only from the principal stress difference. After determining the boundary conditions for the local finite element model, the stresses can be obtained by finite element direct analysis. The effectiveness is demonstrated by applying the proposed method to a perforated plate under tension and contact problems. Results show that the boundary conditions of the local finite element model can be determined from the photoelastic fringes and then the individual stresses can be obtained by the proposed method. Furthermore, the stresses can be evaluated even if the boundary condition is complicated such as at the contact surface. It is expected that the proposed method can be powerful tool for stress analysis.     

The free surface flow of Newtonian and non-Newtonian fluids trapped by surface tension

In this paper the position of the free surface of a swirling fluid held in by surface tension is calculated by the finite element method. A new locally mass-conserving quadratic velocity, linear pressure triangular element is used to overcome non-physical solutions produced by the well known Taylor-Hood element.     

A finite element model for free surface flows on fixed meshes

In this paper, we present a finite element model for free surface flows on fixed meshes. The main novelty of the approach, compared with typical fixed mesh finite element models for such flows, is that we take advantage of the particularities of free surface flow, instead of considering it a particular case of two‐phase flow. The fact that a given free surface implies a known boundary condition on the interface, allows us to solve the Navier–Stokes equations on the fluid domain uncoupled from the solution on the rest of the finite element mesh. This, together with the use of enhanced integration allows us to model low Froude number flows accurately, something that is not possible with typical two‐phase flow models applied to free surface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

A nonlinear,transient finite element method for coupled solvent diffusion and large deformation of hydrogels

Hydrogels are capable of coupled mass transport and large deformation in response to external stimuli. In this paper, a nonlinear, transient finite element formulation is presented for initial boundary value problems associated with swelling and deformation of hydrogels, based on a nonlinear continuum theory that is consistent with classical theory of linear poroelasticity. A mixed finite element method is implemented with implicit time integration. The incompressible or nearly incompressible behavior at the initial stage imposes a constraint to the finite element discretization in order to satisfy the Ladyzhenskaya–Babuska–Brezzi (LBB) condition for stability of the mixed method, similar to linear poroelasticity as well as incompressible elasticity and Stokes flow; failure to choose an appropriate discretization would result in locking and numerical oscillations in transient analysis. To demonstrate the numerical method, two problems of practical interests are considered: constrained swelling and flat-punch indentation of hydrogel layers. Constrained swelling may lead to instantaneous surface instability for a soft hydrogel in a good solvent, which can be regulated by assuming a stiff surface layer. Indentation relaxation of hydrogels is simulated beyond the linear regime under plane strain conditions, in comparison with two elastic limits for the instantaneous and equilibrium states. The effects of Poisson’s ratio and loading rate are discussed. It is concluded that the present finite element method is robust and can be extended to study other transient phenomena in hydrogels.     

各向异性弹塑性中厚度板壳的有限元分析

本文在文献［２，３］的基础上，提出了一个解各向异性弹塑性中厚度板壳问题的有限元方法。考虑材料各向异性的特点，采用了Ｈｉｌｌ推广的Ｈｕｂｅｒ－Ｍｉｓｅｓ屈服准则；借用Ｏｗｅｎ的剪切修正系数，正确计及了叠层复合材料壳体的横向剪切效应；为了避免“自锁”现象，文中采用了９节点的Ｈｅｔｅｒｏｓｉｓ二次壳单元；特别是本文利用插值外推的思想，提出了一个带预测的弧长增量控制法，显著提高了确定变形路径的计算效率。几个数值算例表明本文给出的有限元方法对于各向异性中厚度板壳的弹塑性分析有较好的精度，尤其是对具有复杂变形路径的结构计算，收敛速度提高更快。