FRIC
關注創建者：deep 創建時間：2018-06-26
FRIC的視頻教程
Abaqus子程序開發系列課程-Fric_coef和Fric子程序
本次課程講解了Fric_coef和Fric子程序的編寫及調用過程，分別用兩個子程序實現了摩擦系數隨滑移率變化的仿真（新版本的子程序和老版本子程序接口發生了些許變化，需要注意）。 01_Fric_coef子程序開發教程 02_01_Fric子程序開發教程 03_摩擦子程序在磨損方面的簡單應用（相比umeshmotion簡單，磨損量比較大的情況下會失真。
免費 1小時42分鐘 805播放
查看
應用Fric子程序模擬粘結滑移
本課程主要講解了Fric的子程序的格式，參數輸入，以及編寫過程，并以實際例子講解子程序的計算過程，希望讀者能夠有所收獲，并能自己在本子程序的基礎上改寫子程序的，得到自己想要的結果，希望讀者看有所值！！
￥160 27分鐘 1609播放
查看
FRIC的實例教程
在調用FRIC之前，應力和滑動分量已經旋轉了這個量。該矩陣作為單位矩陣傳入二維和軸對稱接觸問題。
TEMP（2）
從節點和對應的主表面，此時的溫度
PREDEF（2，NPRED）
一個數組，包含當前增量步結束時的所有用戶指定場變量的一對值(分析開始時的初始值和分析期間的當前值)。從接觸對調用FRIC，一對值中的第一個值對應于從節點，第二個值對應于主表面上最近的點。如果從一個大滑動接觸單元調用FRIC，則PREDEF(1,NPRED)對應于該單元積分點處的值，PFREDEF(2,NPRED)對應于對應表面上的最近點。如果從一個小滑動接觸單元調用FRIC，則PREDEF(1,NPRED)對應于第一側積分點處的值，PFREDEF(2,NPRED)對應于該單元對應面上積分點處的值。
NFDIR
摩擦方向的數量
MCRD
接觸點處坐標方向的數量
NPRED
預定義場變量的數量
NSTATV
用戶定義狀態變量的數量
CHRLNGTH
接觸表面的面特征尺寸，可以用于定義最大允許的彈性滑動
PROPS（NPROPS）
用戶定義的屬性值數組，用來定義接觸表面間的摩擦行為
NPROPS
用戶定義的，與摩擦模型相關的屬性值的數量
展開 答案就是fric_coef子程序
二。fric_coef子程序介紹
該子程序定義接觸面的摩擦屬性，其標準格式如下所示：
參考上述標準文件格式和幫助文檔的案例，即可完成摩擦系數的定義。
文獻中給出的結果如下：
附件文獻：
30-Al-Qadi、汪浩的胎路摩擦論文.pdf
最后，大家有相關需求可以關注“320科技工作室”的微信公眾號，更多干貨等你來撩~~
" STATIC_FRIC_COEF="0.0048" STATIC_FRIC_LOAD="4.0" DYNAMIC_FRIC_COEF="0.0024" DYNAMIC_FRIC_LOAD="2.0"/>
<RESTRAINT.JOINT ID="35" NAME="AnkleR_joi" JOINT="AnkleR_jnt" STATIC_FRIC_COEF="0.0048" STATIC_FRIC_LOAD="4.0" DYNAMIC_FRIC_COEF="0.0024" DYNAMIC_FRIC_LOAD="2.0"/>
<INITIAL.JOINT_POS Q5="0.14" Q6="0.0" Q7="0.135" JOINT="Dummy_jnt" ORIENT="Dummy_ori"/>
<INITIAL.JOINT_POS Q1="0.0" JOINT="NeckBracket_jnt"/>
<INITIAL.JOINT_POS JOINT="NeckPivot1_jnt" ORIENT="NeckPivot1_ori"/>
<INITIAL.JOINT_POS JOINT="NeckPivot2_jnt" ORIENT="NeckPivot2_ori"/>
<INITIAL.JOINT_POS JOINT="NeckPivot3_jnt" ORIENT="NeckPivot3_ori"/>
<INITIAL.JOINT_POS JOINT="NeckPivot4_jnt" ORIENT="NeckPivot4_ori"/>
<INITIAL.JOINT_POS Q1="-0.11" Q2="0.25708" JOINT="ShoulderL_jnt"/>
<INITIAL.JOINT_POS Q1="-0.11" Q2="
展開 1 model linear property kn 2e7 ks 1e1 fric 0.5 ...
step 1000
save unaxial.sav
; output the results to file
his write -2 vs 1 file test1.dat
;his write -2 vs 1 file test2.dat
ret
幾何模型及計算過程中的壓縮曲線如圖1所示：
圖1 幾何模型及應力應變曲線
完成以下兩步以進行測試2：
(1) 注釋掉語句：pro den 2500 bulk _bk she _sh co 2e6 fric 45 ten 1e6 并uncomment以下語句：
pro den 2500 young 2e8 poisson 0.25 co 2e6 fric 45 ten 1e6
(2) 注釋掉語句：his write -2 vs 1 file test1.dat 并uncomment以下語句：his write -2 vs 1 file test2.dat
兩個實驗都計算了1000時步，對比計算結果如圖2所示：
圖2 實驗1與實驗2計算結果對比圖
曲線完全吻合。說明用K和G進行彈性參數的賦值，和用E和v進行賦值的效果是一樣的。但是就像手冊里說的那樣，用戶在平時使用時還是推薦通過K和G進行賦值，可以保證在包括熱力耦合或流固耦合時都不出現問題。
展開 
FRIC的相關專題、標簽、搜索
FRIC的最新內容
color: rgb(100, 106, 115);">Plain Text</span></p><p><span style="color: rgb(100, 106, 115);"> </span>*SURFACE PROPERTY, NAME=<表面磨損屬性名稱></p><p>*WEAR SURFACE PROPERTIES, UNITLESS WEAR COEF=YES/NO, FRIC
AUTO
2、創建功能與Patran界面功能相同的文本文件AUTO_CONTACT.bdf，通過INCLUDE 方式被主文件應用，示例如下：
上圖定義：自動創建接觸的屬性、單元等幾何等，常用接觸定義參數如下：
a) 接觸類型定義：ctype=touch，pglue定義接觸類型是否為粘接等
b) 可能接觸單元搜索距離：distance=1.0
c) 接觸對中摩擦系數：FRIC
K0固結后
model restore 'P0d20f0d3'
model calm
ball fix spin
[do_xSevro = true]
[do_ySevro = true]
[k0 = 0.43]
;[K0 = 1.0]
[tyy = -100.0e3]
[txx = tyy*k0]
[ballFrictionC = 0.3]
ball property 'fric
property pb_coh 10e6 pb_ten 10e6 pb_fa 50 fric 0.1 lin_mode 1cmat add 1 model linear method deform emod 10e9 kratio 1.5 ...
property pb_coh 10e6 pb_ten 10e6 pb_fa 50 fric 0.1 lin_mode 1cmat add 1 model linear method deform emod 10e9 kratio 1.5 ...
domain extent -2 2def RdPar rdmin=0.08 rdmax=0.12 poro=0.3end@RdParcmat default type ball-facet model linear property kn 1e8 cmat default type ball-ball model linear property kn 1e7 ks 1e7 fric
"kuaishi.stl"
rblock template create "kuaishi" from-geometry "kuaishi"
rblock replicate "kuaishi" position -121 331 80contact cmat default model linear method deform emod 10e7 kratio 1.5 property fric
property pb_coh [pb_coh] pb_ten [pb_coh*ten_coh] pb_fa 45.7 fric 0.1 contact cmat add 1 model linearpbond method deformability ...
model restore "rblock_slip"rblock export to-geometry "allgeo" cmat default type ball-facet model linear method deform emod 100e6 kratio 1.5 property fric 0.5cmat default type ball-ball model linear method
geometry generate cylinder base @x_pos1 @y_pos1 0 axis @vec height @height radius [jing_he*0.8] endcommandend@ball_range
cmat default model rrlinear method deformability emod 100e6 kratio 1.5 property fric
