2021.3 Matlab与SAP2000静力学与动力学计算结果对比

一、静力学：平面梁单元（剖分为4段）

模型参数：

平面4m跨固接梁，C30截面120x100mm（便于计算惯性矩等），跨中荷载10KN，剖分4段。

模型简图：

注：1.静力学对比方法为：在matlab中建立刚度矩阵K和质量矩阵M，输入边界条件和外荷载进行内力位移计算，与sap2000对比。

2. 做动力学对比，则需要从SAP2000输出K及M，因为SAP2000使用的梁单元考虑了剪切变形，而在matlab中使用的是纯弯曲梁单元。

1.位移对比

 2.部分代码

 

 

二、动力学：二层框架

三个对比：1.周期对比，2.标准模态对比，3.模态参与系数对比，4.部分代码

注：

1.此计算未在matlab中进行刚度矩阵K和质量矩阵M组装，均由SAP2000导出模型K及M，由matlab计算。

2.此部分重点为二阶微分方程求解，及模态归一化理论。

1.周期对比（二层模型）

2.标准模态对比（二层模型）

注：以一个节点为例

模态1：

模态2

 

模态3 

其余模态同，不再列出 

3.模态参与系数对比（另一单层模型）

 注：par Mod1 u3 代表 参与系数-模态1-自由度3

4.部分代码