Study on bending fretting fatigue of LZ50 axle steel considering ratchetting by finite element method

Based on the finite element simulation, the crack nucleation location and failure life of bending fretting fatigue are analysed and predicted for LZ50 axle steel, respectively. In the simulation, a simplified two‐dimensional finite element model with an equivalent normal force is proposed and the role of ratchetting is also considered by using a new cyclic elastoplastic constitutive model. The crack nucleation location and fatigue life are finally predicted by Smith–Watson–Topper (SWT) critical plane criteria referring to the finite element calculation. It is shown that the predicted results are in fairly good agreement with the experimental ones.     

预应力连续箱梁桥的基准动力有限元模型研究z

介绍了作为预应力混凝土连续梁桥的通扬运河特大桥的环境激励振动试验,并利用基     

Contact finite element analysis of three- and four-point short-beam bending of unidirectional composites

Detailed finite element stress analysis of both three-point and four-point short-beam bending specimens of a unidirectional glass-fibre/epoxy composite has been carried out. Contact elements are used to calculate the actual contact stresses and contact length at the roller. This gives more realistic stress distributions than those based on assumed load distributions. If the material is assumed to be linear, a parabolic assumption for the shear stress distribution is good for large parts of the beam between the loading and supporting rollers. However, higher shear stresses arise locally in the vicinity of the rollers. Taking account of the real nonlinear shear response of the material significantly reduces the shear stresses in the specimens.     

New finite element model for analysis of Kirchhoff plate

A new kind of approach to formulate an isotropic thin plate bending element is presented. The strain energy of the element is decomposed into two parts: an integral concerning the first strain invariant and a line integral around the elemental boundary. The former can be discretized by quasi-conforming technique¹ and the latter can be directly calculated using the interpolation of the deflection and its normal slope along the element boundary. By this method, an assumed first strain invariant quadrilateral element (AFSIQ) is derived. The procedure of formulating the element and the numerical examples show that the new kind of element not only simplifies the formulation considerably but also has excellent accuracy.     

A unified 3D‐based technique for plate bending analysis using scaled boundary finite element method

This paper presents a unified technique for solving the plate bending problems by extending the scaled boundary finite element method. The formulation is based on the three‐dimensional governing equation without enforcing the kinematics of plate theory. Only the in‐plane dimensions are discretised into finite elements. Any two‐dimensional displacement‐based elements can be employed. The solution along the thickness is expressed analytically by using a matrix function. The proposed technique is consistent with the three‐dimensional theory and applicable to both thick and thin plates without exhibiting the numerical locking phenomenon. Moreover, the use of higher order spectral elements allows the proposed technique to better represent curved boundaries and to achieve high accuracy and fast convergence. Numerical examples of various plate structures with different thickness‐to‐length ratios demonstrate the applicability and accuracy of the proposed technique. Copyright © 2012 John Wiley & Sons, Ltd.     

Finite element numerical simulation on thermo-mechanical behavior of steel billet in continuous casting mold

The thermo-mechanical behavior of a thin, growing shell during the early stages of solidification in a continuous casting mold is very important to the ultimate quality of the final billet. A two-dimensional, transient finite element model has been developed to treat the heat flow and deformation of the solidifying shell in the continuous casting billet mold as a coupled phenomena. The major application of the model is to predict the extent of the gap between the mold and the shell and focus on the influence of mold taper on the thermo-mechanical behavior of the steel billet to help to understand the formation of off-corner cracks and break-outs in the solidifying shell. The calculations indicate that the gap is initially formed at the corner of the billet, where heat transfer is greatly reduced. Insufficient mold taper contributes to a hot spot in the off-corner region, which corresponds to the lowest shell thickness. At the same time, the solidifying front on the diagonal of the billet is subjected to an excessive mechanical strain, which causes the off-corner cracks and even the break-outs.     

A 3D cylindrical finite element model for thick curved beam stress analysis

A finite element model is proposed to perform stress analysis for thick curved beams and panels subjected to various types of loadings. The model has 18 nodes in a three‐dimensional cylindrical co‐ordinates system. Three stress components on radial surface (σ_rr, τ_rθ, and τ_rz) and three displacement components (u_r, u_θ, and u_z) are used as nodal degrees of freedom. Therefore, the continuity condition for both stresses and displacements is achieved in the radial direction. Formulation of nodal shape functions and equilibrium equations are based on three‐dimensional elasticity theory and a minimum potential energy method. The accuracy of the method is verified with the standard test problems and exact solutions from the theory of elasticity. The model shows no locking phenomena. Convergence is investigated and the application to layered composite panel is illustrated. Copyright © 2003 John Wiley & Sons, Ltd.     

Finite element analysis of creep fracture initiation in a model superalloy material

Detailed finite element analyses were performed for a single edge-cracked specimen geometry under both plane stress and plane strain constraint for a superalloy material that obeys a power-law creep relationship. The objectives of these analyses were to elucidate the stationary creep crack-tip fields and to provide guidance for the experimental measurement of crack-tip deformations. New results demonstrate that, for both plane stress and plane strain, the angular variations in the creep strain fields do not agree with HRR-type predictions, although the radial variations are in agreement with HRR-type creep strain field predictions in a zone very near the crack tip. Thus, the use of experimental measurement of surface displacement and/or strain data for the location of HRR-type fields may not be possible, unless modifications to the existing HRR-type theory are made. It is also noted that the size of the stress-based HRR-dominance zone is only a fraction of the creep zone size in plane stress, and is very small (especially along =0°) compared to the creep zone size in plane strain. Furthermore, the dominance of the singular strain fields are at least two orders of magnitude smaller than the corresponding stress dominance zones. As such, unless the microstructural features of the material are smaller than the dimensions of the dominance zones, the basis for using stress or strain-based fracture parameters derived from the HRR-type fields for prediction of creep fracture initiation is unclear.     

基于COMSOL连铸坯壳测厚横波EMAT优化

为了研究电磁超声传感器(Electromagnetic Acoustic Transducer,EMAT)横波测量连铸坯壳厚度的机理及横波在连铸坯中的传播情况,选取坯壳厚度为10~50 mm的Q235小方坯为研究对象,利用有限元软件COMSOL建立脉冲电磁铁和螺旋线圈的电磁超声模型,分析在不同EMAT结构参数下,连铸坯中电磁场、力场、声场的分布规律。研究结果表明:脉冲电磁铁和螺旋线圈组成的EMAT能够在连铸坯壳集肤层激发出超声波横波。脉冲电磁铁空心螺线管线圈匝数、内半径、线圈导线半径对换能效率的影响依次减小,且当脉冲电磁铁内半径尺寸大于螺旋线圈尺寸时,产生横波的效率最高。坯壳厚度越小,螺旋线圈最优激励频率越大,测量精度越高,信号衰减越快。因此,坯壳厚度为10~50 mm的Q235小方坯选择1.1 MHz为最佳激励频率。     

有限元简化模型应用于超声换能器模拟分析

依据有限元方法的基本物理思想,在某些不需要计算辐射声场的准确声学参数和波束特性的工程应用方面,对流体模型进行充分简化,提出了简化模型处理的有效方法,利用该方法对超声换能器进行模拟分析,并进行了样品的制作和测试,实测结果与模型简化分析处理的结果基本一致。可以证明,用该方法进行换能器的优化设计是可行和高效的。     

A transition element for uniform strain hexahedral and tetrahedral finite elements

A transition element is presented for meshes containing uniform strain hexahedral and tetrahedral finite elements. It is shown that the volume of the standard uniform strain hexahedron is identical to that of a polyhedron with 14 vertices and 24 triangular faces. Based on this equivalence, a transition element is developed as a simple modification of the uniform strain hexahedron. The transition element makes use of a general method for hourglass control and satisfies first‐order patch tests. Example problems in linear elasticity are included to demonstrate the application of the element. Copyright © 1999 John Wiley & Sons, Ltd. This paper was produced under the auspices of the U.S. Government and it is therefore not subject to copyright in the U.S.     

A hierarchic 3D finite element for laminated composites

In this paper a new 3D multilayer element is presented for analysis of thick‐walled laminated composites. This element uses two steps to calculate the full stress tensor. In the first step the in‐plane stresses are computed from the material law using a displacement approximation, and then the transverse stresses are calculated from the 3D equilibrium equations. Since the 3D equilibrium equations require high‐order interpolation functions, a hierarchic interpolation of displacements is used. The new element is compared with existing ones, e.g. from MSC.MARC. Copyright © 2004 John Wiley & Sons, Ltd.     

An integrated,finite element‐based process model for the analysis of flow,heat transfer,and solidification in a continuous slab caster

An integrated, finite element‐based process model is presented for the prediction of full three‐dimensional flow, heat transfer, and solidification occurring in a continuous caster. Described in detail are the basic models for the analysis of turbulent flow and heat transfer in the liquid steel zone, in the zone of mixture of the liquid steel and solidified steel, and in the solidified zone. Then, the models are integrated to form a process model which can take into account the strong interdependence between the heat transfer behaviour and the flow behaviour. The capability of the process model to reveal the detailed aspects of turbulent flow, heat transfer, and solidification occurring in a continuous caster is demonstrated through a series of process simulations. Copyright © 2003 John Wiley & Sons, Ltd.     

铝合金轮毂弯曲疲劳试验的有限元模型

在铝合金轮毂产品开发阶段,针对试压铸小批量样件来进行的弯曲疲劳试验,是必须通过的台架试验之一．为避免试验的盲目性,减少试验次数,从而降低试验成本,提高试验的可靠性,有限元数值分析技术被采用,并已成为铝合金轮毂开发的先进设计技术手段．为正确应用该技术进行虚拟台架试验,必须建立有效的计算模型．基于理论分析与应用实践,考虑螺栓预紧力影响的非线性接触,以及加载杆与轮毂的材料异同,建立了轮毂弯曲疲劳台架试验的3种有限元力学模型：加载杆与轮毂为同种材料的整体线弹性分析模型、加载杆与轮毂为不同材料的线弹性分析模型、考虑各元件间的接触关系及螺栓预紧力影响的非线性分析模型．通过分析与比较,以及典型轮毂的试验验证,研究了以上3种有限元力学模型的优缺点、应用范围及计算的可靠程度,并以典型轮毂为例,用试验验证了模型的准确性．研究表明：加载非线性模型较为准确;线弹性模型可以用于轮辐处易发生疲劳破坏的分析;加载杆与轮毂材料异同对最大等效应力影响不大;对螺栓孔附近易破坏的轮毂,应采用考虑预紧力影响及非线性接触的模型．     

Statically admissible finite element solution of the thin plate bending problem with a posteriori error estimation

Two triangular elements of class C⁰ developed on the basis of the principle of complementary work are applied in the static analysis of a thin plate. Some techniques to widen the versatility of the equilibrium approach for the finite element method are presented. Plates of various shapes subjected to diverse types of loading are considered. The results are compared with outcomes obtained by use of the displacement-based finite element method. By use of these two dual types of solutions, the error of the approximate solution is calculated. The lower and upper bounds for the strain energy are found.     

A generalized phase field multiscale finite element method for brittle fracture

A generalized multiscale finite element method is introduced to address the computationally taxing problem of elastic fracture across scales. Crack propagation is accounted for at the microscale utilizing phase field theory. Both the displacement-based equilibrium equations and phase field state equations at the microscale are mapped on a coarser scale. The latter is defined by a set of multinode coarse elements, where solution of the governing equations is performed. Mapping is achieved by employing a set of numerically derived multiscale shape functions. A set of representative benchmark tests is used to verify the proposed procedure and assess its performance in terms of accuracy and efficiency compared with the standard phase field finite element implementation.     

A finite element model for liquid phase electroepitaxy

A finite element numerical simulation model for the liquid phase electroepitaxial growth process of gallium arsenide is presented. The basic equations obtained from the fundamental principles of electrodynamics of continua, the constitutive equations for the liquid and solid phases derived from a rational thermodynamic theory, and the associated interface and boundary conditions are presented for a two-dimensional axisymmetric growth cell configuration. The field equations are solved numerically by an adaptive finite element procedure. The effect of moving interfaces is taken into account. Numerical simulations are carried out for different convection levels by changing the value of the gravitational constant. Results show that convection has significant effect on the growth process under normal gravity conditions and results in thickness non-uniformity of the grown layers. The thickness non-uniformity leads to curved interfaces of growth and dissolution, which enhance convection.     

Mixed-hybrid scheme of the finite element method for solving the problems on bending,free vibration,and stability of plates

To solve a problem on bending, vibration, and stability of plates, a hybrid finite element has been constructed on the basis of Zienkiewicz’s triangle. A mixed approximation is used for the plate deflection and turning angles. It is shown that with a decrease in the triangle dimensions the mixed approach ensures convergence both for the plate deflection and the bending moments, which is practically independent of the way the plate is split into triangular elements. In the problems on free vibrations and stability of plates, the mixed approach yields more exact values of the eigenfrequencies and critical loads as compared to a classical Zienkiewicz’s triangle. The results of the numerical analysis of the convergence and accuracy of the solutions to a number of test problems on bending, free vibration, and stability of a square plate are presented. __________ Translated from Problemy Prochnosti, No. 4, pp. 108–122, July–August, 2008.