2017-东南大学计算力学习题讲解
东南大学结构力学习题(附答案)
M /8
M图
11、 X = ql (←) （有侧支座水平反力） 1 28 13、
P
12、 M CB
= 2.06kN ⋅ m （上侧受拉）
X1
X2
—— 8 ——
15、
17、
P l
X1
Pl 2Pl
Pl
Pl Pl 2 MP 图 Pl 2 M图
M1 图
δ 11 =
2l Pl , Δ1P = − , 3EI EI
Pa
52、
53、
N=4.5qa 117 81 36 1.5 qa 24 9 0 M 图 0 5 qa2 36 _ N= 9qa 4.5qa2 4.5qa
54、
4 54 54
55、
4 N =－28.8 kN 36 36 4 M 图 ( kN . m ) 4
4
56、
2 qa 2 2 qa 2
M 图
57、
第二部分 静定结构内力计算（参考答案）
1、 （O） 6、 （O） 11、 （O） 16、 （O） 19、
3 Pa+1.5 m 0 4 3m 4 0
2、 （X） 7、 （X） 12、 （O） 17、 （X）
3、 （O） 8、 （X） 13、 （O） 18、 （O）
4、 （O） 9、 （O） 14、 （X） 20、
3 30、 Δ DV = 8 Pa / EI + 125 Pa / 4 EA( ↓ )
4 2 31、 Δ DV = 11qa / 24 EI + 15qa / 8 EA( ↓ )
32、 φ SR = 33、 Δ DV
3Pa 2 EI
（
）
东南大学结构力学练习题(附答案)详解
第一部分 平面体系的几何组成分析一、判断题：1、在任意荷载下，仅用静力平衡方程即可确定全部反力和内力的体系是几何不变体系。
2、图中链杆1和2的交点O 可视为虚铰。
O二、分析题：对下列平面体系进行几何组成分析。
3、 4、ACDBAC DB5、 6、A CDBEABCDE7、 8、ABCD GE FA BCDEFGHK9、 10、11、 12、1234513、 14、15、 16、17、 18、19、 20、1245321、 22、124567831234523、 24、12345625、 26、27、 28、29、 30、31、 32、33、ACBDEF三、在下列体系中添加支承链杆，使之成为无多余约束的几何不变体系。
34、35、第二部分 静定结构内力计算一、判断题：1、静定结构的全部内力及反力，只根据平衡条件求得，且解答是唯一的。
2、静定结构受外界因素影响均产生内力，内力大小与杆件截面尺寸无关。
3、静定结构的几何特征是几何不变且无多余约束。
4、图(a)所示结构||M C =0。
(a)BCa aA ϕ2a2 (b)5、图(b)所示结构支座A 转动ϕ角，M AB = 0，R C = 0。
6、荷载作用在静定多跨梁的附属部分时，基本部分一般内力不为零。
7、图(c)所示静定结构，在竖向荷载作用下，AB 是基本部分，BC 是附属部分。
ABC(c)8、图(d)所示结构B 支座反力等于P /2(。
)↑9、图(e)所示结构中，当改变B 点链杆的方向（不通过A 铰）时，对该梁的影响是轴力有变化。
ＡＢ(e)10、在相同跨度及竖向荷载下，拱脚等高的三铰拱，水平推力随矢高减小而减小。
—— 1 ——11、图(f)所示桁架有9根零杆。
(f)a a a a(g)12、图(g)所示桁架有：=== 0。
N 1N 2N 313、图(h)所示桁架DE 杆的内力为零。
a a(h)14、图(i)所示对称桁架在对称荷载作用下，其零杆共有三根。
东南大学-土力学课后习题答案(第三版)
东南大学土力学第二章2-2、有一饱和的原状土样切满于容积为21.7cm 3的环刀内，称得总质量为72.49g，经105℃烘干至恒重为61.28g，已知环刀质量为32.54g，土粒比重为2.74，试求该土样的湿密度、含水量、干密度及孔隙比（要求汇出土的三相比例示意图，按三相比例指标的定义求解）。
解：3/84.17.2154.3249.72cm g V m =-==ρ %3954.3228.6128.6149.72=--==S W m m ω 3/32.17.2154.3228.61cm g V m S d =-==ρ 069.149.1021.11===S V V V e 2-3、某原状土样的密度为1.85g/cm 3，含水量为34%，土粒相对密度为2.71，试求该土样的饱和密度、有效密度和有效重度（先推导公式然后求解）。
解：（1）VV m WV s sat ρρ⋅+=W S m m m += S W m m =ω 设1=S m ρω+=∴1V W S S S V m d ρ=WS W S S S d d m V ρρ⋅=⋅=∴1()()()()()()3W S S W S S W W satcm /87g .1171.20.341171.285.1d 11d 11d 111d 11111=+⨯+-⨯=++-=+++⎪⎪⎭⎫ ⎝⎛-=+-++=+⎪⎪⎭⎫ ⎝⎛⋅-++=∴ρωρωρωρωρρωρρωρρρωρW S d 有（2）()3'/87.0187.1cm g VV V V V V V m V V m W sat WV Ssat W V W V W S S W S S =-=-=+-=-+-=-=ρρρρρρρρρ （3）3''/7.81087.0cm kN g =⨯=⋅=ργ或3'3/7.8107.18/7.181087.1cmkN cm kN g W sat sat sat =-=-==⨯=⋅=γγγργ2-4、某砂土土样的密度为1.77g/cm 3，含水量9.8%，土粒相对密度为2.67，烘干后测定最小孔隙比为0.461，最大孔隙比为0.943，试求孔隙比e 和相对密实度Dr，并评定该砂土的密实度。
(完整版)计算力学复习题答案
计算力学试题答案1. 有限单元法和经典Ritz 法的主要区别是什么？答：经典Ritz 法是在整个区域内假设未知函数，适用于边界几何形状简单的情形；有限单元法是将整个区域离散，分散成若干个单元，在单元上假设未知函数。
有限单元法是单元一级的Ritz 法。
2、单元刚度矩阵和整体刚度矩阵各有什么特征？刚度矩阵[K ]奇异有何物理意义？在 求解问题时如何消除奇异性？答：单元刚度矩阵的特征：⑴对称性⑵奇异性⑶主元恒正⑷平面图形相似、弹性矩阵D 、厚度t 相同的单元，e K 相同⑸e K 的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两列，其位置与结点位置对应。
整体刚度矩阵的特征：⑴对称性⑵奇异性⑶主元恒正⑷稀疏性⑸非零元素呈带状分布。
[]K 的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足力和力矩的平衡。
为消除[]K 的奇异性，需要引入边界条件，至少需给出能限制刚体位移的约束条件。
3. 列式说明乘大数法引入给定位移边界条件的原理？答：设：j j a a =，则将 jj jj k k α=j jj j P k a α=即：修改后的第j 个方程为112222j j jj j j n n jj j k a k a k a k a k a αα+++++=由于得 jj j jj j k a k a αα≈ 所以 j j a a ≈对于多个给定位移()12,,,l j c c c =时，则按序将每个给定位移都作上述修正，得到全部进行修正后的K 和P ，然后解方程即可得到包括给定位移在内的全部结点位移值。
4. 何为等参数单元？为什么要引入等参数单元？答：等参变换是对单元的几何形状和单元内的场函数采用相同数目的结点参数及相同的插值函数进行变换，采用等参变换的单元称之为等参数单元。
借助于等参数单元可以对于一般的任意几何形状的工程问题和物理问题方便地进行有限元离散，其优点有：对单元形状的适应性强；单元特性矩阵的积分求解方便（积分限标准化）；便于编制通用化程序。
东南大学计算力学课件(研究生课程) 第6章 曲边等参单元和数值积分
关键在于，首先要解决进行任意形状单元和规则单元之间的几何变换。
1 坐标系的映射
m
x Ni'xi
m
y Ni' yi
m
z Ni'zi
说明:根据前面收敛性说明可知,对任一连续函数u的唯一性 要求都可以转化为对(x,y,z)坐标的唯一性要求.当相连单 元在相同节点处具有相同的坐标时,则意味着单元是连 续的.
1 坐标系的映射
映射
前面的变换 不会产生这
种情况
不产生缝隙与插值 函数（对应场变量） 要求的协调性是类
似的。
不协调
很容易保证
2 单元的映射变换
其中
(6-2)
N1
1 4
(1 )(1)
N3
1 4
(1 )(1)
N2
1 4
(1 )(1)
N4
1 (1 )(1)
4
(6-3)
2 单元的映射变换
1 坐标系的映射
J 0的条件： d 0 or d 0
or sin(d ,d) 0
η
4
3
η
3 4
ξ
ξ
η
1, 4
d 0
1
正常
2
1, 2
d 0
J 0 4
η
3
ξ
2
sin(d , d) 0
1
即J 0
ξ
η
3
ξ
2
1 坐标系的映射
计算力学大作
以下确定 Q(y)，和 R(x)； 将 u，V 代入 exy 表达式得
PLx PLx 2 dQ( y) vPy 2 dR( x) (1 v) P h 2 ( y2 ) EJ 2 EJ dy 2 EJ dx EJ 4
再分别对 x 和 y 求导并积分有
Q( y )
vPy 3 (1 v) Py 3 (1 v) Pyh 2 Ay B 6 EJ 3EJ 4 EJ PLx 2 Px3 R( x) Ax C 2 EJ 6 EJ
{
联立解得：
c1 2P h3
c2 0
c3
3P 2h
c4
P 2
故： 有上式求解应力：
A
P y y3 ( L x)(1 3 4 3 ) 2 h h
X
2 A P( L x) y y 2 J
{ 利用广义 Hooked 定律
eij
2 A y 0 x 2 2 A P h2 xy ( y2 ) yx 2J 4
目录
计算力学大作业 ......................................................................................................... 1 一． 解析部分 ................................................................................................... 1 二． 数值计算结果 ............................................................................................ 4 第一组： ..................................................................................................... 4 第二组 ......................................................................................................... 8 第三组 ......................................................................................................... 9 第四组 ....................................................................................................... 10 三． COMSOL 计算 .......................................................................................... 11 1x16 网格................................................................................................... 11 1x20 网格................................................................................................... 13 Finer 网格（由 COMSOL 自动生成） ....................................................... 14 较短梁讨论 ....................................................................................................... 15
力学计算2017版(含答案)
力学计算专题知识点：1、速度：ts v = 2、密度：ρρρmV V m V m =⇒=⇒=（单位 1g/cm 3=1000kg/m 3） 3、重力：G=mg （g 一般取10N/kg ） 4、压强：SFp =（液体和气体还可以用gh p ρ=） 5、浮力：)(::::只适合漂浮或悬浮二力平衡法阿基米德原理称重法压力差法浮排排浮示浮向下向上浮G F gV G F F G F F F F ===-=-=ρ6、功：GhW Fs W ==::竖直方向水平方向7、功率：Fv tWP ==8、机械效率：%100⨯=总有用W W η9、杠杆原理：阻阻动动L F L F =常规题：1、正常情况下，冰的密度为33/109.0m kg ⨯，有一块体积为34.0m 的冰块，请求出： （1）这块冰的质量；（2）如果这块冰全部化为水，水的体积是多少3m2有一瓶子，最多能装4kg 煤油，它最多能装多少千克的水？ 还有一个瓶子能装4kg 水，它最多能装多少千克煤油？ （ρ煤油=0.8×103kg/m 3 ρ水=1.0×103kg/m 3)3、刘忠源的公司准备购进20t 煤油，根据规格，已知每个油桶的容积为318.0m ，他又上网查出了煤油的密度为33108.0-⋅⨯m kg ，请问他应准备这种油桶多少个?4、小明在班里坐在最后面，最近老感觉看不太清楚黑板上的粉笔字，估计自己可能是近视了，于是准备去眼镜店配一副眼镜。
根据八年级学到的光学知识，要矫正近视眼，需要配凹透镜。
出于好奇，他上网查到了组成眼镜主要材料的部分技术指标如下表： （1）求一块体积为35cm 的玻璃镜片的质量？（2）一副铜合金镜架的质量20g ，请求出这副镜架的体积？ （3）若以钛合金代替上一问的铜合金，可以减轻多少质量？5、一个质量为9千克的球，体积是1.4分米3，它是空心的还是实心的?如果是空的，空心部分体积多大?(ρ＝7.2×103千克／米3)6、一个边长为5cm 的正方体，重为10N 。
东南大学计算力学课件(研究生课程) 第10章动力学问题
2、动态分析模型的确定
3) 模型应正确反映结构的实际特性。一个具体结构
的动态特性，主要取决于质量、刚度的大小与分布，取决 于结构边界条件与阻尼特性。因此，模型应尽量保持整个 结构、甚至各构件的质量、质心位置不变，保证构造的刚 度特性和传力路线基本不变。真实反映分析对象的边界条 件与阻尼特性。因此，选取的单元应保持其几何形状、受 力特性、变形特性与实际结构的几何形状特点、受力传力 特点、变形特性相一致。
N T
N dV
(t
e
)
M
e
(t
e
)
V
建立此质量矩阵[M]所用的位移插值函数与刚度矩阵一致 ，故称为一致质量矩阵。一致质量矩阵总是正定满秩的 。
10
2、动态分析模型的确定
单元质量矩阵
[M ]e N T N dV V
平面常应变三角形单元的一致质量阵为：
负号表示惯性力与加速度相反。
显然，整个单元上惯性力即为上式的积分。如何将这个作用于单元上 的惯性力移置到单元节点上，通常有两种方法：
1）虚功原理法——求得一致质量矩阵
2）直接分配法——即按重心不变原则分配，求得集中质量矩阵。
9
2、动态分析模型的确定
1）虚功原理法
设单元中发生虚位移为 f *(t) *(t) e
， (t) e 、 (t) e 则单元节点内任一点的位移
f (t) N (t)e
[N]为形函数，与时间t无关，为x、y、z的函数，它与静力分析 中一样；由于[N]与时间无关，则单元应变矩阵，应力矩阵仍与 静力分析完全相同：
(t) B (t)e (t) D (t) DB (t)e
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
8 l
9
可以使用线性单元的 简化方法
东南大学土木工程学院
q
1 2 3
单元①荷载作用在1-2边上，故等 效结点力只与1、2号结点有关， 形函数
①
4 5
②
6
N1 (1 )
N 2
7 l
8 l
9
0 线性分布面力 q q
{FS } [ N ] {qS }hds
l
η
1 η=1 2
ξ=1
1 4 5
ξ
东南大学土木工程学院
单元②形函数
2 2 s 3
s N2 l
N3 1
s l
N6 0
6
s 0 1 l T Fs N qs ds s l l 0 1 l
0 q s s q l
7
3
三角形单元族插值 函数构造及变节点 数法联合运用
东南大学土木工程学院
1
1 x y xy x2y x3y x2y2 xy2 xy3 y2
完全一次 完全二次
9
4 5 10
x2
8
x3
3
y3 完全三次 y4
2
x4
6
7
•Pascal三角形
u 1 2 x 3 y 4 x 2 5 xy 6 y 2 7 x 3 8 x 2 y 9 xy 2 10 y 3 2 2 3 2 2 3 v x y x xy y x x y xy y 11 12 13 14 15 16 17 18 19 20
1
2
3
4
5
6
7 l
8 l
9
平衡所用到的的量均要属于节点的量，如单元位移、单元力。 载荷亦应如此，必须将体积力、表面力转化到节点上去，成 为等效荷载（载荷）。
东南大学土木工程学院
把总等效荷载 R e分解成体积力、表面力和集中力的等效荷 载之和，有 {R}e {FV } {FS } { pc } FV——单元上体积力的等效荷载 FS——单元上表面力的等效荷载 pC——单元上节点上的集中力 1、体积力的等效荷载
N5 4 ( 1), N6 4 ( 1)(1 )
b
q
6
8
N7 4 ( 1)(1 ), N8 4 ( 1)
x
3
7
4
东南大学土木工程学院
（2）荷载作用在2-3边上，故等效结点力只与2、6、3号结点 有关
N 2 ( 1)(2 2 +1), N 6 4 ( 1)(1 ), N 3 (1 )( 1)(2 2 1)
y a 2 5 1
q
6
8
b
3
7
4
x
东南大学土木工程学院
y a 2 5 1
1 x x2 x3
8
b
y xy y2 xy2 x2y2 xy3
完全一次 完全二次
q
6
x2y x3y
y3 完全三次 y4
x4
3 7 4 x
•Pascal三角形
2 2 2 2 u x y x xy y x y xy 1 2 3 4 5 6 7 8 2 2 2 2 v x y x xy y x y xy 9 10 11 12 13 14 15 16
东南大学土木工程学院
2 2 2 2 u 1 2 x 3 y 4 x 5 xy 6 y 7 x y 8 xy 2 2 2 2 v x y x xy y x y xy 9 10 11 12 13 14 15 16
y a 2 5 1
在 0边上计算
q
6
8
3
7
4
x
Ni N 6 N 3 N 2 4 1, 4(1 2 ), 4 3 x 0, N b N N y 0 3 6 b 2 2, 2
l T
η
1 η=1 2
1 4 5
ds ld
ξ
ξ=1
东南大学土木工程学院
q
1 2 3
在 1边上
②
N1 1
N2
1 2
①
4 5
6
ql P2 y N 2 q y ds ql d 3 l 0
7 l
8 l
9
P2 x N 2 qx ds 0
东南大学土木工程学院
x y 局部坐标 , a b
N1 (2 +2 3), N 2 ( 1)(2 2 +1) N 3 (1 )( 1)(2 2 1), N 4 (2 -2 +1)( 1)
y a 2 5 1
s l 0 0 0 0 s 0 0 l
T
0 T ql 1 2 s ds 0 0 0 0 q 2 3 3 l
0 R ql 2 故结点2的等效荷载列阵为 2 3
东南大学土木工程学院
3，图中两个三角形单元组成平行四边形，已知单元①按局部 编码i , j , m 的单元刚度矩阵和应力矩阵是
{FV } [ N ]T {qV }hdxdy
A
2、表面力的等效荷载
{FS } [ N ] {qS }hds
l
T
东南大学土木工程学院
q
3，某平面结构采用四结点矩形 单元和三结点三角形单元建立有 限元计算模型，其如图所示。试 求结点2的等效荷载列阵{R2}。
1
2
3
①
4 5
②
6
7 l
x1 N1 N 2 N 3 N 4 x2 J , , , 1 x3 x4
1 1 1 1
1 2 (8 12 1) 6 4 (3 2 1) 1 3 1 2 [k ] E (8 12 1) 4 6 1 (3 2 1) 3 1 2 (8 12 -1) 6
东南大学土木工程学院
（1）构造不考虑边结点和内部结点的角结点的插值函数：
ˆ L ,N ˆ L ,N ˆ L N 1 1 2 2 3 3
N1 N 2 N 3 N 4
1 (8 12 2 1), 6 4 (3 2 1), 3 4 (3 2 1), 3 1 (8 12 2 1) 6
东南大学土木工程学院
(b)一维问题中，单元刚度矩阵
[k ] [ B]T [ D][ B] J d E [ B]T [ B] J d
4 (3 2 1) 3
4 (3 2 1) 3
1 (8 12 2 -1) 1d 6
整体坐标非自然坐标结果如何？
东南大学土木工程学院
1
8，试利用变结点数法构造 插值函数，构造出如图所 示的三次三角形单元的形 函数及相应的位移函数。
2
4 5 10
9 8
6
6 2 6 8 0 6 16 6 12 6 4 13.5 9 7.5 3 13.5 3 1.5 对 称 9.5 3 5.5
i i 1 j m
m 2 j
K (1)
S (1)
0 3 0 0 0 3 0 4 0 3 0 1 2 0 1.5 1.5 0.5 1.5
计算力学习题讲解
东南大学土木工程学院 2018.1 东南大学土木工程学院
1，试述弹性力学中按位移求解与有限元单元法中按位移求解 之间的异同点。
弹性力学
物理模型 基本方程 解法 连续体 几何方程 物理方程 平衡微分方程 解微分方程
有限单元法
离散化结构 几何方程 物理方程 结点平衡方程 解代数方程
解答形式
b c
a
东南大学土木工程学院
（1）考虑每个结点有两个自由度，半带宽d=（相邻结点 码的最大差值+1）*节点自由度数
bb
5 (4) 4 (2) 2 (1) 1 (3) 3 (7) 6 (11) (8) (5) 7 (9) (6)
cc
8
(10)
11 (12) 10
aa
9
故结点编号如图所示可使单元内结点编码相差3，使 得带宽d=12(按每节点3自由度假定）。
x
1
1
q
6
8
b
3
7
4
P2 x t N 2 qx ds 2qt N 2 d
s 0
1
qt 3
东南大学土木工程学院
5，如图所示刚架 （1）如何进行结点编号使整体刚度矩阵[K]的带宽最小？ （2）刚架的整体刚度矩阵中a结点的总体刚度矩阵Kaa和总 体刚度矩阵Kbc各由哪些分块矩阵叠加组成（自行确定单元 局部坐标方向）
按图中所示单元②的局部 编码写出单元刚度矩阵和 应力矩阵。
东南大学土木工程学院
东南大学土木工程学院
平面问题的单元刚度矩阵[k]不随单元（或坐标轴）的平行 移动或作n角度（n为整数）的转动而改变。S阵变号
东南大学土木工程=b=2，单元厚度为t。 （1）求该单元的位移函数和形函数，并检验其是否满足收 敛性条件。 （2）求在2-6-3边作用均布水平荷载q时的等效结点荷载。
东南大学土木工程学院
s
-1 1 -1/2 2 1/2 3 1 4
(a)由拉格朗日插值函数可知：