 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 837 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問題，作者挑戰(zhàn)平凡將贈送您 30 個金幣 | ||
挑戰(zhàn)平凡銅蟲 (小有名氣)
|
[求助]
有限梁單元剛度矩陣存儲的問題
|
|
|
最近需要編寫一個梁單元的計算程序，單元剛度已經(jīng)得到，但是組集成總剛搞不定。我直接看的書是“有限單元法基礎(chǔ)及matlab編程”，4.5剛度矩陣存儲這一節(jié)，他先找主對角元素的地址，然后組集了總剛剩下的一個下三角陣。前兩個 形成結(jié)點位移編碼矩陣和計算結(jié)點位移編碼我還能理解（但不確定理解的對不對），后面存儲對角線地址向量我就很難理解。。。如果直接手寫的話我是能寫出來矩陣如何匯總。請教各位。。。 %形成結(jié)點位移編碼矩陣 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 %計算結(jié)點位移編碼 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 |
回復(fù)此樓» 猜你喜歡
德國亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生
已經(jīng)有0人回復(fù)
-
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級好用，推薦試試
已經(jīng)有11人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有97人回復(fù)
-
2026-CJ開始申報了
已經(jīng)有1人回復(fù)
-
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有2人回復(fù)
挑戰(zhàn)平凡
銅蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 187.2
- 散金: 18
- 紅花: 1
- 帖子: 113
- 在線: 25.9小時
- 蟲號: 6689883
- 注冊: 2017-06-05
- 專業(yè): 結(jié)構(gòu)工程
2樓2022-02-27 20:16:33
3樓2022-02-27 23:02:47
|
自己百度搜索里面的代碼，書絕版了，淘寶可以買，徐榮橋老師的matlab有限元編程 發(fā)自小木蟲Android客戶端 |
» 本帖已獲得的紅花（最新10朵）
挑戰(zhàn)平凡4樓2022-02-27 23:04:12
挑戰(zhàn)平凡
銅蟲 (小有名氣)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 187.2
- 散金: 18
- 紅花: 1
- 帖子: 113
- 在線: 25.9小時
- 蟲號: 6689883
- 注冊: 2017-06-05
- 專業(yè): 結(jié)構(gòu)工程
5樓2022-03-02 21:21:20
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 挑戰(zhàn)平凡 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 環(huán)境工程學(xué)碩288求助調(diào)劑 +7 | 多吃億口芝士 2026-03-02 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工 298 調(diào)劑 +3 | 小西笑嘻嘻 2026-03-03 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學(xué)，化學(xué)專業(yè)344分，求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
|
[考研] 一志愿中科大能動297求調(diào)劑，本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
|
[考研] 278求調(diào)劑 +3 | 滿天星11_22 2026-03-02 | 3/150 |
|
|
|
[考研] 環(huán)境工程專碩307求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
|
[考研] 281電子信息求調(diào)劑 +3 | jhtfeybgj 2026-03-02 | 6/300 |
|
|
|
[考研] 求調(diào)劑院校 +6 | 云朵452 2026-03-02 | 8/400 |
|
|
|
[考博] 博士自薦 +4 | kkluvs 2026-02-28 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-02 | 3/150 |
|
|
|
[考研] 哈工大計算機(jī)劉劼團(tuán)隊招生 +4 | hit_aiot 2026-03-01 | 6/300 |
|
|
|
[基金申請] 成果系統(tǒng)訪問量大，請一小時后再嘗試。---NSFC啥時候好哦，已經(jīng)兩天這樣了 +4 | NSFC2026我來了 2026-02-28 | 4/200 |
|
|
|
[碩博家園] 博士自薦 +7 | 科研狗111 2026-02-26 | 11/550 |
|
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|
|
|
[論文投稿]
求助coordination chemistry reviews 的寫作模板
 10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|