Introduction to Mesh Based Generated Lumped Parameter Models for Electromagnetic Problems using Triangular Elements

This paper is an introduction to mesh based generated reluctance network modeling using triangular elements.Many contributions on mesh based generated reluctance networks using rectangular shaped elements have been published,but very few on those generated from a mesh using triangular elements.The use of triangular elements is aimed at extending the application of the approach to any shape of modeled devices.Basic concepts of the approach are presented in the case of electromagnetic devices.The procedure for coding the approach in the case of a flat linear permanent magnet machine is presented.Codes developed under MATLAB environment are also included.

Development and Application of a High-Performance Triangular Shell Element and an Explicit Algorithm in OpenSees for Strongly Nonlinear Analysis

The open-source finite element software,OpenSees,is widely used in the earthquake engineering community.However,the shell elements and explicit algorithm in OpenSees still require further improvements.Therefore,in this work,a triangular shell element,NLDKGT,and an explicit algorithm are proposed and implemented in OpenSees.Specifically,based on the generalized conforming theory and the updated Lagrangian formulation,the proposed NLDKGT element is suitable for problems with complicated boundary conditions and strong nonlinearity.The accuracy and reliability of the NLDKGT element are validated through typical cases.Furthermore,by adopting the leapfrog integration method,an explicit algorithm in OpenSees and a modal damping model are developed.Finally,the stability and efficiency of the proposed shell element and explicit algorithm are validated through the nonlinear time-history analysis of a highrise building.

An explicit finite volume element method for solving characteristic level set equation on triangular grids

Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.

A 3-Node Co-Rotational Triangular Finite Element for Non-Smooth,Folded and Multi-Shell Laminated Composite Structures

Based on the first-order shear deformation theory,a 3-node co-rotational triangular finite element formulation is developed for large deformation modeling of non-smooth,folded and multi-shell laminated composite structures.The two smaller components of the mid-surface normal vector of shell at a node are defined as nodal rotational variables in the co-rotational local coordinate system.In the global coordinate system,two smaller components of one vector,together with the smallest or second smallest component of another vector,of an orthogonal triad at a node on a non-smooth intersection of plates and/or shells are defined as rotational variables,whereas the two smaller components of the mid-surface normal vector at a node on the smooth part of the plate or shell(away from non-smooth intersections)are defined as rotational variables.All these vectorial rotational variables can be updated in an additive manner during an incremental solution procedure,and thus improve the computational efficiency in the nonlinear solution of these composite shell structures.Due to the commutativity of all nodal variables in calculating of the second derivatives of the local nodal variables with respect to global nodal variables,and the second derivatives of the strain energy functional with respect to local nodal variables,symmetric tangent stiffness matrices in local and global coordinate systems are obtained.To overcome shear locking,the assumed transverse shear strains obtained from the line-integration approach are employed.The reliability and computational accuracy of the present 3-node triangular shell finite element are verified through modeling two patch tests,several smooth and non-smooth laminated composite shells undergoing large displacements and large rotations.

Finite Element Analysis of Magnetohydrodynamic Mixed Convection in a Lid-Driven Trapezoidal Enclosure Having Heated Triangular Block

A numerical research on magnetohydrodynamic mixed convection flow in a lid-driven trapezoidal enclosure at non-uniform heating of bottom wall has been studied numerically. The enclosure consists of insulated top wall and cold side walls, too. It also contains a heated triangular block (Rot = 0° - 90°) located somewhere inside the enclosure. The boundary top wall of the enclosure is moving through uniform speed U0. The geometry of the model has been represented mathematically by coupled governing equations in accordance with proper boundary conditions and then a two-dimensional Galerkin finite element based numerical approach has been adopted to solve this paper. The numerical computations have been carried out for the wide range of parameters Prandtl number (0.5 ≤ Pr ≤ 2), Reynolds number (60 ≤ Re ≤ 120), Rayleigh number (Ra = 103) and Hartmann number (Ha = 20) taking with different rotations of heated triangular block. The results have been shown in the form of streamlines, temperature patterns or isotherms, average Nusselt number and average bulk temperature of the fluid in the enclosure at non-uniform heating of bottom wall. It is also indicated that both the streamlines, isotherm patterns strongly depend on the aforesaid governing parameters and location of the triangular block but the thermal conductivity of the triangular block has a noteworthy role on the isotherm pattern lines. Moreover, the variation of Nuav of hot bottom wall and θav in the enclosure is demonstrated here to show the characteristics of heat transfer in the enclosure.

TRIANGULAR ELEMENTS FOR REISSNER-MINDLIN PLATE

A general method to construct locking free Reissner-Mindlin plate elements is presented. According to this method the shear strain is replaced by its proper interpolation polynomial, which corresponds to the Kirchoff conditions at the interpolation points as the thickness of plate tends to zero, so the element is locking free. We construct two triangular elements by this method - a 3-node element and a 6-node element. The numerical results are provided.

The Numerical Integration of Discrete Functions on a Triangular Element

With the application of Hammer integral formulas of a continuous function on a triangular element, the numerical integral formulas of some discrete functions on the element are derived by means of decomposition and recombination of base functions. Hammer integral formulas are the special examples of those of the paper.

The Triangular Hermite Finite Element Complementing the Bogner-Fox-Schmit Rectangle

The Bogner-Fox-Schmit rectangular element is one of the simplest elements that provide continuous differentiability of an approximate solution in the framework of the finite element method. However, it can be applied only on a simple domain composed of rectangles or parallelograms whose sides are parallel to two different straight lines. We propose a new triangular Hermite element with 13 degrees of freedom. It is used in combination with the Bogner-Fox-Schmit element near the boundary of an arbitrary polygonal domain and provides continuous differentiability of an approximate solution in the whole domain up to the boundary.

18-DOF Triangular Quasi-Conforming Element for Couple Stress Theory

The basic idea of quasi-conforming method is that the strain-dis-placement equations are weakened as well as the equilibrium equations.In this paper,an 18-DOF triangular element for couple stress theory is proposed within the framework of quasi-conforming technique.The formulation starts from truncated Taylor expansion of strains and appropriate interpolation functions are chosen to calculate strain integration.This element satisfies C0 continuity with second order accuracy and weak C1 continuity simultaneously.Numerical examples demonstrate that the proposed model can pass the C0??1 patch test and has high accuracy.The element does not exhibit extra zero energy modes and can capture the scale effects of microstructure.

Analysis of Linear Triangular Elements for Convection-diffusion Problems by Streamline Diffusion Finite Element Methods

这份报纸被奉献给学习 superconvergence 优化散开为传送对流散开问题的有限元素方法。在里面[8 ] ，在条件下面那 h2 最佳的有限元素错误估计在 L2 标准被获得。然而，在现在的纸，一样的错误估计结果在更弱的条件下面被获得那 h。

Finite layer and triangular prism element method to subsidence prediction and stress analysis in underground mining

The application of the finite layer & triangular prism element method to the 3D ground subsidence and stress analysis caused by mining is presented. The layer elements and the triangular prism elements have been alternatively used in the numerical simulation system, the displacement pattern, strain matrix, elastic matrix, stiffness matrix, load matrix and the stress matrix of the layer element and triangular prism element have been presented. By means of the Fortran90 programming language, a numerical simulation system based on finite layer & triangular prism element have been built up, and this system is suitable for subsidence prediction and stress analysis of all mining condition and mining methods. Comparing with the infinite element method, this approach dramatically reduces the size of the set of equations that need to be solved, and greatly reduces the amount of data preparation required. It not only saves the internal storage, and the computation time, but also decreases the cost.

Finite Element Analysis of Magnetohydrodynamic Natural Convection within Semi-Circular Top Enclosure with Triangular Obstacles

The phenomena of magneto-hydrodynamic natural convection in a two-dimensional semicircular top enclosure with triangular obstacle in the rectangular cavity were studied numerically. The governing differential equations are solved by using the most important method which is finite element method (weighted-residual method). The top wall is placed at cold Tc and bottom wall is heated Th. Here the sidewalls of the cavity assumed adiabatic. Also all the wall are occupied to be no-slip condition. A heated triangular obstacle is located at the center of the cavity. The study accomplished for Prandtl number Pr = 0.71;the Rayleigh number Ra = 103, 105, 5 × 105, 106 and for Hartmann number Ha = 0, 20, 50, 100. The results represent the streamlines, isotherms, velocity and temperature fields as well as local Nusselt number.

Robust Falk-Neilan Finite Elements for the Reissner-Mindlin Plate

The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conforming or consistent in the sense that no projection/reduction is introduced.Theoretical proof and numerical confirmation are presented.

Four-Order Superconvergent Weak Galerkin Methods for the Biharmonic Equation on Triangular Meshes

A stabilizer-free weak Galerkin(SFWG)finite element method was introduced and analyzed in Ye and Zhang(SIAM J.Numer.Anal.58:2572–2588,2020)for the biharmonic equation,which has an ultra simple finite element formulation.This work is a continuation of our investigation of the SFWG method for the biharmonic equation.The new SFWG method is highly accurate with a convergence rate of four orders higher than the optimal order of convergence in both the energy norm and the L^(2)norm on triangular grids.This new method also keeps the formulation that is symmetric,positive definite,and stabilizer-free.Four-order superconvergence error estimates are proved for the corresponding SFWG finite element solutions in a discrete H^(2)norm.Superconvergence of four orders in the L^(2)norm is also derived for k≥3,where k is the degree of the approximation polynomial.The postprocessing is proved to lift a P_(k)SFWG solution to a P_(k+4)solution elementwise which converges at the optimal order.Numerical examples are tested to verify the theor ies.

仿生平面三角形单元积分声场的波函数构造

波叠加法在求解声源外部辐射声场时,需要对所有离散单元进行数值积分计算导致计算效率较低,而等效源法由于过度简化单元始终存在较大的积分近似误差.针对以上缺陷,利用Helmholtz方程在球坐标系下的解构造了一种与单元积分等效且无需积分的波函数.受仿生复合材料三角形缝合结构启发,以适用范围最广的平面三角形单元为例,构造了波函数的一般形式和内推形式.最后,通过数值仿真对比了两种波函数与直接积分的计算声场的精度和效率.结果表明,两种波函数与直接积分的计算误差低于0.5%,且内推波函数的计算效率约为直接积分的6倍.

接触问题的三角形载荷离散FFT加速算法

接触问题控制方程的有效求解,往往涉及到复杂的数学理论知识,而在实际工程应用中,接触应力分布又具有高度随机性。为高效快速求解任意载荷分布下固体的接触响应,基于三角形载荷离散单元,嵌入离散卷积快速傅里叶变换(DC-FFT)算法,提供了一种高精度、高可靠度的计算方法。相比于通常采用的分段均布载荷离散方法,三角形单元的解析求解略显复杂,但能更好地模拟接触载荷任意分布的特性,对于接触边缘处载荷由零递增或递减为零的情况,也可以予以充分考虑。为优化三角形载荷离散单元的求解方法,根据接触影响系数矩阵的“激励-响应”特性,推导了三角形载荷单元和均布载荷单元作用下的应力分量解析解。通过构造包含影响系数矩阵的离散卷积形式应力解,将某一目标节点在所有载荷单元作用下,重复度极高的矩阵运算叠加过程,采用DC-FFT算法来简化加速计算。通过程序编程计算,分析验证了所提出算法的精确度和高效性。

静电成形薄膜反射面天线的裁剪设计

该文通过裁剪和拼接平面薄膜片来对静电成形薄膜反射面天线的表面进行设计制作。由于抛物面型薄膜反射面是一个不可展开表面,通常由许多平面薄膜片制作而成,将不可避免地降低反射面表面的精度。因此,需要详细地分析裁剪和拼接过程对反射面表面精度的影响,从而获得设计需要的表面精度。该文基于旋转抛物面的无矩理论建立了工作状态的反射面模型,采用常应变三角形单元离散反射面,建立了反射面的展平和拼接有限元模型。将同口径反射面沿经线划分若干组不同数量子曲面,通过子曲面进行展平和拼接数值仿真分析,对薄膜反射面进行裁剪设计。最后,制作了0.55 m口径的薄膜反射器,并进行了摄影测量实验。验证了该方法的有效性,并讨论了制造误差的影响。

基于板-梁理论的双轴对称三角形管翼缘工字形梁组合扭转的理论研究及有限元验证

基于板-梁理论推导双轴对称三角形管翼缘工字形悬臂梁的总应变能方程,得到梁的翘曲刚度和自由扭转刚度,然后通过外力势能和总应变能得到总势能方程。再依据能量变分模型、微分方程模型及边界条件得出悬臂梁端部最大扭转角公式,将求出的最大扭转角与用ANSYS建立的有限元模型计算结果进行对比。结果表明,悬臂梁在端扭矩作用下的扭转角理论值与有限元值较为接近,误差最小的绝对值为2.40%,最大的绝对值为4.44%,可为工程设计提供参考。

混凝土脱空修复对角钢连接件力学性能影响

为探索钢壳混凝土角钢连接件根部混凝土脱空对连接件抗剪性能的不利影响,提出合理的脱空修复方法,基于深中通道沉管隧道角钢连接件的局部构造,以脱空形状、脱空修复和混凝土强度为变化参数,设计8组足尺推出试件,并开展抗剪性能试验和有限元模拟分析。分析显示:1)普通混凝土试件先在角钢内部发生混凝土压溃破坏,试件破坏模式受脱空形式及脱空修复影响不明显;2)脱空对角钢连接件的抗剪承载力和抗剪刚度均有削弱影响,三棱柱整体脱空时角钢连接件的抗剪承载力和抗剪刚度可折减32.8%和76.9%;3)脱空修复对角钢连接件的抗剪承载力和抗剪刚度有提升作用,但由于修复料与脱空界面间难以协同受力,性能提升效果并不显著;4)考虑脱空处未修复,用超高性能混凝土替换普通混凝土时,推出试件的抗剪承载力较未脱空时提高77%,抗剪刚度较脱空时折减8.9%;5)超高性能混凝土的优良力学性能使角钢连接件的抗剪承载力得到明显提高,但由于脱空导致角钢连接件出现较大变形,因此无法有效提升脱空时的抗剪刚度;6)提高修复料力学性能对改善三棱柱脱空修复后的承载性能影响不明显。

采棉机三角履带行走装置设计与仿真

为了适应长江、黄河流域棉区的棉花机械化采收需求,设计了一款采棉机三角履带行走装置。分析棉田工况、整车接地比压和车速等参数,确定了三角履带行走装置的主要技术参数;通过有限元分析,对基架、驱动轮和摆动销轴等关键部件进行了强度校核计算。结果表明:所设计结构与材料选用均能满足使用要求,单轮承载最高可达8 t。该装置可与原车轮胎快速互换,对于推进长江、黄河流域棉花收获机械化具有一定意义。