 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 840 | 回復(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)有271人回復(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ù)
|
自己百度搜索里面的代碼，書絕版了，淘寶可以買，徐榮橋老師的matlab有限元編程 發(fā)自小木蟲Android客戶端 |
挑戰(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)工程
2樓2022-02-27 20:16:33
3樓2022-02-27 23:02:47
挑戰(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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854總分272 +5 | 打小就是老實人 2026-03-02 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
|
[考研] 085700資環(huán)求調(diào)劑，初始279，六級已過，英語能力強 +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
|
[考研] 085602化學(xué)工程350，調(diào)劑，有沒有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
|
[考研] 一志愿華中科技大學(xué)，化學(xué)專業(yè)344分，求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
|
[考研] 環(huán)境工程專碩307求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
|
[考研] 化學(xué)0703求調(diào)劑 學(xué)碩 理/工科均可 總分279 +3 | 1一11 2026-03-03 | 5/250 |
|
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
|
[考博] 26申博 +4 | north， 2026-02-28 | 4/200 |
|
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
|
[考研] 282求調(diào)劑 +4 | 2103240126 2026-03-02 | 7/350 |
|
|
|
[考研] 295求調(diào)劑。一志愿報考鄭州大學(xué)化學(xué)工藝學(xué)碩，總分295分 +8 | yl1 2026-03-02 | 9/450 |
|
|
|
[考研] 085600材料工程一志愿中科大總分312求調(diào)劑 +9 | 吃宵夜1 2026-02-28 | 11/550 |
|
|
|
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學(xué)校有行學(xué)上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +3 | 熬夜的貓頭鷹 2026-03-02 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +6 | 愛擦汗的可樂冰 2026-02-28 | 7/350 |
|
|
|
[考研] 275求調(diào)劑 +3 | L-xin? 2026-03-01 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|