有限元分析课后习题.docx
- 文档编号:12976035
- 上传时间:2023-06-09
- 格式:DOCX
- 页数:24
- 大小:527.46KB
有限元分析课后习题.docx
《有限元分析课后习题.docx》由会员分享,可在线阅读,更多相关《有限元分析课后习题.docx(24页珍藏版)》请在冰点文库上搜索。
有限元分析课后习题
3.1“强”形式相关的场变量要求强的连续性。
定义这些场变屋的所有函数必须可微,而可微的次数必须等于存在于强形式的系统方程中的偏微分方程的次数。
“弱”形式通常是积分形式,且对场变量要求较弱的连续性,弱形式通常能得到更精确的解。
3.2⑻协调性方程
(b)本质边界条件或运动边界条件
(c)在初始刻和末时刻的条件
3.3(a)域的离散
(b)位移插值
(O构造形函数
(d)坐标变换
(e)整体有限元方程的组装
(f)位移约束的施加
(g)求解整体有限元方程
3.4理论上不用必须离散所求解问题的区域。
把问题划分成单元的目的是更容易地假设位移场的模式。
3.5证明:
(1)方程的左边为
=(&J+—(Sal)/~+-(Sa2)/5
方程的右边为
=3[aJ+-aj~+-a2P]
=[(网)/+(旳)尹+(Sa2)-P]
很显然方程的左右两边相等。
(2)方程的左边为
5"⑴=3(a.+2a,x)cLv
方程的右边为
很显然方程的左右两边相等。
3.6再生性和连续性
形函数是线性无关的
德尔塔函数性质
单位分解性
线性场再生性
3.7答案很明显
3.8为了把所有的单元方程组合起来构成整体的系统方程,必须对每个单元进行坐标变换。
3.9组装的过程就是把与某个节点相连的所有单元的贡献相加。
4.1桁架构件通过销钉或较链(而不是焊接)连接在一起,因此构件之间值传递力(而不是力矩)。
因此,一个桁架结构仅仅有轴向变形,我们只分析轴向变形,我们可以沿杆单元的轴向取为局部坐标的x轴,在单元的每个节点处只有一个自由度,即轴向位移。
4.2局部坐标系D0F=2,整体坐标系D0F=4.
在局部坐标系中,桁架单元仅仅考虑轴向变形,因此一个节点仅有一个自由度。
整体坐标系用于描述桁架结构的所有单元,不能保证桁架结构的所有单元坐标轴总是沿着轴向变形的方向,因此,一个节点的自由度需要两个位移来描述,因为一个节点能在两个方向上有位移,所以在整体坐标系中,一个节点有两个自由度。
4.3结果类似于4.2.
4.4
4.5
_o.oip/rc+2ie+i]
6|_-十13/+2J
4.6很明显这个结构师对称的,因此,我们可以仅仅分析结构的一半如下图所示,
简化后,各单元的尺寸参数如下表所示。
单元号
横截面积,Aem2
长度lem
杨氏模量EN/m2
1
0.005
0.5
69x109
2
0.01
1.118
69x109
第一步:
计算各单元的的方向余弦
知道了整体坐标系中的节点坐标后,第一步就是考虑单元对于整体坐标系的方向,由于该问题是一个平面问题,因此只需要计算出lij和mij,各节点在整体坐标系中的坐标和个单元的方向余弦如下表所示:
单元编号
对应的整体节点
整体坐标系中的坐标
方向余弦
1•••■
lijmij
局部节
点1(i)
局部节点2(j)
Xi,Yi
Xj「j
1
1
2
0,-0.5
0,0
01
2
2
3
0,
0.5,1.0
0.44720.8945
由于节点1和2仅沿垂直方向运且外力也是竖直的,所以节点2相对于节点2的位移为
杆1的应力为
£=2£0_=10xl05N/m-
1A0.005
第二步:
计算整体坐标系中的单元矩阵
计算出方向余弦后,就可以求出整体坐标系中的单元矩阵,注意到该问题是一个静力学问题,因此不需要计算单元质量矩阵,只需要求刚度矩阵,局部坐标系中的单元刚度矩阵是2x2的矩阵,因为每个单元总的自由度为2。
但是在变换到整体坐标系时,每个单元的自由度数就变为4,因此,整体坐标系中的刚度矩阵就变成了4x4矩阵,刚度矩阵如卞所示:
0
(0.005)(69xl09)
o3
巧.
■0000
6.90-6.9
X
00
4.938
第三步:
整体有限元矩阵的组装
求出了单元矩阵后,下一步就是把单元矩阵组装为整体有限元矩阵,因为结构的整体自由度为6,所以整体刚度矩阵为6x6矩阵,吧与某一节点相连的所有单元对该节点的贡献加起来就可以完成对该节点的组装,最后得到整体刚度矩阵为
第四步,施加边界条件
施加边界条件通常能见效整体矩阵的大小,在这里,如果约束DI,D2,D3和D5,于是有
DI=D3=D5=D6=0这就是说,第一、第二、第三和第五行实际上对求解矩阵方程无影响,因此,我们只需删除相应的行和列,有
D]D2£>3DqD5
L|
n
)0(
6.9(
0(
-6.9(
(
(
D,
rx
L=10sx
1.2
34~~2/169=T7
11.838-2.
Z34~=2?
469-4.!
)38
「2
T6
1.2
TD5
S
*■
J
TD6
凝聚了的整体矩阵就变成了2x2矩阵,如下所示
x!
0sNni2
6.9-6.9
K=
-6.911.838
很容易证明上述凝聚的刚度矩阵是对称正定的。
施加约束后的有限元方程为KD=F
式中
D?
=[2D4]
且力向量F给出为
[-500]
F斗
0
注意到节点1只作用一个与D2方向相反的向下的力,式(6-3)实际上等价于包含了两个未知量D2和D4的两个联立方程如下
[6.9D,-6.9Z)4]xl08=-500
[-6.9r>2+11.838D4]xlOs=O
第五步:
求解有限元矩阵方程
为了得到D2和D4的解,求解方程组,得到
D2=1.0121x10-6m
D4=1.7367x10-6m
我们也能得到
=2—2=7・246xl0Tm
计算单元应力:
ZTT9xW贡贡
—1.0126x10“
5.1同桁架单元推导公式一样
5.2在划分网格时,一般应该把受力的点划为节点,如果受力点不是节点,一般把力看成一个沿着单元分布的力。
节点力p,能这样描述
q(x)=P5(x_*)
因此,利用公式得到力向量
f严jNYdU+jNHdS,
VSf
0
0
2)
0
0
皿呵
-p/r
"心0)
-P//8
N«=0)
-P/2
N心0).
P//8
0
5.3为了得到力向量,将该结构分成两个单元,分析作用在每一个单元上的外力和力矩,如
F=<
0
0
>=V
0
0
-p
0
P/2
+(-
0
P/2
)
0
0
0
0
卜•图所示
6.1框架单元具有桁架单元和两单元的性质,框架单元不仅有桁架单元轴向的变形,还具有梁单元在x-y平面内的横向变形和绕着z轴的旋转。
由于变形比较小,所以轴向变形不会影响横向变形,这时利用叠加方法就可以构造框架单元的方程。
当轴向变形和横向变形产生耦合时,即变形比较人时,叠加方法将失效变形。
当变形比较人时,轴向变形能引起横向变行,反过来横向变形也能引起轴向变形,因此叠加方法将会失效。
6.2不用,因为对于平面框架,局部坐标系和整体坐标系Z轴同向。
6.3问题分析
该结构明显是轴对称结构,因此我们仅用一半结构如图所示
0.5m0.5m
各参数如下表所示:
.
单元
横截面积,
长度
杨氏模量
横截面积对Z轴
编号
Aem2
lem
EN/m2
的二次矩,lzm4
1
0.005
0.5
69x109
1.989x10-6
2
0.01
1.118
69x109
7.958x10-6
7rr4_A2
44兀
第一步:
计算单元的方向余弦
知道了节点在整体坐标系中的坐标后,第一步就是要得到局部坐标系与整体坐标系之间的关系,这一关系能通过求方向余弦得到,各节点坐标及方向余弦如下表所示:
单元
对应的整体节点
整体坐标系中的坐标
方向余弦
编号
局部节
点l(i)
局部节
点2(j)
Xi,Yi
Xj,Yj
lx
mx
ly
my
1
1
2
0,05
0,0
0
1
-1
0
2
2
3
0,0
0.5,1.0
0.447
0.894
-0.894
0.447
由于节点1和2仅沿垂直方向运且外力也是竖直的,所以节点2相对于节点2的位移为
fIP(0.5)(500)7
d「=——=—―r__—=7・246xlOTm
-EA69xl09x0.005
杆1的应力为
£=卫_=]0x]05N/m2
1A0.005
第二步:
计算整体坐标系中的单元矩阵
计算出方向余弦后,就可以求出整体坐标系中的单元矩阵,注意到该问题是一个静力学问题,因此不需要计算单元质量矩阵,只需要求刚度矩阵,局部坐标系中的单元刚度矩阵是6x6
0
0
-6.17
0
0
0.0471
0.0263
0
-0.0471
0.0263
0.0196
0
-0.0263
0.0098
6.17
0
0
巧.
0.0471
-0.0263
0.0196
6.17
xlOsNnr2
利用坐标变换,得到整体坐标系中的刚度矩阵
0
-0.0329
-0.1320
0
-0.0329
6.9
0
0
-6.9
0
0.0109
0.0329
0
0.0055
0.1320
0
0.0329
sy.
6.9
0
0.0109
0.1320
xlOsNin2
第三步:
组装整体矩阵
将三个节点的自由度组合在一起,形成9x9矩阵,如下所示:
0.132
0
-0.0329
-0.132
0
-0.0329
0
0
0
0
6.9
0
0
一6.9
0
0
0
0
一0.0329
0
0.0109
0.0329
0
0.0055
0
0
0
-0.132
0
0.0329
0.1320+1.2705
0+2.4468
0.0329-0.0235
-1.2705
-2.4468
-0.0235
0
一6.9
0
0+2.4468
6.9+4.9407
0+0.0118
-2.4468
-4.9407
0.0118
xlOsNm-:
一0.0329
0
0.0055
0.0329-0.0235
0+0.0118
0.0109+0.0196
0.0235
-0.0118
0.0980
0
0
0
-1.2705
-2.4468
0.0235
1.2705
2.4468
0.0235
0
0
0
-2.4468
-4.9407
-0.0118
2.4468
4.9407
-0.0118
0
0
0
-0.0235
0.0118
0.0980
0.0235
-0.0118
0.0196
第四步:
施加边界条件
施加边界条件后,这种情况下,
DI=D3=D4=D6=D7=D8=D9=0
因此对应的行和列对方程的解没有任何影响,所以我们去除相应的行和列。
压缩后的矩阵如下所示:
:
我们也能得到
=D2-D4=7.246x10-7m
7.1虽然单元的厚度是不均匀的,但是在2D平面内,厚度是线性的。
这样我们就能利用节点位置处的厚度和线性形函数的性质插补得到任何点处的厚度。
利用形函数的德尔塔函数性质,形函数在自身节点处应取单位值1,而在其他节点处为0.这是德尔塔函数的性质,作为插入函数,该形函数的应用也是可用的,因此单元厚度能写成:
h(x,y)=NJh+NJs+弘饨=[“N2h2>
.復
这里Ni和hi分别为形函数和节点处的节点厚度。
刚度矩阵能写成
kf=JBrcBdV=j(£dz)B7cBdA=J/?
B7cBdA
带入,得到
kc=B?
cBJ、NR+N2h2+NjgcL4
形函数能够表示成面坐标的形式。
由于面坐标也具有德尔塔函数性质,利用数学方法
(1-7)
加!
川“!
2A
(in+n+p+2)\
得到
所以
k=B%B(丛+丛+丛'
eI333丿
=JiABcB
7.2由于单元的厚度变化是线性的,厚度能利用形函数插值方法得到,因此
力(x,y)=NJh+N2h2+他人=[MN\N3]?
「
这里Ni和hi分别为形函数和节点处的节点厚度。
单元质量矩阵能够写成
m(=jpNrNdV=j「d¥pNrNdA=j/?
pNrNcL4
匕人4
带入形函数得到
利用公式
得到
所以
7.3有刚度矩阵的计算公式,我们看到被积函数为hBTcB.应变矩阵B是g和7]的线性函数,厚度能够利用线性形函数和节点处的厚度值得到。
因此,在每一个方向上,被枳函数是一个立方函数,所以在每个方向上有两个高斯积分点就足以计算出含有最高次数为3次的多项式刚度矩阵。
因此,一个矩形单元,4个高斯节点就足够了。
有质量矩阵的计算公式,我们看到被枳函数为hNTNo厚度能够利用线性形函数和节点处的厚度值得到。
因此,在每一个方向上,被枳函数是一个立方函数,所以在每个方向上有两个高斯枳分点就足以计算出含有最高次数为3次的多项式质量矩阵。
因此,一个矩形单元,
4个高斯节点就足够了。
7.4线性四边形单元的质量矩阵能够写成
mr=JpNyNdV=jJ:
drpN'NdA=J/?
pNyNcL4
V4A
=fi丁:
%N"Ndet|j|dgd〃
厚度函数能够写成
h(x,刃二NR+Nh+NJq+N扎=[MN=N,N」{「丨
厚度是的g和q线性函数,NTN的次数为2,雅克比行列式的次数为2,因此,整个被积函数的次数为5,我们能用3x3个高斯点得到准确解。
单元体积能够表示成
匕=Jj?
dA=j*JI力det川dgd〃
V=/?
cL4=J:
]*:
力det|j|dgd〃
所以被积函数次数为3,利用2x2个高斯点就可得到准确解
形函数是g和q的双线性函数,,将双线性函数进行微分并除以其元素也就是双线性函数饿雅可比矩阵就可以求出应变矩阵中的元素,因此被枳函数/?
RTcBdetlJl是不可能用多项式表示的分式函数,也就是说,利用高斯枳分法不可能准确的计算出刚度矩阵。
7.5利用拉线法构造,线不一定全是直线,也可以是曲线。
8.1对于一个厚板的能量方程表达为
/=扛[:
代必出+扛匚:
刊必丈
如果板不是匀质的而是分层的,沿Z方向的积分将会考虑
例如,一个四层对称等厚度的板,我们应该有
/=£上:
/曲反+扛匸:
VydAdz
=-f\eocL4d<+-f[tvdAdz
+-f2(''e7ocL4dz+-f2fZ,T7ycL4dz
2J®」加42以
假设该板为瑞斯那板,得到刚度矩阵为
专虑[BWBW+扛體[BTgu
+£足酹]7。
必+扛冷BjjBOdA
同样的得到质量矩阵,能量为
Te=(沪+,+诉Mz+2gp2(沪+沪+妒)dAdz
=-\
2九
因此质量矩阵为
P\h
2
「pji
~2
Pji'
96
(1・5)
7pJ?
3
96
P»$
~96~
7亦
96
N’NcL4+fNtI2NdA
JAr
(1-6)
9.1可以,3D实体单元可以用于解决二维实体单元平面应力和平面应变的问题,因为2D单元是3D单元的特例,一些特殊的3D单元我们可以简化为2D单元。
例如,如果结构一个方向厚度远小于其他两个方向的厚度,可以简化为平面应力问题。
当然3D单元也能解决考虑所有的3D问题
9.2没有区别。
这是因为形函数和节点总数完全一样。
如果六边形用的是线性函数,结呆可能不同因为四边形有恒定的应变矩阵,也就是它不如六面体单元准确。
9.3五边形单元如下所示,为了解该单元的方程,我们可以通过从多项式构建形函数。
拎一
个简单的方法是该单元类型由两个四边形单元组成。
11.1结构的几何形状、材料性质和支撑条件都是对称的。
外载荷不一定是轴对称的。
11.2Answer:
letL二a+b
Clampedsymmetricbeamstructure
Deflection=0Rotation二Free
Clampedanti-symmetricbeamstructure
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 有限元分析 课后 习题