推 gamer: FEM的contact stess會跟接觸的節點數量有關,你上下模型 11/05 21:25
→ gamer: 的節點有沒有都對在一起? 11/05 21:25
→ j6zoz04: 回g大:感謝提醒,這點我沒有注意到,但不知道如何確認, 11/05 21:51
→ j6zoz04: 我只是依序在接觸的位置分佈目標和接觸元素 11/05 21:52
→ j6zoz04: 還是我應該用contact wizard? 但我記得裡面也沒有相關設 11/05 21:54
→ j6zoz04: 定能確認兩個元素是否有對到 11/05 21:54
推 gamer: 節點只能手動對齊而已 11/06 00:33
→ j6zoz04: 兩邊接觸區域附近的網格設置都相同,所以應該是有對齊, 11/06 17:44
→ j6zoz04: 我之後回去再確認,此外還有其它的建議嗎? 11/06 17:44
推 RASSIS: 建模的時候 share topology或者form new part(? 11/06 23:37
推 gamer: 還有就是接觸節點的層數,可以試試看只讓一層節點接觸 11/06 23:53
推 shin750926: 截圖一下局部接觸應力分佈,這樣比較清楚 11/07 02:06
※ 編輯: j6zoz04 (163.13.136.41 臺灣), 11/08/2021 14:35:29
→ j6zoz04: 回R大:不好意思,不了解什麼是share topology和form new 11/08 14:37
→ j6zoz04: part 11/08 14:37
→ j6zoz04: 回g大:想請問層數的意思,已把接觸對的圖片補上,我這樣 11/08 14:38
→ j6zoz04: 算是一層嗎? 11/08 14:38
→ j6zoz04: 回s大:已補上圖片 11/08 14:39
推 gamer: 我這樣看你就有7層的元素有接觸到 11/08 21:44
→ gamer: 我個人的猜測,範例可以求出接近解析解的結果應該是元素切 11/08 21:47
→ gamer: 比較粗,所以有接觸的節點層數比較少的關係,你可以試試看 11/08 21:47
→ gamer: 放大元素大小 11/08 21:47
→ j6zoz04: 不好意思,還是不太懂您的"層數"是指有初始接觸到的元素 11/09 01:18
→ j6zoz04: 數量嗎? 11/09 01:18
→ j6zoz04: 一般來說不是應該縮小元素大小才會接近解析解嗎? 11/09 01:19
→ j6zoz04: 剛剛發現我的接觸元素節點的contact pressure和y stress 11/09 17:51
→ j6zoz04: 居然不一致.. 11/09 17:51
推 gamer: 層數是指接觸到的節點數量(切線方向) 11/09 20:55
→ gamer: 由於我不知道你列的哪個解析解是怎麼解的,所以以下是我的 11/09 20:56
→ gamer: 猜測,理論上兩個圓的接觸面是一個無限小的點(以圓柱看就是 11/09 20:57
→ gamer: 一條線),所以此理論接觸壓力會是無窮大,我猜測這個解析解 11/09 20:57
→ gamer: 大概是用一個很小的面積去計算出來的。但你把圓用有限元素 11/09 20:58
→ gamer: 離散化後,這個面積大小會是你用了多少元素而定,當元素越 11/09 20:59
→ gamer: 細的時候接觸的節點會越多,而接觸壓力的計算方式是把節點 11/09 21:00
→ gamer: 上力平衡的結果,理論上同樣的力量如果越多節點接觸,這個 11/09 21:00
→ gamer: 壓力的極值會變小(因為力量平均的關係)。 11/09 21:02
→ j6zoz04: 了解,所以我應該把它這裡的接觸壓力理解為 外力/觸點數 11/09 21:28
→ j6zoz04: 而非外力/面積 11/09 21:29
→ j6zoz04: 我試著將元素減少後應力值會更小於理論解 11/09 21:32
→ j6zoz04: 元素由粗到細的應力趨勢是由小→大來接近理論解的 11/09 21:35
推 gamer: 等一下,你看的是應力還是接觸壓力? 11/09 21:37
→ gamer: 另外就是你看的是應力上的值,還是平均後的值 11/09 21:37
→ gamer: 節點上的值(上面打錯字) 11/09 21:38
→ j6zoz04: 我看的應力是節點上的應力(Y-Component Stress) 11/09 21:39
推 gamer: 應力跟contact pressure不一樣喔 11/09 21:41
→ gamer: 應力是二階張量,contact pressure是向量(一階張量) 11/09 21:43
→ gamer: 後處理裡面應該有另一個值是contact pressure 11/09 21:43
→ j6zoz04: 我看赫氏理論在接觸區域的contact pressure與 11/09 21:46
→ j6zoz04: y方向應力(surface stress)是一樣的 11/09 21:46
→ j6zoz04: 在後處理中,同個接觸節點的contact pressure與sy不相同 11/09 21:49
推 samuely: 物體不是剛體的話,接觸面積不會是無窮小。請問你知道這 11/13 17:16
→ samuely: 解析解是有沒有在小變形的假設下嗎? 11/13 17:16
推 samuely: 另外你這個是建2D模型嗎?是的話是plane stress 還plane 11/13 17:21
→ samuely: strain 11/13 17:21
→ j6zoz04: 回s大,此理論是以有小變形為前提的 11/13 18:03
→ j6zoz04: 我選擇使用plane strain 11/13 18:03
→ j6zoz04: 之後做了一些不同嘗試,我在想是不是與我網格劃分不過漂 11/13 18:08
→ j6zoz04: 亮有關 11/13 18:08
→ j6zoz04: 對,是2d模型 11/13 20:25
推 samuely: 這感覺是3D理論解,因為我看輸入的地方有長度可填,你建 11/13 22:25
→ samuely: 個3Dmodel或改plane stress看看,雖然我覺得可能結果更小 11/13 22:25
→ samuely: ...。如果想從mesh下手,這要做成幅射扇形的比較好 11/13 22:25
→ j6zoz04: 他輸入長度好像只是要換算成pressure而已,我看書上也是 11/14 02:43
→ j6zoz04: 使用2d來推導公式 11/14 02:43
→ j6zoz04: 我也曾想改成plane stress 但這例子我覺得plane strain才 11/14 02:44
→ j6zoz04: 對 11/14 02:45
→ j6zoz04: 目前希望是能用2d來跑,因為教育版的網格限制 11/14 02:46
→ j6zoz04: 想請問扇形網格的部分,目前只會用線段分割數或網格大小 11/14 02:56
→ j6zoz04: 來劃分網格。 11/14 02:56
→ j6zoz04: 附上plane stress的結果: 11/14 17:04
推 samuely: 看來改成plane stress就解了 11/26 00:58
→ j6zoz04: 但為何要用plane stress這點我不能理解就是 11/26 14:12