作者caron0225 (淯仔)
看板Cad_Cae
標題[請益] Abaqus 非元素節點 接觸 插值
時間Thu Feb 22 14:03:25 2018
各位大大好:
小弟在abaqus上有些問題想請教~
我想要在接觸面上，輸出一些點的CSTRESS
但是這些點並不在元素的節點上
在abaqus中是否有內建的功能能這樣作呢??
起初，小弟有嘗試過建立與那些點配對的離散剛體並輸出其RF來代替，
但是如果我的點數量要是很多且密，相對離散剛體的元素量也很多且小，
導致計算成本增加，且因為過小的主面造成數值上的誤差增加。不是好辦法...
而解析剛體只能耦合一個點，只能輸出整體的RF。
剩下辦法應該就是利用元素節點的cstress值去內插指定點的值了。
雖然我使用的是一階元素，但是從ODB的結果來看非節點的值不是單純的線性內插??
還是應該要用積分點去外插??但是好像不能輸出積分點的CSTRESS??
還是有其他更好辦法??
希望能給小弟一些建議或參考!!感謝!!
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 60.248.190.51
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Cad_Cae/M.1519279407.A.026.html
1F：推 gamer: 可以輸出surface上的contact force再平均到元素表面上。 02/22 17:53
2F：→ caron0225: 如果一個元素的面有數個點，這樣值都一樣，我的模型的 02/24 15:51
3F：→ caron0225: 網格相對於點的密度大多了，我也不希望細化模型網格 02/24 15:52
4F：推 gamer: 因為contact force的解是用牛頓法從節點上的力平衡求解來的 02/25 04:10
5F：→ gamer: 所以我猜CSTRESS只是把這個力平均到兩節點的Surface上而已 02/25 04:11
6F：推 gamer: 那ABAQUS本身應該沒有輸出任意點上這個平均值解的預設功能 02/25 04:13
7F：→ gamer: 除非你自己寫plug-in做 02/25 04:13
8F：→ gamer: 至於積分點沒值是因為contact求解不是從K矩陣來的，所以解 02/25 04:15
9F：→ gamer: 非來自積分點上。 02/25 04:15
的確如G大所說的，積分點是在元素內部，並不在表面，無法計算CANTACT FORCE，用積分
點做插值並不合理。
其實我想計算的是一個面壓在一個壓力毯上，壓力毯是由數個單位矩形組成的，並輸出每
個單位矩形上的反力來表示壓力。
但是用離散剛體作類似這樣的壓力毯相對於我的模型來說，元素數量實在太多且密，計算
成本高且收斂精度低。
而僅僅輸出各個模型節點上的CSTRESS又不太夠，似乎只能自己寫plug-in作插值
※ 編輯: caron0225 (60.248.190.51), 02/26/2018 08:30:45
10F：推 gamer: 好奇問一下你要怎麼內插，如果是用相鄰節點的值做線性內插 02/27 23:59
11F：→ gamer: 那可以利用X-Y DATA建立分佈圖，然後勾選interpolate betwe 02/28 00:02
12F：→ gamer: en point選項，就可以用probe選出你要的位置，ABAQUS會幫你 02/28 00:03
13F：→ gamer: 計算兩點之間的線性內插值 02/28 00:03
14F：→ caron0225: 原本是打算用element的形狀函數去內插~ 02/28 12:43
15F：→ caron0225: G大的這個方法似乎更方便!!我嘗試一下!! 02/28 12:43