[請益] Abaqus 非元素節點 接觸 插值

看板Cad_Cae (電腦輔助設計)作者 (淯仔)時間6年前 (2018/02/22 14:03), 6年前編輯推噓4(4011)
留言15則, 2人參與, 6年前最新討論串1/1
各位大大好: 小弟在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

02/22 17:53, 6年前 , 1F
可以輸出surface上的contact force再平均到元素表面上。
02/22 17:53, 1F

02/24 15:51, 6年前 , 2F
如果一個元素的面有數個點,這樣值都一樣,我的模型的
02/24 15:51, 2F

02/24 15:52, 6年前 , 3F
網格相對於點的密度大多了,我也不希望細化模型網格
02/24 15:52, 3F

02/25 04:10, 6年前 , 4F
因為contact force的解是用牛頓法從節點上的力平衡求解來的
02/25 04:10, 4F

02/25 04:11, 6年前 , 5F
所以我猜CSTRESS只是把這個力平均到兩節點的Surface上而已
02/25 04:11, 5F

02/25 04:13, 6年前 , 6F
那ABAQUS本身應該沒有輸出任意點上這個平均值解的預設功能
02/25 04:13, 6F

02/25 04:13, 6年前 , 7F
除非你自己寫plug-in做
02/25 04:13, 7F

02/25 04:15, 6年前 , 8F
至於積分點沒值是因為contact求解不是從K矩陣來的,所以解
02/25 04:15, 8F

02/25 04:15, 6年前 , 9F
非來自積分點上。
02/25 04:15, 9F
的確如G大所說的,積分點是在元素內部,並不在表面,無法計算CANTACT FORCE,用積分 點做插值並不合理。 其實我想計算的是一個面壓在一個壓力毯上,壓力毯是由數個單位矩形組成的,並輸出每 個單位矩形上的反力來表示壓力。 但是用離散剛體作類似這樣的壓力毯相對於我的模型來說,元素數量實在太多且密,計算 成本高且收斂精度低。 而僅僅輸出各個模型節點上的CSTRESS又不太夠,似乎只能自己寫plug-in作插值 ※ 編輯: caron0225 (60.248.190.51), 02/26/2018 08:30:45

02/27 23:59, 6年前 , 10F
好奇問一下你要怎麼內插,如果是用相鄰節點的值做線性內插
02/27 23:59, 10F

02/28 00:02, 6年前 , 11F
那可以利用X-Y DATA建立分佈圖,然後勾選interpolate betwe
02/28 00:02, 11F

02/28 00:03, 6年前 , 12F
en point選項,就可以用probe選出你要的位置,ABAQUS會幫你
02/28 00:03, 12F

02/28 00:03, 6年前 , 13F
計算兩點之間的線性內插值
02/28 00:03, 13F

02/28 12:43, 6年前 , 14F
原本是打算用element的形狀函數去內插~
02/28 12:43, 14F

02/28 12:43, 6年前 , 15F
G大的這個方法似乎更方便!!我嘗試一下!!
02/28 12:43, 15F
文章代碼(AID): #1QZbql0c (Cad_Cae)
文章代碼(AID): #1QZbql0c (Cad_Cae)