|
前言
1 N% h: |8 n- Y0 E2 c6 |# E学习有限元分析有两年时间了,非常热爱这个方向,借助此版块记录下自己一些学习体会小结,同时希望能与论坛内热爱有限元的朋友共同交流、相互促进,使自己对有限元分析有更深刻的认识。- ~" b8 W# P6 ?: P1 ]5 p/ _
* C* {6 V: m, {& h一、什么是有限元法(FEM)?
; M' I( S" ?$ M$ H( k有限元法是建立在弹性力学的基础上发展而来的一门学科,是用来求解微分方程组近似解的一种方法。注意,是近似解。3 t; U' r% e; B) T$ w) c7 j
在工程中很多微分方程组只靠纯粹的数值求解是无法得到答案的,而有限元法的出现解决了这一问题,特别是随着计算机的快速发展有限元法得到快速的推广。- g, \6 J5 k) q9 ?2 ^. ]
. r9 B3 v3 f G: ~- G
二、有限元的基本思想是什么?
0 o E9 S6 y4 c( Y% Z2 J8 h+ p其思想为:化整为零,积零为整。
4 K8 E8 N. n7 O6 r1 l& t对这个八个字的解释是:一个连续的物体具有无限的自由度,通过网格划分把连续的物体划分成有限个单元,每个单元通过边缘节点连接到一起,组装成一个整体,这样就把无限自由度问题转化为有限自由度问题。每一个单元内都有一个假设物理场(例如位移场),利用边缘节点数值相同这一条件及变分原理求得我们所需要的物理量。
2 t2 ~$ R: ~) |6 o4 I+ c基本物理量有三个:位移(displacement)、应力(stress)、应变(strain)。/ a* w" B$ F$ c$ s& U2 _: l
位移(displacement):u、v、w分别表示在笛卡尔坐标系下X、Y、Z三个方向的位移量。+ X/ \7 T; D u, ?6 @' K
应力(stress):应力的物理含义为为单位面积上所受的内力。在结构中任意一点的应力状态沿着不同界面都不一样。但学习弹性力学可知,从微观角度取一个微元段进行研究,只需要六个应力即可完全确定一点的应力状态,分别为 σx 、 σy 、 σz 、 τ xy 、 τxz 、τyz。注意在弹性力学中τ xy=τ yx
" ^. y2 J) F! N/ U& [/ [8 D这是根据切应力互等定理得到(与材料力学中的切应力互等定理不同)。在有限元中由这个六个基本应力(考虑到τ xy=τ yx)构成了二维张量。
7 P; a* X! K5 Z' e6 C# ~应变(strain):对应着应力应变也有六个。# l# N& s; c7 I5 [% \
4 V# v- P2 e& u- p三、有限元的求解步骤是什么?
* l; t% L5 B# D5 _9 B1.单元划分及单元节点编号。
- F' R! H- I/ Y$ W- j2.求解单元刚度矩阵,每一个单元都一个单元刚度矩阵。单元类型的不同也就体现在单元刚度矩阵不同。它反映了自身的单元特性。例如,壳单元中的一阶三角形单元是一个常应力单元,即没有应力梯度;而四边形单元就好多了,能够反映结构的应力梯度。+ I9 }5 i+ i( g, M4 g+ Q- I1 E
3.组装成整体单元刚度矩阵。一般为稀疏矩阵或称为带状矩阵。
, c# C' A( w8 Y4.边界条件处理,包括约束与载荷。5 p7 s6 ^3 D0 l7 F0 F- A$ E e; o
5.求解运算。" {& W0 p7 I" i: l- L! X; @" d
6.后处理。# V6 m( D: ~. {3 _% X7 \
以上为理论分析步骤,对应着软件分析步骤:1、2、3、4为软件操作中的前处理;5.对应着计算机的求解,此过程相当于一个黑匣子;6对应后处理。4 Z: g; l, N$ ~; j5 b; [
0 O1 {1 o c7 H. F6 W8 J
前处理常用软件:Hypermesh、ansa我擅长使用Hypermesh,非常方便。6 R3 _' @- G# ~% P% @0 {. m$ h+ Q
求解器:Nastran结构分析中的行业标准、abaqus非线性分析老大、ansys多场耦合。当然这些软件也有自己的前后处理器。
9 F) z9 E7 v( P, h |
|