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