PTT數位生活區 / Cad_Cae (電腦輔助設計)

[問題] CAE固體力學之穩態分析，無法收斂?

看板Cad_Cae (電腦輔助設計)作者 (阿泰斯)時間10年前 (2014/05/02 18:09), 編輯推噓17(17推 0噓 102→)
留言119則, 4人參與, 最新討論串1/1
我用的是Comsol 4.3b 是選固體力學模組 分析一個鐵塊，其中一面受10Pa壓力，然後研究選"穩態"進行分析 發現不管怎麼分析，都會不收斂，網格切越細，誤差就越大... 請問怎麼會這樣? -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 175.180.80.71 ※ 文章網址: http://www.ptt.cc/bbs/Cad_Cae/M.1399025386.A.A4C.html
→
05/02 18:13, , 1F
http://ppt.cc/o~HA 這是上面說的的Comsol4.3b存檔
05/02 18:13, 1F
→
05/02 18:29, , 2F
但是找特徵頻率或掃頻卻沒問題，可以分析。
05/02 18:29, 2F
→
05/02 18:45, , 3F
請問穩態不收斂，去分析特徵頻率或掃頻可嗎？
05/02 18:45, 3F
→
05/03 01:01, , 4F
奇異點?
05/03 01:01, 4F
推
05/03 01:26, , 5F
你邊界條件怎麼設定?
05/03 01:26, 5F
→
05/03 01:28, , 6F
如果是有一面整個固定住，哪在固定面的邊上會是奇異點，應力
05/03 01:28, 6F
→
05/03 01:28, , 7F
為無窮大，你切越細，解就越逼近無窮大。
05/03 01:28, 7F
→
05/03 01:29, , 8F
由於只是應力奇異點，所以做模態並不會影響正確性。
05/03 01:29, 8F
→
05/03 01:42, , 9F
模型是一個長方體，最大的兩個面垂直y軸，其餘的面分別
05/03 01:42, 9F
→
05/03 01:42, , 10F
垂直x和z軸，而邊界條件只有給最大的那兩個面，一個
05/03 01:42, 10F
→
05/03 01:42, , 11F
電極，一個接地。
05/03 01:42, 11F
→
05/03 01:43, , 12F
抱歉，邊界條件說錯了，是最大的那兩個面的其中一個面
05/03 01:43, 12F
→
05/03 01:44, , 13F
，給10Pa的壓力。
05/03 01:44, 13F
→
05/03 01:45, , 14F
可是我只給一個面10Pa的壓力，就只有這個邊界條件而已
05/03 01:45, 14F
→
05/03 01:45, , 15F
沒給任何固定端。
05/03 01:45, 15F
推
05/03 09:26, , 16F
邊界條件這樣設定不對吧
05/03 09:26, 16F
→
05/03 13:04, , 17F
請問要怎設定？
05/03 13:04, 17F
→
05/03 13:38, , 18F
忘了說，comsol內建都會幫我們設定"自由邊界條件"，當
05/03 13:38, 18F
→
05/03 13:39, , 19F
我把其中一面設定壓力，那一個面的自由邊界條件就會
05/03 13:39, 19F
→
05/03 13:39, , 20F
失效。 所以嚴格來講，我的模型有一個面是壓力，其它都
05/03 13:39, 20F
→
05/03 13:40, , 21F
自由。
05/03 13:40, 21F
推
05/03 13:54, , 22F
所以才不對 靠地面的面要設固定面
05/03 13:54, 22F
→
05/03 14:03, , 23F
感謝！但是我用理論下去算的邊界都是自由耶，只有一面
05/03 14:03, 23F
→
05/03 14:03, , 24F
是壓力。
05/03 14:03, 24F
推
05/03 14:06, , 25F
你怎麼算的 有計算式的過程嗎
05/03 14:06, 25F
→
05/03 14:09, , 26F
我是算特徵頻率，就很簡單的偏微分方程式而已。
05/03 14:09, 26F
→
05/03 14:10, , 27F
另外，理論會假設其他面都是自由邊界的原因是實體產品
05/03 14:10, 27F
→
05/03 14:11, , 28F
的設計也都是以"盡量懸空"的方式設計，整個實體幾乎是
05/03 14:11, 28F
→
05/03 14:11, , 29F
懸空，只有一兩個固定的點，那兩個點都在實體邊緣，
05/03 14:11, 29F
→
05/03 14:11, , 30F
佔整體的面積相當的少。
05/03 14:11, 30F
推
05/03 14:12, , 31F
不是少的問題 有固定點的邊界條件就是要放進去
05/03 14:12, 31F
→
05/03 14:12, , 32F
試想 如果沒有那兩個固定點 產品不就掉下去了嗎
05/03 14:12, 32F
→
05/03 14:13, , 33F
可是理論就是給它自由邊界阿...，若CAE給固定，就無法
05/03 14:13, 33F
→
05/03 14:13, , 34F
跟理論比對了...
05/03 14:13, 34F
→
05/03 14:14, , 35F
如果真實情況是有固定點 為何理論要忽略固定點?
05/03 14:14, 35F
→
05/03 14:14, , 36F
因為總總原因...，理論無法計算到完全真實的狀況，但是
05/03 14:14, 36F
→
05/03 14:15, , 37F
會先用半真實的狀況和CAE比對，確定CAE設計的模型是對
05/03 14:15, 37F
→
05/03 14:15, , 38F
的。
05/03 14:15, 38F
→
05/03 14:16, , 39F
而且現在比較奇怪的是，"穩態分析"會發散出錯，
05/03 14:16, 39F
還有 40 則推文
→
05/04 01:10, , 80F
狀況。至於設不設固定邊界條件也不是重點，重點是你模型不能
05/04 01:10, 80F
→
05/04 01:11, , 81F
有鋼體運動，不然就變成動態問題，用固力的方式是找不出解的
05/04 01:11, 81F
→
05/04 01:12, , 82F
至於動態還是靜力分析那個才是正確的，這就得看你的問題到底
05/04 01:12, 82F
→
05/04 01:12, , 83F
是什麼，不先把自己要求的問題弄清楚，而只是把一堆參數輸入
05/04 01:12, 83F
→
05/04 01:13, , 84F
軟體，這叫做garbage in, garbage out.
05/04 01:13, 84F
→
05/04 01:20, , 85F
感謝!!!但是怪怪的...，我用理論解特徵頻率，要用到
05/04 01:20, 85F
→
05/04 01:20, , 86F
governing eq耶...
05/04 01:20, 86F
→
05/04 01:22, , 87F
governing eq解特徵值與特徵向量，再給定邊界條件，
05/04 01:22, 87F
→
05/04 01:23, , 88F
可以計算出一些掃頻的資訊。
05/04 01:23, 88F
→
05/04 01:26, , 89F
另外，您說的"動態問題"是指"給外力後它會振動"還是
05/04 01:26, 89F
→
05/04 01:26, , 90F
"分析原本就會轉動的齒輪"?
05/04 01:26, 90F
推
05/04 14:59, , 91F
解G.E.跟解出G.E.的特徵根是兩回事，簡化一點說，解微分方程
05/04 14:59, 91F
→
05/04 15:00, , 92F
需要積分得到方程式的形式，解特徵根只要解特徵問題即可。
05/04 15:00, 92F
→
05/04 15:01, , 93F
而有沒有固定邊界條件，基本上都可以解出特徵根和特徵向量
05/04 15:01, 93F
→
05/04 15:02, , 94F
但是當問題為靜力分析問題時，模型不固定的情況，會變成平衡
05/04 15:02, 94F
→
05/04 15:02, , 95F
方程式沒辦法滿足，缺少平衡方程式，等於缺少獨立方程式，解
05/04 15:02, 95F
→
05/04 15:03, , 96F
不出獨立變數。
05/04 15:03, 96F
→
05/04 15:04, , 97F
動態問題泛指的所有非靜力平衡問題，以你設的條件來說，一塊
05/04 15:04, 97F
→
05/04 15:05, , 98F
無拘束方塊受到壓力會怎麼樣?牛頓力學告訴我們他會往壓力的
05/04 15:05, 98F
→
05/04 15:06, , 99F
方向移動，以F=ma的加速度。所以你這樣的條件設定其實是在解
05/04 15:06, 99F
→
05/04 15:06, , 100F
動力學問題，那你用固體力學的方法去解，當然解不出來。
05/04 15:06, 100F
→
05/04 16:00, , 101F
感謝!!! 那我要解的是動態問題，就是F=ma，
05/04 16:00, 101F
→
05/04 16:00, , 102F
a=u對時間微分2次，然後最後求的是u。
05/04 16:00, 102F
→
05/04 16:01, , 103F
但是comsol裡面好像沒有分動態跟靜態@@
05/04 16:01, 103F
→
05/04 16:01, , 104F
ansys好像有分"靜態結構"跟"動態"這兩種分析
05/04 16:01, 104F
→
05/04 16:04, , 105F
另外，對於comsol來說，固體力學應該跟"靜態問題"無關
05/04 16:04, 105F
→
05/04 16:05, , 106F
固體力學就只是"固體"的"力學"而已@@
05/04 16:05, 106F
推
05/04 20:47, , 107F
你現在設定的條件是 非穩態 非力平衡
05/04 20:47, 107F
→
05/04 20:47, , 108F
你用力平衡的方法去解 當然解不出來
05/04 20:47, 108F
→
05/04 20:50, , 109F
如果你要解的是動態問題 那你就不能用穩態的方式去解
05/04 20:50, 109F
→
05/04 20:50, , 110F
要用暫態的方式去解
05/04 20:50, 110F
→
05/04 20:52, , 111F
暫態比穩態要花的時間會很多 穩態只要一次收斂就收工
05/04 20:52, 111F
→
05/04 20:53, , 112F
暫態每一個 step 都要收斂 才能進行下一步
05/04 20:53, 112F
推
05/04 21:02, , 113F
comsol 一定會有穩態跟暫態兩種方法
05/04 21:02, 113F
→
05/04 22:12, , 114F
comsol3.5有暫態可以選，4.0以後好像就沒那個選項了
05/04 22:12, 114F
→
05/04 22:13, , 115F
4.0以後的暫態應該就是"時間相依"，分析一段時間，看
05/04 22:13, 115F
→
05/04 22:14, , 116F
每秒的暫態資訊。
05/04 22:14, 116F
→
05/04 22:23, , 117F
非常感謝各位高手幫助!!!
05/04 22:23, 117F
→
05/04 22:24, , 118F
特徵頻率跟掃頻跟時間無關，故模型沒有靜力平衡應該沒
05/04 22:24, 118F
→
05/04 22:25, , 119F
差。
05/04 22:25, 119F
‣ 更多 candy88257 的文章
文章代碼(AID): #1JOsxgfC (Cad_Cae)
Cad_Cae 近期熱門文章
更多 近期熱門文章 >>
PTT數位生活區 即時熱門文章
更多 即時熱門文章 >>
‣ 更多 candy88257 的文章
文章代碼(AID): #1JOsxgfC (Cad_Cae)