曲面薄壳有限元离散误差分析和修正方法研究

基于曲面薄壳几何/物理方程和有限元离散原理推导,揭示了壳单元离散方法特点、刚度矩阵误差产生机理和范数表示方法。以航空涡轮发动机测压探针为例,深入研究了利用协调平板壳单元、非协调平板壳单元、一次和二次实体退化壳单元等不同类型单元离散曲面薄壳时,应力计算误差影响因素及变化规律,综述了单元的适用条件。基于Boltzmann模型提出了应力计算结果修正通式,通过FORTRAN语言编程修正应力仿真结果和ANSYS二次开发语言APDL命令实现应力修正结果的读取和后处理,应力修正结果与理论解相对误差小于5%,满足曲面薄壳设计和稳定性分析应用要求。

Abstract：

Based on the geometric physical equations and deduction of discrete process of finite element for curved thin-walled structure,characteristics of shell elements for discretion and discrete error for element stiffness matrix were analyzed theoretically,and the norm expression of discrete error for element stiffness matrix was put forward.Then,as an illustration,a probe for detecting the pressure inside turboprops was meshed by curved thin-walled shell element,compatible flat shell element,incompatible flat shell element and 3D degenerated shell element,respectively,and the stress of the probe was numerically calculated by FEM.On the basis of solved results,the stress error was analyzed thoroughly,the influence factors on the error and the applying situations for different elements were determined.Finally,the general form of correctional function was put forward according to the Boltzmann model.The program developed with FORTRAN was used to correct the numerical results and the APDL commands inside ANSYS software were utilized to read and plot the correctional stress results.Through the correction,the discrete error can be reduced within 5% that satisfies the requirement of application for structural design and stability analysis.

Discrete Error Analysis and Correction Method of Finite Element for Curved Thin-walled Structures

Based on the geometric physical equations and deduction of discrete process of finite element for curved thin-walled structure,characteristics of shell elements for discretion and discrete error for element stiffness matrix were analyzed theoretically,and the norm expression of discrete error for element stiffness matrix was put forward.Then,as an illustration,a probe for detecting the pressure inside turboprops was meshed by curved thin-walled shell element,compatible flat shell element,incompatible flat shell element and 3D degenerated shell element,respectively,and the stress of the probe was numerically calculated by FEM.On the basis of solved results,the stress error was analyzed thoroughly,the influence factors on the error and the applying situations for different elements were determined.Finally,the general form of correctional function was put forward according to the Boltzmann model.The program developed with FORTRAN was used to correct the numerical results and the APDL commands inside ANSYS software were utilized to read and plot the correctional stress results.Through the correction,the discrete error can be reduced within 5% that satisfies the requirement of application for structural design and stability analysis.