PFC3D建模流程和代码解释(2)

; fname: cmat2.p2dat ;

; Illustrate simple usage of the CMAT

; (Default slots filled with the different data)

;============================================================================== new

set random 10001

domain extent -1.5e-2 1.5e-2 condition destroy wall generate box -1.0e-2 1.0e-2 onewall

ball generate number 200 radius 0.4e-3 0.6e-3 box -1.0e-2 1.0e-2

ball attribute density 2500.0;生成了和cmat1一样的模型，只是在接触模型上有了变化。

; set different contact models for the ball-ball and ball-facet contact types cmat default type ball-ball ...

model hertz ...;赫兹接触模型，颗粒与颗粒之间

property hz_shear 30e9 hz_poiss 0.3 ...;赫兹切向模量30的9次方，赫兹泊松比为0.3

fric 0.25 ...;摩擦系数0.25 dp_nratio 0.2;法向临界阻尼比为0.2

cmat default type ball-facet ... model linear ...

property kn 1e6 dp_nratio 0.2

set gravity 9.81 solve

save cmat2

return;对两种接触模型分别定义

;============================================================================== ; eof: cmat2.p2dat 例子三

; fname: cmat3.p2dat ;

; Illustrate advanced usage of the CMAT

; (Default and Optional slots filled with different data)

;============================================================================== new

set random 10001

domain extent -1.5e-2 1.5e-2 condition destroy;生成了两组颗粒，在ball generate当中用了

group命令

wall generate box -1.0e-2 1.0e-2 onewall

ball generate number 100 radius 0.4e-3 0.6e-3 box -1.0e-2 1.0e-2 group glass ;玻璃球 ball generate number 100 radius 0.4e-3 0.6e-3 box -1.0e-2 1.0e-2 group steel ;钢球

ball attribute density 2500.0 range group glass;定义密度 ball attribute density 7500.0 range group steel

;contact group left range x -1.5e-2 0 ;contact group right range x 0 1.5e-2

contact groupbehavior and ;contact groupbehavior后面有三个关键词，and，or ，contact。and关键词指的是颗粒与颗粒之间的接触是cmat add 定义的模型作用的范围只限于该组内部的接触，而对于另外组的接触则不定义。如果说是选的关键词是or定义的范围只要颗粒属于这个组，则定义的模型范围是所有与这个组相接触。contact关键词是用来根据颗粒的分组来判断接触的。通过对接触分组来对某个范围内的颗粒定义接触。

;不论cmat add 1 model hertz ，cmat add 2 model hertz ，cmat default type ball-ball，cmat default type ball-facet 这四个的顺序怎么样，总是按照这个顺序来执行 cmat add 1 model hertz ...

property hz_shear 30e9 hz_poiss 0.3 ...;30Gpa的切向模量，0.3的泊松比 fric 0.25 ...;0.25的摩擦系数

dp_nratio 0.2 ...;0.2的法向临界阻尼比 range group glass;区域范围为玻璃球组合

cmat add 2 model hertz ...;添加赫兹模型2 property hz_shear 70e9 hz_poiss 0.3 ... fric 0.05 ... dp_nratio 0.2 ... range group steel;区域范围为钢球组合

cmat default type ball-ball ...

model hertz ... property hz_shear 50e9 hz_poiss 0.3 ... fric 0.05 ... dp_nratio 0.2

cmat default type ball-facet ... model linear ... property kn 1e6 dp_nratio 0.2

set gravity 9.81 solve

save cmat3

return

;============================================================================== ; eof: cmat3.p2dat

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库PFC3D建模流程和代码解释(2)在线全文阅读。