| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 846 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者挑戰(zhàn)平凡將贈(zèng)送您 30 個(gè)金幣 | ||
挑戰(zhàn)平凡銅蟲(chóng) (小有名氣)
|
[求助]
有限梁?jiǎn)卧獎(jiǎng)偠染仃嚧鎯?chǔ)的問(wèn)題
|
|
|
最近需要編寫(xiě)一個(gè)梁?jiǎn)卧挠?jì)算程序,單元?jiǎng)偠纫呀?jīng)得到,但是組集成總剛搞不定。我直接看的書(shū)是“有限單元法基礎(chǔ)及matlab編程”,4.5剛度矩陣存儲(chǔ)這一節(jié),他先找主對(duì)角元素的地址,然后組集了總剛剩下的一個(gè)下三角陣。前兩個(gè) 形成結(jié)點(diǎn)位移編碼矩陣和計(jì)算結(jié)點(diǎn)位移編碼我還能理解(但不確定理解的對(duì)不對(duì)),后面存儲(chǔ)對(duì)角線地址向量我就很難理解。。。如果直接手寫(xiě)的話我是能寫(xiě)出來(lái)矩陣如何匯總。請(qǐng)教各位。。。 %形成結(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ì)算存儲(chǔ)主對(duì)角元素地址 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 |
銅蟲(chóng) (小有名氣)
|
用matlab直接用滿陣就可以了,沒(méi)必要費(fèi)那個(gè)力氣去用半帶寬。看徐榮橋那本書(shū),可以直接用他的代碼。 發(fā)自小木蟲(chóng)Android客戶端 |
|
自己百度搜索里面的代碼,書(shū)絕版了,淘寶可以買(mǎi),徐榮橋老師的matlab有限元編程 發(fā)自小木蟲(chóng)Android客戶端 |
銅蟲(chóng) (小有名氣)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料與化工353分求調(diào)劑 +4 | NIFFFfff 2026-03-09 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +4 | 上班不著吉 2026-03-09 | 4/200 |
|
|
[考博] 2026博士申請(qǐng) +5 | 起泡酒 2026-03-08 | 5/250 |
|
|
[考研] 考研調(diào)劑,一志愿山東大學(xué)材料與化工,328分,政治51 +5 | 關(guān)你西紅柿929 2026-03-08 | 6/300 |
|
|
[考研] 337求調(diào)劑 +3 | 睡醒,。 2026-03-09 | 3/150 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
[基金申請(qǐng)] 進(jìn)入個(gè)人成果庫(kù)好難,一下午都沒(méi)進(jìn)去 +6 | mi_dilee 2026-03-05 | 6/300 |
|
|
[考研] 085701環(huán)境工程專(zhuān)業(yè),初試305,均過(guò)國(guó)家A區(qū)線 +7 | 卡卡來(lái)了@ 2026-03-07 | 8/400 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | G212 2026-03-03 | 6/300 |
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類(lèi)323分 誠(chéng)求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
[考研] 化學(xué)專(zhuān)業(yè)調(diào)劑 +5 | 好好好1233 2026-03-04 | 6/300 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財(cái) 2026-03-05 | 5/250 |
|
|
[考研] 0856材料與化工求調(diào)劑! +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
[考研] 271求調(diào)劑 +7 | 月色c 2026-03-05 | 8/400 |
|
|
[考研] 求調(diào)劑 +3 | 泡了個(gè)椒 2026-03-04 | 4/200 |
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
[論文投稿]
100+4
|
Stray2021 2026-03-03 | 4/200 |
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|