东南大学计算力学习题及答案版讲解.docx
《东南大学计算力学习题及答案版讲解.docx》由会员分享,可在线阅读,更多相关《东南大学计算力学习题及答案版讲解.docx(36页珍藏版)》请在冰点文库上搜索。
东南大学计算力学习题及答案版讲解
第三章
1
如图所示一三角形钢板,两个结点固定,对第三个结点施以单位水平位移,测出所施加的力,从而得出相应的刚度系数。
其他点依此类推,这样测得的刚度系数所组成的刚度矩阵,是否与按照常规三角形单元刚度矩阵计算公式所得结果一样?
用这样实测所得的刚度矩阵能否进行有限元分析?
为什么?
解:
不一样。
单元刚度矩阵中每个元素的物理意义:
kij表示单元第j个自由度产生单位位移,其它自由度固定时,
第i个自由度产生的节点力。
单元刚度矩阵是在单元处于平衡状态的前提下得出的,单元作为分离体看待,作用在它上面的外力(单元力)必是平衡力系,然而研究单元平衡时没有引入约束承受平衡力系作用的无约束单元,其变形是确定,但位移是不能确定的,即单元可发生任意的刚体位移。
不能。
因为与有限元中单元与单元之间的约束情况不一样,不能进行有限元分析。
2以位移为基本未知量的有限元法其解具有下限性质,试证明之。
解:
系统总位能的离散形式p1aTKaaTP
p2
将求解的方程KaP带入可得
p1aTKaaTKa1aTKaU
p22在平衡情况下,系统总位能等于负的应变能。
在有限元解中,由于假定的近似位移模式一般来说总与精确解有差别的。
设近似解为p、U、[K]、a、KaP,真实解为p、U、[K]、a、KaP
且根据最小势能原理,得到的系统的总位能总会比真正的总位能要大,故pp则UU
aTKaaTKaaTPaTP则近似解的位移总体上小于精确解的位移
解释如下:
单元原是连续体的一部分,具有无限多个自由度,在假定了单元的位移函数后,自由度限制为只有以结点位移表示的有限自由度,引入了更多的约束和限制,使得单元刚度较实际连续体加强了,连续体的整体刚度随之增加,所以有限元解整体上较真实解偏小。
3请分别阐述单元刚度矩阵和整体刚度矩阵中任一元素的物理意义。
解:
在单刚Ke中,kiej表示单元第j个位移产生一单位位移,其它位移为零时,第i个位移方向上引起的节点力。
在整体刚度中,Kij表示第j个自由度产生一单位位移,其它自由度为零时,第i个自由度上引起的节点力。
4简述虚功原理,且使用虚功原理导出外荷载与节点荷载的等效关系式。
解:
虚功原理:
变形体中任意满足平衡的力系在任意满足协调条件的变形状态上作的虚功等于零,即体系外力的虚功与内力的虚功之和等于零。
设q为外荷载(此处为体力),p为节点荷载,w为单元内位移场,为结点位移场
eTeeTe
根据虚功原理pwqdV
V
由于we[N]e故weqedVe[N]TqedVe[N]TqedV
VVV
则epee[N]TqedVpe[N]TqedV
VV
5
试述弹性力学中按位移求解与有限单元法中按位移求解之间的异同点。
解:
弹性力学
有限单元法
物理模型
连续体
离散化结构
基本方程
几何方程
几何方程
物理方程
物理方程
平衡微分方程
结点平衡方程
解法
解微分方程
解代数方程
解答形式
用函数表示
用数值表示
解答精度
精确解
近似解
6如果三节点三角形单元绕其中某一个节点作小的刚体转动,其转角为,证明单元内所有的应力均为零。
解:
在三角形单元中DB
则可知节点位移向量0,0,yj,xj,ym,xm
由于弹性矩阵D为常量矩阵,应变向量
为零向量,故D为零向量,即单元内所有的应力为零。
7
刚度矩阵相同。
在平面内旋转时,刚度矩阵也相同。
二维单元在x,y坐标内平面平移到不同位置,单元刚度矩阵相同吗?
在平面内旋转时又怎样?
试证明之。
解:
二维单元在x,y坐标内平面移到不同位置时,
单元平移或旋转时,bi,ci不变,故单元刚度矩阵不变。
8判断有限元网格离散合理性
a)对图1(a)所示的有限元网格,评论网格的优劣性,指出模型中的错误,并加以改正。
b)评论图1(b)的网格划分合理吗?
为什么?
请加以改正。
图1
解:
(a)网格划分不合理。
1)无过渡单元
2)无边界条件
3)夹角区应力集中,应适当加密风格
4)对称结构网格应对称划分
(b)不合理。
1)左部网格应适当加密
2)由于三角形单元会造成局部精度不够,过渡区可采用其它单元划分
3)右部单元的长宽比较大,就进行适当调整。
9如图2所示,平面三角形构件以x-y坐标系表示的刚度矩阵方程如下
10
2.5
1.83
2.5
ux1
Px1
4
1.83
2.5
5.0
2.5
vy1
Py1
104
Px2
2.5
4.5
2.5
2.5
ux2
2.5
2.5
2.5
2.5
vy2
Py2
ux1
ux1
解:
用坐标变换
'T则
vy1
vy1
'y1
ux2
ux2cos
vy2
u'x2sin
1
0
0
0
其中T其中T
0
1
0
0
0
0
cos
0
0
0
sin
0
试建立以ux1,uy1,ux2,(与图中Px2,同向的位移)及Px1,Py1,Px2'
来表示的刚度矩阵方程。
由K'PKTP
10
1.83
2.5
2.5
1.83
5.0
2.5
2.5
KT
2.5
4.5
2.5
2.5
2.5
2.5
2.5
2.5
1000
0100
102.52.9640
4
1.832.52.50
000
5
2.54.50.50
3
2.52.50.50
000
5
10
4
1041.83
0
3所示。
试求结点2的等效
2.52.964ux1Px1
2.52.5vy1Py1
20u'x2Px2'
10某平面结构采用四节点矩形单元和三节点三角形单元建立有限元计算模型,其如图荷载列阵R2。
11
解:
单元①,荷载作用于12边上,故等效节点力只与1、2号节点有关
则P1yN1qydsql2dql
l03
单元③,
形函数N1
1l,N
ls,N
0
在1-2边上,qss
sql
FS[N]
l
ss
T
0
10
0
0
0
ll
s
01s0l
sl
0
0
qls
qSds
ds
l
0T
ql01
23
230
故节点2的等效荷载列阵
0
R2ql2
3
1,且设Ke表示第e个单元的单元
12试求如图4所示的有限元网格的整体刚度矩阵,假设每个节点的自由度数为刚度矩阵(注意:
结果应该用kiej表示)。
图5
,K
(2)
k5(35)
K(3)k7(35)
k8(35)
解:
单元刚度矩阵K
(1)
k(3)k(4)k(4)
k58k55k56
k7(38),K(4)k6(45)k6(64)
k(3)k(4)k(4)
888586
k1(11)
k2(11)
0
k1(21)k2(12)k2(22)
k3(22)
0k2(23)k3(23)
K
k4(11)
k4(12)
0
k5(11)
k5(12)k5(22)
k5(23)
0
k6(22)
k6(23)
0
0
0
0
0
0
整体刚度矩阵:
K①和应力矩阵S①是
8
0
6
6
2
6
16
6
12
6
4
0
0
-3
0
3
0
K
(1)
对
13.5
9
7.5
3
S
(1)
0
4
0
-3
0
-1
13.5
3
1.5
2
0
1.5
-1.5
-0.5
1.5
9.5
3
称
5.5
12图5中两个三角形单元组成平行四边形,已知单元①按局部编码
i,j,m的单元刚度矩阵
按图5示单元②的局部编码写出K②,S②。
解:
由图可知m
(1)i
(2),i
(1)j
(2),j
(1)m
(2)
3
0
0
0
-3
0
S
(2)
0
-1
0
4
0
-3
-0.5
1.5
2
0
-1.5
-1.5
13如图6所示8结点矩形单元(每边中点为结点),3点为坐标原点,元厚为t。
①求该单元的位移函数和形函数和并检验其是否满足收敛性条件。
②求在2-6-3边作用均布水平荷载q时的等效结点荷载。
解:
(1)位移函数:
u12x3y4x25xy6y27x2y8xy2
2222
v910x11y12x213xy14y215x2y16xy2
则n2时,10,21,10,21l11,l2,p11,p2则角节点的形函数为
1111
N14
(1)
(1),N24
(1)
(1)
(1)
2222
N34(12)
(1)(12)
(1),N44(21)(12)
(1)
边中节点的形函数为
N54
(1),N64
(1)
(1),N74
(1)
(1),N84
(1)
证明收敛性:
位移函数中
2222
u12x3y4x5xy6y7xy8xy
2222
v910x11y12x13xy14y15xy16xy
1,9表示刚体位移,2,3和9,10表示常应变,故位移函数具有完备性设相邻单元公共边界上的直线方程是yb(或xa),代入位移函数中
222
u13b6b2(25b8b2)x(47b)x2
222
v911b14b2(1013b)16b2x(1215b)x2
为x(或y)的2次函数,而边界上三点确定的位移函数为也为二次曲线,故单元在公共边界连续,
故位移函数收敛
(2)荷载作用在23边上,故等效节点力只与2,6,3号节点有关
1111
N24()()
(1),N64
(1)
(1),N34()
(1)()
(1)
2
222
N2N6N3
N241,N64(12),N343
x0,y0N3bN6
bN22
ds
(x)2(y)2d2d
1
P3xtN3qxds2qtN3d
s0
qt
3
4
P6xtN6qxds2qtN6d3qtP2xtN2qxds
3
2qtN2dqt
03
第四章
1经典梁理论和Timoshenko梁理论有哪些相同点和哪些不同点?
基于以上两种理论的梁单元各有何特性?
解:
经典梁理论
Timoshenko梁理论
相同点
Kirchhoff假设
C1型单元
C0型单元
不同点
弯曲梁单元
考虑剪切变形影响
截面转动是挠度w的一阶导数,只有挠度w是独立的
挠度w和截面转动各自独立插值
采用Hermite插值
采用拉格朗日插值
特性
梁的高度远小于跨度
梁很薄时,会造成剪切锁死现象
2写出杆件的应变能计算公式,并给出推导过程。
解:
将只考虑轴向变形的杆件划分成n个单元,节点坐标为x0,x1,,xi,xi1,,xn
单元的位移函数u(x)12x(xi1xxi)
用形函数近似位移函数得u(x)Nie1(x)ui1Nieui,其中Nie1(x)i,Nie(x)i1
xi1xixixi1
1[11]uie[B]uie
xixi1
单元的应力EE[B]uie
xi1
其中Kie[B]TEA[B]dx
xi
3
还有什么方法可以实现以上条件,并比较这几种方
在杆系系统中,除了采用凝聚自由度的方法实现铰接端条件,
法的优缺点。
解:
优点
缺点
凝聚自由度法
8
4利用最小势能原理,推导图1所示弹性基础上梁单元方程,其中该梁的势能为:
2
L
0wvdx
图1
解:
根据最小势能原理可知p0
lll
故有p(EIv'')v''dxkfvvdxwvdx0
000
l
对第一项分部积分(EIv'')v''dx(EIv'')v'l0(EIv'')'v'(EIv'')v'l0(EIv'')'vl0(EIv'')''v0
[(EIv'')''kfvw]vdx(EIv'')v'l0(EIv'')'vl000
引入强制边界条件和自然边界条件使(EIv'')v'l0(EIv'')'v00
由于v的任意性故控制微分方程为(EIv'')''kfvw0
此梁的位移函数vh(x)N1e(x)v1N2e(x)1N3e(x)v3N4e(x)4[N]d,则v[N]d
由于物理关系可知v''(x)Bd则v''[B]d
lll
由(EIv'')v''dxkfvvdxwvdx0得000
5图2所示刚架
L12Lkfv
p012EI(v")2dx0f2dx
1)如何进行节点编号使整体刚度矩阵[K]的带宽最小?
2)刚架的整体刚度矩阵中a节点的总刚度矩阵Kaa和的总刚度矩阵Kbc各由哪些单元的哪些分块矩阵叠加组成(自行确定单元局部坐标方向)
3)试按照二维等带宽存储和一维变带宽存储方式确定Kaa中对角元素的在相应存储数组中的位置。
图2有铰点的刚架
解:
1)考虑每个节点有两个自由度
由于半带宽d=(相邻结点码的最大差值+1)*2故节点编号如图所示可使单元内节点编码相差3,使得带宽d=8
2)KaaK2(22)K2(32)K1(14)K1(51)KbcK1(62)
3)考虑单元①节点1自由度的凝聚
可知Kaa中对角线元素在原整体刚度矩阵中第6行第7列和第7行第8列
则用二维等带宽存储后在矩阵中的第6行第2列和第7行第2列用一维变带宽存储后在Kaa中对角线元素在数组中的位置为9和10
(1)
(1)
K
K
0
0
0
0
0
0
0
0
0
11
12
K
(1)
K
(1)K
(2)
0
K
(2)
0
0
0
0
0
0
0
21
2211
12
(3)
(3)
0
0
K
K
0
0
0
0
0
0
0
11
12
K
(2)K(3)
0
K
(2)
K(3)
2222
K(4)
0
K(5)
0
0
0
0
21
21
K(4)K(5)
12
12
1111
0
0
0
K(4)
K(4)K(6)
0
0
K(6)
0
0
0
21
2211
12
(7)
(7)
0
0
0
0
0
K
K
0
0
0
0
11
12
K(5)K(7)
0
0
0
K(5)
0
K(7)
2222
K(8)
0
0
0
21
21
K(8)K(9)
12
1111
0
0
0
0
K(6)
0
K(8)
K(6)
K(8)K(10)
0
0
K(10)
21
21
22
2211
12
(11)
(11)
0
0
0
0
0
0
0
0
K
K
0
11
12
0
0
0
0
0
0
K(9)
0
K(11)
K(9)
K(11)K(12)
K(12)
21
21
22
2211
12
(10)
(12)
(10)(12)
0
0
0
0
0
0
0
K
0
K
KK
21
21
2222
10
第五章
1根据以下形函数表达式
132313223
N13(2x33x2LL3)N23(x3L2x2L2xL3)
1321322N3L3(2x33x2L)N4L3(x3Lx2L2)
画出形函数N1和N3以及导数(dN2/dx)和(dN4/dx),它们代表梁单元整个长度上形状变化2、对于图1中给出的四节点二次应变一维等参单元,试确定:
a)形函数N1,N2,N3,N4;
b)单元刚度矩阵[k]。
图1
x1
N2N3N4x2
x33x4
16(81221)
41
[k]3E
31
4(321)
3
4(321)
31(81221)
6
1(81221)
4(32
1)4(321)
16(81221)d
11
3试利用变节点数法构造插值函数的,构造出图2所示的三次三角形单元的形函数及相应的位移函数。
解:
位移函数:
223223
u12x3y4x5xy6y7x8xy9xy10y
v1112x13y14x215xy16y217x318x2y19xy220y3
1)构造不考虑边节点和内部节点的角节点的插值函数:
4)修正角节点的插值函数:
211
N1N?
132(N4N9)13(N5N8)13N10
333
2991991
L1[L1L2(3L11)L1L3(3L11)][L1L2(3L21)L1L3(3L31)]27L1L2L3
13212121313212221333123
1
2(3L11)(3L12)L1
211
N2N?
23(N5N6)3(N4N7)3N10
333
2991991
L23[2L1L2(3L21)2L2L3(3L21)]3[2L1L2(3L11)2L2L3(3L31)]327L1L2L3
3223223
1
12(3L21)(3L22)L2
12
211N3N?
33(N7N8)3(N6N9)3N10
333
2991991
L332[92L2L3(3L31)92L1L3(3L31)]31[29L2L3(3L21)29L1L3(3L11)]1327L1L2L3
1
2(3L31)(3L32)L3
4
试构造如图3所示的15结点三棱柱体单元的插值函数,并判断其构造的位移函数是否收敛。
(4)修正角节点的插值函数:
N1N?
11(N10N12)1N71L1
(1)1[2L1L2
(1)2L1L3
(1)]1L1(12)
22222
1
L1
(1)(2L22L3)
2
13
111112
N2N?
2(N10N11)N8L2
(1)[2L1L2
(1)2L2L3
(1)]L2(12)
22222
1
L2
(1)(2L12L3)
2
N3N?
31(N11N12)1N91L3
(1)1[2L2L3
(1)2L1L3
(1)]1L3(12)
22222
1
L3
(1)(2L12L2)
2
N4N?
41(N13N15)1N7
22
1112
12L1
(1)21[2L1L2
(1)2L1L3
(1)]12L1(12)
1
L1
(1)(2L22L3)
2
111112
N5N?
5(N13N14)N8L2
(1)[2L1L2
(1)2L2L3
(1)]L2(12)
22222
1
L2
(1)(2L12L3)
2
111112
N6N?
621(N14N15)12N921L3
(1)12[2L2L3
(1)2L1L3
(1)]21L3(12)
1
L3
(1)(2L12L2)
2
1