3第三章 单元类型及单元刚度矩阵.ppt

《3第三章 单元类型及单元刚度矩阵.ppt》由会员分享，可在线阅读，更多相关《3第三章 单元类型及单元刚度矩阵.ppt（107页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、3第三章单元类型及单元刚度矩阵第三章第三章 单元类型及单元刚度矩阵单元类型及单元刚度矩阵 有限元法的基本原理是将结构划分成单元，在单有限元法的基本原理是将结构划分成单元，在单元内用较简单的函数描述单元位移，即元内用较简单的函数描述单元位移，即 这是对单元位移这是对单元位移u(x)u(x)的近似。在前面两章的介绍的近似。在前面两章的介绍中，我们讲过，是用单元的节点位移来描述单元内中，我们讲过，是用单元的节点位移来描述单元内点位移，这里所用的变量点位移，这里所用的变量q qi i，是节点位移的一种推，是节点位移的一种推广，即一组广义坐标，或称广义节点位移，包括节广，即一组广义坐标，或称广义节点位移

2、，包括节点位移和节点位移导数。点位移和节点位移导数。N Ni i为形状函数。根据单元为形状函数。根据单元广义节点位移的不同，形状函数分两类：广义节点位移的不同，形状函数分两类：LangrangeLangrange和和HermiteHermite型。型。三三维维的的帕帕斯斯卡卡三三角角形形一、形状函数类型及其特征一、形状函数类型及其特征形状函数应该满足以下条件形状函数应该满足以下条件1.1.2.2.3.3.保证所定义位移函数在相邻单元之间的连续保证所定义位移函数在相邻单元之间的连续4.4.保证所定义位移函数反映常应变状态保证所定义位移函数反映常应变状态一、形状函数类型及其特征一、形状函数类型及其

3、特征 工程实际中有一种结构，特征为：存在一个长维，但工程实际中有一种结构，特征为：存在一个长维，但相对而言又不像平面应变那样，长短比略小，且载荷可相对而言又不像平面应变那样，长短比略小，且载荷可以为任意。比较典型的是井架、塔架等框架结构，这类以为任意。比较典型的是井架、塔架等框架结构，这类结构可用有限元中的一维单元来离散，根据问题的不同，结构可用有限元中的一维单元来离散，根据问题的不同，一维单元又可分为杆单元和梁单元。一维单元又可分为杆单元和梁单元。二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵1.1.杆单元杆单元 杆单元受轴向力，在单元端点处无弯矩和扭矩作用，杆单元受轴向力，在单元端点

4、处无弯矩和扭矩作用，将此单元独立出来进行受力分析时为二力杆。根据单元将此单元独立出来进行受力分析时为二力杆。根据单元形状函数的阶次，又可分为一次杆单元和二次杆单元。形状函数的阶次，又可分为一次杆单元和二次杆单元。一次杆单元一次杆单元 单元有两个节点，如图所示，编号为单元有两个节点，如图所示，编号为i i、j,j,采用局部采用局部坐标坐标 ，记，记 ，并取，并取i i为为x x坐标的原点，则有坐标的原点，则有i(1i(1)j(2j(2)l lFFx 根据形状函数的定义，我们知道，形状函数是根据形状函数的定义，我们知道，形状函数是描述或反映单元内点位移与单元节点位移的关系。描述或反映单元内点位移与

5、单元节点位移的关系。对于上述问题，已知节点位移为对于上述问题，已知节点位移为u ui i,u,uj j,而要求节点而要求节点间任一内点的位移，显然可以根据线性插值来计算间任一内点的位移，显然可以根据线性插值来计算（二点一次拉氏插值），即（二点一次拉氏插值），即 1.1.杆单元杆单元 一次杆单元一次杆单元二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵代入代入 ，有，有 令令 得得 所以单元内点位移为所以单元内点位移为 单元应变单元应变 1.1.杆单元杆单元 一次杆单元一次杆单元二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵所以，几何矩阵为所以，几何矩阵为 单元应力为单元应力为 弹性矩

6、阵弹性矩阵 单元刚度矩阵通式为单元刚度矩阵通式为 二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵代入，得代入，得 这是一次杆单元的单刚阵，它对这是一次杆单元的单刚阵，它对称、对角线元素大于零且奇异！称、对角线元素大于零且奇异！1.1.杆单元杆单元 一次杆单元一次杆单元 当上述单元用于描述仅受扭转变形的杆件时，当上述单元用于描述仅受扭转变形的杆件时，其单刚阵类似于一次杆单元的单刚阵，为：其单刚阵类似于一次杆单元的单刚阵，为：i(1i(1)j(2j(2)l lMnxMn二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵1.1.杆单元杆单元 一次杆单元一次杆单元二次杆单元二次杆单元 单元有三

7、个节点，如图所示，端点编号为单元有三个节点，如图所示，端点编号为i i、j,j,三三个节点依次为个节点依次为1 1、3 3、2 2。单元位移可以根据抛物线插。单元位移可以根据抛物线插值（亦称三点两次拉氏插值）获得，即值（亦称三点两次拉氏插值）获得，即i(1i(1)j(2j(2)l lFFx(3)(3)同样令同样令 二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵1.1.杆单元杆单元 令令二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵1.1.杆单元杆单元 二次杆单元二次杆单元所以单元内点位移为所以单元内点位移为 单元应变单元应变 二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵1.

8、1.杆单元杆单元 二次杆单元二次杆单元几何矩阵为几何矩阵为 单元应力为单元应力为 单元刚度矩阵单元刚度矩阵 二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵1.1.杆单元杆单元 二次杆单元二次杆单元元素的计算元素的计算 可以直接应用可以直接应用 元素的计算元素的计算 二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵1.1.杆单元杆单元 二次杆单元二次杆单元元素的计算元素的计算 二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵1.1.杆单元杆单元 二次杆单元二次杆单元其余元素利用对称性可求得其余元素利用对称性可求得 2.2.三次梁单元三次梁单元 梁单元如图所示，仅考虑节点在梁单元如图

9、所示，仅考虑节点在xoyxoy平面内的位移平面内的位移为为v v、，这时一个单元有四个自由度，形状函数为，这时一个单元有四个自由度，形状函数为三次多项式，即使用三次三次多项式，即使用三次HermiteHermite插值多项式。插值多项式。二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵1 12 2l lMz1xMZ2yzQy1Qy12.2.三次梁单元三次梁单元二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵HermiteHermite位移插值多项式位移插值多项式2.2.三次梁单元三次梁单元 其中其中二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵2.2.三次梁单元三次梁单元 根据平

10、面梁弯曲变形公式（忽略剪切变形）根据平面梁弯曲变形公式（忽略剪切变形）其中其中二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵2.2.三次梁单元三次梁单元 二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵同样令同样令 2.2.三次梁单元三次梁单元单元应力为单元应力为 单元刚度矩阵单元刚度矩阵 引入引入 二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵2.2.三次梁单元三次梁单元单单元元刚刚度度矩矩阵阵 元素的计算元素的计算 二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵2.2.三次梁单元三次梁单元二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵元素的计算元素的计算 2.2

11、.三次梁单元三次梁单元二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵元素的计算元素的计算 2.2.三次梁单元三次梁单元二、一维单元及其单元刚度阵二、一维单元及其单元刚度阵元素的计算元素的计算 其余元素利用对称性可求的其余元素利用对称性可求的 二维单元用于分析和解决平面问题和轴对称为二维单元用于分析和解决平面问题和轴对称为题。在第二章中已详细介绍过，而且是在直角坐题。在第二章中已详细介绍过，而且是在直角坐标中推导的。在下面这一节中，我们将介绍两种标中推导的。在下面这一节中，我们将介绍两种平面单元，即三角形单元和四边形单元，包括一平面单元，即三角形单元和四边形单元，包括一次和二次三角形单元以

12、及一次四边形单元。次和二次三角形单元以及一次四边形单元。三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元三角形单元按其位移的阶数分为一、二、三次单元。三角形单元按其位移的阶数分为一、二、三次单元。一次三角形单元一次三角形单元 第二章详细介绍过这种单元，其形状函数是坐标第二章详细介绍过这种单元，其形状函数是坐标的一次多项式，推导采用直角坐标。对于高次三角的一次多项式，推导采用直角坐标。对于高次三角形单元，这类坐标不方便，特此引入面积坐标。形单元，这类坐标不方便，特此引入面积坐标。1.1.三角形单元三角形单元 一次三角形单元一次三角形单元三、二维单元及其单元刚度阵三

13、、二维单元及其单元刚度阵面积坐标面积坐标xyA1A2A3P(x,y)如图所示，在三角形单元如图所示，在三角形单元A A1 1A A2 2A A3 3中，中，有任意一点有任意一点P(x,y)P(x,y)连接连接PAPA1 1、PAPA2 2、PAPA3 3,得到三个小三角形：得到三个小三角形：PAPA2 2A A3 3、PAPA3 3A A1 1、PAPA1 1A A2 2,记面积比为：记面积比为：称称1 1、2 2、3 3为为P P点的面积坐标点的面积坐标三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标面积坐标 由于由于1

14、1+2 2+3 3=1=1，因此该三，因此该三个坐标不独立。其负号的规定为：个坐标不独立。其负号的规定为：分子分母对应的三角形顶点编号分子分母对应的三角形顶点编号转向相同时为正，反之为负，由转向相同时为正，反之为负，由于三角形于三角形A A1 1A A2 2A A3 3的顶点编号一般规的顶点编号一般规定为逆时针，因此，子三角形顶定为逆时针，因此，子三角形顶点编号为逆时针时面积坐标为正，点编号为逆时针时面积坐标为正，反之为负。反之为负。xyA1A2A3P(x,y)三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标面积坐标xyA1

15、A2A3GA6A5A4三角形中一些特殊点的面积坐标三角形中一些特殊点的面积坐标顶顶点点边边中中点点形心形心三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标与直角坐标的关系面积坐标与直角坐标的关系 三角形三个顶点在直角坐标系中的坐标为（三角形三个顶点在直角坐标系中的坐标为（x xi i,y,yi i）,则则A A1 1A A2 2A A3 3的面积为的面积为类似地三个小三角形的面积依次为类似地三个小三角形的面积依次为三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面

16、积坐标与直角坐标的关系面积坐标与直角坐标的关系三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标与直角坐标的关系面积坐标与直角坐标的关系其中其中三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标与直角坐标的关系面积坐标与直角坐标的关系同理同理三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标与直角坐标的关系面积坐标与直角坐标的关系其中其中令令i=1,2,3;ii=1,2,3;i、j j、k k

17、按按1,2,31,2,3轮转轮转三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标与直角坐标的关系面积坐标与直角坐标的关系则则于是有于是有三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标与直角坐标的关系面积坐标与直角坐标的关系对于单元的三个角点，应有对于单元的三个角点，应有A A1 1点点A A2 2点点三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标与直角坐标的关系面积坐标与直角坐标的关

18、系A A3 3点点重写上述三式重写上述三式三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元面积坐标与直角坐标的关系面积坐标与直角坐标的关系所以有所以有于是有于是有三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元形状函数形状函数根据形状函数的定义根据形状函数的定义i=1i=1时，对于时，对于A A1、A A2、A A3 3点点 i=1,2,3i=1,2,3所以所以三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元形

19、状函数形状函数同理，同理，i=2i=2、3 3时，对于时，对于A A1、A A2、A A3 3点点 所以所以于是于是单元应变单元应变 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元位移函数位移函数三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元其中其中 为几何矩阵为几何矩阵单元应变单元应变 利用复合函数求偏导数的公式利用复合函数求偏导数的公式三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元几何矩阵几何矩阵 三

20、、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元可得可得平面应力问题的应力为平面应力问题的应力为单元刚度矩阵单元刚度矩阵 几何矩阵几何矩阵 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 一次三角形单元一次三角形单元单元刚阵单元刚阵二次三角形单元二次三角形单元xy123645o 二次三角形单元，如图所示，二次三角形单元，如图所示，单元共六个节点，单元共六个节点，1212个自由度，个自由度，三个角节点，三个边中点。三个角节点，三个边中点。三、二维单元及其单刚阵三、二维单元及其单刚阵1.1.三角形单元三角

21、形单元 二次三角形单元二次三角形单元形状函数形状函数i=1,2i=1,2，6 6节点面积坐标节点面积坐标xy123645o对于对于A A1、A A2、A A3 3、A A4 4、A A5 5、A A6 6点点 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 二次三角形单元二次三角形单元形状函数形状函数形状函数分量形状函数分量三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 二次三角形单元二次三角形单元形状函数形状函数形状函数分量形状函数分量同样，同样，对于对于A A1、A A2、A A3 3、A A4 4、A A5 5、A A6

22、 6点点 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 二次三角形单元二次三角形单元形状函数形状函数形状函数分量形状函数分量同样，同样，对于对于A A1、A A2、A A3 3、A A4 4、A A5 5、A A6 6点点 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 二次三角形单元二次三角形单元形状函数形状函数形状函数分量形状函数分量同样，同样，对于对于A A1、A A2、A A3 3、A A4 4、A A5 5、A A6 6点点 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 二次三角

23、形单元二次三角形单元形状函数形状函数形状函数分量形状函数分量同样，同样，对于对于A A1、A A2、A A3 3、A A4 4、A A5 5、A A6 6点点 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 二次三角形单元二次三角形单元形状函数形状函数形状函数分量形状函数分量同样，同样，对于对于A A1、A A2、A A3 3、A A4 4、A A5 5、A A6 6点点 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 二次三角形单元二次三角形单元位移函数位移函数其中其中 为几何矩阵为几何矩阵单元应变单元应变 三、二维单元及其

24、单元刚度阵三、二维单元及其单元刚度阵几何矩阵几何矩阵1.1.三角形单元三角形单元 二次三角形单元二次三角形单元三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵几何矩阵几何矩阵平面应力问题的应力为平面应力问题的应力为1.1.三角形单元三角形单元 二次三角形单元二次三角形单元j j、k k根据根据i=4i=4、5 5、6 6依次按依次按1,2,31,2,3轮转轮转三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵1.1.三角形单元三角形单元 二次三角形单元二次三角形单元注意使用积分公式注意使用积分公式单元刚度矩阵单元刚度矩阵 2.2.矩形单元矩形单元三、二维单元及其单元刚度阵三、二维单元及其

25、单元刚度阵 矩形单元如图所示，共四个节点，每个节点两个矩形单元如图所示，共四个节点，每个节点两个自由度，单元共自由度，单元共8 8个节点位移。为计算方便，引入个节点位移。为计算方便，引入新的变量：新的变量：xy123o4令令 一次矩形单元一次矩形单元2.2.矩形单元矩形单元 一次矩形单元一次矩形单元三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵其中其中 四个角点新坐标四个角点新坐标 xy123o4三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵形状函数形状函数利用形状函数的性质，可得利用形状函数的性质，可得 验证验证 2.2.矩形单元矩形单元 一次矩形单元一次矩形单元三、二维单元及其

26、单元刚度阵三、二维单元及其单元刚度阵位移函数位移函数其中其中 为几何矩阵为几何矩阵单元应变单元应变 2.2.矩形单元矩形单元 一次矩形单元一次矩形单元三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵其中其中 几何矩阵几何矩阵2.2.矩形单元矩形单元 一次矩形单元一次矩形单元三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵几何矩阵几何矩阵2.2.矩形单元矩形单元 一次矩形单元一次矩形单元三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵平面应力问题的应力为平面应力问题的应力为单元刚度矩阵单元刚度矩阵 2.2.矩形单元矩形单元 一次矩形单元一次矩形单元三、二维单元及其单元刚度阵三、二维单

27、元及其单元刚度阵其中其中对于平面应变问题对于平面应变问题2.2.矩形单元矩形单元 一次矩形单元一次矩形单元三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵2.2.矩形单元矩形单元 一次矩形单元一次矩形单元三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵2.2.矩形单元矩形单元二次矩形单元二次矩形单元xy123o48765令令 八个节点新坐标八个节点新坐标 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵2.2.矩形单元矩形单元 二次矩形单元二次矩形单元xy123o48765形状函数形状函数利用形状函数的性质，可得利用形状函数的性质，可得 三、二维单元及其单元刚度阵三、二维单元及其单

28、元刚度阵2.2.矩形单元矩形单元 二次矩形单元二次矩形单元三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵2.2.矩形单元矩形单元 二次矩形单元二次矩形单元位移函数位移函数其中其中 为几何矩阵为几何矩阵单元应变单元应变 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵2.2.矩形单元矩形单元 二次矩形单元二次矩形单元几何矩阵几何矩阵平面应力问题的应力为平面应力问题的应力为单元刚度矩阵单元刚度矩阵 三、二维单元及其单元刚度阵三、二维单元及其单元刚度阵2.2.矩形单元矩形单元 二次矩形单元二次矩形单元 工程中的一切问题都对应着空间三维问题，都可以工程中的一切问题都对应着空间三维问题，都可以

29、用三维单元来构成其总体结构。本课程介绍三维单元用三维单元来构成其总体结构。本课程介绍三维单元及其单刚阵，包括六面体、四面体和曲线等参单元。及其单刚阵，包括六面体、四面体和曲线等参单元。四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元 单元如图所示，共单元如图所示，共8 8个节点，个节点，每个节点的位移参数是每个节点的位移参数是u u、v v、w w，在进行单元分析时，同矩形，在进行单元分析时，同矩形单元一样，常用局部坐标单元一样，常用局部坐标 表示，其原点位于六面表示，其原点位于六面体形心，坐标方向同体形心，坐标方向同x x、y y、z z一一致，其相互关系为：

30、致，其相互关系为：yz123o4x8567yz123o4x8567四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元令令 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元形状函数形状函数八个角点新坐标八个角点新坐标 yz123o4x8567四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元位移函数位移函数验证验证 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元单元应变单元应变其中其中 为几何矩阵为几何矩阵四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体

31、单元六面体单元几何矩阵几何矩阵四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元几何矩阵几何矩阵三维问题的应力为三维问题的应力为四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元弹性矩阵弹性矩阵单元刚度矩阵单元刚度矩阵 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元单元刚度矩阵单元刚度矩阵 为为3x33x3的块方阵，的块方阵，i i，j=1j=1，2 2，8 8，且，且 四、三维单元及其单元刚度阵四、三维单元及其单元刚度

32、阵1.1.六面体单元六面体单元单元刚度矩阵单元刚度矩阵 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元单元刚度矩阵单元刚度矩阵 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵1.1.六面体单元六面体单元单元刚度矩阵单元刚度矩阵 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵2.2.四面体单元四面体单元 工程实际中的结构往往比较复杂，仅用形状规则工程实际中的结构往往比较复杂，仅用形状规则的单元难于较好的近似结构的几何边界，下面介绍的单元难于较好的近似结构的几何边界，下面介绍多用于过度单元的多用于过度单元的4 4面体面体4 4节点三维单元。节点三维单元

33、。1234P 采用体积坐标，单元内任意一采用体积坐标，单元内任意一点点P P的位置由的位置由4 4个比值来确定：个比值来确定：2.2.四面体单元四面体单元V V是四面体的体积是四面体的体积 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵2.2.四面体单元四面体单元 其中其中 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵2.2.四面体单元四面体单元其中其中 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵2.2.四面体单元四面体单元形状函数形状函数 从而可以推的位移函数、单元应变、几何矩阵、弹从而可以推的位移函数、单元应变、几

34、何矩阵、弹性矩阵、刚度矩阵。性矩阵、刚度矩阵。积分中用到积分中用到四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元 前面介绍的几种单元几何形状规则，便于进行运前面介绍的几种单元几何形状规则，便于进行运算，但难以适应工程实际的需要，工程实际中零部算，但难以适应工程实际的需要，工程实际中零部件的外形基本上都比较复杂。曲边等参单元可以解件的外形基本上都比较复杂。曲边等参单元可以解决这个矛盾，这种单元可以用曲边单元划分实际结决这个矛盾，这种单元可以用曲边单元划分实际结构，而按直边单元进行计算，中间用坐标变换来转构，而按直边单元进行计算，中间用坐标变换来转换之，把（换

35、之，把（x,y,zx,y,z）转换成（）转换成（）。）。坐标变换坐标变换 用平面单元来说明。子单元是任意四边形，母单元用平面单元来说明。子单元是任意四边形，母单元是正方形。是正方形。四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元坐标变换坐标变换 xy123o41234四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元坐标变换坐标变换 用等分四边形的两族直线划分子单元，以子单用等分四边形的两族直线划分子单元，以子单元形心为新坐标（元形心为新坐标（）的原点，在这个坐标系）的原点，在这个坐标系中，母单元内任意一点的位移为中，母单

36、元内任意一点的位移为 其中其中 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元坐标变换坐标变换 其中的其中的N Ni i与形状函数中的与形状函数中的N Ni i相同。相同。设单元内点坐标与单元节点坐标之间的关系：设单元内点坐标与单元节点坐标之间的关系：四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元坐标变换坐标变换 在上式中令在上式中令 得得 这是直线这是直线2-32-3的方程，同样可以得到直线的方程，同样可以得到直线1-21-2、3-43-4、4-14-1的方程。的方程。由于在坐标变换时使用的形状函数与位移函数使用由于

37、在坐标变换时使用的形状函数与位移函数使用的形状函数相同，所以这类单元称为等参数单元。的形状函数相同，所以这类单元称为等参数单元。四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元等参单元的单元分析等参单元的单元分析 单元分析单元分析 对对n n个节点的空间等参单元，可以类似地写出位个节点的空间等参单元，可以类似地写出位移函数和坐标变换关系：移函数和坐标变换关系：四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元单元分析单元分析 有了位移函数后，就可以的得到几何矩阵，在有了位移函数后，就可以的得到几何矩阵，在几何矩阵中要出现形状

38、函数分量对整体坐标的导几何矩阵中要出现形状函数分量对整体坐标的导数，根据复合函数导数规则：数，根据复合函数导数规则：等参单元的单元分析等参单元的单元分析 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵可表示为可表示为 式中，式中，为雅克比矩阵为雅克比矩阵 3.3.曲边等参单元曲边等参单元单元分析单元分析 等参单元的单元分析等参单元的单元分析 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵可以推得可以推得3.3.曲边等参单元曲边等参单元单元分析单元分析 等参单元的单元分析等参单元的单元分析 将将代入代入 ，得：，得：四、三维单元及其

39、单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元单元分析单元分析 等参单元的单元分析等参单元的单元分析 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵有有求出单元的形状函数后，如求出单元的形状函数后，如3.3.曲边等参单元曲边等参单元单元分析单元分析 等参单元的单元分析等参单元的单元分析 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元单元分析单元分析 等参单元的单元分析等参单元的单元分析 对局部坐标求导对局部坐标求导数，并和节点坐标数，并和节点坐标一块代入一块代入 中并中并求出求出 的逆阵，的逆阵，从而可以计算出形从而可以计算出形

40、状函数对整体坐标状函数对整体坐标的导数，这样就可的导数，这样就可以获得以获得 。单元刚度矩阵单元刚度矩阵 其分块矩阵为其分块矩阵为 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元单元分析单元分析 等参单元的单元分析等参单元的单元分析 四、三维单元及其单元刚度阵四、三维单元及其单元刚度阵3.3.曲边等参单元曲边等参单元等参单元的单元分析等参单元的单元分析 单元刚度矩阵单元刚度矩阵 上式是在整体坐标中进行，且有上式是在整体坐标中进行，且有 是是 的行列式的行列式 于是有于是有高斯积分高斯积分此此课件下件下载可自行可自行编辑修改，修改，仅供参考！供参考！感感谢您的支持，我您的支持，我们努力做得更好！努力做得更好！谢谢!