DOI厲倥rkivS||振動與沖擊物料沖擊破碎過程的一種非線性力模型秦志英趙月靜侯書軍（河北科技大學(xué)機械電子工程學(xué)院，石家莊050054）壓力。

　　A狀態(tài)―B狀態(tài)為壓縮階段，B狀態(tài)―C狀態(tài)為膨脹階段。h2〈h.，h1〈fc即在，沖擊破碎過程的接觸力可以簡化為一個分段線性的接觸力模型。

　　AB為壓縮階段（X>0），knl為壓縮剛度，為初始接觸條件，接觸力為為膨脹階段（fC0），kn2為膨脹剛度，心為破碎后脫離接觸條件，接觸力為0而B為壓縮階段和膨脹階段的分界（x=0），為最大接觸變開形Fm為最大接觸力。

　　假設(shè)A點剛體的沖擊速度為va則在壓縮階段剛體損失的動能，等于物料層吸收的能量，即jmva 1假設(shè)C點剛體的沖擊速度為Vc則在膨脹階段剛體獲得的動能，等于物料層釋放的能量，即-1mv2為了宏觀地表示剛體在整個沖擊破碎過程中所損失的能量，也即用于物料破碎的能量，參照正碰撞過程引入恢復(fù)系數(shù)e=-VcXa再根據(jù)B點的力相等條則膨脹剛度為2=4kn1，破碎后脫因此對于這個分段線性接觸力模型，關(guān)鍵是確定壓縮剛度knl和恢復(fù)系數(shù)e=01的取值，并在大量實驗的基礎(chǔ)上，建立料層性態(tài)與模型參數(shù)之間的相互對應(yīng)關(guān)系。而初始接觸條件和最大接觸變形主要表示剛體施力的幾何條件、及剛體施力的大小，破碎后脫離接觸條件心則是一個非獨立變量，由其他參數(shù)推出。

　　同理，也可以擴展為一個分段非線性接觸力模型只要壓縮和膨脹階段采用相同n次冪的多項式，且引入恢復(fù)系數(shù)e就可以很方便地得到下面的關(guān)系：膨脹剛度為kn2=~ifcl，破碎后脫離接觸條件仍保持為心e =、-e2Um-1）。這樣增加了一個模型參數(shù)，能模擬更復(fù)雜的物料層性態(tài)，分段線性模型只是n= 1的特例。而在大量，針對振動破磨過程中一個最簡的單球沖擊破碎過程，將其簡化為一個垂直方向的單自由度系統(tǒng)，建立其DEM計算模型。

　　斷，t時刻小球的動力學(xué)方程為mx（t）+f（x二0其中f（x）為小球和固定面之間的接觸力，力模型與大小取決于兩者之間待破碎物料層的性態(tài)。

　　法求解運動方程，得t +V時刻的位移、速度根據(jù)待破碎物料層的性態(tài)，確定接觸力模型和模型參數(shù)。

　　DEM計算中常用的線性彈簧阻尼模型為），其中線性阻尼為cn而阻尼系數(shù)和恢復(fù)系數(shù)之間的關(guān)系為接觸剛度5*105Nm不同恢復(fù)系數(shù)能量損失率的誤差（％）線性模型分段線性模型Q101728%-Q0192%Q302453%Q不同接觸剛度能量損失率的誤差（％）線性模型分段線性模型005*105N/m63196%-Q0060%q球質(zhì)量m=026kg小球初始高度h=50nm初始接觸條件為（1）在給定接觸剛度kn1 =5X1Nm恢夏系數(shù)e=05時，比較沖擊破碎的微觀過程。由圖可見，整個沖擊破碎過程的時間很短，大約只有幾毫秒，因此如何用。可以看出，在接觸參數(shù)相同時，分段線性模型的計算精度高于線性模型。

　　表1接觸參數(shù)對計算誤差的影響-1），去除了沖擊速度的影響。設(shè)沖擊前的速度為v沖擊后的速度為u= 4結(jié)論沖擊破碎過程是一個剛散耦合動力學(xué)系統(tǒng)，物料的性態(tài)是系統(tǒng)非線性的根本來源。

　　分段線性的接觸力模型可以反映沖擊破碎過程的本質(zhì)，模型參數(shù)具有明確的物理意義，尤其是恢復(fù)系數(shù)的引入，可以從宏觀上反映剛體在沖擊破碎過程中的能量損失。

　　通過計算機模擬，從微觀上來看接觸力在整個過程中保持正壓力，且由于物料的破碎產(chǎn)生了不可恢復(fù)的變形，符合物料破碎的實際情況。