计算结构力学-有限元.pptVIP

计算结构力学;第2章有限单元法;2.1概述;2.1概述;2.1.2有限单元法概念;平面问题的常用单元：;2)单元分析

假设位移模式，分析单元力学特性：

{F}e=[k]{?}e

体力、面力——?等效结点荷载

3)整体分析——建立{F}=[K]{?}

[K]{?}={R}

4)再次单元分析

求出各单元的应变和应力。;2.2.1两类平面问题

1)平面应力问题

?z=0,?zx=0,?zy=0

?x=?x(x,y)

?y=?y(x,y)

?xy=?xy(x,y)

例如：深梁、剪力墙(受竖向力);2)平面应变问题

w=0,?z=0,?zx=0,?zy=0

?x=?x(x,y)

?y=?y(x,y)

?xy=?xy(x,y)

?x=?x(x,y)

?y=?y(x,y)

?xy=?xy(x,y)

例如：挡土墙、重力坝;1)根本量;几何方程：

物理方程：

(平面应力问题)

[D]:弹性矩阵;

E:弹性模量;

?:泊松比。;;平衡方程弱形式——能量原理;在集中力{F}作用下,;单元局部编号：i,j,m

单元结点位移向量:{?}e=[uiviujvjumvm]T

单元结点力向量:{F}e=[UiViUjVjUmVm]T

;2.3.2单元位移模式;将位移模式写成结点位移的显式：

u=Niui+Njuj+Nmum

v=Nivi+Njvj+Nmvm

Ni、Nj、Nm：形函数(插值函数);ai、bi、ci：分母行列式第1行元素代数余子式：

ai=xjym-xmyj

bi=yj-ym〔i,j,m轮换〕

ci=-xj+xm;形函数的性质:

(i)(Ni)i=1，(Ni)j=0，(Ni)m=0

(ii)单元内任一点：Ni+Nj+Nm=1;2.3.3单元应变矩阵、应力矩阵;{?}=[B]{?}e;2.3.4有限元方程的建立;令：{?}e=[G]{?}，[G]——6?2n位置矩阵;令：;2.3.5单元刚度矩阵;[k]元素的物理意义

kpq：第q个结点位移分量为单位位移〔其它结点位移=0〕，所引起的第p个结点力分量。

如:k26。

[k]的性质：

(1)对称性：kpq=kqp

(2)奇异性；

(3)每行〔列〕元素之和为零。

(4)[k]取决于单元的形状、方位和弹性常数，与所在位置〔即平移或np转动〕无关。;3)单刚列式推导的另一方法;虚功原理;2.3.6单元等效结点荷载;2.3.7整体刚度矩阵;总刚[K]的集成

例1:E,t,n=常数

1)离散

4个单元，6个结点

单元编号：(1)~(4)；

结点编号：1~6

2)建立局部-整体编码关系;;4)换码，对号入座，形成总刚;对一总刚元素有奉献的单刚：;例2:E,t=常数，?=0；l12=l13=l35=1m，形成[K]。;(1)

(3);(1)

(3);整体刚度元素Kpq的物理意义：

结构第q个结点位移为单位位移〔其它结点位移=0〕时，所引起的第p个结点力。

[K]的性质：

(1)对称性(Kpq=Kqp),主对角元素必为正；

(2)稀疏性，且一般为带状分布；

平面问题最大半带宽=

2?(单元结点号之差最大值+1)

(3)引入约束条件后为正定矩阵。;利用对称性、带状稀疏性，

[K]可用半带宽存储。;2.3.8位移边界条件的引入;(1)

(3);(1)

(3);2.对角元素改1法〔零位移边界〕

?i=0Kii=1,Kij=0,Kji=0〔j≠i〕,Ri=0

;3.乘大数法

?i=?iKii=?(大数)，Ri=?Kii?i;2.3.9位移模式与解答的收敛准那么;三结点三角形单元的