可视化 | MATLAB画杆系结构变形图_MATLAB-仿真秀干货文章

可视化 | MATLAB画杆系结构变形图

一图胜千言。将计算结果用图表达出来定是极好的！

调用MATLAB中的line函数可以画直线。例如line([1,2],[3,4])将画出（1,3）到(2,4)的一条直线，而不是(1,2)到(3,4)。

当调用line(A,B)中，a b是相同大小的矩阵时，将会在对应的每一列做一条直线。当调用line(A,B,C)时，相应地会在三维图中画一条线。此时的A,B,C 均为2xN矩阵。总结为line([起点横坐标，终点横坐标],[起点纵坐标，终点纵坐标])，注意line与plot两个函数的区别。

现在用一个简单平面桁架的变形来实现

结点2有一个向下的集中荷载，在计算出结点位移后就可以画变形图了。画图时调用的函数为

为便于比较，将变形前后的图画在一起

来源：数值分析与有限元编程

矩阵方程

对于矩阵 A(n,n) 和 B(n,m) 组成的矩阵方程[A][X] = [B]记 X(n,m) 的第i列向量为 Xi(i = 1,2...m)，矩阵B的第i列向量为 Bi(i = 1,2...m), 则上述方程等价为但在实际计算时不应该分别求解，如果是这样的话就造成计算机资源极大浪费，而应该是对所有向量一次选主元消去，然后分别回代。即可以得到方程的解矩阵X。具体做法是将矩阵A(n,n)和B(n,m)组成增广矩阵[AB],通过选主元消去将AB的第1列至第n列变成上三角矩阵，用解上（下）三角方程组的回带方法解方程组 [Aup][Xi] = [Bi] (i = 1,2...m) 。再把Xi拼接起来就得到 [X] 了。以下是模块代码：用以下的矩阵方程来验证输出结果为来源：数值分析与有限元编程