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

计算结构力学 课 件（二） 第2章 有限单元法 2.1 概述 2.2 弹性力学平面问题的矩阵描述 2.3 三结点三角形单元分析 2.4 高精度三角形单元、矩形单元 2.5 C0连续型单元形函数的构造 2.6 平面等参数单元分析 2.7 有限元程序实现 2.8 平面杆件结构有限元 2.9 板弯矩的有限元 2.1 概述 2.1.1 发展概况 Courant——1943年应用三角形分片插值函数; Turner, Clough等——1956年推广直接刚度法; Clough——1960年提出“有限单元法”名称； Zienkiewicz等——编写第一本有限元方面专著； Melosh ——证明有限元位移法是里兹法另一形式； 冯康——独立证明了有限元法； Wilson——第一个编写通用有限元软件：SAP； 1960-70年代理论基础研究 1960至今：实际工程应用、复杂问题理论研究 2.1 概述 1960年代前：起步阶段 1960年代-70初：理论体系建立，快速发展期； 1970初-80中：巩固期； 1980中之后：推广、综合应用； 2.1.2 有限单元法概念 1) 离散化 划分为有限数目的单元； 单元间在指定点连接——结点。 ——单元形状、连接方式可不同。 平面问题的常用单元： 2) 单元分析 假设位移模式，分析单元力学特性： {F}e=[k]{?}e 体力、面力 ——? 等效结点荷载 3) 整体分析—— 建立{F}=[K]{?} [K]{?}= {R} 4) 再次单元分析 求出各单元的应变和应力。 2.2 弹性力学平面问题的矩阵描述 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) 例如：挡土墙、重力坝 2.2.2 基本量及基本方程的矩阵表示 1) 基本量 2) 基本方程 几何方程： 物理方程： (平面应力问题) [D]: 弹性矩阵; E: 弹性模量; ? : 泊松比。 平面应力 ? 平面应变 平衡方程弱形式—— 能量原理 在集中力{F}作用下, {F}=[U1 V1 U2 V2 … Un Vn]T {?}=[u1 v1 u2 v2 … un vn]T 2.3 三结点三角形单元分析 单元局部编号：i, j, m 单元结点位移向量: {?}e=[ui vi uj vj um vm]T 单元结点力向量: {F}e=[Ui Vi Uj Vj Um Vm]T 2.3.2 单元位移模式 1) 什么是位移模式（位移函数） 利用单元的结点位移将整个单元的位移分量表示为坐标的函数。 将位移模式写成结点位移的显式： 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]= [Bi Bj Bm ], [S]= [Si Sj Sm ]，对平面应力问题： 2.3.4 有限元方程的建立 2.3.5 单元刚度矩阵 单刚列式 [k]元素的物理意义 kpq：第q个结点位移分量为单位位移（其它结点位移=0），所引起的第p个结点力分量。 如: k26。 [k]的性质： (1) 对称性