形状记忆合金管接头空间轴对称有限元分析

本文采用形状记忆合金（ＳＭＡ）的三维本构方程和有限变形理论，考虑拉、压不同应力状态对相变点移动的规律，编制了ＳＭＡ轴对称大变形的有限元程序，与单向拉伸下解析所得的应力、应变曲线相比，证实程序的正确性．文末计算一ＳＭＡ管接头，并指出按空间轴对称计算的必要性．     

轴对称体扭转问题的轴对称有限元模拟解

轴对称体的轴对称问题与扭转问题一向被认为是两个互相不能模拟的问题.前者的未知量与方程多于后者,形式也不相同.本文提出一种退化模拟方法.能够把扭转问题模拟为轴对称问题的一类特殊情况来解.一般的结构分析程序都能够分析轴对称问题,但轴对称体的扭转问题通常作为三维问题处理.按本文提出的方法,可用结构分析程序的轴对称分析功能模拟扭转分析.本文还给出模拟计算的算例.计算结果表明与理论解完全一致.本文对退化模拟的材料本构关系进行了研究,建议在数值计算时以各向异性材料的轴对称问题模拟任何材料的扭转问题.当限定用各向同性材料的轴对称问题来模拟时,采用了罚系数法.     

含热传导的冲击动力学有限元程序的研究与应用

从有限变形的基本框架出发 ,建立了含热传导的冲击动力学基本控制方程和合适的初边条件。应用变分原理和伽辽金方法得出了控制方程相应的有限元列式 ,并探讨了数值计算中的几个关键算法 ,主要包括构型转换、旋转张量的算法、本构算法以及计算流程和程序框图等。利用自行研制的程序 ,对长脉冲激光辐照下靶目标的变形和破坏、冲击压缩变形及变形局部化等问题进行了数值模拟和研究 ,所得结果与实验和理论分析相吻合。     

轴对称冲击有限元一致质量矩阵迭代解

给出高速冲击动力有限元一致质量矩阵解的迭代过程，即把集总质量矩阵解作为初值进行有限次迭代得到满足工程需要的一致质量矩阵近似解．实际算例说明，一致质量近似解较集总质量解改善了对应力波传播过程的分析，而且在高速冲击计算中能给出与实验接近的计算结果．     

弹性力学轴对称问题的有限元线法

给出了解弹性力学空间轴对称问题的有限元线法的基本理论。该法包括了２－４条结线的等参数单元，沿结线方向的两点边值问题采用插值矩阵法解之。算例表明，本法具有良好的收敛性和较高的计算精度。     

压电材料轴对称有限元分析

给出了横观各向同性压电材料轴对称问题有限元法列式,叙述了压电材料有限元法分析中的一些特殊处理步骤,包括无量纲化和特征值问题计算中对电势自由度的凝聚。对压电材料轴对称问题的一些经典静力问题进行有限元计算,同时对压电圆板的轴对称自由振动进行了计算,和解析结果作了比较,二者符合良好。     

关于离散热传导物理模型的探讨

本文从热传导离散物理模型上论证了集中质量热容矩阵模型可以在离散点上满足热量守恒定律，采用它就可以避免许多热传导时间积分中的不合理现象。几个算例表明了该模型具有良好的精度。     

非均质材料热传导问题的扩展有限元法

针对非均质材料,提出了以导热系数为基本参数的热传导扩展有限元法。划分网格时不需要考虑材料界面的存在,因此网格的形成可以大大地简化,且可以获得高质量的网格。不含材料界面的单元,其温度场函数将退化为常规有限元的函数。含材料界面的单元,采用基于水平集的加强函数加强常规温度的近似,加强函数用于模拟界面。数值算例结果体现了该方法...     

近场动力学与有限元方法耦合求解热传导问题

求解含裂纹等不连续问题一直是计算力学的重点研究课题之一,以偏微分方程为基础的连续介质力学方法处理不连续问题时面临很大的困难. 近场动力学方法是一种基于积分方程的非局部理论,在处理不连续问题时有很大的优越性. 本文提出了求解含裂纹热传导问题的一种新的近场动力学与有限元法的耦合方法. 结合近场动力学方法处理不连续问题的优势以及有限元方法计算效率高的优势,将求解区域划分为两个区域,近场动力学区域和有限元区域. 包含裂纹的区域采用近场动力学方法建模,其他区域采用有限元方法建模. 本文提出的耦合方案实施简单方便,近场动力学区域与有限元区域之间不需要设置重叠区域. 耦合方法通过近场动力学粒子与其域内所有粒子(包括近场动力学粒子和有限元节点)以非局部方式连接,有限元节点与其周围的所有粒子以有限元方式相互作用. 将有限元热传导矩阵和近场动力学粒子相互作用矩阵写入同一整体热传导矩阵中,并采用Guyan缩聚法进一步减小计算量. 分别采用连续介质力学方法和近场动力学方法对一维以及二维温度场算例进行模拟,结果表明,本文的耦合方法具有较高的计算精度和计算效率. 该耦合方案可以进一步拓展到热力耦合条件下含裂纹材料和结构的裂纹扩展问题.      

Finite element analysis of convective heat transfer in porous media

The finite element method is used to analyse convective heat transfer in a porous medium. Convection past a vertical surface embedded in the medium and convection in a confined porous medium enclosure are analysed using the above method. The results are compared with those available in the literature and the agreement is found to be good. The method is applicable for two-dimensional analysis in a porous body of any arbitrary shape. The restriction of the boundary layer assumption is relaxed.     

Hybrid finite element/volume method for 3D incompressible flows with heat transfer

An implicit hybrid finite element (FE)/volume solver has been extended to incompressible flows coupled with the energy equation. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centred finite volume (FV) method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centres and the auxiliary variable at vertices. The Generalized Minimal Residual (GMRES) matrix-free strategy is adapted to solve the governing equations in both FE and FV methods. The presented 2D and 3D numerical examples show the robustness and accuracy of the numerical method.     

Streamline upwind finite element method for conjugate heat transfer problems

This paper presents a combined finite element method for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow. The streamline upwind finite element method is used for the analysis of thermal viscous flow in the fluid region, whereas the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the proposed method is to consistently couple heat transfer along the fluid-solid interface. Three test cases, i.e. conjugate Couette flow problem in parallel plate channel, counter-flow in heat exchanger, and conjugate natural convection in a square cavity with a conducting wall, are selected to evaluate the efficiency of the present method. The English text was polished byYunming Chen.     

Finite element solutions of laminar flow and heat transfer of air in a staggered and an in-line tube bank

A finite element method is used to solve the full Navier-Stokes and energy equations for the problems of laminar flow and heat transfer characteristics of air around three isothermal heated horizontal cylinders in a staggered tube bank and around four isothermal heated horizontal cylinders in an in line tube bank. The variations of surface shear stress, pressure and Nusselt number are obtained over the entire cylinder surface, including the zone beyond the separation point. The predicted values of total drag, pressure drag and friction drag coefficients, average Nusselt number, and the plots of velocity flow fields and isotherms are also presented.     

固体非傅立叶温度场的时域间断Galerkin有限元法

运用时域间断Galerkin有限元法[1],对高频非傅立叶热波动问题[2-3]进行分析。其主要特点是:取温度及温度的时间导数为基本未知量,对其分别进行3次Hermite插值和线性插值。在保证节点温度自动保持连续的基础上,温度的时间导数在离散时域存在间断。数值结果表明所提出的方法能够滤掉虚假的数值震荡,能够良好地模拟固体中的非傅立叶热波动行为。     

A comparison of boundary element and finite element methods for modeling axisymmetric polymeric drop deformation

A modified boundary element method (BEM) and the DEVSS‐G finite element method (FEM) are applied to model the deformation of a polymeric drop suspended in another fluid subjected to start‐up uniaxial extensional flow. The effects of viscoelasticity, via the Oldroyd‐B differential model, are considered for the drop phase using both FEM and BEM and for both the drop and matrix phases using FEM. Where possible, results are compared with the linear deformation theory. Consistent predictions are obtained among the BEM, FEM, and linear theory for purely Newtonian systems and between FEM and linear theory for fully viscoelastic systems. FEM and BEM predictions for viscoelastic drops in a Newtonian matrix agree very well at short times but differ at longer times, with worst agreement occurring as critical flow strength is approached. This suggests that the dominant computational advantages held by the BEM over the FEM for this and similar problems may diminish or even disappear when the issue of accuracy is appropriately considered. Fully viscoelastic problems, which are only feasible using the FEM formulation, shed new insight on the role of viscoelasticity of the matrix fluid in drop deformation. Copyright © 2001 John Wiley & Sons, Ltd.     

MHD axisymmetric flow of third grade fluid between porous disks with heat transfer

The magnetohydrodynamic(MHD) flow of the third grade fluid between two permeable disks with heat transfer is investigated.The governing partial differential equations are converted into the ordinary differential equations by suitable transformations.The transformed equations are solved by the homotopy analysis method(HAM).The expressions for square residual errors are defined,and the optimal values of convergencecontrol parameters are selected.The dimensionless velocity and temperature fields are examined for various dimensionless parameters.The skin friction coefficient and the Nusselt number are tabulated to analyze the effects of dimensionless parameters.     

轴对称非协调元分片检验检验函数的结构

导出了轴对称非协调元分片检验中检验函数的形式,其位移函数是一种不完全的线性函数。从而给出了轴对称非协调调元分析检验更准确,完全的描述。     

计及热传导影响对长杆弹侵彻陶瓷靶的数值分析

采用有限元方法离散瞬态热传导方程,编写成侵彻过程热传导计算模块,并将之嵌入已有的冲击动力学程序中,然后运用于长杆弹在900~1 800 m/s着速范围内侵彻AD95陶瓷靶的数值分析,得到了符合物理事实的计算图像,所得的计算结果比采用传统的绝热模型得到的计算结果更符合实验结果。探讨了计及热传导效应对长杆弹侵彻AD95陶瓷靶数值模拟的影响:着速在900~1 350 m/s范围内时,计及热传导的数值计算所得侵深小于绝热模型计算结果;着速在在1 350~1 450 m/s范围内时,两种模型计算侵深接近;着速在在1 450~1 800 m/s范围内时,热传导模型计算侵深大于绝热模型计算结果。