[問題] 有限元素法 任意四邊形板元素

看板Mechanical作者 (阿輝)時間9年前 (2015/06/11 22:52), 9年前編輯推噓7(7推 0噓 13→)
留言20則, 7人參與, 最新討論串1/1
各位前輩好 小弟我研究所做的研究是有關FEM及振動學 我作的分析是要自行寫程式(MATLAB) 求出剛性矩陣與質量矩陣 再用F=Kq 解出位移q 目前的Plate bending element板元素 有導出矩形元素及三角形元素 與ANSYS的振動頻率結果比較差異不大 因為我題目是做圓板及橢圓板 必須用三角形或任意四邊形切割 三角形元素做出來的結果與ANSYS差異不大 但我推導出的任意四邊形做出來的結果 在與ANSYS同樣Mesh下 答案比較起來有明顯的差距 任意四邊形是用自然座標系導的 Jocobian矩陣與一般四節點元素相同 高斯積分也取很多點了 看起來應該是收斂 找了很多書板的任意四邊形元素都很少提到 期刊跟論文找過很多 正在努力看懂@@ 有沒有前輩做FEM 希望可以指導一下> < PS. 用我推導的任意四邊形元素做矩形板的振動模態 MESH都切割成正方形 答案會與ANSYS接近 -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 1.34.51.102 ※ 文章網址: https://www.ptt.cc/bbs/Mechanical/M.1434034349.A.AD2.html ※ ahuie:轉錄至看板 Cad_Cae 06/11 22:53
推
06/11 23:00, , 1F
你是把元素切到多細?
06/11 23:00, 1F
→
06/11 23:11, , 2F
1000多個元素 再增加看起來答案也收斂了@@
06/11 23:11, 2F
→
06/11 23:13, , 3F
1200
06/11 23:13, 3F
推
06/11 23:23, , 4F
喔喔！元素量不多~ 那就跟初始設定有關了(大問題惹)
06/11 23:23, 4F
→
06/11 23:29, , 5F
QQ 我在想應該是形狀函數的問題
06/11 23:29, 5F
※ 編輯: ahuie (1.34.51.102), 06/11/2015 23:31:42
推
06/11 23:31, , 6F
任意四邊形沒啥涉獵的說 抱歉啦～
06/11 23:31, 6F
→
06/11 23:32, , 7F
沒關係QQ 還有很多時間可以做...
06/11 23:32, 7F
推
06/12 05:37, , 8F
有專版Cad_Cae
06/12 05:37, 8F
推
06/12 07:57, , 9F
你的元素編號是哪個？幾個節點？
06/12 07:57, 9F
→
06/12 07:58, , 10F
喔喔看到了4節點，你可以改成五節點看看
06/12 07:58, 10F
推
06/12 10:21, , 11F
cadcae很多高手
06/12 10:21, 11F
推
06/12 21:39, , 12F
其實簡單想你用四個點去舖圓形，你的邊界不會好看到哪
06/12 21:39, 12F
→
06/12 21:40, , 13F
三角點去舖面，自然與原來形狀近似度最高
06/12 21:40, 13F
→
06/12 22:45, , 14F
我一開始想法也跟樓上一樣 三角形本來就適合斜邊與曲面
06/12 22:45, 14F
→
06/12 22:46, , 15F
邊界分割 但是我們教授叫我試看看用四邊形割
06/12 22:46, 15F
→
06/12 22:46, , 16F
ANSYS用四邊形做答案也滿準 只是不知道他的方法...
06/12 22:46, 16F
→
06/14 23:32, , 17F
我不懂任意四邊形元素，不過標準元素三角形的收斂性遠差於
06/14 23:32, 17F
→
06/14 23:32, , 18F
四邊形，至於為什麼，修過有限元素法你就會知道了。和數值
06/14 23:32, 18F
→
06/14 23:33, , 19F
積分時轉換到高斯積分域的形狀函數有關，四邊形元素在轉換
06/14 23:33, 19F
→
06/14 23:33, , 20F
時的誤差會比三角形小。
06/14 23:33, 20F
‣ 更多 ahuie 的文章
文章代碼(AID): #1LUQ2jhI (Mechanical)