初始参数只有这些n = 3; rho = 2.7 * 1e3;%这个变量表示材料的密度。 S = 0.1 * 0.01;%这个变量表示横截面积。 E = 7.2 * 1e10;%这个变量表示杨氏模量。I = 0.1^3*0.01/12;L=3
时间: 2024-03-27 22:40:46 浏览: 21
好的,根据提供的初始参数,我们可以先计算出单元刚度矩阵和质量矩阵,然后将其组装成6*6的全局矩阵。这里假设每个单元有2个节点,每个节点有3个自由度(x、y方向的位移和转角),因此单元刚度矩阵和质量矩阵的大小均为6*6。我们使用类似上面的循环遍历方法,将单元矩阵组装成6*6的全局矩阵。
```matlab
% 根据给定的参数计算单元的材料和几何参数
L = 3; % 单元长度
E = 7.2e10; % 弹性模量
rho = 2.7e3; % 密度
S = 0.1 * 0.01; % 横截面积
I = 0.1^3 * 0.01 / 12; % 截面惯性矩
% 计算单元刚度矩阵和质量矩阵
Ke = (E*I/L^3) * [12, 6*L, -12, 6*L; 6*L, 4*L^2, -6*L, 2*L^2; -12, -6*L, 12, -6*L; 6*L, 2*L^2, -6*L, 4*L^2];
Me = (rho*S*L/420) * [156, 22*L, 54, -13*L; 22*L, 4*L^2, 13*L, -3*L^2; 54, 13*L, 156, -22*L; -13*L, -3*L^2, -22*L, 4*L^2];
% 定义全局刚度矩阵和质量矩阵的大小
n = 3; % 假设有3个节点
K = zeros(3*n, 3*n); % 初始化全局刚度矩阵
M = zeros(3*n, 3*n); % 初始化全局质量矩阵
% 循环遍历每个单元,将其贡献添加到全局刚度矩阵和质量矩阵中
for i = 1:n-1
K_i = Ke; % 第i个单元的刚度矩阵
M_i = Me; % 第i个单元的质量矩阵
% 将第i个单元的刚度矩阵和质量矩阵添加到全局刚度矩阵和质量矩阵中
K(3*i-2:3*i+3, 3*i-2:3*i+3) = K(3*i-2:3*i+3, 3*i-2:3*i+3) + K_i;
M(3*i-2:3*i+3, 3*i-2:3*i+3) = M(3*i-2:3*i+3, 3*i-2:3*i+3) + M_i;
end
```
在上面的代码中,我们首先根据给定的参数计算出单元的材料和几何参数,然后计算出单元的刚度矩阵和质量矩阵。接着,我们定义了全局刚度矩阵和质量矩阵的大小,并初始化为零矩阵。最后,我们循环遍历每个单元,将其贡献添加到全局刚度矩阵和质量矩阵中。这样,我们就可以得到组装后的6*6全局刚度矩阵和质量矩阵了。
相关推荐
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)