高等有限元课后题答案 (1)
有限元课后习题答案
有限元课后习题答案1.1有限元法的基本思想和基本步骤是什么首先，将表示结构的连续离散为若干个子域，单元之间通过其边界上的节点连接成组合体。
其次，用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。
步骤：结构的离散化，单元分析，单元集成，引入约束条件，求解线性方程组，得出节点位移。
1.2有限元法有哪些优点和缺点优点：有限元法可以模拟各种几何形状复杂的结构，得出其近似解；通过计算机程序，可以广泛地应用于各种场合；可以从其他CAD软件中导入建好的模型；数学处理比较方便，对复杂形状的结构也能适用；有限元法和优化设计方法相结合，以便发挥各自的优点。
缺点：有限元计算，尤其是复杂问题的分析计算，所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。
对无限求解域问题没有较好的处理办法。
1.3有限元法在机械工程中有哪些具体的应用静力学分析模态分析动力学分析热应力分析其他分析2.1杆件结构划分单元的原则是什么？1）杆件的交点一定要取为节点2）阶梯形杆截面变化处一定要取为节点3）支撑点和自由端要取为节点4）集中载荷作用处要取为节点5）欲求位移的点要取为节点6）单元长度不要相差太多2.2简述单元刚度矩阵的性质。
单元刚度矩阵是描述单元节点力与节点位移之间关系的矩阵。
2.3有限元法基本方程中每一项的意义是什么？{Q}---整个结构的节点载荷列阵（包括外载荷、约束力）;{}---整个结构的节点位移列阵；[K]---结构的整体刚度矩阵，又称总刚度矩阵。
2.4简述整体刚度矩阵的性质和特点。
对称性奇异性稀疏性主对角上的元素恒为正2.5位移边界条件和载荷边界条件的意义是什么由于刚度矩阵的线性相关性不能得到解，从而引入边界条件。
2.6写出平面刚架问题中单元刚度矩阵的坐标变换式2.7推导平面刚架局部坐标系下的单元刚度矩阵。
2.8简述整体坐标的概念。
单元刚度矩阵的坐标变换式把平面刚架的所有单元在局部坐标系X’O’Y’下的单元刚度矩阵变换到一个统一的坐标系xOy下，这个统一的坐标系xOy称为整体坐标系。
有限元习题及答案ppt课件
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
病原体侵入机体，消弱机体防御机能 ，破坏 机体内 环境的 相对稳 定性， 且在一 定部位 生长繁 殖，引 起不同 程度的 病理生 理过程
中国科学院大学有限元试题及答案
(1)引入边界条件： v1 0,1 0, v2 0, M 3 m, M 2 0, Y3 0 由后三个方程可求得 2、v3、 3 ，然后把 2、v3、 3 代入前三个方程，求得 Y1、M 1、Y2 。
例１：已知：p，l，EA。求： u 2 , v 2
解：方法１：１）划分单元，给节点编号 ２）单元分析 ①单元： 0, cos 1, sin 0
3
p
10
9
7
y
8 5
1
1
解：
6
9
8
x
6
3
7
5
2
2
4
3
题3 图
4
题3图. 三角形结构网 格
（2） d 4,
M B 2(d1 v4 0
4
4
7
15 10
11
3
1
2
6
13 15
题3图
5
9 12 14
答： （2） d=4 , B=2(d+1)=10 (3) u1 u15 v1 v15 0
p 作用。杆件沿 y 轴方向，长为 a 1 m ，截面积 A 0.01m 2 ，
E2 E0 。载荷及约束信息如图示，自重不计。试采用图示的
1个三角形常应变元和1个平面杆元求： （1）结构整体的等效结点力列阵； （2）采用划行划列法引入已知结 点位移，计算出结点1和2的 a 位移； （3）杆件中内力。 i j m 单元2： 1 3 2 单元1： 2 4
答： 在有限单元法中，采用低阶多项式拟合振型。结构的低阶振 型曲线与低阶多项式比较通配，结构的高阶振型曲线与低阶 多项式曲线有着显著的差异。因而，有限元法中求出的低阶 频率和振型是可信的，而所求出的高阶频率和振型误差较大 ，甚至无效。
有限元试题及答案
有限元试题及答案 有限元试题及答案 一 判断题（20分）（×）1. 节点的位置依赖于形态，而并不依赖于载荷的位置（√）2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元 （×）3. 不能把梁单元、壳单元和实体单元混合在一起作成模型 （√）4. 四边形的平面单元尽可能作成接近正方形形状的单元（×）5. 平面应变单元也好，平面应力单元也好，如果以单位厚来作模型化 处理的话会得到一样的答案（×）6. 用有限元法不可以对运动的物体的结构进行静力分析 （√）7. 一般应力变化大的地方单元尺寸要划的小才好（×）8. 所谓全约束只要将位移自由度约束住，而不必约束转动自由度 （√）9. 同一载荷作用下的结构，所给材料的弹性模量越大则变形值越小 （√）10一维变带宽存储通常比二维等带宽存储更节省存储量。
二、填空（20分）1．平面应力问题与薄板弯曲问题的弹性体几何形状都是 薄板 ，但前者受力特点是： 平行于板面且沿厚度均布载荷作用 ，变形发生在板面内；后者受力特点是： 垂直于板面 的力的作用，板将变成有弯有扭的曲面。
2．平面应力问题与平面应变问题都具有三个独立的应力分量： σx ，σy ，τxy ，三个独立的应变分量：εx ，εy ，γxy ，但对应的弹性体几何形状前者为 薄板 ，后者为 长柱体 。
3．位移模式需反映 刚体位移 ，反映 常变形 ，满足 单元边界上位移连续 。
4．单元刚度矩阵的特点有：对称性 ， 奇异性 ，还可按节点分块。
5．轴对称问题单元形状为：三角形或四边形截面的空间环形单元 ，由于轴对称的特性，任意一点变形只发生在子午面上，因此可以作为 二 维问题处理。
6．等参数单元指的是：描述位移和描述坐标采用相同的形函数形式。
等参数单元优点是：可以采用高阶次位移模式，能够模拟复杂几何边界，方便单元刚度矩阵和等效节点载荷的积分运算。
7．有限单元法首先求出的解是 节点位移 ，单元应力可由它求得，其计算公式为{}{}[][]eD B σδ=。
2012级有限元1标准答案
一1，数值分析2，离散3，节点节点点二1，三个基本方程①平衡方程②本构方程平面应力平面应变应变协调方程③几何方程二类边界条件①力的边界条件②位移边界条件如今给定的位移边界为，则有（在），，其中分别为边界上x，y方向上的位移分量2，①步骤：⑴将结构离散化⑵单元分析，求得单元节点位移与节点力的关系，计算单元刚度矩阵⑶以节点为隔离体，建立平衡方程⑷施加荷载⑸引入边界条件⑹求解方程，求得节点位移⑺对每一单元循环，由单元节点位移通过单元刚度矩阵求得单元应力或杆件内力②表达式：位移模式几何矩阵[B]=弹性矩阵应力矩阵3，加权余量法：当n有限时，定解方程存在偏差（余量），取权函数，强迫余量在某种平均意义上均为采用使余量的加权积分为0的等效积分以“弱”形式来求得微分方程近似解的方法。
半解析法：离散与解析相结合的方法，减少计算工作量，降低费用。
样条有限元法：具有紧凑型及良好的光滑性，明确的表达式的优点，所得到的结果均在单元节点上，在数据的后处理方面更为方便和精确。
边界单元法：将所研究问题的偏微分方程，设法转换为在边界上定义的边界积分方程，然后将边界积分方程离散化为只含有边界结点未知量的代数方程组，解此方程组可得边界节点上的未知量并可由此进一步求得所研究区域中的未知量，它除了能处理有限元方法所适应的大部分问题外，还能处理有限元法不易解决的无限域问题。
4，⑴整体刚度矩阵是对称矩阵⑵整体刚度矩阵的主对角线上的元素总是正的⑶整体刚度矩阵是一个稀疏阵⑷施加荷载没有支承的整体刚度矩阵是一个奇异阵5，不能。
因为不满足完备性，缺少表示刚体位移的常数项和表示应变是位移一阶导数的常应变项不能保证解的收敛性。
6，维数、单元格划分，节点数目。
7，①②理由：单元刚度矩阵不随单元（或坐标轴）的平行移动或作n元（n为整数）角度的移动而改变8，有限元解位移大于解析解的原因是单元为非完全协调单元。
挠度w是弯曲问题中的基本未知函数且由于忽略了Z方向的变化，因此它只是x,y的函数：，若w已知，则唯一、内力、应力均可按上述相应公式求出。
有限元习题及答案
有限元习题及答案有限元习题及答案有限元方法是一种常用的数值计算方法，用于求解各种工程和科学问题。
在学习有限元方法的过程中，练习习题是非常重要的，可以帮助学生巩固所学的知识，并提高解决实际问题的能力。
本文将介绍一些有限元习题及其答案，希望对学习有限元方法的同学有所帮助。
习题一：一维热传导问题考虑一个长度为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方向上的分量。
边界条件根据具体情况给定。
通过数值方法，如有限元方法，可以求解位移场的近似解。
习题三：三维流体力学问题考虑一个三维流体力学问题，流体在一个封闭容器内流动，容器的形状为一个长方体，已知流体的速度场和压力场的初始条件，求解流体的运动状态。
有限元方法例题解答
2023《有限元技术》习题一参考答案1、用欧拉方程求泛函()1022[()]'2(0)0,(1)0J y x y y xy dx y y ⎧=--⎪⎨⎪==⎩⎰的极值曲线。
解：22'2F y y xy =--，代入欧拉方程'0y y dF F dx-=， 得：''++0y y x =，解微分方程得通解：12sin cos y C x C x x =+-，代入边界条件(0)0,(1)0y y ==，解得sin sin1xy x =-。
2、如图所示，一长度为L 质量为M 的项链悬挂在跨度为2a 的A 和B 两点，项链在重力场中自然下垂，试求该链悬在稳定状态时的曲线方程。
（重力加速度为g ）解： （方法一）将原坐标系(),x y 向下平移1C 个单位(),x y ，拟采用1cosh y C t =代换求解 在新坐标系中，悬链在稳定状态时能量处于最小值。
悬链线质量密度MLλ=， 长度为dl 的势能为：()Mgy dW dmgy dl gy dl L λ====，悬链总势能泛函：(a a a a Mg W dW dx dx L --===⎰⎰⎰，约束条件为：悬链线长度aL -=⎰，泛函的被积函数：(),F y y '=，势能泛函取极小值时的欧拉方程为：'1'y F y F C -=， 即：21C -=，化简得：y C =于是：dx =x =，令1cosh y C t =（在新坐标系下才能作此代换），得：1sinh sinh dy C tdt t =⎧=，代入x =，得112x C dt C t C ==+⎰所以，21x C t C -=，21cosh cosh x C t C ⎛⎫-= ⎪⎝⎭回代1cosh y C t =得：211cosh x C y C C ⎛⎫-= ⎪⎝⎭，曲线关于y 轴对称得20C =，1C由悬链线长度112sinhaaL C C -==⎰给出， 故新坐标系下所求曲线方程为11cosh x y C C ⎛⎫=⎪⎝⎭， 1C 由11sinh 2L aC C =确定。
有限元法基础习题答案
有限元法基础习题答案有限元法是一种常用的工程分析方法，广泛应用于结构力学、热传导、流体力学等领域。
它通过将复杂的物理问题离散化为一系列简单的子问题，并利用数值方法求解这些子问题，从而得到整体问题的近似解。
在学习有限元法的过程中，习题是必不可少的一环。
本文将给出一些有限元法基础习题的答案，希望能够帮助读者更好地理解和掌握这一方法。
习题一：一维线性弹性力学问题考虑一根长度为L的弹性杆，杆的截面积为A，杨氏模量为E。
在杆的一端施加一个沿杆轴向的拉力F，另一端固定。
假设杆轴向变形u(x)满足以下方程：EAu''(x) = -F，0 < x < Lu(0) = 0, u(L) = 0其中，u''(x)表示u(x)对x的二阶导数。
解答：根据上述方程，我们可以得到杆的位移函数u(x)的表达式。
首先，对方程两边进行积分，得到：EAu'(x) = -Fx + C1其中，C1为积分常数。
再次对方程两边进行积分，得到：EAu(x) = -F/2*x^2 + C1*x + C2其中，C2为积分常数。
根据边界条件u(0) = 0，可得C2 = 0。
代入边界条件u(L) = 0，可得：EAu(L) = -F/2*L^2 + C1*L = 0由此可得C1 = F/2*L。
将C1代入上式，可得：EAu(x) = -F/2*x^2 + F/2*L*x最终得到杆的位移函数u(x)的表达式为：u(x) = (-F/2*E)*(x^2 - L*x)，0 < x < L习题二：二维平面弹性力学问题考虑一个正方形薄板，边长为L，板的厚度为h。
假设薄板的杨氏模量为E，泊松比为ν。
在薄板的一侧施加一个沿法向的均匀表面压力P，另一侧固定。
求薄板的位移和应力分布。
解答：根据平面弹性力学理论，我们可以得到薄板的位移和应力分布。
首先，根据杨氏模量E、泊松比ν和薄板的厚度h，可以计算出薄板的弹性模量D：D = E*h^3 / (12*(1-ν^2))接下来，根据薄板的边界条件和平衡方程，可以得到薄板的位移和应力分布。
- 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 iiT e 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 hq R ⎭⎬⎫⎩⎨⎧=⎭⎬⎫⎩⎨⎧=0320310020360001自重为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法只能用于给定零位移的情形，而对角元素充大数法可以适用任意指定位移的情形。