pfc3d中文講解例題

1、通過例子學(xué)習(xí)幾種常見命令通過例子學(xué)習(xí)幾種常見命令例二：PFC3D目錄下的GuideStartfooting.dat，程序如下（注意：本程序與原程序不同，特加了plotsetrotation(24.0352.0340.0)這句以使得模型轉(zhuǎn)過一定角度，讓用戶看的更清楚，另外還將原程序一分為二，主要是為了更好地學(xué)習(xí)每個(gè)命令的作用），以下是程序代碼。fname:footing.DAT(tutialexamplefPFC3D)newclearpr

2、ogramstatetobeginnewproblemsetromresetromnumbergenerattitleTutialExamplewallid=1face(000)(1000)(1005)(005)wallid=2face(000)(005)(025)(020)wallid=3face(1005)(1000)(1020)(1025)wallid=4face(1025)(1020)(020)(025)wallid=5face

3、(005)(1005)(1025)(025)wallid=6face(000)(020)(1020)(1000)genid=1750rad0.150.20x=010y=02z=50plotcreateFootingplotsettitletextBasiccollectionofparticlesinaboxplotsetrotation(24.0352.0340.0)plotaddballyellowplotaddwallwhitep

4、lotaddaxesbrownplotshowpausewallid=1kn=1e8ks=1e8wallid=2kn=1e8ks=1e8wallid=3kn=1e8ks=1e8wallid=4kn=1e8ks=1e8wallid=5kn=1e8ks=1e8wallid=6kn=1e8ks=1e8propdensity2000kn1e8ks1e8iniradmul1.47sethist_rep=5histballzvel310histdi

5、agnosticmufsetdtdscalecycle3000以下是程序中逐條命令的解釋：1、我們首先定義了模型的邊界，使用WALL命令，如下：wallid=5kn=1e8ks=1e8wallid=6kn=1e8ks=1e8其中kn表示法向剛度，ks表示切向剛度，單位為Nm5、定義球的密度，法向剛度和切向剛度propdensity2000kn1e8ks1e8使用命令property，可簡(jiǎn)寫為prop，定義密度density為2000kg

6、m3，法向剛度和切向剛度均為108Nm6、改變小球半徑iniradmul1.47ini應(yīng)該是initial的簡(jiǎn)寫，用于初始化小球半徑（有待考證）rad表示半徑mul是multiply的簡(jiǎn)寫，是倍數(shù)的意思。1.47是倍數(shù)因子，表示現(xiàn)在半徑是原來半徑的1.47倍。7、設(shè)定監(jiān)測(cè)變量和采樣間隔時(shí)間sethist_rep=5histballzvel310histdiagnosticmufhist是histy歷史記錄是簡(jiǎn)寫，sethist_rep=

7、5表示每隔5個(gè)運(yùn)算步記錄一次變量的值（默認(rèn)的是10個(gè)運(yùn)算步記錄一次）；本程序設(shè)定了兩個(gè)監(jiān)測(cè)對(duì)象，一個(gè)是histballzvel310，表示的是中心坐標(biāo)最接近x=3，y=1，z=0的球的z方向的速度被監(jiān)測(cè)和記錄（在這個(gè)例子中，282號(hào)球最接近這個(gè)坐標(biāo)，因此作為監(jiān)測(cè)球）；另一個(gè)是histdiagnosticmuf，是用來監(jiān)測(cè)系統(tǒng)的不平衡力，主要用于靜力學(xué)分析，當(dāng)這個(gè)值趨于很小時(shí)，表示系統(tǒng)達(dá)到穩(wěn)定狀態(tài)。8、設(shè)定求解選項(xiàng)setdtdscale因

8、為這是一個(gè)靜態(tài)分析，因此我們開啟densityscaling選項(xiàng)（簡(jiǎn)寫為dscale）來優(yōu)化求解效率，dt應(yīng)該是這個(gè)選項(xiàng)中的某個(gè)模式（有待考證）。9、設(shè)定求解步數(shù)cycle3000當(dāng)寫入這句程序后，計(jì)算就開始了，計(jì)算的步數(shù)，時(shí)間步，總時(shí)間，不平衡力平均值和最大不平衡力這些信息均顯示在屏幕上，并且每隔5個(gè)計(jì)算步更新一次（這是sethist_rep=5的作用），圖像動(dòng)畫每隔20個(gè)計(jì)算步更新一次（默認(rèn)），這個(gè)更新時(shí)間可用命令setpinter