常用单元的刚度矩阵.docx
- 文档编号:29640632
- 上传时间:2023-07-25
- 格式:DOCX
- 页数:11
- 大小:138.99KB
常用单元的刚度矩阵.docx
《常用单元的刚度矩阵.docx》由会员分享,可在线阅读,更多相关《常用单元的刚度矩阵.docx(11页珍藏版)》请在冰豆网上搜索。
常用单元的刚度矩阵
2龙(/・+“)一2加・_u
2岔r
由于各点在圆周方向上无位移,因而剪应变%和怙均为
根据上式,可推导出几何方程{^}=[b]M
3.弹性方程和弹性矩阵[D]
依照广义虎克定律,同样可以写出在轴对称中应力和应变之间的弹性方程,其形式为
6胡勺一心+空)]
J=云&一“(6+刃)]
2(1+//)
N_—E—找
所以弹性方程为匕}=[切{耳
式中应力矩阵{cr}={br(T0a.rj
4.单元刚度矩阵[灯“
与平面问题相同,仍用虚功原理来建立单元刚度矩阵,其积
分式为
[屮訂[町[卬抄
在柱面坐标系中,dV=iTttlrdz
将dV=2加加7z代入\k]ey=J[BY[b[B}lV,则[k]e}=2叩[j?
]r\p[B\-drdz
V
即为轴对称问题求单元刚度矩阵的积分式。
与弹性力学平面问题的三角形单元不同,在轴对称问题中,几何矩阵[B]有的元素(如业勺等)是坐标r.z的函r
数,不是常量。
因此,乘积[bY[dIb]不能简单地从式=2叩[町[D\B}Mz的积分号中提出。
如果对该乘积逐项求积分,将是一个繁重的工作。
一般采用近似的方法:
用三角形形心的坐标值代替几何矩阵[B]的r和Z的值。
用屈表示在形心(展)处计算出的矩阵[B]。
其中
-a+5+乙)-(Zi+zj+zk)
,=3—3
只要单元尺寸不太大,经过这样处理引起的误差也不大。
被积函数又成为常数,可以提出到积分号外面:
[k]e)=2/可[Q][可[Jrdrdz.=2/可[功[亦△式中△——三角形的面积。
由式旧⑺=2兀两[功広j]rdrdz.=2龙两[功厨込可以看出,两轴对称的三角形单元,当形状、大小及方位完全相同而位置不同时,其刚度矩阵也不相同。
距离主轴线越远的单元,其刚度越大。
这与平面问题不一样。
二、等参数的刚度矩阵
对一些由曲线轮廓的复杂结构,如果釆用直角边单元进行离散,由于用直线代替了曲线,除非网格划分得很细,否则不能获得较高的精度;对另一些应力随坐标急剧变化的结构,釆用简单的常应力单元离散时,也必须划分成大量的微小单元,以保证足够的精度。
为此引入一种髙精度的单元一—等参数单元。
它既能简化复杂单元划分的工作,又能在满足同样精度的要求时,大大减少使用的单元数。
目前流行的大程序中较常用,它成功地解决了许多二维和三维的弹性力学问题。
为导出等参数单元的刚度矩阵,首先要建立根据每个单元的形状确定的自然坐标系,然后将位移模式和形状函数都写成自然坐标的函数。
一个单元在自然坐标系的点余元整体坐标系的点成一一对应的关系。
通过映射,可以将整体坐标系中的图形转化为自然坐标系中的相应徒刑。
例如可以将整体坐标系中的一个任意四边形(实际单元)映射到自然坐标系中成为一个正方形(基本单元)。
同样也可以将任意四面体、六面体(包括直边和曲边的)分别映射成正四面体和正六面体。
这里只介绍较简单的一种平面问题的情况,将整体坐标系中的一个任意四边形映射成自然坐标系中的一个正方体,并导出单元刚度矩阵。
其它种单元的映射,可依次原理进行。
不再叙述。
1.位移模式和形状函数
图4-2中的任意四边形单元上,作连接对边中点的直线,取其交点为原点,这两条直线分别为纟和”轴,并令四条边上的纟和〃值分别为±1,建立一个新的坐标系,称之为该单元的自然坐标系。
原坐标系XOY称为整体坐标系。
在整体坐标系中,自然坐标系非正交,它由任意四边形的形状所确定。
图4-19
如果将自然坐标系改画成直角坐标系,那么图4-19(a)中的任意四边形单元就成为图4-19(b)所示的正方形。
上述两个四边形的点(包括顶点)一一对应,即它们之间相互映射。
因此,需要写出整体坐标X、Y和自然坐标承〃之间的坐标转换式,即
X=④+勺歹+切+也勿*
Y=a5+a^+如〃+禺勿
四边形四个顶点的坐标值在XOY坐标系中分别为(x1,yJ(x2,y2),(X3,y3),(x4,Kl):
在蚀坐标系中相应为(7-1),(1,-1),(1,1),(71)。
将有关数据代入*中的第一式,则有
X]=«)-(z2-a3+of4,X2=at+a2
Xi=al+a1+ay+(z4,X4=a(-a2+a3-(z4
求解上述方程组得:
x,+x.+x.+x4-x.+x.+x.-x4
务=—亍—s=—r-—
-X.-X.+X^+X.E—X,+(—Xq
5=——T-—‘久=—冷——1
坐标变换方程*成为
x+(1+§-§77)X2+(1+§+〃+纫)X3+(1--
同理
y=#[(i—§一"+勿加+(1+"〃一切必+(】+§+〃+切必+(i-初必]当引入函数Ni(M后,坐标变换方程成为
f-1
y=(刖比
r-I
式中mG,〃)=]i+£§Xi+7")
变量§、〃的正负号由相应节点的坐标值知”,决定。
例如当i二4时,勺=-1,〃厂1,因此,n4(的)=上犁也。
下面再来研究函数n,D的特性。
对节点1(X「YJ,相应的自然坐标值为(-1,-l)o从式mG,“)=土(1+/X1+M)中很容易看出,除叫二1外,n2=n3=n4=o0对其余各节点也一样。
总而言之,对节点i(1=1,2,3,4),除Ni二1外,其余三个N值均为零。
同时,不难看出
”佐〃)+“2(§,〃)+弘低〃)+汕低77)=1,即四个节点的Ni函数之
和等于lo
函数n,(m)具备上章所介绍的形状函数应满足的条件,可作为本单元的形状函数。
采用做形状函数,其位移模式为
"=工7]兀」=工v(艮“加
在这种实际单元(任意四边形)中,坐标变换式和位移模
式不仅采用了相同的形状函数NS,而且具有相同的数学模型。
这种性质的实际单元称为等参数单元。
对用节点位移值Ui(或Vi等)求单元某一点位移量U(或V等)的插值公式(釦7”,只要将11(或V等)换成X(或Y
J-1
等),便成为利用节点值Xi(或Yi等)求相应点坐标X(或Y等)的插值公式。
相反也是这样。
2.几何矩阵[B]
由于几何矩阵[B]通过对位移求偏导数而得出,所以首先
必须利用复合函数求导的规则得出下述公式
du
[ondXdudY]
du
du
i
oX6$dY
>或打成<
苗
•=[4
6X
du
oudXdudY+
du
du
cXdtj*dYdt]
6Xar
式中[小勇第,此式称为雅可比矩阵。
dr]d/]
为了将几何矩阵[B]写成变量品〃的函数,必须将式
从表示单元各点位移与其应变关系的几何方程可知:
0
_100o'
J
d
u
<
>=
0001
dududv6v;
dY
V
施dYdXdY
d
0110
•■丿
ax.
d
dX
对单元(e),任意一点的位移u,v对自然坐标乡〃的偏
导数可利用上式求出,写成矩阵形式为:
[色色色色[沙)
1%讪g期JL
式中{肖}心={"JV,M1vtMtV,Vjy
[0]
[o]
[0]
K\
Ml
=LM
kip]
[o]
K]
[0]
[心]
[0]
[ml
在整体坐标系中位移和应变关系的几何方程:
⑷⑷=冈%严
式中的几何矩阵[B]是自然坐标弓"的函数:
r.7X=YNi(^fl)Xi
也可利用k]=簧}求得的[砒以及;■和
1-1
3.单元刚度矩阵旧⑷
设单元板厚为t,根据虚功方程有:
叶)胡[3]⑪][30,
A
此式中几何矩阵[B]和弹性矩阵[D]都已求出。
因为几何矩阵[B]中的变量是自然坐标乡〃,所以也要用自然坐标表示微分面积dAo
在实际单元中任取一点P,其整体坐标位X、Y,其相应的自然坐标为乡“。
过P点做乡〃的等值线,同时做
§+妬,〃+切的等值线,围成一小块微分面积dA,如图4-20
(a)所示。
为便于分析,将四边形pqrs放大,如图4-20
(b)所示。
实际上,始切取得很小,因此该四边形可视为平行四边形。
若相邻的两边用向量{必紛表示,则两者的乘积等于该平行四边形的面积dA。
图4-20
dA=”制sinO=\ax/?
|
若ci=5i+ciyj,b=bj+byj,则
dA=腕丿+ayj\bxi+byj]=:
;(
•'*
为了求出代,心,》,伏的值,要先写出a和b两端节点p、q、s的坐标值。
点p:
Xp=X(§,讥丫厂Y(§,小
点q:
x厂xG+姑,〃比=«§+〃§,〃)
点s:
XS=XG,“+〃讥Y、=心"+dri)
濮亦OX一期
利用泰勒技术展开并略去高阶项,可得
x(§+d的)=x(的)+
xG,"+d“)=xG,“)+
对吃,也可写出相应的展开式。
利用式
可得:
畑伽)=如)+『xG,〃+〃〃)=xG,〃)+学切
J=Xq-Xp=養站心=Yq—Yp=鲁旳
»=X$-X卩==YS-Yp=—drj
将此式代入式"A=|(dJ+a:
j\bxi+化j]="',得到:
dyDy
单元刚度矩阵为:
[田二jj[町[D间叶歟,这个积分可以采用"数值方法”,用
-1-1
高斯求积分公式很方便的求出,在此不作介绍。
例:
求如图所示四边形的雅可比矩阵。
解:
求雅可比矩阵可在整体坐标系中进行,也可以在实际单元的局部坐标系中进行。
为便于计算,本例在局部坐标系中进行。
对单元
(1):
将四个节点的自然坐标值(-1,-1)、(1,-1)、(1,1)、(-1,1)代入下式:
n@〃)=扌(1+爲Xi+〃刃),再将所得到的m(§,司值及四个节点
实际单元在局部坐标系的坐标值(-3,-2)、(3,-2)、(3,2)、(-3,2)代入下式计算:
dXdY
雅克比矩阵为[小馬就
df]di]
对单元
(2):
四个节点在局部坐标系中的坐标值分别为(-1,-3/4)、
(1,—3/4)、(1,5/4)、(-1,l/4)o因此,
x=才{(1-§Xi-"X-i)+(i+歹X】-祕1)+(1+纟)0+祕1)+(1-欽1+祕-1)}
因而雅可比矩阵为
"_丄「4(1+初
[匕[o(3+効
也可利用[丿卜购“In“In”In钟]][:
:
:
:
彳求雅可比矩阵,其结果与上相同,同学们可自行验证。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 常用 单元 刚度 矩阵