斗輪機挖掘取料切削阻力時間歷程仿真模型的創(chuàng)建方法
【技術領域】
[0001] 本發(fā)明涉及設計分析新型仿真建模技術領域，特別地，涉及一種斗輪機挖掘取料 切削阻力時間歷程仿真模型的創(chuàng)建方法。
【背景技術】
[0002] 斗輪機的設計制造以及理論研究已有100多年歷史。這期間，德國學者研究了斗 輪機的切削阻力系數(shù)與巖石性質(zhì)之間的關系；前蘇聯(lián)學者提出動載荷的經(jīng)驗公式、挖掘物 料黏性與切割力的關系式；我國斗輪機設計通常采用切割阻力系數(shù)計算平均切削阻力。這 些相關的阻力計算模型建立的方法與公式都與時間變量無關，無法得到時間-載荷歷程關 系。
[0003] 人們試圖描述斗輪機取料切削阻力時間歷程，用以指導設計。困難在于：斗輪挖掘 散料時，料堆受力復雜多變，與鏟斗的形狀、斗齒的布置、切削位置、切削深度、物料性質(zhì)、切 割速度、斗輪的傾角以及工況等等都有關系。由于技術手段的局限，至今尚未見到與工程實 際相符的切削阻力仿真模型構建方法。

【發(fā)明內(nèi)容】

[0004] 本發(fā)明目的在于提供一種斗輪機挖掘取料切削阻力時間歷程仿真模型的創(chuàng)建方 法，以解決斗輪挖掘散料過程無法使用計算機模型進行仿真的技術問題。
[0005] 為實現(xiàn)上述目的，本發(fā)明提供了一種斗輪機挖掘取料切削阻力時間歷程仿真模型 的創(chuàng)建方法，包括以下步驟：
[0006]A、建造料堆模型：
[0007]A1、根據(jù)實際料堆的堆積角，確定料堆模型的觀堆積體，所述觀堆積體由標準幾何 體組合而成，所述標準幾何體的單位為圓球形顆粒單元；
[0008]A2、將料堆模型納入空間直角坐標系，X方向為堅直面的水平方向，Y方向為堅直 面的堅直方向，Z方向為垂直堅直面的深度方向，以料堆模型的X、Y、Z坐標方向的尺寸分別 為Lx、Ly與Lz，圓球形顆粒單元的半徑為R;
[0009] 則X方向的顆粒單元的個數(shù)為
【主權項】
1. 一種斗輪機挖掘取料切削阻力時間歷程仿真模型的創(chuàng)建方法，其特征在于，包括以 下步驟： A、 建造料堆模型： A1、根據(jù)實際料堆的堆積角，確定料堆模型的觀堆積體，所述觀堆積體由標準幾何體組 合而成，所述標準幾何體的單位為圓球形顆粒單元； A2、將料堆模型納入空間直角坐標系，X方向為堅直面的水平方向，Y方向為堅直面的 堅直方向，Z方向為垂直堅直面的深度方向，以料堆模型的X、Y、Z坐標方向的尺寸分別為 Lx、Ly與Lz，圓球形顆粒單元的半徑為R ;
則模型顆?？倲?shù)N = Offn+Nx*Ny*Nz_Nw ;第i顆粒的坐標Xi = Offx+Xc，Yi = Offy+Yc, Zi = Offz+Zc，i = 1、2、3. ? ? N ;其中 0ffx、0ffy、0ffz 為顆粒坐標偏移值，Xc、Yc、 Zc為顆粒中心的坐標值，Nw為料堆邊界之外的顆粒數(shù)，Offn為顆粒數(shù)量偏移值； B、 建造鏟斗有限單元模型： B1、取鏟斗切削過程中的一個取料時刻作為建模狀態(tài)； B2、提取有關擋料板、斗輪體、鏟斗形狀與斗齒的布置實際尺寸信息，建立鏟斗模型； C、 將料堆模型與鏟斗有限單元模型組合： C1、建立一個定位坐標原點，依據(jù)斗輪機工作時鏟斗切削物料的技術要求，確定料堆模 型的圓球顆粒與鏟斗模型運動的圓周方向、法線方向與側向方向的初始位置及邊界條件， 使得顆粒中心點與單元面接觸耦合； C2、以斗輪軸為轉動軸定義一旋轉定位矢量，確保鏟斗轉動軌跡正確； C3、計算顆粒單元間的摩擦系數(shù)yy = tan((t)，$為摩擦角，$ =堆積角*(0. 65- 0.75);物料顆粒單元與纟產(chǎn)斗單元面的滑動摩擦系數(shù)取y的1-2倍，滾動摩擦系數(shù)取iig= U 八5-15)的 1-2 倍。
2. 根據(jù)權利要求1所述的創(chuàng)建方法，其特征在于，步驟A2所述將料堆模型納入空間直 角坐標系為： a、 按照Y、X、Z方向的順序依次掃描料堆模型， 建立顆粒單元中心坐標值（Xc，Yc，Zc);而掃描得到的顆粒單元數(shù)量M = i*j*k，當A2 中的Offn = 0時，M大于或等于模型顆?？倲?shù)N ;當Offn大于0時，M小于N ; b、 利用標準幾何體空間邊界面數(shù)學方程的約束表達式去除物料堆之外的顆粒。
3. 根據(jù)權利要求1所述的創(chuàng)建方法，其特征在于，步驟C1所述顆粒中心點與單元面接 觸耦合為： 設顆粒G與鏟斗的一單元面Sm接觸，其顆粒中心點法向接觸力為&與顆粒的位置、單 元面Sm的幾何形狀和材料常數(shù)有關，耦合表達式為：
式中Km，Lm，Am分別為單元面Sm的體積彈性模量、對角線長度與表面積；1為顆粒G中心 到單元面Sm接觸點連線矢量；nm為單元面Sm接觸點的法向矢量。
4. 根據(jù)權利要求3所述的創(chuàng)建方法，其特征在于，單元面Sm包含的4個節(jié)點與顆粒G 的接觸力耦合的計算式為： f1= 0.25(1+1 c) (i+nc)fG f2= 0.25(1-1 c)(i+nc)fG f3= 0.25(1-l c)(i-nc)fG f4= 0.25(1+1 c) (l-nc)fG €。，n。為接觸點的局部坐標值。
5. 根據(jù)權利要求1所述的創(chuàng)建方法，其特征在于，步驟A1所述的標準幾何體包括四棱 錐、四棱臺、球體、半球體、圓柱體、圓錐體和方形六面體。
【專利摘要】本發(fā)明提供了一種斗輪機挖掘取料切削阻力時間歷程仿真模型的創(chuàng)建方法，包括建造料堆模型、建造鏟斗有限單元模型、將料堆模型與鏟斗有限單元模型組合的步驟，本發(fā)明首次采用球體顆粒單元與有限單元相耦合創(chuàng)建仿真模型，提供仿真削阻力時間歷程所需的模型數(shù)據(jù)。
【IPC分類】G06F17-50
【公開號】CN104615833
【申請?zhí)枴緾N201510073903
【發(fā)明人】孫小衛(wèi), 王德森, 何路茵, 丁寧, 賀浩
【申請人】湖南長重機器股份有限公司
【公開日】2015年5月13日
【申請日】2015年2月12日