工程弹塑性力学课件.ppt

工程弹塑性力学(有限元、塑性力学部分) 演示稿 第0章 平面问题的有限单元法 0.1 概述、基本量及基本方程的矩阵表示 0.2 有限单元法的概念 0.3 位移模式与解答的收敛性 0.4 单元刚度矩阵 0.5 等效结点荷载 0.6 整体刚度矩阵 0.7 单元划分应注意的问题 0.1 概述、基本量及基本方程的矩阵表示 弹性力学问题的求解方法： 解析方法：函数解、级数解——少数简单问题 数值方法： 有限单元法的发展概况 1956年提出 1960-70年代理论基础研究 1960至今：实际工程应用、复杂问题理论研究 通用有限元软件：SAP、ADINA、NASTRAN ANSYS、ABAQUS等 基本量及基本方程的矩阵表示 1. 基本量： 2. 基本方程： 3. 虚功方程： 外力虚功=内力虚功 在集中力{F}作用下，虚功方程简化为 {F}=[U1 V1 U2 V2 … Un Vn]T 为结点力向量； {d}=[u1 v1 u2 v2 … un vn]T 为结点位移向量。 0.2 有限单元法的概念 1. 离散化 划分为有限数目、有限大小的单元。 平面问题的常用单元： 2. 单元分析 建立： {F}e=[k]{d}e [k]：单元刚度矩阵 单元结点力向量{F}e=[Ui Vi Uj Vj Um Vm]T 单元结点位移向量{d}e=[ui vi uj vj um vm]T 体力、面力 ——? 等效结点荷载 3. 整体分析 建立： {F}=[K]{d}， [K]：整体刚度矩阵 由各结点平衡{F}={R}，得有限元方程： [K]{d}= {R} 0.3 位移模式与解答的收敛性 1. 什么是位移模式（位移函数） 利用单元的结点位移将整个单元的位移分量表 示为坐标的函数。 2. 三结点三角形单元的位移模式 设：u=a1+a2 x+a3 y v=a4+a5 x+a6 y 系数a1~a6由结点位移 ui , vi , uj , vj , um , vm确定。 将位移模式写成结点位移的显式： u= Niui+ Njuj +Nmum v= Nivi+ Njvj +Nmvm Ni、Nj、Nm：形函数 (插值函数) 形函数的性质 (1) (Ni )i=1，(Ni )j=0，(Ni )m=0 (2) 单元内任一点：Ni+Nj+Nm=1 3. 位移模式的矩阵表示 4. 位移模式应满足以下条件，才能保证有限元解答收敛： (1) 位移模式必须能反映单元的刚体位移 (2) 位移模式必须能反映单元的常量应变 (3) 位移模式应尽可能反映位移的连续性 三结点三角形单元的完备性和连续性： (1) 反映刚体位移： (2) 反映常量应变：ex=a2， ey=a6，gxy=a2+a3 三结点三角形单元的完备性和连续性： (3) 位移连续性： ▲单元内：单值连续； ▲相邻单元之间：uij(1)=uij(2)？vij(1)=vij(2) ？ ij边的方程：y=ax+b，则 uij=a1+a2 x+a3(ax+b)= cx+d uij(1)、uij(2)均为坐标的线性函数，故可由i、j两点的结点位移唯一确定。 0.4 单元刚度矩阵 建立： {F}e=[k]{d}e 单元刚度矩阵： 结点位移 位移 应变 应力 结点力 {d}e ——?{f} ——?{e} ——?{s} ——? {F}e [B]：应变矩阵； [S]：应力矩阵； 三结点三角形单元中，[B]、 [S] 的元素均为常数，故这种单元又称 常应变单元，或常应力单元。 [k]元素的物理意义 kpq：第q个结点位移分量为单位位移（其它结点位移=0），所引起的第p个结点力分量。 如 k25： [k]的性质： (1) 对称性： kpq= kqp (2) 奇异性； (3) 每行（列）元素之和为零。 (4) [k]取决于单元的形状、方位和弹性常数，与所在位置（即平移或np 转动）无关。 0.5 等效结点荷载 非结点荷载需等效移置 到结点上。 移置原则：静力等效