一、 课程设计题目 一）矩阵方程1. 利用全选主元的高斯约当（Gauss-Joadan ）消去法求解如下方程组，并给出详细的程序注解和说明：⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∙⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡1536353424543214019753910862781071567554321x x x x x 2. 利用追赶法求解如下方程组，并给出详细的程序注解和说明。

⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧-=⎪⎪⎪⎭⎪⎪⎪⎬⎫⎪⎪⎪⎩⎪⎪⎪⎨⎧∙⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡862031234567891011121354321x x x x x 3. 利用全选主元的高斯约当（Gauss-Joadan ）消去法如下求解大型稀疏矩阵的大型方程组，并给出详细注解及说明。

⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧----=⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎨⎧∙⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡4292728642-01-0100001-0402003-0001050006000102-00034-0002000006-00060020001-0087654321x x x x x x x x 二） 结构力学1. 试求解图示平面桁架各杆之轴力图，已知各材料性能及截面面积相同，27.90,210cm A Gpa E ==。

（注：在有限元分析中，桁架杆的模拟只能选择Ansys 的Link 单元）。

2. 试求解图示平面刚架内力图（轴力图、剪力图和弯矩图），已知各材料性能及截面面积相同，Gpa E 210=，泊松比0=μ。

3. 试求解图示平面刚架内力图（轴力图、剪力图和弯矩图），已知各杆432104.5,18.0,30m I m A Mpa E -⨯===。

KN101=二、 工作计划在老师给我们开完动员会以后，我们小组就进行了对要做那题进行讨论和分工，我做的是结构力学部分的第三题用位移法解出这道题。

（1）我在分工以后我就专一对结构力学部分第三题进行分析，因为有一个学期的时间没有接触结构力学，我在看完题后我就查找教材和网上的资料，看明白位移法的解题思路和特点。

（2）在明白位移法的解题思路后我就查找关于结构力学位移法经典例题的解题过程和老师所给课件上的解题过程。

从而对解决我所选择这题有个明了的思路。

（3）我开始老师所给课件的例题和查找的位移法经典例题解题过程解决这道题，在做这道题的过程中也遇到不理解的地方，最后在老师的帮助下我终于解出了这道题的结果，我把我做出的结果和我组用力法和ANSYS分析法所作出的结果进行比较，得出的结果基本相似，因为每种方法所用的条件不同，故答案有所误差，可以认为我们所做的答案是正确的。

课程负责人签名：指导教师签名：年月日年月日目录一：位移法介绍 (1)二：位移法主要内容 (1)三：元素刚度矩阵 (1)四：位移法解题步骤 (3)五：位移法结构分析 (5)六：结构力学题解题过程 (5)七：与其他方法结果对比和分析 (12)八：体会心得......................................... .14 九：参考文献......................................... .15结构力学课程设计一：位移法介绍位移法是以矩阵运算作为数学工具来处理结构位移计算的。

在结构力学的计算中，通过采用对结点位移作为基本未知量，进而通过矩阵的形式对各基本参数进行组织，编排，求出未知量的方法。

位移法解题思路为：首先对结构进行分析，确定独立的结点角位移和线位移的数量。

然后在每个独立的角位移上加上附加刚臂，每个独立的线位移上加上附加支座链杆，得到位移法的基本体系，再由基本体系与原体系等价根据平衡条件得到位移法典型方程，解位移法方程得到结点位移最后利用叠加原理得到结构的弯矩图。

二：位移法主要内容矩阵位移法主要内容包括两个部分：（1）单元分析，即将结构分解为有限个较小的单元，进行所谓的离散化。

对杆系结构，一般一根杆件或杆件的一段作为一个单元，分析单元的内力与位移的关系，建立单元刚度矩阵。

（2）整体分析，即将各单元又集成原来结构，要求各结构满足原结构的变形协调条件和平衡条件，从而建立整个结构的刚度方程，以求解原来的结构的位移和内力。

在杆系结构中，若单元只受轴力作用，则称为杆元素，如桁架；若单元不仅受轴力，还受剪力和弯矩的作用，则称为梁元素，如梁，刚架等。

三：元素刚度矩阵一般的平面杆件单元，每一杆端有三个杆端位移，即两个线位移和一个角位移；与此对应有三个杆端力，即两个集中力和一个集中力矩。

在局部坐标系（或称单元坐标系）轴与单元的的轴线重合，其正方向由单元始端i指向末端j。

杆端位移和杆端力统一以沿坐标正向为正，转角和力矩的正方向按照右手法则的方向，在图示坐标系中即为顺时针方向。

将单元（e）两个杆端的各杆端位移和杆端力分别按顺序组成列向量，可得到局部坐标系中的杆端位移向量和杆端力向量如下：单元（e）的杆端力和杆端位移之间的关系方程就是该单元的刚度方程，即其中，为局部坐标系中的单元刚度矩阵。

在矩阵位移法中，单元分析完成后，还要进行结构的整体分析。

整体分析是为了形成结构的整体刚度矩阵，建立整体结构的结点力和结点位移之间的关系，即整体刚度方程；然后根据整体结构各结点处的结点力与结点荷载之间的平衡条件建立位移法基本方程，由该方程解出结点位移。

结构的整体刚度方程可写成{F}=[K] {△}式中，{F} 为结构的整体结点力向量{△}为整体结点位移向量，[K]为整体刚度矩阵。

[K]可直接由各单元刚度矩阵集成得到。

集成的方法是根据各单元的局部和整体结点位移编码之间的关系，先将各单元刚度元素的下标局部编码换成整体编码，再将该元素叠加到整体刚度矩阵中与整体码对应的行和列的位置上，即所谓的“对号入座”。

该集成方法反映了整体结构的变形协调条件和平衡条件（未考虑结点荷载）。

集成时可采用先处理法将位移边界考虑进去，从而使整体刚度矩阵成为非奇异矩阵。

四：位移法解题步骤用矩阵位移法分析平面结构的一般步骤为：（1）划分单元，并对单元和结点位移进行编号，选取整坐标和局部坐标系。

（2）建立局部坐标系中的单元刚度矩阵，并经坐标变换得到整体坐标系中的单元刚度矩阵。

（3）按“对号入座”集成整体刚度矩阵。

（4）计算各单元的固端力，形成局部坐标系中的单元等效结点载荷向量，并变换成整体坐标系中的单元等效结点载荷向量，再集成结构的整体等效结点载荷向量。

（5）如整体刚度矩阵采用后处理法形成，则引入支撑条件，修改结构原始刚度矩阵。

（6）解位移法基本方程，求出结点位移。

（7）计算结构各杆（单元）的杆端力，最后杆端力应等于等效节点载荷引起的杆端力与单元固端力之和。

对于平面钢架用矩阵位移法进行分析并编制相应程序计算，可按下面的程序框图进行编程。

除了上述利用位移法基本体系在附加约束处的平衡条件来建立位移法的典型方程的方法，还可以直接由原结构的结点和截面平衡来建立位移法的方程。

具体步骤为：（1）结点位移分析 对结构进行分析，确定独立的角位移和线位移的数量，即确定位移法的基本未知量。

（2）确定基本体系 在每个独立的角位移上加上附加刚臂，每个独立的线位移上加上附加支座链杆得到位移法基本体系。

（3）列位移法方程 根据基本体系与原体系等价，由平衡条件列出位移法方程。

（4）作弯矩图 为计算附加约束上的反力，分别作出位移法基本结、构在荷载单位位移作用下的弯矩图。

（5）计算反力与反力系数 根据平衡条件，求出单位位移作用下附加约束上的力，系数，及单位位移作用下附加约束上的反力（自由项）。

（6）求结点位移 解位移法方程得到结点位移。

（7）作最终弯矩图 利用基本结构的荷载弯矩图和单位位移弯矩图，根据叠加原理求出杆端弯矩，然后用简支梁法作每个杆段弯矩图。

（8）作剪力图和轴力图 由杆件平衡利用杆端弯矩可计算出杆端剪 力并作图，再根据结点平衡计算出各杆的轴力并作图。

五：位移法结构分析位移法结构分析的步骤为：（1）选定结构总体坐标系，进行结点及元素编号；（2）按照编号顺序列出结点位移向量列阵及结点力向量列阵； （3）建立元素在总体坐标下的刚度矩阵； （4）用集合方法建立总体刚度矩阵； （5）利用给定位移边界条件，求得结点位移； （6）求得结构的支反力，结点力，内力等。

六：结构力学题解题过程3. 试求解图示平面刚架内力图（轴力图、剪力图和弯矩图），已知各杆432104.5,18.0,30m I m A Mpa E -⨯===。

KN101=解：1F134（1） 结构离散化将结构划分为4个结点，3个单元，截面积A=0.18m 2，惯性矩I=5.4×10-3m 4。

（2）求结点载荷首先需求局部坐标系中固定端内力{F 0}e 。

7.812510<a>单元1<b>单元3====18.75KN单元1===7.8125KN •m====10KN单元二====10KN •m在局部坐标系下单元载荷列向量单元1 单元2 单元3{}0.F 1= ⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-8125.775.1808125.775.180 {}0.F 2= ⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡000000 {}0.F 3= ⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-1010010100为求出在整体坐标下的载荷列向量，先求单元得坐标转换矩阵[]T 。

单元1,2 α=00[]T 1=[]T 2=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--10cos sin 0000sin cos 0000001000000cos sin 0000sin cos αααααααα =⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡100000010000001000000100000010000001单元3 α=900[]T 3=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--1000000cos sin 0000sin cos 0000001000000cos sin 0000sin cos αααααααα=⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--100001000010000000100000001000010求个单元在整体坐标下的等效结点载荷{}TT P 10-=⎭⎬⎫⎩⎨⎧-{}{}{}{}1020110108125.775.1808125.775.180⎭⎬⎫⎩⎨⎧=⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧--=-=P P F F{}2P []{}{}{}20302202000000⎭⎬⎫⎩⎨⎧=⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=-=P P F T T{}[]{}{}{}30204303301001010010⎭⎬⎫⎩⎨⎧=⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧-=-=P P F T P T求刚架的等效结点载荷{}0p{}0p {}{}{}=++=32010p p p ⎪⎪⎪⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎪⎨⎧---0000008125.7175.1808125.775.180⎪⎪⎪⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎪⎨⎧+000000000000+⎪⎪⎪⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎪⎨⎧-1001000010010000⎪⎪⎪⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎪⎪⎪⎬⎫⎪⎪⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎪⎪⎨⎧----=100100008125.1775.18108125.775.180{}[]Tp 1001001008125.4775.18108125.775.180-----=(3)求单元刚度矩阵表达式由于单元1、2的尺寸相同，材料弹性模量相同，故[][]21K K =[]=e K ⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--------l EI l EI lEI l EI l EI l EI l EI l EI l EA l EAl EI l EI l EI l EI l EI l EI l EI l EI l EA l EA 460260612061200000260460612061200000222323222323 则[][]==21K K ⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--------2592001555201296001555200155520124416015552012441600021600000021600001296001555200259200155520015552012441601555201244160002160000002160000[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡---------=162000607508100060750060750303750607503037500013500000013500008100060750016200060750060750303750607503037500013500000013500003K（4）、求整体坐标系中的[]eK单元1、[][][][][]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--------=⎥⎥⎦⎤⎢⎢⎣⎡===25920015552001296001555200155520124416015552012441600021600000021600012960015552002592001555200155520124416015552012441600216000000216000012212111211111111K K K K K T K T K T 单元2[][][][][]123323222322222222_K K K K K K T K T K T =⎥⎥⎦⎤⎢⎢⎣⎡===⎥⎥⎦⎤⎢⎢⎣⎡ 单元3[][][][]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--------=⎥⎥⎦⎤⎢⎢⎣⎡==162000060750810006075001350000001350000060750030375607500303758100006075016200006075001350000001350000060750030375607500303753223243443423333K K K K T K T K T(5)、求结构刚度矩阵[]K利用刚度集成法[]⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡++=3443422332323242233222221221121121110000K K K K K K K K K K K K K(6)、建立原始平衡方程式⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡++432143213443422332323242233222221221121121110000P P P P K K K K K K K K K K K K δδδδ (7）、引入约束条件解方程组由于1为固定支座，4为固定端，则0u 44411=====A u v v .则建立平衡方程为⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡-----------129600155520012960015552000155520124416015552012441600002160000002160000012960015552006804000607501296001555201244160015988320155520002160000607500435037500001296001555200259200⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡3332221θθθv u v u ⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=3332221W Q P W Q P W ⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡---=01008125.4775.18108125.7解得：⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡3332221θθθv u v u =⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--0214.057615.000256.0001217.002037.0002472.001956.0 ⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡W Q P Q P 44411=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡----0162000060750000001350000000050750030375000015552012441601555200000021600000=⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡3332221θθθv u v u =⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡--742.15357.27765.14072.19427.5 求得弯矩图，剪力图，轴力图，如下图所示七：与其他方法结果对比和分析力法结果：ANSYS 分析结果：最后得到位移和受力以及弯矩力：PRINT F SUMMED NODAL LOADS***** POST1 SUMMED TOTAL NODAL LOADS LISTING *****LOAD STEP= 1 SUBSTEP= 1 TIME= 1.0000 LOAD CASE= 0THE FOLLOWING X,Y,Z SOLUTIONS ARE IN THE GLOBAL COORDINATE SYSTEMNODE FX FY FZ1 14465. -27522.17 -20000.32 5534.8 -19978.62 10000.TOTAL VALUESVALUE -0.85493E-09 -37500. 0.0000位移：PRINT U NODAL SOLUTION PER NODE***** POST1 NODAL DEGREE OF FREEDOM LISTING *****LOAD STEP= 1 SUBSTEP= 1TIME= 1.0000 LOAD CASE= 0THE FOLLOWING DEGREE OF FREEDOM RESULTS ARE IN THE GLOBAL COORDINATESYSTEMNODE UX UY UZ USUM1 0.0000 0.0000 0.0000 0.00002 0.25624E-02 -0.20387E-01 0.0000 0.20547E-0117 0.32225 -0.10194E-01 0.0000 0.3224132 0.0000 0.0000 0.0000 0.000062 0.25624E-02 -0.57615 0.0000 0.57616弯矩：NODE MX MY MZ1 0.15024E-12 -0.72235E-12 15930.上图是我用位移法和另外两个小组成员，用力法和ANSYS分析方法得出的结果，进行的柱状图的对比。