Adaptive nodeless variable finite elements with flux-based formulation for thermal–structural analysis

A nodeless variable element method with the fluxbased formulation is developed to analyze two-dimensional thermal-structural problems. The nodeless variable formula- tion provides accurate temperature distributions to yield more accurate thermal stress solutions. The flux-based formulation is used to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method. The solution accuracy is further improved by implementing an adaptive meshing technique to generate finite element meshes that can adapt and move along with the transient solution behavior. A version of a nearly optimal element size determination is proposed to provide high convergence rate of the predicted solutions. The combined procedure is evaluated by solving several thermal, structural, and thermal stress problems.     

Hybrid-Trefftz有限元法的研究进展

系统概述了Ｈｙｂｒｉｄ Ｔｒｅｆｆｔｚ有限元法及其在弹性力学中的应用．该单元模型由于在插值函数上的灵活选择性使其比普通有限元能更有效地处理局部效应问题,如孔洞,集中荷载等．通过适当选择单元插值函数可构造出高精度的ｐ－扩展元和一系列满足特殊条件的新单元,以在同等条件下提高计算精度      

A Finite Element Analysis of a Viscoplastic Plate Subjected to Rapid Heating

Abstract

This paper models the response of a thin metallic plate that is subjected to a rapid heat input. In order to accurately model plate response, both the dynamic mechanical and transient heat transfer problems must be solved. The solution is complicated by nonlinearities due to radiation boundary conditions and material inelasticity. Furthermore, the viscoplastic constitutive equations that model the mechanical material behavior are numerically stiff. Nonlinear finite element algorithms are developed for both heat transfer and mechanical analyses. The algorithms are both stable and efficient for solving the problems considered herein. Example problems presented in the paper demonstrate the importance of including material nonlinearity in the model     

Combined adaptive meshing technique and characteristic-based split algorithm for viscous incompressible flow analysis

A combined characteristic-based split algorithm and an adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented.The method uses the three-node triangular element with equal-order interpolation functions for all variables of the velocity components and pressure.The main advantage of the combined method is that it improves the solution accuracy by coupling an error estima- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change in solution gradients,and at the same time,larger elements in the other regions.The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder,for their transient and steady-state flow behaviors.     

Thermal analysis of composite superconductors subjected to time-dependent disturbances

A one-dimensional heat conduction equation with time- and temperature-dependent heat sources was employed to study the steady-state and transient response of a composite superconductor subjected to a thermal disturbance. An integral formulation was used to solve the steady-state problem of current redistribution and heat generation. The results of the integral formulation are compared with those of an analytical solution. The two solutions agree with each other except when the analytical solution fails as the temperature in the superconductor begins to exceed the critical temperature. Transient solutions were obtained by the finite-difference technique and the results are compared with a known analytical solution. Results of numerical calculations of the transient response of a composite superconductor subjected to an initial pulsed disturbance are presented. It is demonstrated that the superconductor can switch between the superconducting and the current-sharing state. The transient response and the stability of the composite conductor depend on the magnitude and duration of the disturbance, the dimensionless temperature θ^*, and the dimensionless parameter φ. Received on 18 November 1996     

Transient heat conduction analysis using the NURBS-enhanced scaled boundary finite element method and modified precise integration method

The Non-uniform rational B-spline(NURBS)enhanced scaled boundary finite element method in combination with the modified precise integration method is proposed for the transient heat conduction problems in this paper.The scaled boundary finite element method is a semi-analytical technique,which weakens the governing differential equations along the circumferential direction and solves those analytically in the radial direction.In this method,only the boundary is discretized in the finite element sense leading to a reduction of the spatial dimension by one with no fundamental solution required.Nevertheless,in case of the complex geometry,a huge number of elements are generally required to properly approximate the exact shape of the domain and distorted meshes are often unavoidable in the conventional finite element approach,which leads to huge computational efforts and loss of accuracy.NURBS are the most popular mathematical tool in CAD industry due to its flexibility to fit any free-form shape.In the proposed methodology,the arbitrary curved boundary of problem domain is exactly represented with NURBS basis functions,while the straight part of the boundary is discretized by the conventional Lagrange shape functions.Both the concepts of isogeometric analysis and scaled boundary finite element method are combined to form the governing equations of transient heat conduction analysis and the solution is obtained using the modified precise integration method.The stiffness matrix is obtained from a standard quadratic eigenvalue problem and the mass matrix is determined from the low-frequency expansion.Finally the governing equations become a system of first-order ordinary differential equations and the time domain response is solved numerically by the modified precise integration method.The accuracy and stability of the proposed method to deal with the transient heat conduction problems are demonstrated by numerical examples.     

PRECISE INTEGRAL ALGORITHM BASED SOLUTION FOR TRANSIENT INVERSE HEAT CONDUCTION PROBLEMS WITH MULTI-VARIABLES

IntroductionIHCPs (InverseHeatConductionProblems)arecloselyassociatedwithmanyengineeringaspects,andwelldocumentedintheliteratures,coveringtheidentificationsofthermalparameters[1,2 ],boundaryshapes[3],boundaryconditions[4 ]andsource_relatedterms[5 ,6 ]etc .Howeveritseemsthatonlylittleworkisdirectlyconcernedwithmulti_variablesidentificationsbyauthors’knowledge.Tsengetal.presentedanapproachtodeterminingtwokindsofvariables[7],butonlygavefewnumericalexamplestodeterminethemsimultaneously .Thesol…     

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.     

基于特征分裂有限元准隐格式的共轭传热整体耦合数值模拟方法

共轭传热现象在科学和工程领域中大量存在.随着计算能力的发展,对共轭传热现象进行准确有效的数值模拟,成为科学研究和工程设计上的重要挑战.共轭传热数值模拟的方法可以分为两大类:分区耦合和整体耦合.本文采用有限元法对共轭传热问题进行整体耦合模拟.固体传热求解采用标准的伽辽金有限元方法.流动求解采用基于特征分裂的有限元方法.该...     

Thermomechanical formulation of ductile damage coupled to nonlinear isotropic hardening and multiplicative viscoplasticity

In this paper, we present a thermomechanical framework which makes use of the internal variable theory of thermodynamics for damage-coupled finite viscoplasticity with nonlinear isotropic hardening. Damage evolution, being an irreversible process, generates heat. In addition to its direct effect on material's strength and stiffness, it causes deterioration of the heat conduction. The formulation, following the footsteps of Simó and Miehe (1992), introduces inelastic entropy as an additional state variable. Given a temperature dependent damage dissipation potential, we show that the evolution of inelastic entropy assumes a split form relating to plastic and damage parts, respectively. The solution of the thermomechanical problem is based on the so-called isothermal split. This allows the use of the model in 2D and 3D example problems involving geometrical imperfection triggered necking in an axisymmetric bar and thermally triggered necking of a 3D rectangular bar.     

Heat conduction between confocal elliptical surfaces using the theory of a Cosserat shell

This paper is concerned with steady-state heat conduction in rigid shell-like interphase regions. By analogy this work may provide insight into related problems of electric, dielectric and magnetic behavior. Although the field equations for three-dimensional linear Fourier heat condition are rather simple, the solution of problems in shell regions is significantly complicated when the shell has a general geometry and variable thickness. Here, the problem of heat conduction between confocal elliptical surfaces is solved within the context of the theory of a Cosserat shell. This problem is of particular interest because the Cosserat solution can be compared with an exact solution and the influences of variable shell thickness and strong variations of the temperature field through the shell’s thickness can be explored independently. The results show that the Cosserat approach is reasonably accurate even for moderately thick shells, moderate ellipticity, and moderately strong variation of the temperature through the shell’s thickness.     

Absolute nodal coordinate formulation of large-deformation piezoelectric laminated plates

The Absolute Nodal Coordinate Formulation (ANCF) has been initiated in 1996 by Shabana (Computational Continuum Mechanics, 3rd edn., Cambridge: Cambridge University Press, 2008). It introduces large displacements of planar and spatial finite elements relative to the global reference frame without using any local frame. A sub-family of beam, plate and cable finite elements with large deformations are proposed and employed the 3D theory of continuum mechanics. In the ANCF, the nodal coordinates consist of absolute position coordinates and gradients that can be used to define a unique rotation and deformation fields within the element. In contrast to other large deformation formulations, the equations of motion contain constant mass matrices as well as zero centrifugal and Coriolis inertia forces. The only nonlinear term is a vector of elastic forces. This investigation concerns a way to generate new finite element in the ANCF for laminated composite plates. This formulation utilizes the assumption that the bonds between the laminae are thin and shear is non-deformable. Consequently, the Equivalent Single Layer, ESL model, is implemented. In the ESL models, the laminate is assumed to deform as a single layer, assuming a smooth variation of the displacement field across the thickness. In this paper, the coupled electromechanical effect of Piezoelectric Laminated Plate is imposed within the ANCF thin plate element, in such a way as to achieve the continuity of the gradients at the nodal points, and obtain a formulation that automatically satisfies the principle of work and energy. Convergence and accuracy of the finite-element ANCF Piezoelectric Laminated Plate is demonstrated in geometrically nonlinear static and dynamic test problems, as well as in linear analysis of natural frequencies. The computer implementation and several numerical examples are presented in order to demonstrate the use of the formulation developed in this paper. A comparison with the commercial finite element package COMSOL MULTIPHYSICS () is carried out with an excellent agreement.     

Mixed-mode thermal stress intensity factors from the finite element discretized symplectic method

A finite element discretized symplectic method is introduced to find the thermal stress intensity factors (TSIFs) under steady-state thermal loading by symplectic expansion. The cracked body is modeled by the conventional finite elements and divided into two regions: near and far fields. In the near field, Hamiltonian systems are established for the heat conduction and thermoelasticity problems respectively. Closed form temperature and displacement functions are expressed by symplectic eigen-solutions in polar coordinates. Combined with the analytic symplectic series and the classical finite elements for arbitrary boundary conditions, the main unknowns are no longer the nodal temperature and displacements but are the coefficients of the symplectic series after matrix transformation. The TSIFs, temperatures, displacements and stresses at the singular region are obtained simultaneously without any post-processing. A number of numerical examples as well as convergence studies are given and are found to be in good agreement with the existing solutions.     

基于高阶剪切变形理论的四边形求积元板单元及其应用

近年来由各类新型复合材料或功能梯度材料构成的板结构在工程领域得到了广泛应用,其显著特点是材料性能沿板厚变化.为合理考虑横向剪切应变,许多学者基于Reddy高阶剪切变形理论,构建了不同的有限元单元对该类板结构进行分析,但其中满足$C^{1}$连续条件的单元相对较少.本文基于Reddy高阶剪切变形理论,采用求积元方法,建立了$C^{1}$连续的四边形板单元.利用该单元对均质材料、复合材料、功能梯度材料构成的等厚度矩形板、变厚度矩形板及等厚度斜板的线弹性弯曲和自由振动问题进行了计算分析,并与现有文献中的相应计算结果进行了对比.研究表明:基于高阶剪切变形理论的四边形求积元板单元具有较高的计算效率和良好的适应性,文中各类材料构成的等变厚度矩形板及等厚度斜板均只需1个单元即可得到理想的计算结果.对于等/变厚度矩形板,可仅使用9$\times$9个积分点,而对于等厚度斜板,随着斜角的增大,所需积分点的数目逐渐增多至15$\times $15.该四边形求积元板单元可进一步用于新型复合材料板的非线性分析.      

Winkler地基上变厚度圆板的轴对称弯曲

本文提出了Winkler地基上变厚度圆板轴对称弯曲的传递矩阵算法。首先,根据贝塞尔函数理论获得了等厚度圆板和环板单元在任意荷载作用下轴对称弯曲的解析解,这些解均由通解和特解两部分组成。基于这些解析解,导出了等厚度圆板和环板单元的传递矩阵。然后沿径向将变厚度圆板划分成一个等厚度圆板单元和一系列等厚度环板单元,应用传递矩阵算法原理获得了变厚度圆板的整体传递矩阵。引入圆板的边界条件,给出了该板每条节线上的挠度、径向转角、径向弯矩和径向剪力。最后,讨论了受均布荷载作用的简支线性变厚度圆板的弯曲,将本文数值解与解析解进行比较,证实了本文方法的有效性,并简要地讨论了地基参数对板挠度和径向弯矩的影响。     

金属结构航天器陨落过程三维瞬态传热有限元算法研究

构建金属桁架结构航天器陨落再入气动热环境有限元传热模型,是准确预测在轨服役期满大型航天器陨落再入解体过程温度分布的关键。本文采用四节点四面体单元对空间进行离散,依据泛函理论,将传热控制方程离散为代数方程组;利用有限单元法总体合成得到具有对称正定、高度稀疏和非0元素分布的规则性刚度矩阵,发展一维变带宽压缩存贮技术,有效解决大型稀疏矩阵的数据存贮问题;为有效抑制求解过程出现的温度在时间和空间上的振荡问题,发展集中热容矩阵系数处理方法,将热容矩阵的同行或同列元素相加代替对角线元素,使非对角线元素化为0,构造求解三维瞬态温度场的两点向后差分格式、Crank-Nicolson格式和Galerkin格式。通过对正方体瞬态传热计算验证分析,在相同条件下,采用以上三种格式均可获得一致稳定的温度解,并得到与现有ANSYS有限元软件较为吻合的计算结果,验证了所建立三维瞬态传热有限元计算模型的准确可靠性。在此基础上,对铝合金低轨航天器薄壳结构进行了传热计算,给出了类天宫飞行器两舱体陨落飞行107.5 km~90 km不同高度的瞬态温度分布,为这类寿命末期航天器陨落再入解体预报提供理论支撑与可计算模型。     

Ground thermal response to heat conduction in a power transmission tower foundation

An analytical formulation is developed to predict transient heat conduction in a semi-infinite medium with a vertical finite line heat source, which represents a buried tower of a power transmission line foundation. Unlike past studies with a constant line heat source, the current model develops a time-dependent variable heating strength, as well as a time varying surface temperature of the ground. An approximate VHS model (variable heating strength) is developed for sinusoidal variations of the line source strength and surface temperature, in order to simulate seasonal variations of ground temperatures. The VHS model reduces computational time and exhibits good accuracy, when compared against a full exact solution. The model is applied to heat conduction in a tower foundation, with time-varying ground surface temperatures. Effects of ground thermal conductivity and diffusivity, as well as variations of the line source strength, are investigated in this article.     

Cavity detection in a heat conductor using linear sampling method

In this paper, solution of inverse problems in heat conduction transient fields is investigated. For this purpose, a new time-domain version of linear sampling method (TDLSM) is developed for cavity detection in a heat conductor. The linear sampling method (LSM) is an effective approach to image the geometrical features of unknown targets. Although this method has been used in the context of inverse scattering problems such as solid, acoustics, and electromagnetism, there is no specific attempt to apply this method to identification of cavities in heat conductors. This study emphasizes the implementation of the LSM in the time-domain fields using finite element method. A set of numerical simulations on two-dimensional transient heat conduction problems is presented to highlight many effective features of the proposed TDLSM fast qualitative identification method.     

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin (SUPG) method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid are performed by the Galerkin method. The second-order semiimplicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are linearized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effiectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.