 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 850 | 回復(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ù)此樓» 猜你喜歡
幫我的英語口語老師找學(xué)生
已經(jīng)有4人回復(fù)
-
計算三階力常數(shù)用thirdorder擴胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有172人回復(fù)
-
2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有1人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
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)工程
2樓2022-02-27 20:16:33
|
自己百度搜索里面的代碼，書絕版了，淘寶可以買，徐榮橋老師的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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[文學(xué)芳草園] 伙伴們，祝我生日快樂吧 +4 | myrtle 2026-03-10 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 334求調(diào)劑 +3 | mjqzz 2026-03-06 | 3/150 |
|
|
|
[考研] 一志愿中科大材料與化工。353分能調(diào)劑到哪？ +5 | MayUxw1 2026-03-08 | 7/350 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-09 | 11/550 |
|
|
|
[考研] 求調(diào)劑！ +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
|
[考研] 求調(diào)劑 一志愿蘇州大學(xué)，0856化工323分 | 本科應(yīng)化 | 有專利/競賽/科研助手經(jīng)歷 | +7 | 橙子cyx 2026-03-06 | 9/450 |
|
|
|
[考研]
 調(diào)劑的同學(xué)，走過路過，不要錯過
+11
|
likeihood 2026-03-06 | 16/800 |
|
|
|
[考研] 085701環(huán)境工程295求調(diào)劑，一志愿A區(qū)211 +8 | 等春來， 2026-03-06 | 9/450 |
|
|
|
[考研] 數(shù)二英二309分請求調(diào)劑 +3 | dtdxzxx 2026-03-09 | 4/200 |
|
|
|
[考研] 320求調(diào)劑 +4 | 魏zy 2026-03-08 | 4/200 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 很老實人 2026-03-09 | 4/200 |
|
|
|
[考研] 290 材料與化工求調(diào)劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
|
[考研] 269求調(diào)劑 +3 | 朔朔話 2026-03-08 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程317分求調(diào)劑 +9 | 6汆尼9 2026-03-07 | 9/450 |
|
|
|
[考研] 一志愿鄭大071000分數(shù)282求調(diào)劑 +3 | 研研顏 2026-03-05 | 7/350 |
|
|
|
[考研] 一志愿211 化學(xué)305分求調(diào)劑 +3 | 0703楊悅305分 2026-03-05 | 3/150 |
|
|
|
[考研] 一志愿武漢理工085601，初試301分，請問能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
|
[考研] 282求調(diào)劑 +7 | 夕～日 2026-03-05 | 8/400 |
|
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|