| 查看: 770 | 回复: 4 | ||
| 【悬赏金币】回答本帖问题，作者挑战平凡将赠送您 30 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
有限梁单元刚度矩阵存储的问题
|
||
|
最近需要编写一个梁单元的计算程序，单元刚度已经得到，但是组集成总刚搞不定。我直接看的书是“有限单元法基础及matlab编程”，4.5刚度矩阵存储这一节，他先找主对角元素的地址，然后组集了总刚剩下的一个下三角阵。前两个 形成结点位移编码矩阵和计算结点位移编码我还能理解（但不确定理解的对不对），后面存储对角线地址向量我就很难理解。。。如果直接手写的话我是能写出来矩阵如何汇总。请教各位。。。 %形成结点位移编码矩阵 function nf = formnf (nf) [nRow , nCol] = size(nf); m=0; for j=1:nCol for i=1:nRow if nf(i,j)~=0 m= m + 1; nf(i,j) = m; end end end %计算结点位移编码 function g= numtog (num,nf) nod = length(num); nodof = size(nf,1); for i = 1:nod k = i*nodof; g(k-nodof+1:k) = nf(:,num(i)); end end %计算存储主对角元素地址 function kdiag = fkdiag (kdiag,g) idof = size(g,2); for i= 1:idof iwp1 = 1; if g(i) ~= 0 for j = 1:idof if g(j) ~= 0 im = g(i) - g(j) +1; if im > iwp1 iwp1 = im; end end end k = g(i); if iwp1 > kdiag(k) kdiag(k) = iwp1; end end end end %集装总体刚度矩阵 function kv = fsparv(kv,km,g,kdiag) idof = length(g); for i = 1:idof k=g(i); if k~= 0 for j =1:idof if g(j)~= 0 iw = k- g(i); if iw >= 0 ival = kdiag(k) - iw kv(ival) = kv(ival) + km(i,j); end end end end end end |
回复此楼» 猜你喜欢
求助VISSIM破解版软件
已经有0人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有155人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有17人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有0人回复
送红花一朵
5楼2022-03-02 21:21:20
2楼2022-02-27 20:16:33
3楼2022-02-27 23:02:47
4楼2022-02-27 23:04:12