FRIC的案例
Abaqus子程序系列：FRIC（定義接觸表面的摩擦行為）
在調(diào)用FRIC之前，應(yīng)力和滑動(dòng)分量已經(jīng)旋轉(zhuǎn)了這個(gè)量。該矩陣作為單位矩陣傳入二維和軸對(duì)稱接觸問(wèn)題。
TEMP（2）
從節(jié)點(diǎn)和對(duì)應(yīng)的主表面，此時(shí)的溫度
PREDEF（2，NPRED）
一個(gè)數(shù)組，包含當(dāng)前增量步結(jié)束時(shí)的所有用戶指定場(chǎng)變量的一對(duì)值(分析開(kāi)始時(shí)的初始值和分析期間的當(dāng)前值)。從接觸對(duì)調(diào)用FRIC，一對(duì)值中的第一個(gè)值對(duì)應(yīng)于從節(jié)點(diǎn)，第二個(gè)值對(duì)應(yīng)于主表面上最近的點(diǎn)。如果從一個(gè)大滑動(dòng)接觸單元調(diào)用FRIC，則PREDEF(1,NPRED)對(duì)應(yīng)于該單元積分點(diǎn)處的值，PFREDEF(2,NPRED)對(duì)應(yīng)于對(duì)應(yīng)表面上的最近點(diǎn)。如果從一個(gè)小滑動(dòng)接觸單元調(diào)用FRIC，則PREDEF(1,NPRED)對(duì)應(yīng)于第一側(cè)積分點(diǎn)處的值，PFREDEF(2,NPRED)對(duì)應(yīng)于該單元對(duì)應(yīng)面上積分點(diǎn)處的值。
NFDIR
摩擦方向的數(shù)量
MCRD
接觸點(diǎn)處坐標(biāo)方向的數(shù)量
NPRED
預(yù)定義場(chǎng)變量的數(shù)量
NSTATV
用戶定義狀態(tài)變量的數(shù)量
CHRLNGTH
接觸表面的面特征尺寸，可以用于定義最大允許的彈性滑動(dòng)
PROPS（NPROPS）
用戶定義的屬性值數(shù)組，用來(lái)定義接觸表面間的摩擦行為
NPROPS
用戶定義的，與摩擦模型相關(guān)的屬性值的數(shù)量
展開(kāi) Abaqus中利用fric_coef子程序定義摩擦系數(shù)與速度的關(guān)系
答案就是fric_coef子程序
二。fric_coef子程序介紹
該子程序定義接觸面的摩擦屬性，其標(biāo)準(zhǔn)格式如下所示：
參考上述標(biāo)準(zhǔn)文件格式和幫助文檔的案例，即可完成摩擦系數(shù)的定義。
文獻(xiàn)中給出的結(jié)果如下：
附件文獻(xiàn)：
30-Al-Qadi、汪浩的胎路摩擦論文.pdf
最后，大家有相關(guān)需求可以關(guān)注“320科技工作室”的微信公眾號(hào)，更多干貨等你來(lái)撩~~
安全帶預(yù)模擬(相對(duì)簡(jiǎn)單 自己編的)
" 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="
展開(kāi) Range fish的使用
1 model linear property kn 2e7 ks 1e1 fric 0.5 ...
FLAC3D中K/G 與 E/v之間的關(guān)系及其它
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
幾何模型及計(jì)算過(guò)程中的壓縮曲線如圖1所示：
圖1 幾何模型及應(yīng)力應(yīng)變曲線
完成以下兩步以進(jìn)行測(cè)試2：
(1) 注釋掉語(yǔ)句：pro den 2500 bulk _bk she _sh co 2e6 fric 45 ten 1e6 并uncomment以下語(yǔ)句：
pro den 2500 young 2e8 poisson 0.25 co 2e6 fric 45 ten 1e6
(2) 注釋掉語(yǔ)句：his write -2 vs 1 file test1.dat 并uncomment以下語(yǔ)句：his write -2 vs 1 file test2.dat
兩個(gè)實(shí)驗(yàn)都計(jì)算了1000時(shí)步，對(duì)比計(jì)算結(jié)果如圖2所示：
圖2 實(shí)驗(yàn)1與實(shí)驗(yàn)2計(jì)算結(jié)果對(duì)比圖
曲線完全吻合。說(shuō)明用K和G進(jìn)行彈性參數(shù)的賦值，和用E和v進(jìn)行賦值的效果是一樣的。但是就像手冊(cè)里說(shuō)的那樣，用戶在平時(shí)使用時(shí)還是推薦通過(guò)K和G進(jìn)行賦值，可以保證在包括熱力耦合或流固耦合時(shí)都不出現(xiàn)問(wèn)題。
展開(kāi) 【PFC6.0】三維Cluster模擬混合料混凝土
property pb_coh 20e6 pb_ten 20e6 pb_fa 50 fric 0.1 lin_mode 1 range group "Rock3" match 2cmat add 5 model linearpbond method deform emod 8e9 kratio 1.5 pb_deform emod 8e9 kratio 1.5 ...
參賽| 煤礦開(kāi)采
change mat 2 range 0,200 15,17 ;底板老底粉砂巖
change mat 3 range 0,200 17,19 ;直接底泥巖
change mat 4 range 0,200 19,24 ;13-1煤
change mat 3 range 0,200 24,26 ;頂板直接頂泥巖
change mat 2 range 0,200 26,32 ;頂板老頂粉砂巖
change mat 2 range 0,200 32,37 ;粉砂巖
change mat 5 range 0,200 37,42 ;細(xì)砂巖
change mat 6 range 0,200 42,48 ;砂質(zhì)泥巖
change mat 5 range 0,200 48,54 ;細(xì)砂巖
change mat 6 range 0,200 54,90 ;砂質(zhì)泥巖
prop mat 1 d=2550 bulk=26.23e9 shear=9.4e9 fric=34 coh=23.13e6 ten=0
prop mat 2 d=2450 bulk=12.2e9 shear=9.92e9 fric=31 coh=17.64e6 ten=0
prop mat 3 d=2500 bulk=14.17e9 shear=5.43e9 fric=35 coh=8.25e6 ten=0
prop mat 4 d=13500 bulk=4.31e9 shear=0.73e9 fric=36 coh=2.96e6 ten=0
prop mat 5 d=2500 bulk=16.67e9 shear=7.69e9 fric=35 coh=20.8e6 ten=0
prop
展開(kāi) PFC中成樣方法簡(jiǎn)單介紹四種
new
domain extent -10 10
wall generate box -5 5 expand 1.5ball generate number 2291 radius 0.05 box -5 5
ball attribute radius 0.1ball attribute density 2.7e3 damp 0.7
cmat default model linear method deformability emod 100e6 kratio 1.5 property fric 0.5
cycle 2000 calm 50
solve
save sampleFromPengZhang
3、壓縮法
這里一樣是用generate，這里不是把半徑減小，而是將區(qū)域擴(kuò)大，然后壓縮的方法。
new
domain extent -20 20
wall generate box -5 5 -5 15 expand 1.5ball generate number 2291 radius 0.1 box -5 5 -5 15
ball attribute density 2.7e3 damp 0.7
cmat default model linear method deformability emod 100e6 kratio 1.5 property fric 0.5
wall attribute yvelocity -2 range id 3
solve time 5wall attribute yvelocity 0 range id 3
solve
save sampleFromYasuo
可以看到下部的顆粒分布并不是很好。
展開(kāi) [轉(zhuǎn)帖]流體與動(dòng)力的耦合計(jì)算-大壩
;接觸面參數(shù)
interface 2 prop coh=0 fric 60. dil 0. kn=1e8 ks=1e8 ten 0.
model fl_iso;各項(xiàng)同性流體模型,計(jì)算流體必須的.
; mechanical boundary and initial conditions物理邊界和初始條件
fix z range z=-.1 .1;固定z=0的面
fix x range x=-.1 .1
fix x range x=39.9 40.1;固定x=40的面
fix y range y=-.1 .1
fix y range y=9.9 10.1
;初始應(yīng)力,垂直應(yīng)力為水平的2倍,在z方向有梯度變化
ini szz -3.6e5 grad 0 0 1.8e4
ini sxx -1.8e5 grad 0 0 0.9e4
ini syy -1.8e5 grad 0 0 0.9e4
set grav 0 0 -10;設(shè)置重力加速度
;記錄監(jiān)測(cè)數(shù)據(jù)
hist unbal;不平衡力
hist gp zdisp 16,5,20;點(diǎn)(16,5,20)的z方向位移
hist gp zdisp 30,5,20
;求解
solve
save block1.sav;保存
; assign realistic strength properties設(shè)置現(xiàn)實(shí)的強(qiáng)度參數(shù)
prop cohes=0 tens 0 fric 35 range group 'soil embankment'
prop cohes=0 tens 0 fric 40 range group 'dense soil foundation'
solve
save block1.sav
;流體分析
rest block1.sav;調(diào)用保存的文件
; specify fluid properties設(shè)定流體參數(shù)
prop perm
展開(kāi) 【PFC6.0.30】三維Cluster模擬GBM礦物晶粒巖石單軸
property pb_coh 20e6 pb_ten 20e6 pb_fa 50 fric 0.1 lin_mode 1 range fish @isRock3In cmat add 5 model linearpbond method deform emod 8e9 kratio 1.5 pb_deform emod 8e9 kratio 1.5 ...
PFC休止角（堆積角）計(jì)算的滾筒方法——砂土顆粒
plane o 0 0 dip 0 belowcmat default model linear method deformability emod 100e6 kratio 1.74 property fric 0.5
ball attribute density 2.7e3 damp 0.7
cycle 2000 calm 50
set gravity 9.8
solve
save sample
以上為得到初始狀態(tài)所需的代碼。
然后簡(jiǎn)單的給圓形容器一個(gè)滾動(dòng)速度便可以了，注意墻和球之間的摩擦效率要很大，這樣模擬效率高。
『分享』預(yù)應(yīng)力錨桿、錨索的模擬！
這兩個(gè)參數(shù)的確定，是應(yīng)該按抗拔試驗(yàn)來(lái)確定的，將試驗(yàn)錨桿的最大抗拔力/錨固長(zhǎng)度作為y軸，圍巖壓力*錨固體直徑為x軸，截距就是C，軟件中是gr_coh，傾角就是gr_fric。實(shí)際應(yīng)用中，軟件可以采用gr_coh等于D（鋼筋拔出）或者D+2t（錨固體拔出）乘于pi再乘于max（圍巖和錨固砂漿的單軸抗壓強(qiáng)度）/2來(lái)計(jì)算；如果不考慮圍壓的影響，gr_fric可以設(shè)置為0。
預(yù)應(yīng)力的施加，俺認(rèn)為使用sel cable pretension較好，直接加在錨固段上，錨頭嘛，你應(yīng)該和梁連在一起．
PFC休止角（堆積角）滾筒計(jì)算方法——碎石
;plane o 0 0 dip 0 belowclump attribute density 3e3 damp 0.7
cmat default model linear method deformability emod 100e6 kratio 1.74 property fric 0.5cycle 2000 calm 50
set gravity 9.8
solvesave sample
clump attribute displacement multiply 0
cmat add model linear method deformability emod 100e6 kratio 1.74 property fric 1.5 range contact type pebble-facetset mech age 0
[time_record=mech.age]
[baocunpinlv=2]
def savefile
if mech.age-time_record > baocunpinlv then
filename=string.build("jieguo%1",count)
command
save @filename
endcommand
time_record=mech.age
count +=1
endif
end
set fish callback -1.0 @savefilewall attribute centrotation 0 0 spin 50
solve time 70
這里將稻谷和此處的輪廓圖附上，我這里是用adobe illustrate畫(huà)的二維圖
展開(kāi) PFC砂土柔性雙軸
restore ball_sample[fric_shiyang=0.5][rrfric_shiyang=0.2][emod_shiyang=100e6]ball group shiyang cmat add 1 model rrlinear method deformability emod @emod_shiyang kratio 1.5 property fric @fric_shiyang rr_fric @rrfric_shiyang ...
PFC休止角（堆積角）滾筒計(jì)算——稻谷
;plane o 0 0 dip 0 belowclump attribute density 3e3 damp 0.7
cmat default model linear method deformability emod 100e6 kratio 1.74 property fric 0.5cycle 2000 calm 50
set gravity 9.8
solvesave sample
clump attribute displacement multiply 0
cmat add model linear method deformability emod 100e6 kratio 1.74 property fric 1.5 range contact type pebble-facetset mech age 0
[time_record=mech.age]
[baocunpinlv=2]
def savefile
if mech.age-time_record > baocunpinlv then
filename=string.build("jieguo%1",count)
command
save @filename
endcommand
time_record=mech.age
count +=1
endif
end
set fish callback -1.0 @savefilewall attribute centrotation 0 0 spin 50
clump delete range y -10 -1
solve time 300
注意clump計(jì)算是非常吃計(jì)算速度的
展開(kāi) 