黏弹夹芯板的有限元建模及实验研究

从有限元分析和数值模拟及实验验证的角度研究了黏弹夹芯板的频率依赖振动特性。夹芯板中间层为黏弹性材料,其刚度和阻尼的频率依赖性行为直接影响系统的模态频率和阻尼,并导致振动模式求解的复杂化。采用三阶七参数Biot模型描述黏弹性材料频率相关的黏弹性行为。开发了三层四节点28自由度的夹芯板单元,基于经典板理论和哈密顿原理建立了黏弹夹芯板的有限元动力学方程。通过引入辅助耗散坐标,将Biot模型和黏弹夹芯板的有限元动力学模型结合起来,并将其转化为常规二阶线性系统形式,极大简化了求解非线性振动特性的过程。对一边固定、另三边自由的黏弹夹芯板进行了前三阶固有频率和损耗因子的预测,并与实验结果对比。数值模拟结果和实验结果吻合良好,说明所提有限元方法是正确有效的。     

分数阶导数粘弹性模型的有限元法

在分析分数阶导数三元件模型理论的基础上,把分数阶导数三元件模型引入有限元模型中,推导出具有分数阶导数三元件本构关系的粘弹性结构动力学有限元格式。同时,应用分数阶导数型粘弹性结构动力学方程的数值算法求解了该有限元格式的数值解。并以二维沥青路面结构为例进行了路面动态粘弹性响应分析。算例分析表明,该方法能够正确有效地进行路面动态粘弹性分析。     

Finite element modeling of tube hydroforming of polycrystalline aluminum alloy extrusions

Finite element modeling of tube hydroforming requires information about the anisotropy of the extruded aluminum tube. Unlike sheet metals, the complex geometry of extruded tubes makes it difficult, except in extrusion direction, to directly measure material properties. Therefore, polycrystalline models provide a good alternative for calculating the anisotropy of the tube in all directions and under various loading conditions. Using a rate-independent single crystal yield surface and rigid plasticity, a Taylor-type polycrystalline model was developed and implemented into ABAQUS/Explicit finite element (FE) code using VUMAT. The constitutive model was then used to calculate the crystallographic texture evolution during the hydroforming of an extruded aluminum tube. Initial crystallographic texture measured using orientation imaging microscopy (OIM) and uniaxial tensile test data obtained along the extrusion direction were input to this FEA model. In order to efficiently and practically simulate the tube hydroforming process using the polycrystalline model, sensitivity to the number of grain orientation, total simulation time, and number of finite elements were studied. Predicted results agreed very well with experimentally measured strain obtained from tube hydroforming process.     

Finite element modeling of contact conditions in multibody system dynamics

In this paper a new method is developed for the dynamic analysis of contact conditions in flexible multibody systems undergoing a rolling type of motion. The relative motion between the two contacting bodies is treated as a constraint condition describing their kinematic and geometric relations. Equations of motion of the system are presented in a matrix form making use of Kane's equations and finite element method. The method developed has been implemented in a general purpose program called DARS and applied to the simulation and analysis of a rotating wheel on a track. Both the bodies are assumed flexible and discretized using a three dimensional 8-noded isoparametric elements. The time variant constraint conditions are imposed on the nodal points located at the peripheral surfaces of the bodies under consideration. The simulation is carried out under two different boundary conditions describing the support of the track. The subsequent constraint forces associated with the generalized coordinates of the system are computed and plotted. The effects of friction are also discussed.     

Finite element solution to passive scalar transport behind line sources under neutral and unstable stratification

This study employed a direct numerical simulation (DNS) technique to contrast the plume behaviours and mixing of passive scalar emitted from line sources (aligned with the spanwise direction) in neutrally and unstably stratified open‐channel flows. The DNS model was developed using the Galerkin finite element method (FEM) employing trilinear brick elements with equal‐order interpolating polynomials that solved the momentum and continuity equations, together with conservation of energy and mass equations in incompressible flow. The second‐order accurate fractional‐step method was used to handle the implicit velocity–pressure coupling in incompressible flow. It also segregated the solution to the advection and diffusion terms, which were then integrated in time, respectively, by the explicit third‐order accurate Runge–Kutta method and the implicit second‐order accurate Crank–Nicolson method. The buoyancy term under unstable stratification was integrated in time explicitly by the first‐order accurate Euler method. The DNS FEM model calculated the scalar‐plume development and the mean plume path. In particular, it calculated the plume meandering in the wall‐normal direction under unstable stratification that agreed well with the laboratory and field measurements, as well as previous modelling results available in literature. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite element assembly strategies on multi‐core and many‐core architectures

We demonstrate that radically differing implementations of finite element methods (FEMs) are needed on multi‐core (CPU) and many‐core (GPU) architectures, if their respective performance potential is to be realised. Our numerical investigations using a finite element advection–diffusion solver show that increased performance on each architecture can only be achieved by committing to specific and diverse algorithmic choices that cut across the high‐level structure of the implementation. Making these commitments to achieve high performance for a single architecture leads to a loss of performance portability. Data structures that include redundant data but enable coalesced memory accesses are faster on many‐core architectures, whereas redundancy‐free data structures that are accessed indirectly are faster on multi‐core architectures. The Addto algorithm for global assembly is optimal on multi‐core architectures, whereas the Local Matrix Approach is optimal on many‐core architectures despite requiring more computation than the Addto algorithm. These results demonstrate the value in making the correct choice of algorithm and data structure when implementing FEMs, spectral element methods and low‐order discontinuous Galerkin methods on modern high‐performance architectures. Copyright © 2012 John Wiley & Sons, Ltd.     

Finite element analysis of forming limit in bore expanding of aluminium alloy sheets

Summary A recently proposed method to predict the forming limit of sheet metals is applied to the problem of bore expanding. Axisymmetric bore-expanding processes of various aluminium alloy sheets are simulated by the rigid-plastic finite element method. From the calculated histories of stress and strain, the forming limit, i.e. the fracture initiation, is predicted by means of the ductile fracture criterion. The comparison with the experimental results shows that the fracture initiation site and the critical punch stroke are successfully predicted by the present approach. Received 23 February 1998; accepted for publication 12 May 1998     

摄动有限元法在结构动力模型修改中的应用

本文将摄动理论与有限元法相结合,提出了用于小变参数结构分析的摄动有限元法(P-FEM),导出了在结构参数发生小变化的情况下,结构摄动单元矩阵的一般公式及结构的动特性随结构参数变化的二阶渐近展开式,并将这一结果运用于结构的动力模型修改中,提出了一种新的适合于工程应用的结构动力模型修改方法,把这一方法应用于实际的复杂结构动力模型修改中,获得了十分满意的结果.     

Finite element modeling of tire/terrain interaction: Application to predicting soil compaction and tire mobility

Tire/terrain interaction has been an important research topic in terramechanics. For off-road vehicle design, good tire mobility and little compaction on terrain are always strongly desired. These two issues were always investigated based on empirical approaches or testing methods. Finite element modeling of tire/terrain interaction seems a good approach, but the capability of the finite element has not well demonstrated. In this paper, the fundamental formulations on modeling soil compaction and tire mobility issues are further introduced. The Drucker-Prager/Cap model implemented in ABAQUS is used to model the soil compaction. A user subroutine for finite strain hyperelasticity model is developed to model nearly incompressible rubber material for tire. In order to predict transient spatial density, large deformation finite element formulation is used to capture the configuration change, which combines with soil elastoplastic model to calculate the transient spatial density due to tire compaction on terrain. Representative simulations are provided to demonstrate how the tire/terrain interaction model can be used to predict soil compaction and tire mobility in the field of terramechanics.     

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.     

近场动力学与有限元的混合建模方法

综述了近场动力学与有限元混合建模方法的研究进展,阐明了各种混合建模方法的基本原理与特点,并重点介绍本课题组在近场动力学与有限元方法混合建模方面的研究工作。现有近场动力学与有限元混合建模方法包括位移协调约束、力耦合、混合函数方法以及子模型方法等,除子模型方法外,都可归结为并行式多尺度分析方法,其基本思想是将计算结构划分为近场动力学子域、有限元子域以及两者的交界区域(或重叠区域、或界面单元、或过渡区域)。子模型方法可归结为显-显分析方法,先采用显式有限元进行整体分析,后采用近场动力学方法对重点区域进行分析。混合建模方法需要着重提高交界区域的计算精度,并且消除虚假力和虚假应力波问题。提出了通过力耦合的近场动力学与有限元混合建模的隐式分析方法,该方法不再设置重叠区,通过杆单元连接近场动力学子域与有限元子域,其中界面上的有限元结点不仅与其所在单元的其他结点发生作用,还通过杆单元与以其为圆心、一定半径的圆域内的其他物质点相互作用。研究表明,本文提出的混合模型和求解方法既能有效解决裂纹扩展等不连续问题,又可提高计算效率,为工程结构破坏问题的计算分析提供一种有效方法。     

Micromechanics of breakage in sharp-edge particles using combined DEM and FEM

By combining DEM （Discrete Element Method） and FEM （Finite Element Method）, a model is established to simulate the breakage of twodimensional sharp-edge particles, in which the simulated particles are assumed to have no cracks. Particles can, however, crush during different stages of the numerical analysis, if stress-based breakage criteria are fulfilled inside the particles. With this model, it is possible to study the influence of particle breakage on macro- and micro-mechanical behavior of simulated angular materials. Two series of tests, with and without breakable particles, are simulated under different confining pressures based on conditions of biaxial tests. The results, presented in terms of micromechanical behavior for different confining pressures, are compared with macroparameters. The influence of particle breakage on microstructure of sharp-edge materials is discussed and the related confining pressure effects are investigated. Breakage of particles in rockfill materials are shown to reduce the anisotropy coefficients of the samples and therefore their strength and dilation behaviors.     

Micromechanics of breakage in sharp-edge particles using combined DEM and FEM

By combining DEM (Discrete Element Method) and FEM (Finite Element Method),a model is established to simulate the breakage of two-dimensional sharp-edge particles,in which the simulated particles are assumed to have no cracks.Particles can,however,crush during different stages of the numerical analysis,if stress-based breakage criteria are fulfilled inside the particles.With this model,it is possible to study the influence of particle breakage on macro- and micro-mechanical behavior of simulated angular materials.Two series of tests,with and without breakable particles,are simulated under different confining pressures based on conditions of biaxial tests.The results,presented in terms of micromechanical behavior for different confining pressures,are compared with macroparameters.The influence of particle breakage on microstructure of sharp-edge materials is discussed and the related confining pressure effects are investigated.Breakage of particles in rockfill materials are shown to reduce the anisotropy coefficients of the samples and therefore their strength and dilation behaviors.     

基于SGMD及LWOA-ELM的有限元模型修正

为得到待修正参数与结构响应之间的关系,提高模型修正的效率和精度,提出了一种基于辛几何模态分解(SGMD)和Lévy飞行鲸鱼优化算法(LWOA)优化极限学习机(ELM)的有限元模型修正(FEMU)方法。首先,对加速度频响函数(AFRF)进行SGMD分解,采用能量熵增量法确定重组辛几何分量(SGC)构成SGC矩阵。然后,利用LWOA对ELM的权值和阈值进行优化,提高ELM模型的预测效率,以LWOA-ELM为代理模型映射出待修正参数与SGC矩阵之间的关系。最后,以试验频响函数SGC矩阵与LWOA-ELM模型输出所得矩阵差值的F-范数最小为目标函数,结合LWOA求解待修正参数。算例分析表明,提出的方法用于有限元模型修正有较好的可行性和有效性。以SGC矩阵表征AFRF的修正方法,有较好的噪声鲁棒性;LWOA-ELM作为代理模型预测精度高,泛化能力强。     

Discrete element modeling of acoustic emission in rock fracture

The acoustic emission (AE) features in rock fracture are simulated numerically with discrete element model (DEM). The specimen is constructed by using spherical particles bonded via the parallel bond model. As a result of the heterogeneity in rock specimen, the failure criterion of bonded particle is coupled by the shear and tensile strengths, which follow a normal probability distribution. The Kaiser effect is simulated in the fracture process, for a cubic rock specimen under uniaxial compression with a constant rate. The AE number is estimated with breakages of bonded particles using a pair of parameters, in the temporal and spatial scale, respectively. It is found that the AE numbers and the elastic energy release curves coincide. The range for the Kaiser effect from the AE number and the elastic energy release are the same. Furthermore, the frequency-magnitude relation of the AE number shows that the value of B determined with DEM is consistent with the experimental data.     

Finite element and experimental method for analyzing the effects of martensite morphologies on the formability of DP steels

Abstract

In this article, we investigated the effect of martensite morphology on the mechanical properties and formability of dual phase steels. At first, three heat treatment cycles were subjected to a low-carbon steel to produce ferrite–martensite microstructure with martensite morphology of blocky-shaped, continuous, and fibrous. Tensile tests were then carried out so as to study mechanical properties, particularly the strength and strain hardening behavior of dual phase steels. In order to study the formability of dual phase samples, Forming Limit Diagram was obtained experimentally and numerically. Experimental forming limit diagram was obtained using Nakazima forming test, while Finite Element Method was utilized to numerically predict the forming limit diagram. The results indicated that the dual phase samples with fibrous martensite morphology had the highest tensile properties and strain rate hardening out of the three different microstructures. Blocky-shaped martensite morphology, on the other hand, had the worst mechanical properties. The study of the strain hardening behavior of dual phase sample by Kocks–Mecking-type plots, evinced two stages of strain hardening for all specimens with different microstructures: stages III and IV. The forming limit diagram of dual phase steels also proved that samples with fibrous martensite morphology had the best formability compared to other two microstructures. The simulated forming limit diagram manifested that there is a good agreement between experimental results and those obtained by FEM.     

三维Cosserat连续体模型与微结构尺寸相关效应的有限元分析

分析了三维Cosserat连续体理论中的应力应变特征,推导了三维Cosserat连续体的有限元方程,基于ABAQUS计算软件提供的用户单元子程序(UEL)接口编写了弹性Cosserat连续体三维20节点有限元程序,并分析了微悬臂梁自由端的挠度问题和微杆扭转问题。通过与基于经典连续体理论的解析解及有限元数值计算结果进行比较,表明所发展的三维Cosserat连续体有限元能有效地模拟微结构尺寸相关效应问题,即随着微结构尺寸与材料内部长度参数的接近,基于Cosserat连续体有限元分析得到的微梁的挠度以及微杆的转角与经典连续体的解析解及有限元解相比越来越小;反之,Cosserat连续体有限元的计算结果与经典连续体的解析解及有限元数值解较为一致。     

Finite element study of free surface turbulent flow in rotating channels

Under certain conditions of liquid flow through rotating channels, the Coriolis force can induce a free surface to be formed. This problem is of practical importance in a Coriolis wear tester, which is used for determining the sliding wear coefficient of wear materials in slurry handling equipment. A deforming Galerkin finite element method is presented for predicting two‐dimensional turbulent free surface mean flow in rotating channels. Reynolds‐averaged Navier–Stokes (RANS) equations are cast into weak(algebraic) form using primitive variables (velocity and pressure). Eddy viscosity is determined via a mixing length model. Velocity is interpolated biquadratically, while pressure is interpolated bilinearly. The kinematic condition is used to form the Galerkin residual for the free surface. The free surface is represented by Hermite polynomials of zeroeth order for continuity of position and slope. Combined Newton's iteration is used to simultaneously solve for the free surface and the field variables. Results of velocity and pressure fields, as well as the free surface are shown to converge with mesh‐size refinement. There is excellent respect for mass conservation. Results are presented for various values of Rossby number (Ro) and height‐based Reynolds number (Re_H). Parameter continuation in Ro and Re_H space is used to compute solutions at higher values of flow rate and angular velocity. Copyright © 2002 John Wiley & Sons, Ltd.     

光滑拉伸试件中不同初始形状孔洞长大的有限元模拟

本文对具有不同硬化指数（ｎ＝０．０５，ｎ＝０．１，ｎ＝０．２）的幂硬化材料的光滑拉伸试样的拉伸变形过程进行了有限元模拟，通过有限元计算和经Ｂｒｉｄｇｍａｎ修正分别得到了试样变形过程中心部应力三维度随应变的变化情况；在此基础上运用控制体胞宏观应力三维度的方法，对含不同初始形状孔洞的体胞模型进行了有限元分析，计算结果表明：（１）孔洞初始形状和材料硬化特性对试样的拉伸破坏过程有重要影响；（２）Ｂｒｉｄｇ