常用单元的刚度矩阵Word格式.docx
- 文档编号:3230350
- 上传时间:2023-05-01
- 格式:DOCX
- 页数:16
- 大小:22.06KB
常用单元的刚度矩阵Word格式.docx
《常用单元的刚度矩阵Word格式.docx》由会员分享,可在线阅读,更多相关《常用单元的刚度矩阵Word格式.docx(16页珍藏版)》请在冰点文库上搜索。
u(r
)
T
zrz
所以弹性方程为式中应力矩阵
弹性矩阵D
(1~~)(12)
12
2
4.单元刚度矩阵k⑹
与平面问题相同,仍用虚功原理来建立单元刚度矩阵,其积分式为
k(e)BTDBdV
V
在柱面坐标系中,dV2drdz
将dV2drdz代入k(e)BTDBdV,贝k⑹2BTDBrdrdz
即为轴对称问题求单元刚度矩阵的积分式。
与弹性力学平面问题的三角形单元不同,在轴对称问题
中,几何矩阵[B]内有的元素(如丛1田等)是坐标r、z的函r
数,不是常量。
因此,乘积btdb不能简单地从式k(e)2BTDBrdrdz的积分号中提出。
如果对该乘积逐项求积分,将是一个繁重的工作。
一般采用近似的方法:
用三角形形心的坐标值代替几何矩阵[B]内的r和z的值。
用B表示
在形心(&
)处计算出的矩阵[B]。
其中
(rirjrk)(乙ZjzQ
r,z
33
只要单元尺寸不太大,经过这样处理引起的误差也不
大。
被积函数又成为常数,可以提出到积分号外面:
k(e)2BTDBrdrdz2BTDBr式中三角形的面积。
由式k⑹2BTDBrdrdz2BTDBF可以看出,两轴对称的三角形单元,当形状、大小及方位完全相同而位置不同时,其刚度矩阵也不相同。
距离主轴线越远的单元,其刚度越大。
这与平面问题不一样。
二、等参数的刚度矩阵
对一些由曲线轮廓的复杂结构,如果采用直角边单元进行离散,由于用直线代替了曲线,除非网格划分得很细,否则不能获得较高的精度;
对另一些应力随坐标急剧变化的结构,采用简单的常应力单元离散时,也必须划分成大量的微小单兀,以保证足够的精度。
为此引入一种高精度的单兀一—等参数单元。
它既能简化复杂单元划分的工作,又能在满足同样精度的要求时,大大减少使用的单元数。
目前流行的大程序中较常用,它成功地解决了许多二维和三维的弹性力学问题。
为导出等参数单元的刚度矩阵,首先要建立根据每个单元的形状确定的自然坐标系,然后将位移模式和形状函数都写成自然坐标的函数。
一个单元在自然坐标系内的点余元整体坐标系内的点成对应的关系。
通过映射,可以将整体坐标系中的图形转化为自然坐标系中的相应徒刑。
例如可以将整体坐标系中
的一个任意四边形(实际单元)映射到自然坐标系中成为一个正方形(基本单元)。
同样也可以将任意四面体、六面体(包括直边和曲边的)分别映射成正四面体和正六面体。
这里只介绍较简单的一种平面问题的情况,将整体坐标系中的一个任意四边形映射成自然坐标系中的一个正方体,并导出单元刚度矩阵。
其它种单元的映射,可依次原理进行不再叙述。
1.位移模式和形状函数
图4-2中的任意四边形单元上,作连接对边中点的直线,取其交点为原点,这两条直线分别为和轴,并令四条边上的和值分别为1,建立一个新的坐标系,称之为该单元的自然坐标系。
原坐标系XOY称为整体坐标系。
在整体坐标系中,自然坐标系非正交,它由任意四边形的形状所确定
图4-19
如果将自然坐标系改画成直角坐标系,那么图4-19(a)中的任意四边形单元就成为图4-19(b)所示的正方形。
上述两个四边形的点(包括顶点)一一对应,即它们之
间相互映射。
因此,需要写出整体坐标X、Y和自然坐标
之间的坐标转换式,即
12
56
34
78
四边形四个顶点的坐标值在
XOY坐标系中分别为
Xi,Yi,X2,Y2,X3,Y3,X4,Y4
坐标系中相应
1,1,1,1,1,1,1,1。
将有关数据代入
中的第一式,则有
X11234,X2
X31234,X4
求解上述方程组得:
X1X2X3X4
14,2
34,
X1X2
4
X3X4
坐标变换方程
X11
*成为
X11
X21
X31
X4
同理
Y11
丫21
丫4
当引入函数Ni
XNi
i1
YNi
后,
Xi
Yi
式中Ni
41
变量、
当i=4时,
坐标变换方程成为
的正负号由相应节点的坐标值
41,41,因此,N4
i、i决定。
例如
下面再来研究函数Ni,的特性。
对节点1X1,Y1,相应的自然坐标值为(-1,-1)。
从式
Ni,11i1i中很容易看出,除N1=1外,
N2=N3=N4=0。
对其余各节点也一样。
总而言之,对节点i(i=1,2,3,4),除Ni=1外,其余三个N值均为零。
同时,不难看出
N1,N2,N3,N4,1,即四个节点的Ni函数之和等于1。
函数Ni,具备上章所介绍的形状函数应满足的条件,
可作为本单元的形状函数。
采用Ni,做形状函数,
U
i
Ni,Ui,V
1
Ni
i1
XNi,
对比
和U
丫Ni,
丫
Vi
Ui,VNi
Vi可以看出:
其位移模式为
在这种实际单元(任意四边形)中,坐标变换式和位移模式不仅采用了相同的形状函数Ni,,而且具有相同的数学
模型。
这种性质的实际单元称为等参数单元。
对用节点位移值Ui(或Vi等)求单元内某一点位移量u(或v
等)的插值公式uNi,Ui,只要将u(或V等)换成X(或丫
等),便成为利用节点值Xi(或Yi等)求相应点坐标X(或丫等)的插值公式。
相反也是这样。
2.几何矩阵[B]
由于几何矩阵[B]通过对位移求偏导数而得出,所以首先
必须利用复合函数求导的规则得出下述公式
uuXuYuu
XuX
丫、/或写成jx
uYuu
yy
X
Y
式中J
,此式称为雅可比矩阵。
为了将几何矩阵[B]写成变量、的函数,必须将式
jX改写成
y
v
J1u,同理,
J1
从表示单元内各点位移与其应变关系的几何方程可知:
对单元(e),任意一点的位移
u,v对自然坐标
的偏
10
‘u
vv
00
XY
01
将式
J
1u和
合并,则
vJ
uvv
Y0
导数可利用上式求出,写成矩阵形式为:
Np
式中
U1
v1u1
V|u1
V1
N1p
N2p
N3p
N4p
对于
i=1,2,3,4
Nip
则可得出表示
在整体坐标系中位移和应变关系的几何方程:
式中的几何矩阵[B]是自然坐标,的函数:
01Np
也可利用
Nip以及
和
求出
o
Y4
X1X2X3
N1pN2pN3pN4pyyy
Y1Y2y3
3.单元刚度矩阵k
设单元板厚为t,根据虚功方程有:
keBTDBtdA,
A
此式中几何矩阵[B]和弹性矩阵[D]都已求出。
因为几何矩阵
[B]中的变量是自然坐标,,所以也要用自然坐标表示微分面积dA。
在实际单元中任取一点p,其整体坐标位X、丫,其相应的自然坐标为,。
过p点做,的等值线,同时做d,d的等值线,围成一小块微分面积dA,如图4-20
(a)所示。
为便于分析,将四边形pqrs放大,如图4-20
(b)所示。
实际上,d,d取得很小,因此该四边形可视为
平行四边形。
若相邻的两边用向量a,b表示,则两者的乘积
等于该平行四边形的面积dA。
图4-20
dAa|bsinab
若aaxiayj,bbxibyj,则
axbx
dAaxiayjbxibyj匚
ayby
为了求出ax,ay,bx,by的值,要先写出a和b两端节点p、q、
s的坐标值。
点p:
XpX
,YpY,
点q:
XqX
d,,YqY
d,
点s:
XsX,
d,丫sY
d
利用泰勒技术展开并略去高阶项,可得
Xd
d,也可写出相应的展开式。
利用式
-Ad
X可得:
—d
ax
XqXp
Yd
bx
XsXp
by乂Yp
将此式代入式dA
axiayjbxibyj
ay
by
得到:
dA
axbxayby
XYdd,简写为dAJdd
单元刚度矩阵为:
11
keBTDBtJdd,这个积分可以采用“数值方法”,用
高斯求积分公式很方便的求出,在此不作介绍。
例:
求如图所示四边形的雅可比矩阵。
解:
求雅可比矩阵可在整体坐标系中进行,也可以在实际单元的局部坐标系中进行。
为便于计算,本例在局部坐标系中进行。
对单元
(1):
将四个节点的自然坐标值(-1,-1)、(1,-1)、(1,1)、
(-1,1)代入下式:
M「1i,再将所得到的2,值及四个节点实际单元在局部坐标系的坐标值(-3,-2)、(3,-2)、(3,
2)、(-3,2)代入下式计算:
XNi,Xi
,则X3,Y2
YM,Y
雅克比矩阵为
30
JXYo2
对单元
(2):
四个节点在局部坐标系中的坐标值分别为(-1,-3/4)、
(1,-3/4)、(1,
X111
3
Y—1
因而雅可比矩阵为
5/4)、(-1,1/4)
11111
11-1
因此,
5
-1
也可利用J
N1pN2pN3pN4pY1Y2Y3Y4求雅可比矩阵,
Y1丫2Y3Y4
其结果与上相同,同学们可自行验证
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 常用 单元 刚度 矩阵