| 5 | 1/1 | 返回列表 |
| 查看: 841 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問題,作者挑戰(zhàn)平凡將贈送您 30 個金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
挑戰(zhàn)平凡銅蟲 (小有名氣)
|
[求助]
有限梁單元剛度矩陣存儲的問題
|
|
|
最近需要編寫一個梁單元的計(jì)算程序,單元剛度已經(jīng)得到,但是組集成總剛搞不定。我直接看的書是“有限單元法基礎(chǔ)及matlab編程”,4.5剛度矩陣存儲這一節(jié),他先找主對角元素的地址,然后組集了總剛剩下的一個下三角陣。前兩個 形成結(jié)點(diǎn)位移編碼矩陣和計(jì)算結(jié)點(diǎn)位移編碼我還能理解(但不確定理解的對不對),后面存儲對角線地址向量我就很難理解。。。如果直接手寫的話我是能寫出來矩陣如何匯總。請教各位。。。 %形成結(jié)點(diǎn)位移編碼矩陣 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 %計(jì)算結(jié)點(diǎn)位移編碼 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 %計(jì)算存儲主對角元素地址 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 |
|
自己百度搜索里面的代碼,書絕版了,淘寶可以買,徐榮橋老師的matlab有限元編程 發(fā)自小木蟲Android客戶端 |
銅蟲 (小有名氣)
銅蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 298 +12 | 小西笑嘻嘻 2026-03-03 | 12/600 |
|
|---|---|---|---|---|
|
[考研]
|
15779376950 2026-03-01 | 7/350 |
|
|
[考研] 085600求調(diào)劑 +4 | LRZZZZZZ 2026-03-02 | 6/300 |
|
|
[考研] 289求調(diào)劑 +8 | yang婷 2026-03-02 | 10/500 |
|
|
[考研] 085601 材料工程 320 +3 | 和樂瑤 2026-03-03 | 3/150 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 中國林科院林化所(南京)2026年招收化學(xué)/材料/環(huán)境工程等背景碩士研究生3名 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
[基金申請] 沒有青基直接申請面上,感覺自己瘋了 +5 | kevin63t 2026-03-02 | 6/300 |
|
|
[考研] 0805總分292,求調(diào)劑 +12 | 幻想之殤 2026-03-01 | 12/600 |
|
|
[考研]
材料270求調(diào)劑
6+6
|
Eiiiio 2026-03-01 | 11/550 |
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +4 | 青科11 2026-03-02 | 4/200 |
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
[考研]
|
glwshine 2026-03-02 | 5/250 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
[考研] 材料類考研調(diào)劑 +6 | gemmgemm 2026-03-01 | 7/350 |
|
|
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5 | awaystay 2026-03-02 | 6/300 |
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
[考研] 一志愿華南理工大學(xué)材料與化工326分,求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
[高分子] 求環(huán)氧樹脂研發(fā)1名 +3 | 孫xc 2026-02-25 | 11/550 |
|