matlab 单元质量矩阵,这样组装平面桁架的质量矩阵对不对呢？

柔情只为你懂 2023-01-15 14:30 118阅读 0赞

%单元刚度矩阵

k1=\[0.7375000E+08  0.0000000E+00 -0.7375000E+08  0.0000000E+00;

0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00;

\-0.7375000E+08  0.0000000E+00  0.7375000E+08  0.0000000E+00;

0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00;\];

k2=\[0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00;

0.0000000E+00  0.9833333E+08  0.0000000E+00 -0.9833333E+08;

0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00;

0.0000000E+00 -0.9833333E+08  0.0000000E+00  0.9833333E+08;\];

k3=\[0.3776000E+08  0.2832000E+08 -0.3776000E+08 -0.2832000E+08;

0.2832000E+08  0.2124000E+08 -0.2832000E+08 -0.2124000E+08;

\-0.3776000E+08 -0.2832000E+08  0.3776000E+08  0.2832000E+08;

\-0.2832000E+08 -0.2124000E+08  0.2832000E+08  0.2124000E+08;\];

k4=\[0.7375000E+08  0.0000000E+00 -0.7375000E+08  0.0000000E+00;

0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00;

\-0.7375000E+08  0.0000000E+00  0.7375000E+08  0.0000000E+00;

0.0000000E+00  0.0000000E+00  0.0000000E+00  0.0000000E+00;\];

%单元质量矩阵

m1=\[0.9466667E-01  0.0000000E+00  0.4733333E-01  0.0000000E+00;

0.0000000E+00  0.9466667E-01  0.0000000E+00  0.4733333E-01;

0.4733333E-01  0.0000000E+00  0.9466667E-01  0.0000000E+00;

0.0000000E+00  0.4733333E-01  0.0000000E+00  0.9466667E-01;\];

m2=\[0.7100000E-01  0.0000000E+00  0.3550000E-01  0.0000000E+00;

0.0000000E+00  0.7100000E-01  0.0000000E+00  0.3550000E-01;

0.3550000E-01  0.0000000E+00  0.7100000E-01  0.0000000E+00;

0.0000000E+00  0.3550000E-01  0.0000000E+00  0.7100000E-01;\];

m3=\[0.1183333E+00  0.0000000E+00  0.5916667E-01  0.0000000E+00;

0.0000000E+00  0.1183333E+00  0.0000000E+00  0.5916667E-01;

0.5916667E-01  0.0000000E+00  0.1183333E+00  0.0000000E+00;

0.0000000E+00  0.5916667E-01  0.0000000E+00  0.1183333E+00;\];

m4=\[0.9466667E-01  0.0000000E+00  0.4733333E-01  0.0000000E+00;

0.0000000E+00  0.9466667E-01  0.0000000E+00  0.4733333E-01;

0.4733333E-01  0.0000000E+00  0.9466667E-01  0.0000000E+00;

0.0000000E+00  0.4733333E-01  0.0000000E+00  0.9466667E-01;\];

%组装刚度矩阵

KK=zeros(8,8);

KK=Bar2D2Node\_Assembly(KK,k1,1,2);

KK=Bar2D2Node\_Assembly(KK,k2,3,2);

KK=Bar2D2Node\_Assembly(KK,k3,1,3);

KK=Bar2D2Node\_Assembly(KK,k4,3,4);

%组装质量矩阵

MM=zeros(8,8);

MM=Bar2D2Node\_Assembly(MM,m1,1,2);

MM=Bar2D2Node\_Assembly(MM,m2,3,2);

MM=Bar2D2Node\_Assembly(MM,m3,1,3);

MM=Bar2D2Node\_Assembly(MM,m4,3,4);

function z = Bar2D2Node\_Assembly(KK,k,i,j)

%该函数进行单元刚度矩阵的组装

%输入单元刚度矩阵 k，单元的节点编号 i、 j

%输出整体刚度矩阵 KK

%--------------------------------------------------------

DOF(1)=2\*i-1;

DOF(2)=2\*i;

DOF(3)=2\*j-1;

DOF(4)=2\*j;

for n1=1:4

for n2=1:4

KK(DOF(n1),DOF(n2))= KK(DOF(n1),DOF(n2))+k(n1,n2);

end

end

z=KK;

![t-9385622-1][]

qq.png

[t-9385622-1]: http://muchong.com/t-9385622-1

发表评论取消回复

表情：

评论列表（有 0 条评论，118人围观）

还没有评论，来说两句吧...

相关阅读

相关 MATLAB 对矩阵的操作

1.常见矩阵函数运算 2.矩阵翻转 3.矩阵抽取矩阵函数运算： inv（A）：对矩阵进行求逆 det（A）：对矩阵求行列式 rank（A）：求矩阵的秩 eig（

古城微笑少年丶/ 2024年02月19日 15:29/ 0 赞/ 76 阅读

相关对保障代码质量至关重要的单元测试，你真的了解吗？

为什么要有这篇文章？个人在不断实践中越发觉得，单元测试对于代码质量的保障真的太有意义了，至少能体现在如下两个方面： ①让你写出更好的代码，可测试的代码一定是优雅的代码

分手后的思念是犯贱/ 2023年09月26日 10:32/ 0 赞/ 16 阅读

相关 matlab 单元质量矩阵,这样组装平面桁架的质量矩阵对不对呢？

%单元刚度矩阵 k1=\[0.7375000E+08 0.0000000E+00 -0.7375000E+08 0.0000000E+00; 0.0000000E+00

柔情只为你懂/ 2023年01月15日 14:30/ 0 赞/ 119 阅读

相关 matlab矩阵建模,MATLAB矩阵与数组

该楼层疑似违规已被系统折叠隐藏此楼查看此楼 MATLAB既然以矩阵实验室命名，就说明该软件在矩阵计算方面具有非常优异的表现。在MATLAB中，一般情况下一个矩阵就是指一个长

本是古典何须时尚/ 2023年01月14日 15:57/ 0 赞/ 198 阅读

相关 matlab 函数矩阵,MATLAB常用矩阵函数

《MATLAB常用矩阵函数》由会员分享，可在线阅读，更多相关《MATLAB常用矩阵函数(2页珍藏版)》请在人人文库网上搜索。 1、1. 矩阵的构造与操作zeros 生成元素全

水深无声/ 2023年01月13日 12:57/ 0 赞/ 179 阅读

相关谈谈对于保障代码质量至关重要的单元测试

为什么要有这篇文章呢？个人在不断实践中越发觉得,单元测试对于代码质量的保障真的太有意义了，至少能体现在如下两个方面: ①让你写出更好的代码，可测试的代码一定是优雅的代

墨蓝/ 2022年11月13日 04:24/ 0 赞/ 177 阅读

相关 matlab-线性代数 det+inv 对n阶矩阵求伴随矩阵

[2019独角兽企业重金招聘Python工程师标准>>> ][2019_Python_] ![hot3.png][] > matlab : R2018a 64bit

向右看齐/ 2022年10月02日 05:52/ 0 赞/ 244 阅读

相关 MATLAB矩阵

目录 MATLAB矩阵 MATLAB引用一个矩阵的元素 MATLAB删除行或列矩阵详细例子 -------------------- MATLAB矩阵矩阵

待我称王封你为后i/ 2022年08月28日 05:44/ 0 赞/ 251 阅读

相关 matlab矩阵

一、矩阵构造简单矩阵的构造：采用方括号\[\]，同行元素用逗号“，”隔开，行与行之间用“；”隔开特殊矩阵的构造函数： ones（n）构建n\n的矩阵，矩阵元素全为1；

以你之姓@/ 2022年06月14日 01:22/ 0 赞/ 345 阅读

相关 MATLAB对矩阵中元素的访问

[MATLAB矩阵的建立与访问][MATLAB] 矩阵访问一律用()，元胞数组用\{\} a(i,j)%第i行第j列队元素 a(:,j)%第j列的所有元素 a(

Bertha 。/ 2022年03月14日 11:50/ 0 赞/ 400 阅读