结构力学-矩阵位移法课件.docxVIP

这学期我在给土木工程专业大三学生讲授结构力学课程时，发现矩阵位移法这一章节总是让同学们感到困惑。上周五下午的课后答疑时间，有七八个学生围着我，说他们虽然理解了基本概念，但在实际计算过程中总是出错。特别是张明同学，他在计算一个两跨连续梁时，刚度矩阵组装出现了问题，导致最终位移结果偏差达到了15%。这让我意识到，现有的课件可能需要更详细的步骤说明和实例演示。

在实际操作中，我们发现矩阵位移法的难点主要集中在刚度矩阵的组装和边界条件的处理上。以我们上周分析的五层两跨钢框架结构为例，该结构共有15个节点，每个节点考虑3个自由度（水平位移、竖向位移和转角），总自由度数为45个。我们采用直接刚度法进行计算，需要建立每个单元的刚度矩阵。

对于标准梁单元，其局部坐标系下的刚度矩阵为12×12阶矩阵，其中弯曲刚度EI=2.1×10?kN·cm2，轴向刚度EA=3.3×10?kN。在转换到整体坐标系时，需要通过坐标转换矩阵进行变换，转换矩阵中的角度θ根据单元的实际倾角确定。例如，对于倾斜30°的斜柱，其转换矩阵中的cosθ=0.6，sinθ=0.5。

在组装整体刚度矩阵时，我们采用对号入座的方法。以第3层第2跨的梁单元为例，该单元连接节点8和节点9，在整体刚度矩阵中对应的自由度位置为第2224行和第2527行。通过MATLAB编程实现时，我们使用双重循环语句，将单元刚度矩阵的元素准确叠加到整体刚度矩阵的相应位置。

边界条件的处理是另一个关键点。对于固定支座，我们采用置1法处理，即将对应自由度的主对角线元素置1，其余行和列的元素置0。以底层柱脚为例，节点13的平动自由度被约束，我们在整体刚度矩阵中将第16行和第16列进行相应处理。

通过这种方法计算得到的结构顶点水平位移为12.3mm，与ANSYS软件计算结果的误差仅为2.1%，完全满足工程精度要求。学生在掌握这个具体案例后，对类似结构的分析就有了明确的参照。

总的来看，下一阶段的重点是让同学们独立完成一个完整的矩阵位移法计算作业。我会在下周三的课堂上发放一份详细的作业指导书，包含一个三层单跨框架结构的全部参数和计算要求。作业将分三个部分提交：第一部分是手算的刚度矩阵组装过程，占30分；第二部分是MATLAB编程实现，占40分；第三部分是结果分析并与ANSYS软件对比，占30分。请同学们务必在11月25日前通过课程平台提交电子版，纸质版请在课堂结束时交到讲台上。对于计算过程中遇到的具体问题，欢迎在每周二下午24点的答疑时间来办公室找我讨论。

王建国

2025年11月17日