有限元课件西安交通大学文档.pptVIP

研究生《有限元法》授课大纲

绪论、弹性力学基本方程、虚位移原理、最小势能原理

2杆系结构有限元

3平面问题有限元分析

4空间与轴对称问题

5板弯曲有限元分析

6壳弯曲有限元分析

7广义变分原理

理论一公式一程序一实际应用

《有限元法》考核方式

1有限元学习报告(打字)

2一个有限元程序(MATLAB语言编写)

3口试

时间:第19周前全部完成

第20周考试结束

第1章有限元法绪论

第1节概述Clough-Thefiniteelementmethod

起源:50年代飞机结构矩阵分析Argyris,Turner,Clough

60年代弹性力学平面问题,目前已涉及众多领域

实质:对力学模型进行近似数值计算的方法

将无限自由度问题变成有限自由度问题

分析过程:结构离散化,确定位移模式,单元特性分析

整体分析,解方程输出计算结果,其他处理

学习方法:与矩阵位移法对比一相同与不同之处

杆系结构

了解基本原理,各种方法的共性与实质

通过自编程序进一步熟悉原理

应用状况:标准通用软件SAP200,ANSYS,

各种专用程序

连续体

第2节弹性力学基本方程

、平衡方程

+X=0

二、几何方程

ayaz.

t

ayaz

avaw

+8+2+Z=0

azay

axaaz

三、本构关系

6=10-(,+0,=21+B)zn

五、边界条件(应力,位移)

6,=Et,0+o)E

应力

62=0:-A(G,+0,)yx=

=gl+tmttn

E

①,=r1+a,m+n

四、协调方程

④:=t2+。m+0

Oyaxandyax(ayazax

ayaz

位移

a

az-ay-ava

azaray)azar

ax2axazaraz(axayaz)axay

续第2节弹性力学基本方程矩阵表示

位移列阵d=

对

体积力列阵E=[Xyz]

D

应力列阵g=a,anzr]

00

应变列阵

00002(1+4)

000

表面力列阵9-[.中,

已知位移列阵d

、平衡方程Aa+F=0

微分算子列阵

、几何方程

A=0

0

本构关系a=DE

08

四、协调方程C≌=0

00

C协调算子矩阵(略

表面外法线方

00m0

五、应力边界条件Φ-Lg

向余弦矩阵

L-0m0In0

位移边界条件d-d=0

第3节虚位移原理为

虛位移原理

弹性体处于平衡状态的必要与充分条件:对于任意的、满足相容条

件的虚位移,外力所做的功等于弹性体所接受的总虚变形功。

6W=6U

外力虚功:

(XSu+YSv+ZSw)dxdydz+l(XSu+YS+ZSw)dS

忌虚变形功:

(od+o,d+o,,+tu,syv+r,+trdxdydz

对于平面问题:

fOXau+Yo)dxdy+(Xou+YS)ds=(o,06,+o,,+t,fray)dxdy

第4节最小势能原理

最小势能原理

在几何可能的一切容许位移和形变中,真正的位移和形变使总势能取

最小值;反之,使总势能取最小值者也必是真正的位移和形变。

总势能:∏I=U+V

形变势能:U=(c+,+E2+7m7+m7n+n)d

外力势能:V=-J(x+1+z)dh-(xm+Y+z)d

形变势能变分

前(c+,B,+:+1例mx++n)hn

即:形变势能的变分表达式与虚变形功的表达式完全相同

外力势能变分

XSu+YSu+Zow)dxdydz-[(XSu+YSv+Zow)dS

即:外力势能的变分表达式与外力虚功负值的表达式完全相同。

第2章杆系结构有限元

第1节等直杆单元分析

1、用结点位移表示单元的位移模式

位移列阵另-a=[an,

设位移模式=a+4v=b+bx+bx2+bx

h(x)=[x]H(x)=[xxx

由结点位移得=4]}-2]

其中:

111213

00

012132

待定参数为

l}=4]=团4]24

结点位移表示的位移模式为

2

=xlrv=lHxl4每形函数矩阵为:n=N]v=四v1v[N.=[n(xI4上N]=xl4

}=lwl}=[,,o=可

第1节等直杆单元分析续f

H(x)

AS]=No

100000

H

010000

001000

.(小-oox0对-b1x0x-700700

2

2、用结点位移表示应变和应力

拉压应变,E弯曲应变

ej

yH”(x)

Ix

00026

}={}=E}=Er

第1节等直杆单元分析续2

3、用虚位移原理导出梁单元的刚度矩阵

单元内各点的虚位移:{}-{NVs}{s}为结点虚位移

单元内的虚应变:{}=[B

单元内的虚功:况}kv=E)1vy

单元结点力:{F}=NgM,NQM

单元外力虚功:a=如k4()=(3)(rkx+

虚位移原理∞=oW

令:状RF=Nr