关于几何非线性有限元计算
1、工程应用场合 S; B! [7 a P& e
几何非线性主要应用计算如下情况:5 Y# P7 l+ q5 L- o
-单体部件(梁和壳体)发生较大的变形;
* I6 `; j( g I( R7 ?-非线性屈曲计算;/ I+ s5 Y/ \* C4 Q0 P
-塑性成形计算9 I! h, `# G) V* b* w; G
-瞬态动力学的计算
4 I" D9 @0 W6 C% C-含有橡胶部件的计算等
& B* @: |* c, W& ~ v
! z3 q/ h0 `' C O8 q" e
( K3 E. F( P3 ~1 D
6 l. x' V1 V) m4 P* C
! v% b8 }4 i5 S6 f/ t5 C这是一个大应变分析的例子,一个轴对称的橡胶密封件受压缩。分析包括接触,当密封件折叠时会发生自身接触。
) X0 Y, U# y3 P7 {( j9 Q
/ [0 h6 I* \- X4 L+ Y8 S* G4 U7 ^" L" w
3 G- Y4 d/ \) j) P3 L
) c5 ^5 t1 @ q& X
此例显示了绕轴线捆扎一根钢条。将金属弯曲成不同的形状是生产中常见的操作。在此例中,应变达到25%,顶端旋转接近270 度
% ~+ [5 X7 m2 W/ q' G' b
' E0 B9 [+ `" A1 z# [, T7 B2、计算原理
* [; A$ Z A- Z+ v" N! a- }几何非线性的本质是把结构的变形对结构的刚度考虑到总体刚度矩阵中。0 O3 k6 A. B; G! f* b! {) M) U
0 B, J4 ?$ T( z
% I$ _3 u; M) O2 M+ \5 _1 j% _
! h; T, ~% Z$ k1 [; k# C
对于几何非线性问题包括以下问题:
8 Q+ @& q9 b" l9 s% L-大变形,大应变;
, D; G p' Q) k- C9 o ] Y-大变形,小应变;
7 G+ y9 z9 Q' u-旋转软化;1 U$ g7 t7 V! C+ C! ^
-应力刚化等3 R0 B3 k1 T8 o; b- D3 W% s
目前在ANSYS Workbench平台中,只需要在分析系统中的分析设置中进行激活即可。
$ I: E7 H7 P; E9 J如果一个单元的形状发生改变(面积、厚度等),它本身的单元矩阵会发生改变。7 o3 K% z! w3 k
9 Z: d: u6 e+ u( B9 g
7 f0 ~$ W6 Y ] Z$ f& u1 H. k* i- Y: M9 K: E% q* r; G X. N" |5 }6 i' T
如果一个单元的取向改变,它的单元刚度向整体刚度的转换矩阵将发生变化。! D/ _- Q4 {! [- ^/ J! E
" z n$ C! s1 X+ K8 Q' t' r2 o
. r" ?/ x; \; V9 H2 D- @( N
3 J( ]% n% {8 Y9 H, O$ ^! R% k
5 k1 L: f. g/ n$ Q" Y如果单元应变产生了明显的面内应力 (膜应力),垂直于面的刚度会明显受影响。
, A8 y- m/ U0 }/ _- Q } 6 L( e. V" h" M h8 Y% k
/ g+ R4 }/ o7 ^3 f* L8 `2 h
. c0 Z+ s) S4 ?' X- [
2 C! o, \- V) b- n- z( P) k/ q/ R. ]2 w5 s2 X1 F5 W
/ x8 I7 k2 r: u
随着垂直位移的增加 (Y),大的膜应力(SX)导致刚化响应。; `& D) m. T0 ^. e
将非线性应变定义推广到一般的三维情况 # X- t2 g% Y: e1 S- K
在二维和三维中, 当一个元件经历大应变变形时, 不仅长度元素改变, 厚度、面积和体积也改变) M! W x2 p7 W8 h- k' d2 N" h/ _
0 w4 i1 M+ M6 `
4 F- ^# k8 L) W
4 h4 E- [: E* ?( ]) C运动和变形" W4 x( y& J0 J3 q8 J0 {
& F! a8 H p" R% Z/ g: d
物体在外载荷的作用下会移动和变形 z" ]3 g- r4 R7 i- y' g* I
: _; c7 [; b0 v- k' D$ Q+ c
v9 ^# z, o T, F8 Y2 Q: p
! k2 S N2 l% }1 W& w如果考察该物体上某一点的运动, 它的初始位置为 {X }, 最终位置为 {x }, 则位移量为 {u }, 图中! c' O/ D/ Q* n. b' t0 H
" J$ f' ^% Q n1 i$ P& q' i2 n
" u: h+ n- F/ q
: d) Q+ W7 W+ [; Z; H' ~4 |1 c变形梯度
8 s! b7 v/ m. p3 P变形梯度是物体变形程度的度量, 定义为, g- ^% ?* `1 n
. {$ q8 H" u+ z3 O( N) o* H" ]
9 M+ c" ?7 b" L% q' x6 ?$ F. {
" c8 e) A" e+ F! k3 y变形梯度 [F ] 包含如下信息:
% e" L, P3 _/ W/ r8 V, n体积改变4 J _4 D" N/ b. S; n7 D ]
转动
b4 z: `4 }) i" j+ v1 n9 |由于应变而引起的形状改变7 m- }/ a6 x0 v0 N D" s& D I; y5 U
• 大应变分析设定应变不再是无限小的,而是有限的或相当大的。
; s, N/ h# g& R• 当应变超过一定百分比及不能忽视几何形状的改变时,可认为% d0 I, b0 G1 I8 ^; k
是大应变。
3 K/ h- r* {; {• 大应变理论考虑了形状的改变(例如厚度,面积等等)及任何大旋转。
; e9 u! {9 ]: I" E- C: J• 当单元旋转“ 大”到明显影响求解精度时,需激活大挠度效应。
8 Y9 N" |' q% c. @5 [+ c• 遗憾的是,没有明确的规定“ 大”到底有多大。“ 小”与“ 大”之间根据问题的不同相差非常大。
+ j, \! B+ ]8 i& I |& H* [' J) P9 h. q/ w/ m* i
3、软件基本设置5 \) @. R( N6 }
' S" ]$ |* H! R, U! n; L- {
- [7 s H( y$ T! h4 \# F
0 [4 S: A# I5 W8 t* @
4 G$ }1 T& n; }9 a0 T- i
2 Q; x7 [3 R- f7 D+ Z
4 Y8 k- ]# z2 H) k& Q4 `0 R% z9 {# Y; l9 Q/ R
其他类型的大变形问题,根据实际情况,进行选择合理的材料模型。
& B z1 u9 u! d$ _! D6 H, p
# g# W% ^0 c- d* V! m# F: S; U* A
. F1 H' `( U3 R5 l& v
2 a' |# z' R2 d
* S2 F3 T% F, p* T) b
对于几何模型,没有特殊要求,用户可以自己在DM中建模也可以从其他CAD软件中导入。
, h* m% P$ x" U) _将定义的材料模型,分配给对于的几何模型。
# C5 l- a0 b+ t7 A* w z+ {1 D( j" W1 l5 c- \' C
+ o2 m4 h* R1 h3 X1 H9 _" Y+ G5 T5 y( d
对于几何非线性问题,在Mesh中的高级选项中,将形状检查选项设置为Aggressive Mechanical
! `4 O4 s- T" A8 \+ a1 g3 C' d在建立大应变、大位移模型时避免使用带中间节点的单元。8 `) d( I( a2 M- F5 d; Q6 R& y
带中间节点的(既高阶)单元在更新几何形状时,单元的中间节点可能穿过单元移动,导致产生负的旋转。3 j: {& z, K1 V9 H/ K+ r3 \
8 J- Q% [. l! A8 k/ j
9 L( A6 M7 h; b- n, H
; i* q: }* y$ W
- k+ N# H+ m: }& G
& d8 c# q+ y1 C4 z B: |5 l
4 K5 o; W9 U' ~4 o4 k) U. ~. ^# N
2 \/ v4 A7 b- w* L: _* i, i
( P3 r; r/ t' s$ P2 C在分析设置中,必须将Large Deflection选项设置为On,才能激活结构计算的大变形。
" @/ F5 Q, t' T- P( [3 ~) X这个选项在存在橡胶模型的问题,必须打开,否则不能计算。7 Z- i/ S' U% s6 l) S* R3 w" y! P% D
5 c) l# q" }9 B$ U8 A4 n/ j
* u! a8 o* @4 e7 F
* d. ?6 ~7 G" Q6 Q+ U+ T( a( u4 ^
当结构经历大位移与旋转时,载荷发生了什么变化。
. ~" h0 F" R- F5 _在许多情况下,载荷将在变形过程中保持一致的方向。 # R6 e% i% X: F# O" x8 j& C+ M
在另一些情况下,承受大旋转时,力将“ 跟随”单元改变方向。
. Y2 m \0 c4 \. i( a: @这两种情况ANSYS都可模拟,取决于施加载荷的类型。
! r+ z: O8 F9 N' ]2 H& p# Y加速度与集中力:5 A' q% J/ [7 i* b3 f
保持它们的初始方向,忽略单元取向。3 f0 T" _* p* |. A+ W
表面压力载荷:. [, E, B8 t' \, x
随单元旋转,因此总是垂直于变形单元的表面(这些是真正的追随力)。
' J& w: M4 \% Z1 D1 P% @% u7 L更新压力表面以计算大应变效应。因此,对于施加的常压力,总压力载荷将随表面积的改变而改变。1 ]* O* ]4 {8 x6 s7 G
- g9 h; N! y5 r
4 ^! H( i/ p; J& Y& ^
r6 A/ A( c5 V2 l4 A' N9 u$ {; |$ ?) `; ^$ Y
约束边界处的所有自有度,由于Poisson 效应将产出非常大的应变。! u1 j V/ c# |
; Z6 ?$ [, |, S+ Q
* \" {7 r/ | K
) |$ u# P$ M0 x" M: @
- X& }+ }# s: y) A4 Q6 ~0 U/ s% L: V在非线性计算中,子步(时间步)的设置非常关键,在静力学计算中时间没有实际意义,本质就是一个计数器和加载比例系数。 $ G7 f1 A7 W# W
) l) h" W! b3 T6 B1 m4 Y
$ Z1 m4 ]* p U/ r
2 ^0 r3 c, J. m; n; r, H8 u
# G- M) I x M2 C% U1 f. ~( m
8 r% a& F4 B5 J* G; l
: [) Q5 J3 B6 R9 R
' e, ~ ]8 z# c: I( V
! ^- d0 R; j: W1 z# @4 @5 g; L. P& J S0 x; Z9 g5 T
! [# \+ R8 J% H2 l* l3 g' S |
|