ansys提取接觸面反力的案例
LS_DYNA接觸力：分面提取與正面提取結果有區別？
自接觸使用*CONTACT_AUTOMATIC_SINGLE_SURFACE，包含所有shell單元，靜/動摩擦系數0.2。
同時建立segment set使用*FORCE_TRANSDUCER_PENALTY提取自接觸力，一個件分面提取，一個件整體提取：
接觸面3
ji
接觸面4
接觸面5
接觸面6
接觸面7整體提取接觸力
模型校核：
讀出來個接觸面接觸力如下
接觸面3、4、5碰撞過程中接觸力均為0，接觸面7的接觸力比接觸面6要大1.1%。
本人理解接觸面6與接觸面7碰撞條件相同，接觸力也應相等，數值差異是否可以看做是誤差。
展開 一點心得，接觸面節點力提取及再施加
這幾天做了個，需要將接觸面上的力提取做被接觸部件的進一步分析，而ABAQUS中沒有自動提取施加接觸力的功能，所以自己動手用matlab+inp的方式，實現了一個很原始的接觸力提取和施加。可能用戶自定義載荷（VDLoad）能夠實現，但是時間不夠沒有嘗試這方面的東西。
首先說一下，abaqus中的一個不是很完整的功能。
就是Analytical Field - Mapped field，Mapped field 是個很有意思的功能，它包含了兩種一個是Point cloud field，即用一組對應空間點的載荷，來表達載荷的空間分布，然后通過插值得到節點上載荷。
第二個叫做.odb mapped field（看到odb你懂的吧），這個就是自動提取odb文件中的載荷分布，這個的用法很簡單，在viewport菜單中新建一個viewport，叫做viewport2然后在這個viewport2里打開載荷存儲的odb文件，并打開你要提取的載荷的Field output，然后在viewport1中進入到.odb mapped field，有按鈕可選提取哪個viewport的分布。很簡單，很智能。
但是，還是要但是，在純力學的范疇內，這兩個功能都只能用在dload就是，pressure里，其余的載荷類型均無法識別。（除非只要接觸正壓強，切向的力完全弄不出來）
所以另辟道路，簡單的說吧，只提個方法好了！就是用編程語言改動后處理輸出的接觸面節點力的prt文件成為inp的格式，然后放入對應代碼模塊即可！！！
很簡單的思路，但是有一些注意的地方：
1. 節點對應，即網格劃分得相同；所以最好copy一個part，然后換成mesh part；
2. 注意提取接觸力時別弄錯了，接觸力有方向的問題，自己注意選擇；
3.
展開 根據命名自動提取并導出支反力_Ansys ACT Python ￥9.9
本文代碼，主要能實現的功能是：
1.根據name selection 定位到相關面，并施加reaction force
2. Evaluate 各個所需時間下的結果，并寫成list
3. 將list寫成一定規律導出到txt文件，以便后續處理
APDL應該也能實現此功能，這里單說說python的事情。
只需要定義面的name、時間、txt文件路徑、提取的數據類型，就可以自動化實現。
二 分析模型
提取支撐柱上和平板接觸的三個面的reaction force。
2.1 建立named selections用于程序參數識別；
2. 復制代碼，運行。
修改代碼main函數中的自定義參數，以及提取信息函數。Txt結果如下，
展開 結構仿真逆向邏輯：深度解析如何在 Ansys 中給定位移并精確提取支反力 ￥2
05 結語
在 Ansys Workbench 中，雖然沒有直接名為“全局方程”的模塊來求解這種“已知位移反求載荷”的問題，但通過 “位移約束 + 探針提取反力” 這一組合，我們可以更直觀地獲得等效結果。
思考拓展：
如果需要模擬彈簧在拉伸 2cm 后，再增加 100N 載荷的情況，僅用靜力學分析是不夠的，需要引入 Multi-Step 分析，即第一步強制位移 2cm，第二步鎖定位移并施加載荷。
