(完整版)有限元第二章课后题答案
机械有限元习题答案——哈工大
机械有限元习题答案——哈⼯⼤第⼆章习题2.1 解释如下的概念：应⼒、应变，⼏何⽅程、物理⽅程、虚位移原理。
解○1应⼒是某截⾯上的应⼒在该处的集度。
○2 应变是指单元体在某⼀个⽅向上有⼀个ΔU 的伸长量，其相对变化量就是应变。
X U Xx ??=ε表⽰在x 轴的⽅向上的正应变，其包括正应变和剪应变。
○3⼏何⽅程是表⽰弹性体内节点的应变分量与位移分量之间的关系，其完整表⽰如下：Txz yz xy z y x x w z u zv y w y u x v z w y vx u x w z u z v y w y u x v z w y v x u ??+++=+????+????+?=??????=γγγεεεε○4物理⽅程：表⽰应⼒和应变关系的⽅程某⼀点应⼒分量与应变分量之间的关系如下：=???????????????????=666564636261565554535251464545434241363534333231262524232221161514131211αααααααααααααααααααααααααααααααααααατττσσσσxz yz xy z y xxz yz xy zz yy xx γγγεεε○5虚位移原理：在弹性有⼀虚位移情况下，由于作⽤在每个质点上的⼒系，在相应的虚位移上虚功总和为零，即为：若弹性体在已知的⾯⼒和体⼒的作⽤下处于平衡状态，那么使弹性体产⽣虚位移，所有作⽤在弹性体上的体⼒在虚位移上所做的⼯就等于弹性体所具有的虚位能。
2.2说明弹性体⼒学中的⼏个基本假设。
○1 连续性假设：就是假定整个物体的体积都被组成该物体的介质所填满，不存在任何间隙。
○2 完全弹性假设：就是假定物体服从虎克定律。
○3 各向同性假设：就是假定整个物体是由同意材料组成的。
○4 ⼩变形和⼩位移假设：就是指物体各点的位移都远远⼩于物体原来的尺⼨，并且其应变和转⾓都⼩于1。
2.3简述线应变与剪应变的⼏何含义。
高等有限元课后题答案
⾼等有限元课后题答案2 弹性⼒学问题的有限单元法思考题2.1 有限元法离散结构时为什么要在应⼒变化复杂的地⽅采⽤较密⽹格，⽽在其他地⽅采⽤较稀疏⽹格？答：在应⼒变化复杂的地⽅每⼀结点与相邻结点的应⼒都变化较⼤，若⽹格划分较稀疏，则在应⼒突变处没有设置结点，⽽使得所求解的误差很⼤，若⽹格划分较密时，则应⼒变化复杂的地⽅可以设置更多的结点，从⽽使得所求解的精度更⾼⼀些。
2.2 因为应⼒边界条件就是边界上的平衡⽅程，所以引⽤虚功原理必然满⾜应⼒边界条件，对吗？答：对。
2.3 为什么有限元只能求解位移边值问题和混合边值问题？弹性⼒学中受内压和外压作⽤的圆环能⽤有限元⽅法求解吗？为什么？答：有限元法是⼀种位移解法，故只能求解位移边值问题和混合边值问题。
⽽应⼒边值问题没有确定的位移约束，不能⽤位移法求解，所以也不能⽤有限元法求解。
2.4 矩形单元旋转⼀个⾓度后还能够保持在单元边界上的位移协调吗？答：能。
矩形单元的插值函数满⾜单元内部和单元边界上的连续性要求，是⼀个协调元。
矩形的插值函数只与坐标差有关，旋转⼀个⾓度后各个结点的坐标差保持不变，所以插值函数保持不变。
因此矩形单元旋转⼀个⾓度后还能够保持在单元边界上的位移协调。
2.5 总体刚度矩阵呈带状分布，与哪些因素有关？如何计算半带宽？答：因素：总体刚度矩阵呈带状分布与单元内最⼤结点号与最⼩结点号的差有关。
计算：设半带宽为B ，每个结点的⾃由度为n ，各单元中结点整体码的最⼤差值为D ，则B=n(D+1)，在平⾯问题中n=2。
2.6为什么单元尺⼨不要相差太⼤，如果这样，会导致什么结果？答：由于实际⼯程是⼀个⼆维或三维的连续体，将其分为具有简单⽽规则的⼏何单元，这样便于⽹格计算，还可以通过增加结点数提⾼单元精度。
在⼏何形状上等于或近似与原来形状，减⼩由于形状差异过⼤带来的误差。
若形状相差过⼤，使结构应⼒分析困难加⼤，误差同时也加⼤。
2.7 剖分⽹格时，在边界出现突变和有集中⼒作⽤的地⽅要设置结点或单元边界，试说明理由。
有限单元法课后习题全部答案_王勖成
∫
∂ 2φ ∂ 2φ ∂φ ∂φ k − ∫ k 2 + k 2 + Q δφ d Ω + ∫ δφ d Ω − ∫ αφ − q − k δφ d Γ Ω Γ − Γ Γ q q ∂y ∂n ∂n ∂x
欧拉方程： k
∂ 2φ ∂ 2φ + +Q = k 0 ∂x 2 ∂y 2
习题 1.2： 在用有限元法求解时，边界条件总是满足的，控制方程的不完全匹配，会产生误差。题中所 ，代入边 给出的近似函数： φ =a0 + a1 x + a2 x + a3 x ，应该满足边界条件，对于情况（1）
2 3
界条件可得 = a0 0, = a3
1 − a1 L − a2 L2 ，从而 L3 x3 x3 x3 2 ) + a ( x − ) + 2 L2 L L3
∫
= =
∑{ A
m k =1 m
T
( N j ( xk )) [ A( N i ( xk )ai ) − f ( xk )]
m
}
( N j )A( N i )ai − ∑ AT ( N j ) f = k 1= k 1
T
∑A
= Ka-P
(写成矩阵形式)
因此， kij =
d 2 w dw d 3 w 0 dx 2 δ dx − dx3 δ w = 0
L
1.5 如有一问题的泛函为 = Π ( w)
∫
L
0
EI d 2 w 2 kw2 + qwdx ，其中 E, I, k 是常数，q 2 + 2 dx 2
最新有限元法课后习题答案
1、有限元是近似求解一般连续场问题的数值方法2、有限元法将连续的求解域离散为若干个子域，得到有限个单元，单元和单元Z间用节点连接3、直梁在外力的作用下，横截面的内力有剪力和弯矩两个.4、平面刚架结构在外力的作用下，横截面上的内力有轴力、剪力、弯矩.5、进行直梁有限元分析，平面刚架单元上每个节点的节点位移为挠度和转角6、平面刚架有限元分析，节点位移有轴向位移、横向位移、转角。
7、在弹性和小变形下，节点力和节点位移关系是线性关系。
8、弹性力学问题的方程个数有15个，未知量个数有15个。
9、弹性力学平面问题方程个数有8,未知数8个。
10、几何方程是研究应变和位移之间关系的方程11、物理方程是描述应力和应变关系的方程12、平衡方程反映了应力和体力之间关系的13、把经过物体内任意一点各个截面上的应力状况叫做一点的应力状态14、9形函数在单元上节点上的值，具有本点为JL_•它点为零的性质，并且在三角形单元的任一节点上，三个行函数之和为_1_15、形函数是—三角形一单元内部坐标的—线性一函数，他反映了单元的—位移—状态16、在进行节点编号时，同一单元的相邻节点的号码差尽量小.17、三角形单元的位移模式为—线性位移模式」18、矩形单元的位移模式为—双线性位移模式_19、在选择多项式位移模式的阶次时，要求—所选的位移模式应该与局部坐标系的方位无关的性质为几何一各向同性20、单元刚度矩阵描述了—节点力_和_节点位移之间的关系21、矩形单元边界上位移是连续变化的1.诉述有限元法的定义答：有限元法是近似求解一般连续场问题的数值方法2.有限元法的基本思想是什么答：首先，将表示结构的连续离散为若干个子域，单元之间通过其边界上的节点连接成组合体。
其次，用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。
3.有限元法的分类和基本步骤有哪些答：分类：位移法、力法、混合法；步骤：结构的离散化，单元分析，单元集成，引入约朿条件，求解线性方程组，得出节点位移。
机械结构有限元分析第二章课后答案 哈工大
εz =
∂w =0 ∂z
γ xy =
∂u ∂v + =0 ∂y ∂x
γ yz =
∂v ∂w + =0 ∂z ∂y
γ zx =
∂u ∂w + = 12 × 10 2 ∂z ∂x
2.9 一具有平面应力场的物体，材料参数为 E、v。有如下位移场
u (x, y ) = ax 3 − bxy 2
v(x, y ) = cx 2 y − dy 3
εz
γ xy
γ yz
γ zx ]
T
式中， D —弹性矩阵，是一个常数矩阵。 虚位移原理：一个弹性体在外力和内力作用下处于平衡状态，则对于任何约束允许的虚位移来 说，外力所做的虚功等于内力的虚功。
2.2 说明弹性力学中的几个基本假设。
答：弹性力学中的几个基本假设有：
（1）连续性假定，指假定整个物体的体积都被组成该物体的介质所填满，不存在任何空隙。 （2）完全弹性假定，指假定物体服从胡克定律，即应变与引起该应变的应力成正比。 （3）均匀性假定，指假定整个物体是由同一材料组成的。
2 ∂ 2ε x ∂ ε y ∂ 3u ∂ 3v + = + = 2c − 2b ∂y 2 ∂x 2 ∂x∂y 2 ∂y∂x 2
∂ 2 γ xy ∂x∂y
= 2c − 2b
2 2 ∂ 2 ε x ∂ ε y ∂ γ xy 所以满足 的相容方程 + = ∂x∂y ∂y 2 ∂x 2
2.10 一具有平面应力场的物体，材料参数为 E、v。有如下位移场
σy =
⎛ ∂v ∂u ⎞ 4v ⎡⎛ 3aE Eb ⎞ 2 ⎛ ⎞ 2⎤ ⎜ ⎜ ∂y + µ ∂x ⎟ ⎟ = 4v − E ⎢⎜ c + 2v − 3a ⎟ x − ⎜ 3d + 2v − b ⎟ y ⎥ ⎠ ⎝ ⎠ ⎦ ⎣⎝ ⎝ ⎠
有限元习题及答案ppt课件
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
有限元习题及答案
有限元习题及答案有限元习题及答案有限元方法是一种常用的数值计算方法，用于求解各种工程和科学问题。
在学习有限元方法的过程中，练习习题是非常重要的，可以帮助学生巩固所学的知识，并提高解决实际问题的能力。
本文将介绍一些有限元习题及其答案，希望对学习有限元方法的同学有所帮助。
习题一：一维热传导问题考虑一个长度为L的一维杆，其两端固定，杆上的温度满足以下热传导方程：∂²T/∂x² = 0，其中T为温度，x为位置。
已知杆的两端温度分别为T1和T2，求解杆上的温度分布。
解答一：根据热传导方程，可以得到温度分布的一般解为T(x) = Ax + B，其中A和B为常数。
根据边界条件，可以得到方程组：T(0) = B = T1T(L) = AL + B = T2解方程组可得A = (T2 - T1) / L，B = T1。
因此，温度分布为T(x) = ((T2 - T1) / L) * x + T1。
习题二：二维弹性问题考虑一个矩形薄板，其长为L，宽为W，材料的弹性模量为E，泊松比为ν。
已知薄板的边界上施加了一定的边界条件，求解薄板上的位移场。
解答二：对于二维弹性问题，可以使用平面应力假设，即假设薄板内部的应力只有两个分量σx和σy，并且与z轴无关。
根据平面应力假设和胡克定律，可以得到位移场的偏微分方程：∂²u/∂x² + ν * (∂²u/∂y²) + (1 - ν) * (∂²v/∂x∂y) = 0∂²v/∂y² + ν * (∂²v/∂x²) + (1 - ν) * (∂²u/∂x∂y) = 0其中u和v分别为位移场在x和y方向上的分量。
边界条件根据具体情况给定。
通过数值方法，如有限元方法，可以求解位移场的近似解。
习题三：三维流体力学问题考虑一个三维流体力学问题，流体在一个封闭容器内流动，容器的形状为一个长方体，已知流体的速度场和压力场的初始条件，求解流体的运动状态。
有限元分析基础课后习题答案ppt课件
2
0
2(1
1
)
1 2 1
0
0
0
1
E 8
4 2 0
2 4 0
0 0 1
4
3
24
4 2 0
D
E 8
2
4
0
0 0 1
1 0 0 0 1 0
B 0 0 0 1 0 1
0 1 1 0 1 1
K e BT DBtA
1 0 0 0 1
0
T
4
2
01
0
0
0
1
0
E 8
0
0
0
1
0
1 1 4 00 0 0 1
ppt精选版19etppt精选版20ppt精选版21ppt精选版22ppt精选版23ppt精选版24ppt精选版25ppt精选版26ppt精选版27ppt精选版28ppt精选版29ppt精选版30
三角形单元i,j,m的j m边作用有 如图所示线形分布面载荷，求 结点载荷向量。 解：面力移置公式：
Re NT ptds
0
v1
0
3
Et
1
20 1 1
0
0
1 1 1 3 0 0
4 0 0 2 1 1
0 4 2 0 1 3
0 2 4 0 3 1
2 0 0 4 1 1
1 1 3 1 4 2
421113uuuvvv432432
0
0
0
P
0
0
Et 4 20 0
0 4
uv33
a
am a * a 0*0 a2,bm 0 a a, cm 0 a a
N
Ni 0
有限单元法部分课后题答案
1.1 有限单元法中“离散”的含义是什么？有限单元法是如何将具有无限自由度的连续介质问题转变成有限自由度问题的？位移有限元法的标准化程式是怎样的？（1）离散的含义即将结构离散化，即用假想的线或面将连续体分割成数目有限的单元，并在其上设定有限个节点；用这些单元组成的单元集合体代替原来的连续体，而场函数的节点值将成为问题的基本未知量。
（2）给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律，即通过插值以单元节点位移表示单元内任意点的位移。
因节点位移个数是有限的，故无限自由度问题被转变成了有限自由度问题。
（3）有限元法的标准化程式：结构或区域离散，单元分析，整体分析，数值求解。
1.3 单元刚度矩阵和整体刚度矩阵各有哪些性质？各自的物理意义是什么?两者有何区别?单元刚度矩阵的性质：对称性、奇异性（单元刚度矩阵的行列式为零）。
整体刚度矩阵的性质：对称性、奇异性、稀疏性。
单元 Kij 物理意义 Kij 即单元节点位移向量中第 j 个自由度发生单位位移而其他位移分量为零时，在第 j 个自由度方向引起的节点力。
整体刚度矩阵 K 中每一列元素的物理意义是：要迫使结构的某节点位移自由度发生单位位移，而其他节点位移都保持为零的变形状态，在所有个节点上需要施加的节点荷载。
2.2 什么叫应变能？什么叫外力势能？试叙述势能变分原理和最小势能原理，并回答下述问题：势能变分原理代表什么控制方程和边界条件？其中附加了哪些条件？(1)在外力作用下，物体内部将产生应力σ和应变ε，外力所做的功将以变形能的形式储存起来，这种能量称为应变能。
(2)外力势能就是外力功的负值。
(3)势能变分原理可叙述如下：在所有满足边界条件的协调位移中，那些满足静力平衡条件的位移使物体势能泛函取驻值，即势能的变分为零δ∏p=δ Uε+δV=0此即变分方程。
对于线性弹性体，势能取最小值，即δ2∏P=δ2Uε+δ2V≥0此时的势能变分原理就是著名的最小势能原理。
完整版有限元法课后习题答案
1、有限元是近似求解一般连续场问题的数值方法2、有限元法将连续的求解域离散为假设干个子域,得到有限个单元,单元和单元之间用节点连接3、直梁在外力的作用下,横截面的内力有剪力和弯矩两个.4、平面刚架结构在外力的作用下横截面上的内力有轴力、剪力、弯矩.5、进行直梁有限元分析,平面刚架单元上每个节点的节点位移为挠度和转角6、平面刚架有限元分析,节点位移有轴向位移、横向位移、转角 .7、在弹性和小变形下,节点力和节点位移关系是线性关系.8、弹性力学问题的方程个数有15个,未知量个数有15个.9、弹性力学平面问题方程个数有8,未知数8个.10、几何方程是研究应变和位移之间关系的方程11、物理方程是描述应力和应变关系的方程12、平衡方程反映了应力和体力之间关系的13、把经过物体内任意一点各个截面上的应力状况叫做一点的应力状态14、9形函数在单元上节点上的值 ,具有本点为_1_.它点为零的性质,并且在三角形单元的任一节点上,三个行函数之和为_1_15、形函数是三角形单元内部坐标的线性函数他反映了单元的位移状态16、在进行节点编号时,同一单元的相邻节点的号差尽量小.17、三角形单元的位移模式为_线性位移模式_-18、矩形单元的位移模式为双线性位移模式19、在选择多项式位移模式的阶次时,要求_所选的位移模式应该与局部坐标系的方位无关的性质为几何各向同性20、单元刚度矩阵描述了节点力和节点位移之间的关系21、矩形单元边界上位移是连续变化的1.诉述有限元法的定义答：有限元法是近似求解一般连续场问题的数值方法2.有限元法的根本思想是什么答：首先,将表示结构的连续离散为假设干个子域,单元之间通过其边界上的节点连接成组合体.其次,用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量.3.有限元法的分类和根本步骤有哪些答：分类：位移法、力法、混合法；步骤：结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移.4.有限元法有哪些优缺点答：优点：有限元法可以模拟各种几何形状复杂的结构,得出其近似解；通过计算机程序,可以广泛地应用于各种场合；可以从其他CAD软件中导入建好的模型；数学处理比较方便, 对复杂形状的结构也能适用；有限元法和优化设计方法相结合,以便发挥各自的优点.缺点：有限元计算,尤其是复杂问题的分析计算, 所消耗的计算时间、内存和磁盘空间等计算资源是相当惊人的. 对无限求解域问题没有较好的处理方法. 尽管现有的有限元软件多数使用了网络自适应技术, 但在具体应用时,采用什么类型的单元、多大的网络密度等都要完全依赖适用者的经验.5.梁单元和平面钢架结构单元的自由度由什么确定答：由每个节点位移分量的总和确定6.简述单元刚度矩阵的性质和矩阵元素的物理意义答：单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵单元刚度矩阵中元素aml的物理意义为单元第L个节点位移分量等于1,其他节点位移分量等于0时,对应的第m个节点力分量.7.有限元法根本方程中的每一项的意义是什么P14答：Q——整个结构的节点载荷列阵〔外载荷、约束力〕；整个结构的节点位移列阵；结构的整体刚度矩阵,又称总刚度矩阵.8.位移边界条件和载荷边界条件的意义是什么答：由于刚度矩阵的线性相关性不能得到解,引入边界条件,使整体刚度矩阵求的唯一解.9.简述整体刚度矩阵的性质和特点P14答：对称性；奇异性；稀疏性；对角线上的元素恒为正.10简述整体坐标的概念P25答：在整体结构上建立的坐标系叫做整体坐标,又叫做统一坐标系.11.简述平面钢架问题有限元法的根本过程答：1〕力学模型确实定,2〕结构的离散化,3〕计算载荷的等效节点力,4〕计算各单元的刚度矩阵,5〕组集整体刚度矩阵,6〕施加边界约束条件,7〕求解降价的有限元根本方程, 8〕求解单元应力,9〕计算结果的输出.12.弹性力学的根本假设是什么.答：连续性假定,弹性假定,均匀性和各向同性假定,小变形假定,无初应力假定.13.弹性力学和材料力学相比,其研究方法和对象有什么不同.答：研究对象：材料力学主要研究杆件,如柱体、梁和轴,在拉压、剪切、弯曲和扭转等作用下的应力、形变和位移.弹性力学研究各种形状的弹性体,除杆件外,还研究平面体、空间体,板和壳等.因此,弹性力学的研究对象要广泛得多.研究方法：弹性力学和材料力学既有相似之外,又有一定区别.弹性力学研究问题,在弹性体区域内必须严格考虑静力学、几何学和物理学三方面条件,在边界上严格考虑受力条件或约束条件,由此建立微分方程和边界条件进行求解,得出较精确的解答.而材料力学虽然也考虑这几方面的条件,但不是十分严格的,材料力学只研究和适用于杆件问题. 14.简述圣维南原理. 答；把物体一小局部上的面力变换为分布不同但静力等效的面力,但影响近处的应力分量, 而不影响远处的应力.“局部影响原理〞15.平面应力问题和平面应变问题的特点和区别各是什么试各举出一个典型平面应力和平面应变的问题的实例.答：平面应力问题的特点：长、宽尺寸远大于厚度,沿板面受有平行板的面力,且沿厚度均匀分布,体力平行于板面且不沿厚度变化,在平板的前后外表上无外力作用平面应变问题的特点：Z向尺寸远大于x、y向尺寸,且与z轴垂直的各个横截面尺寸都相同,受有平行于横截面且不沿z向变化的外载荷,约束条件沿z向也不变,即所有内在因素的外来作用都不沿长度变化.区别：平面应力问题中z方向上应力为零,平面应变问题中z方向上应变为零、应力不为零.举例：平面应力问题等厚度薄板状弹性体,受力方向沿板面方向,荷载不沿板的厚度方向变化,且板的外表无荷载作用.平面应变问题一一水坝用于很长的等截面四柱体,其上作用的载荷均平行于横截面,且沿柱长方向不变法.16.三角形常应变单元的特点是什么矩形单元的特点是什么写出它们的位移模式.答：三角形单元具有适应性强的优点,较容易进行网络划分和逼近边界形状,应用比较灵活.其缺点是它的位移模式是线性函数,单元应力和应变都是常数,精度不够理想.矩形单元的位移模式是双线性函数,单元的应力、应变式线性变化的,具有精度较高, 形状规整,便于实现计算机自动划分等优点,缺点是单元不能适应曲线边界和斜边界,也不能随意改变大小,适用性非常有限.17.写出单元刚度矩阵表达式、并说明单元刚度与哪些因素有关.答：单元刚度矩阵与节点力坐标变换矩阵,局部坐标系下的单元刚度矩阵,节点位移有关的坐标变换矩阵.18.如何由单元刚度矩阵组建整体刚度矩阵〔叠加法〕答：〔1〕把单元刚度矩阵扩展成单元奉献矩阵 ,把单元刚度矩阵中的子块按其在整体刚度矩阵中的位置排列, 空白处用零子块填充.〔2〕把单元的奉献矩阵的对应列的子块相叠加, 即可得出整体刚度矩阵 .19.整体刚度矩阵的性质.答：〔1〕整体刚度矩阵中每一列元素的物理意义为：欲使弹性体的某一节点沿坐标方形发生单位为移,而其他节点都保持为零的变形状态,在各节点上所需要施加的节点力；〔2〕整体刚度矩阵中的主对角元素总是正的；〔3〕整体刚度矩阵是一个对称阵；〔4〕整体刚度矩阵式一个呈带状分布的稀疏性矩阵.〔5〕整体刚度矩阵式一个奇异阵,在排除刚体位移后,他是正定阵.20.简述形函数的概念和性质.答：形函数的性质有：〔1〕形函数单元节点上的值,具有“本点为一、他点为零〞的性质；〔2〕在单元的任一节点上,三角函数之和等于1; 〔3〕三角形单元任一一条边上的形函数,仅与该端点节点坐标有关,而与另外一个节点坐标无关；〔4〕型函数的值在0〜1之间变换.21.结构的网格划分应注意哪些问题 .如何对其进行节点编号.才能使半带宽最小.P50, P8相邻节点的号差最小答：一般首选三角形单元或等参元.对平直边界可选用矩形单元,也可以同时选用两种或两种以上的单元.一般来说,集中力,集中力偶,分布在和强度的突变点,分布载荷与自由边界的分界点,支撑点都应该取为节点,相邻节点的号差尽可能最小才能使半带宽最小22.为了保证解答的收敛性,单元位数模式必须满足什么条件答：〔1〕位移模式必须包含单元刚体位移；〔2〕位移模式必须包含单元的常应变；〔3〕位移模式在单元内要连续,且唯一在相邻单元之间要协调.在有限单元法中,把能够满足条件1和条件2的单元称为完备单元,把满足条件3的单元叫做协调单元或保续单元.23有限元分析求得的位移解收敛于真实解得下界的条件.答：1.位移模式必须包含单元的刚体位移,2.位移模式必须包含单元的常应变,3.位移模式在单元内要连续,且位移在相邻单元之间要协调.24.简述等参数单元的概念.答：坐标变换中采用节点参数的个数等于位移模式中节点参数的个数,这种单元称为等参单元.25.有限元法中等参数单元的主要优点是什么答：1〕应用范围广.在平面或空间连续体,杆系结构和板壳问题中都可应用.2〕将不规那么的单元变化为规那么的单元后,易于构造位移模式.3〕在原结构中可以采用不规那么单元,易于适用边界的形状和改变单元的大小.4〕可以灵活的增减节点,容易构造各种过度单元.5〕推导过程具有通用性.一维,二维三维的推导过程根本相同.26.简述四节点四边形等参数单元的平面问题分析过程.答：〔1〕通过整体坐标系和局部坐标系的映射关系得到四节点四边形等参单元的母单元,并选取单元的唯一模式；〔2〕通过坐标变换和等参元确定平面四节点四边形等参数单元的几何形状和位移模式；〔3〕将四节点四边形等参数单元的位移模式代入平面问题的几何方程,得到单元应变分量的计算式,再将单元应变代入平面问题的物理方程,得到平面四节点等参数单元的应力矩阵〔4〕用虚功原理球的单元刚度矩阵,,最后用高斯积分法计算完成.27.为什么等参数单元要采用自然坐标来表示形函数为什么要引入雅可比矩阵答：简化计算得到形函数的偏导关系.28. ANSYS软件主要包括哪些局部各局部的作用是什么答：1.前处理模块：提供了一个强大的实体建模及网络划分工具,用户可以方便地构造有限元模型.2.分析计算模块：包括结构分析、流体力学分析、磁场分析、声场分析、压电分析以及多种物理场的耦合分析,可以模拟多种物理介质的相互作用,具有灵敏度分析及优化分析水平.3.后处理模块：可将计算后果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示等图形方式显示出来,也可将计算结果以图表、曲线形式显示出来或输出.29. ANSYS软件提供的分析类型有哪些答：结构静力分析、机构动力分析、结构非线性分析、动力学分析、热分析、流体力学分析、电磁场分析、声场分析、压电分析.30.简述ANSYS软件分析静力学问题的根本流程.答：1.前处理器：1〕定义单元类型,2〕定义实常数,3〕定义材料属性,4〕创立实体几何模型,5〕划分网络；2.求解器：1〕定义分析类型,2〕施加载荷和位移约束条件,3〕求解；三角形三节点单元的位移是连续的,应变和应力在单元内是常数,因而其相邻单元将具有不同的应力和应变,即在单元的公共边界上和应变的值将会有突变.矩形单元的边界上,位移是线性变化的,显然,在两个相邻矩形单元的公共边界上,其位移是连续的.节点的选用原那么：一般说,集中力、集中力偶、分布载荷强度的突变点、分布载荷与自由边界的分界点、支承点都能赢取为节点.单元的划分原那么：〔1〕划分单元的数目,视要求的计算精度和计算机的性能而定.〔2〕单元的大小,可根据部位的不同而有所不同.1、试述街节点力和节点载荷的区别.节点力是单元与节点之间的作用力；如果取整个结构为研究对象,节点力为内力,节点载荷是作用在节点上的外载荷.2、试述求整体刚度矩阵的两种方法.分别建立各节点的平衡方程式,写成矩阵形式,可求得整体刚度矩阵；将各单元刚度矩阵按规律叠加,也可得整体刚度矩阵.3、平面问题中划分单元的数目是否越多越好不是越多越好.划分单元的数目,视要求的计算精度和计算机的性能而定.随着单元数目的接连多,有限元解逐步逼近于真实解,但是,单元数目接连加,刚求解的有限元线性方程组的数目接连多, 需要占用更多的计算机内存资源,求解时间接连长,所以,在计算机上进行有限元分析时,还要考虑计算机的性能.单元数过多并不经济.4、写出单元刚度矩阵的表达式,并说明单元刚度与那些因素有关[B]-单元应变矩阵,[D]-弹性矩阵,t-厚度〕单元刚度矩阵取决于单元的大小、方向、和弹性常数,而与单元的位置无关,即不随单元或坐标轴的平移而改变.5、选择多项式为单元的位移模式时,除了要满足单元的完备性和协调性要求,还须考虑什么因素还须考虑两个因素：1、所选的位移模式应该与局部坐标系的方位无关,即几何各向同性. 2、多项式位移模式中的项数必须等于或稍大于单元边界上的外节点的自由度数,通常取多项式的项数与单元的外节点的自由度数想等.。
(完整版)有限元第二章课后题答案
2 弹性力学问题的有限单元法思考题2.1 有限元法离散结构时为什么要在应力变化复杂的地方采用较密网格，而在其他地方采用较稀疏网格？答：在应力变化复杂的地方每一结点与相邻结点的应力都变化较大，若网格划分较稀疏，则在应力突变处没有设置结点，而使得所求解的误差很大，若网格划分较密时，则应力变化复杂的地方可以设置更多的结点，从而使得所求解的精度更高一些。
2.2 因为应力边界条件就是边界上的平衡方程，所以引用虚功原理必然满足应力边界条件，对吗？答：对。
2.3 为什么有限元只能求解位移边值问题和混合边值问题？弹性力学中受内压和外压作用的圆环能用有限元方法求解吗？为什么？答：有限元法是一种位移解法，故只能求解位移边值问题和混合边值问题。
而应力边值问题没有确定的位移约束，不能用位移法求解，所以也不能用有限元法求解。
2.4 矩形单元旋转一个角度后还能够保持在单元边界上的位移协调吗？答：能。
矩形单元的插值函数满足单元内部和单元边界上的连续性要求，是一个协调元。
矩形的插值函数只与坐标差有关，旋转一个角度后各个结点的坐标差保持不变，所以插值函数保持不变。
因此矩形单元旋转一个角度后还能够保持在单元边界上的位移协调。
2.5 总体刚度矩阵呈带状分布，与哪些因素有关？如何计算半带宽？ 答：因素：总体刚度矩阵呈带状分布与单元内最大结点号与最小结点号的差有关。
计算：设半带宽为B ，每个结点的自由度为n ，各单元中结点整体码的最大差值为D ，则B=n(D+1)，在平面问题中n=2。
2.6 为什么单元尺寸不要相差太大，如果这样，会导致什么结果？ 答：由于实际工程是一个二维或三维的连续体，将其分为具有简单而规则的几何单元，这样便于网格计算，还可以通过增加结点数提高单元精度。
在几何形状上等于或近似与原来形状，减小由于形状差异过大带来的误差。
若形状相差过大，使结构应力分析困难加大，误差同时也加大。
2.7 剖分网格时，在边界出现突变和有集中力作用的地方要设置结点或单元边界，试说明理由。
高等有限元课后题答案
2 弹性力学问题的有限单元法思考题2.1 有限元法离散结构时为什么要在应力变化复杂的地方采用较密网格，而在其他地方采用较稀疏网格？答：在应力变化复杂的地方每一结点与相邻结点的应力都变化较大，若网格划分较稀疏，则在应力突变处没有设置结点，而使得所求解的误差很大，若网格划分较密时，则应力变化复杂的地方可以设置更多的结点，从而使得所求解的精度更高一些。
2.2 因为应力边界条件就是边界上的平衡方程，所以引用虚功原理必然满足应力边界条件，对吗？答：对。
2.3 为什么有限元只能求解位移边值问题和混合边值问题？弹性力学中受内压和外压作用的圆环能用有限元方法求解吗？为什么？答：有限元法是一种位移解法，故只能求解位移边值问题和混合边值问题。
而应力边值问题没有确定的位移约束，不能用位移法求解，所以也不能用有限元法求解。
2.4 矩形单元旋转一个角度后还能够保持在单元边界上的位移协调吗？答：能。
矩形单元的插值函数满足单元内部和单元边界上的连续性要求，是一个协调元。
矩形的插值函数只与坐标差有关，旋转一个角度后各个结点的坐标差保持不变，所以插值函数保持不变。
因此矩形单元旋转一个角度后还能够保持在单元边界上的位移协调。
2.5 总体刚度矩阵呈带状分布，与哪些因素有关？如何计算半带宽？答：因素：总体刚度矩阵呈带状分布与单元内最大结点号与最小结点号的差有关。
计算：设半带宽为B ，每个结点的自由度为n ，各单元中结点整体码的最大差值为D ，则B=n(D+1)，在平面问题中n =2。
2.6 为什么单元尺寸不要相差太大，如果这样，会导致什么结果？ 答：由于实际工程是一个二维或三维的连续体，将其分为具有简单而规则的几何单元，这样便于网格计算，还可以通过增加结点数提高单元精度。
《有限元法及其应用》课后习题
《有限元法及其应用》课后习题目录第1章绪论 (3)第2章有限单元法理论基础 (4)第3章杆系结构单元 (5)第4章平面三角形单元 (7)第5章平面四边形等参数单元 (9)第6章常用有限元软件及其在岩土工程中的应用 (10)第1章绪论1-1试说明有限元法解题的基本思路。
1-2试说明用有限元法解题的主要步骤。
1-3有限元法主要有哪些优点？第2章有限单元法理论基础2-1 何为虚功，虚功原理的具体思路是什么？2-2 虚功原理的适用条件有哪些？2-3 位移模式的概念是什么？2-4 如何构造位移模式？2-5 弹性力学问题的求解需要满足哪些条件？第3章 杆系结构单元3-1 推导横截面积为A 的一维桁架结构的单元刚度矩阵。
3-2 图示（见题图3-1）为一平面超静定桁架结构，在载荷P 作用下，求各杆件的轴力。
此结构可看成由14、24、34三个杆单元组成，每个杆单元的两端为杆单元的结点，各结点的水平、铅直位移分别用u 、v 表示。
题图3-1 平面超静定桁架结构a —平面结构；b —单元组成；c —各结点位移3-3 图示（见题图3-2）刚架中，两杆为尺寸相同的等截面杆件，横截面面积为20.5m A =，截面惯性矩为41m 24I =，弹性模量7310kPa E =⨯，求解此结构。
题图3-2 等截面刚架结构第4章平面三角形单元4-1 按位移求解的有限单元法中：（1）应用了哪些弹性力学的基本方程?（2）应力边界条件及位移边界条件是如何反映的？（3）力的平衡条件是如何满足的？（4）变形协调条件是如何满足的？4-2 在有限单元法中，如何应用虚功原理导出单元内的应力和结点力的关系式，并将外荷载静力等效地变换为结点荷载？4-3 为了保证有限单元法解答的收敛性，平面三角形单元位移模式应满足哪些条件？μ=，记杨氏弹性模4-4 题图4-1所示等腰直角三角形单元，设14量为E，厚度为t，求形函数矩阵[]N、应变矩阵[]B、应力矩阵[]S与单元刚度矩阵[]eK。
有限单元法课后习题全部答案-王勖成
Ω
=0
∫ ∫ = Ω w j R dΩ
∂
Ω
∂a
j
AT
(
N i
(
x)ai
)δ
(
x
−
x k
)
[
A(
N i
(
x)ai
)
−
f (x)]dΩ
∫ =
Ω
AT
(N
j (x))δ
(x
−
x k
)
[
A(
Ni
(
x)ai
)
−
f
( x)]d Ω
∑{ } m
=
k =1
AT
(
N
j
(
xk
))
[
A(
N i
(
xk
)ai )
= R( L) 3
0= , R( 2L) 3
0 ，从而可以解出待定系数 a1, a2 。带入（1）式可以得到φ 。
配点法仅考虑了有限个点的局部特性，子域法则要求在有限个子域 Ωi 内残量的积分
∫ R(x)dx = 0 为零，子域的个数仍然取决于未知函数个数，通常选取各子域的并集为整个 Ωi
待求区域，一般情况可以选择各子域大小相同，但对于某些局部变化较复杂的区域，可以缩 小子域的大小，使得子域分布更合理。例如取子域为
−
f
( xk
)]
m
m
∑ ∑ = =k
1= AT (N j )A (N i )ai − k
1
AT
(N
j)f
= Ka-P
(写成矩阵形式)
m
∑ = 因此， kij k= =1 AT (N j )A (N i ) k ji ， 系数矩阵对称，且无需积分。
有限元思考题答案
红字为答疑时老师给的解答第一章思考题1-1“用加权余量法求解微分方程，其权函数V和场函数u的选择没有任何限制”，这种说法对吗？答：不对，有连续性要求。
1-2“加权余量法仅适用为传热学问题建立基本的有限元方程，而基于最小势能原理的虚功原理仅适合为弹性力学问题建立基本的有限元方程”，这种说法对吗？答：不对。
虚位移原理不仅可以应用于弹性力学问题，还可以应用于非线性弹性以及弹塑性等非线性问题，虚功原理可以用来推导各种力学问题的有限元基本方法中的基本方程。
最小势能原理仅适用于弹性力学问题。
加权残值法尤其适用于具有连续场的非力学问题，如声、电、磁学的有限元方程的建立。
1-3现代工程分析中的数值分析方法主要有有限差分法、有限元法和边界元法。
这些方法本质上是将求解区域进行网格离散化，然后求解方程获得数值结果。
是否可以将求解区域离散成结点群，但是没有网格进行求解？答：可以，无网格方法是近年发展起来的一种新的数值计算方法。
与基于网格的方法不同，无网格方法只需要节点的信息，不需要节点的信息而不需要节点之间相互联系的信息。
典型无网格方法有配点法、Galerkin方法、Petrov-Galerkin方法等。
（无网格方法数值求解的基本思想：在每个节点上构建待求物理量近似值的插值函数，并用加权残量法和该近似函数对微分方程进行离散，形成与待求物理量相关的各节点近似值的离散方程，并求解之。
）第二章思考题2-1ANSYS软件有哪些模块？在GUI方式下的六个窗口有何功能特点？主要包括前处理模块，分析计算模块和后处理模块①前处理模块提供了一个强大的试题建模及网格划分工具，用户可以方便地构造有限元②分析计算模块包括结构分析、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力③后处理可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示等图形方式显示出来，也可将计算记过以图表、曲线形式显示或输出。
《有限单元法》1-5章课后习题答案
2
dw dw
′′′
w aL2ax , wa ?2
2
dx dxL
EI k?
2 22 2
∏w 4a+ a x x? L + qax x? L dx
∫?
0
22
25 3
ka L qaL
2
2EILa+60 6
53 2
?∏kaL qL 5qL 40 EILa+→a?
4
?+ a 30 6 120EI kl
2 4
习题1.2 :
在用有限元法求解时,边界条件总是满足的,控制方程的不完全匹配,会产生误差。题中所
23
给出的近似函数:φa+ ax+ a x+ a x ,应该满足边界条件,对于情况(1 ) ,代入边
0 1 2 3
2
1aL a L
12
界条件可得aa 0, ,从而
03
3
L
3 3 3
x xx
2
φax ? + a x?+ (1 )
(注意:只要满足强制边界就可以,怎么判断是强制还是自然?)习题1.7
δ∏φδφφδφ?
k + k Qd δφα φ? q δφdΓ
[ ]
∫∫Γ
q
?xxy ?y?
22φφφ?
k + k + Qd δφ ?+ k δφ dα φ? q δφdΓ
[ ]
∫∫22?∫ΓΓ
q xy ?n22? φφ ?φ ?φ?
12
∫∫
12
可以求出待定系数aa ,。
12
伽辽金法作为加权余量法的特殊形式,权函数选择为插值函数NN , ,
高等有限元课后题答案
本文档如对你有帮助，请帮忙下载支持！2 弹性力学问题的有限单元法思考题2.1 有限元法离散结构时为什么要在应力变化复杂的地方采用较密网格，而在其他地方采用较稀疏网格？答：在应力变化复杂的地方每一结点与相邻结点的应力都变化较大，若网格划分较稀疏，则在应力突变处没有设置结点，而使得所求解的误差很大，若网格划分较密时，则应力变化复杂的地方可以设置更多的结点，从而使得所求解的精度更高一些。
2.2 因为应力边界条件就是边界上的平衡方程，所以引用虚功原理必然满足应力边界条件，对吗？答：对。
2.3 为什么有限元只能求解位移边值问题和混合边值问题？弹性力学中受内压和外压作用的圆环能用有限元方法求解吗？为什么？答：有限元法是一种位移解法，故只能求解位移边值问题和混合边值问题。
而应力边值问题没有确定的位移约束，不能用位移法求解，所以也不能用有限元法求解。
2.4 矩形单元旋转一个角度后还能够保持在单元边界上的位移协调吗？答：能。
矩形单元的插值函数满足单元内部和单元边界上的连续性要求，是一个协调元。
矩形的插值函数只与坐标差有关，旋转一个角度后各个结点的坐标差保持不变，所以插值函数保持不变。
因此矩形单元旋转一个角度后还能够保持在单元边界上的位移协调。
本文档如对你有帮助，请帮忙下载支持！2.5 总体刚度矩阵呈带状分布，与哪些因素有关？如何计算半带宽？答：因素：总体刚度矩阵呈带状分布与单元内最大结点号与最小结点号的差有关。
计算：设半带宽为B，每个结点的自由度为n，各单元中结点整体码的最大差值为D，则B=n(D+1) ，在平面问题中n=2。
2.6 为什么单元尺寸不要相差太大，如果这样，会导致什么结果？答：由于实际工程是一个二维或三维的连续体，将其分为具有简单而规则的几何单元，这样便于网格计算，还可以通过增加结点数提高单元精度。
在几何形状上等于或近似与原来形状，减小由于形状差异过大带来的误差。
若形状相差过大，使结构应力分析困难加大，误差同时也加大。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2 弹性力学问题的有限单元法思考题2.1 有限元法离散结构时为什么要在应力变化复杂的地方采用较密网格，而在其他地方采用较稀疏网格？答：在应力变化复杂的地方每一结点与相邻结点的应力都变化较大，若网格划分较稀疏，则在应力突变处没有设置结点，而使得所求解的误差很大，若网格划分较密时，则应力变化复杂的地方可以设置更多的结点，从而使得所求解的精度更高一些。
2.2 因为应力边界条件就是边界上的平衡方程，所以引用虚功原理必然满足应力边界条件，对吗？答：对。
2.3 为什么有限元只能求解位移边值问题和混合边值问题？弹性力学中受内压和外压作用的圆环能用有限元方法求解吗？为什么？答：有限元法是一种位移解法，故只能求解位移边值问题和混合边值问题。
而应力边值问题没有确定的位移约束，不能用位移法求解，所以也不能用有限元法求解。
2.4 矩形单元旋转一个角度后还能够保持在单元边界上的位移协调吗？答：能。
矩形单元的插值函数满足单元内部和单元边界上的连续性要求，是一个协调元。
矩形的插值函数只与坐标差有关，旋转一个角度后各个结点的坐标差保持不变，所以插值函数保持不变。
因此矩形单元旋转一个角度后还能够保持在单元边界上的位移协调。
2.5 总体刚度矩阵呈带状分布，与哪些因素有关？如何计算半带宽？ 答：因素：总体刚度矩阵呈带状分布与单元内最大结点号与最小结点号的差有关。
计算：设半带宽为B ，每个结点的自由度为n ，各单元中结点整体码的最大差值为D ，则B=n(D+1)，在平面问题中n=2。
2.6 为什么单元尺寸不要相差太大，如果这样，会导致什么结果？ 答：由于实际工程是一个二维或三维的连续体，将其分为具有简单而规则的几何单元，这样便于网格计算，还可以通过增加结点数提高单元精度。
在几何形状上等于或近似与原来形状，减小由于形状差异过大带来的误差。
若形状相差过大，使结构应力分析困难加大，误差同时也加大。
2.7 剖分网格时，在边界出现突变和有集中力作用的地方要设置结点或单元边界，试说明理由。
答：有限元处于弹性力学问题的方法是离散法。
它将一个受外力作用的连续弹性体离散成一定数量的有限小的单元集合体，单元之间只在结点上相互联系，即只有结点才能传递力。
所以在边界出现突变和有集中力作用的地方要设置结点和单元边界。
2.8 为什么说三角形三结点单元是常应变单元，如果在每边中点增加一个结点，那么单元内应力如何分布？答：(1)应变矩阵[B]中的参数m j i m j i c c c b b b 、、、、、由坐标变量x 、y 之差确定。
当单元的坐标差确定之后，这些参数与坐标变量x 、y 无关，因此[B]为常量阵。
当单元的结点位移{a}确定后，由[B]转换求得的单元应变都是常量，也就是说在荷载作用下单元中各点具有统一的xy y x γεε、、值。
因此三结点三角形单元称为常应变单元。
(2)如果在每边中点增加一个结点，单元内的应力为线性分布。
习 题2.1试证明x 、y 与面积坐标的关系 证明:设P 点坐标为（x,y ）j jpijy x yxy x A ii11121=()()()()[]()y c x b a y x x x y y y x y x xy y x y x y x x y y x m m m i j j i i j j i j i j i j i j i ++=-+-+-=---++=212121同理可求得：由面积坐标定义得：()y c x b a AA A L i i i ijmPjm i ++==21()()y c x b a y x y x y x A y c x b a y x y x yxA j j j iim m pmii i i mm j jpjm++==++==21111212111121()y c x b a AA A L j j j ijm Pmi j ++==21()y c x b a AA A L m m m ijmPij m ++==21由此推出坐标y x 、与面积坐标的函数关系：()()22j i i j j i i j i j j i m j j m m j m j m j j m A c L c L a c a c x b c b c A b L b L a b b a y b c b c ⎧-+-=⎪-⎪⎨-+-⎪=⎪-⎩式（2.1）面积：m i i m j m m j i j j i m j i c b c b c b c b c b c b a a a A -=-=-=++=2代入式（2.1）有：mj j m m j j m j m m j i j j i j i i j j i i j c b c b b a b a L b L b y c b c b c a c a L c L c x --+-=--+-=其中形状参数由下式确定：mj mji m j mji j m m j mm jji x x x x c y y y y b y x y x y x y x a +-==-=-=-==1111代入上式（2.1）可转化为：m m j j i i m m j j i i L y L y L y y L x L x L x x ++=++=再加上 m j i L L L ++=1 所以用面积坐标表示直角坐标矩阵形式如下：⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧m j i m jim j i L L L y y y x x x y x 1111 2.2 试证明两相似三角形的单元刚度矩阵相同。
证明：由于两个三角形相似，故设h A A =21， h 为一常数。
三角形：()111121i j ji i c b c b A -=111111i m j m j i y y b y y b -=-=111111i m j m j i x x c x x c +-=+-=参数Λj i j i c c b b 、、、，只与坐标差有关，所以hb b i i 121=hc c i i 121=单元刚度矩阵通式为：hb b b b s r s r 12211=h A A 121= 故[][]21rs rs K K =所以[][][][][][][][][][][][][]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==mm mjmijm jjjiimij ii Te K K K K K K K K K tA B D B K [][]eeK K 21=因此两相似三角形的单元刚度矩阵相同。
2.3 直角三角形固定在刚性基础上，受齐顶的水压力和自重作用，如图2.14所示。
若按一个单元计算，水的容重g γ，三角形平面构件容重g ρ,取泊松比v =1/6，试求顶点位移和固定面上的反力。
解：按逆时针编码，局部编码与整体编码相同：1-2-3建立坐标())0,0(3)3,0(20,21:a a xoy（1） 求形函数矩阵：[]()⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-+-+-+-+=-=s r s r s r s r s r s r s r s r rsb bc c c b b c b c c b c c b b A Et K 21212121142νννννννa a a a 600321===a b b a b 303321-===a c a c c 220321-===图（2.14） 形函数:)(21y c x b a A N i i i i ++=233221a a a A =⨯⨯=所以：ay a x N a yN a xN 32132321--===形函数的矩阵为：[][][][]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----==a y a x a y a x a y a x a y a x N N N N m ji3210302003210302（2） 刚度矩阵[][][][][][][][][][]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333231232221131211K K K K K K K K K K e[]()⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+-+-+-+=-=s r s r s r sr s r s r s r s r rsb bc c c b b c b c c b c c b b A Et K 21212121142ννννννν()125213531416122=-=-==νννa EA Et t可得：[][]⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎦⎤⎢⎢⎣⎡=400353534150093532211EK E K [][]⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=0251035343127273323531233E K E K[]⎥⎥⎦⎤⎢⎢⎣⎡----=215251935313E K []⎥⎥⎦⎤⎢⎢⎣⎡----=41253535323E K（3）位移列向量和右端项由边界条件可确定：{}{}Te u a 000022υ=水压力和构件厚分别为：10==t ghp γ[]⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=43127425215127332135259414001253503525021525025415019109353E K e{}T Te t l q h q h q R ⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧=03203100203600001自重为W 与支座反力：{}Ty x y x eW R R W WR R R ⎭⎬⎫⎩⎨⎧---=330333112所以：{}Ty x y x eW R h q R W h q W R R R ⎭⎬⎫⎩⎨⎧-+--=33363303011由[]{}{}eeeR a K =得到下列矩阵方程组：⎪⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎨⎧-+--=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧33363000030301122W R h q R W h q W R R u y x y x υ 化简得：⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧-=⎭⎬⎫⎩⎨⎧⎥⎥⎦⎤⎢⎢⎣⎡3640035353022W h q u E υ]⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=431274252151273321352594140012535035250215250254150191009353E K e可得：⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧-=⎭⎬⎫⎩⎨⎧E W E h q u 363567022υ将⎭⎬⎫⎩⎨⎧22υu 代入下式： ⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧-+-=⎭⎬⎫⎩⎨⎧⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡----333425135025103533031122W R h q R W R R u E y x y x υ 固定面上的反力：a h ga gh q 330===γγ从而可得支座反力为：4322123412033011h q W R h q W R W h q R WR y x y x -=-=+=-=2.4 试从式（2.69）说明对角线元素改1法只能用于给定零位移的情形，而对角元素充大数法可以适用任意指定位移的情形。