ZKX's LAB

你好,请问一下桁架结构带未知参数的刚度矩阵如何在matlab中实现? MATLAB怎么写总体刚度矩阵

2020-07-24知识15

matlab 求特征值问题 [K-M*W^2]=0 质量矩阵M 和 刚度矩阵K都。 matlab 求特征值问题[K-M*W^2]=0 质量矩阵M 和 刚度矩阵K都.matlab 求特征值问题[K-M*W^2]=0 质量矩阵M 和 刚度矩阵K都为3*3的矩阵(3自由度结构)怎么求W。.求用MATLAB编写空间梁的单元刚度矩阵和整体刚度矩阵的程序最近有一个作业是用MATLAB编写程序计算空间梁,单我不知道如何写求刚度矩阵的程序, matlab 求特征值问题 [K-M*W^2]=0 质量矩阵M 和 刚度矩阵K都为3*3的矩阵(3自由度结构)怎么求W。 你好,请问一下桁架结构带未知参数的刚度矩阵如何在matlab中实现? 你提出的问题我之前刚好做过,使用有限元方法来进行桁架结构分析.Matlab编程实现平面杆单元分析首先,明确Matlab程序要实现的5个重要模块分别为:单元刚度矩阵的求解、单元组装、节点位移的求解、单元应力的求解、节点力的求解.下面给出这5个模块的实现.1.\\x09单元刚度矩阵求解定义函数Bar2D2Node_Stiffness,该函数计算单元的刚度矩阵,输入弹性模量E,横截面积A,两个节点坐标输出单元刚度矩阵k(4X4).具体代码如下:function k=Bar2D2Node_Stiffness(E,A,x1,y1,x2,y2)L=sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));x=acos((x2-x1)/L);C=cos(x);S=sin(x);k=E*A/L*[C*C C*S-C*C-C*S;C*S S*S-C*S-S*S;C*C-C*S C*C C*S;C*S-S*S C*S S*S];2.\\x09单元组装定义函数Bar2D2Node_Assembly,该函数进行单元刚度矩阵的组装,输入单元刚度矩阵k,单元的节点编号i、j.输出整体刚度矩阵KK,具体代码如下:function z=Bar2D2Node_Assembly(KK,k,i,j)DOF(1)=2*i-1;DOF(2)=2*i;DOF(3)=2*j-1;DOF(4)=2*j;for n1=1:4for n2=1:4KK(DOF(n1),DOF(n2))=KK(DOF(n1),DOF(n2))+k(n1,n2);endendz=KK;3.\\x09节点位移的求解定义函数Bar2D2Node_Disp(KK,num,p),该函数输入KK为总体刚度矩阵;num。矩阵迭代幂法编程 已知质量矩阵和刚度矩阵如何用matlab求固有频率和对应的正则振型?求程序 关注者 8 被浏览 2,571 发现 等你来答 。请发送邮件到 jobs@zhihu.com求用MATLAB编写空间梁的单元刚度矩阵和整体刚度矩阵的程序 请参看《MATLAB有限元分析与应用》一书,该书作者P.I.Kattan,由韩来彬翻译,第7章介绍了空间梁单元建立刚度矩阵的MATLAB程序!单元刚度矩阵和整体刚度矩阵有什么特征 单元刚度矩阵特征:1、对称性2 奇异性3 主对角元素恒正4 所有奇数(偶数)行的和为 0结构刚度矩阵的特征:1、对称性2奇异性3主对角元素恒正4稀疏性5非零带状分布单元刚度矩阵集成为整体刚度矩阵怎么用MATLAB编程 %函数定义function K_assembly=Bar1D2Node_Assembly(K_elemnet)计算整体刚度矩阵计算单元数量ii=length(K_elemnet);K_assembly(1:1+i,1:1+i)=0;K_elemneti(1:2,1:2)=0;for j=1:i%依次取出单元刚度矩阵K_elemneti=K_elemnet(j);for n1=1:2for n2=1:2K_assembly(n1+j-1,n2+j-1)=K_assembly(n1+j-1,n2+j-1)+K_elemneti{1,1}(n1,n2);依次将单元矩阵加到整体矩阵中endendendK_assembly;知道刚度矩阵和质量矩阵后用matlab怎么求固有频率和固有振型 特征向量都求出来了,用哪一种归一化还不跟玩儿一样吗?okok已经回答你了,这里再贴一下,主要因为已经消失几天了,出来透口气:复制内容到剪贴板代码:<;/h5>;<;code id=\"code0\">;function ziyouzhendong1(k,m)%k=600*[1-1 0;1 3-2;0-2 5];m=diag([1 1.5 2]);clc[C,B]=eig(inv(m)*k);n=size(m,2);for i=1:nCt=C(:,i);zhenxing(:,i)=Ct/Ct(1);endw=sqrt(diag(B));B1=diag(B)';Eigen=inv(m)*k;for i=1:nt1(:,i)=B1(i)*zhenxing(:,i);t2(:,i)=Eigen*zhenxing(:,i);endB,zhenxing,t1,t2,w<;/code>;<;/div>;其实完全可以不用循环,这是上学初学MATLAB一个月左右时写的程序,目的是帮助同班美女从繁重的结构动力学作业中解脱出来并有充分的时间买化妆品and keep us amused.,当时还没人学MATLAB这玩意儿,全拿计算器算动力学矩阵的问题,显然即使最多只有4×4的K矩阵也足够让人受折磨了.转载求用MATLAB编写空间梁的单元刚度矩阵和整体刚度矩阵的程序最近有一个作业是用MATLAB编写程序计算空间梁,单我不知道如何写求刚度矩阵的程序,哪位高手请帮帮忙啊! 请参看《MATLAB有限元分析与应用》一书,该书作者P.I.Kattan,由韩来彬翻译,第7章介绍了空间梁单元建立刚度矩阵的MATLAB程序!

#刚度#matlab#单元刚度矩阵#矩阵#matlab矩阵

随机阅读

qrcode
访问手机版