自動控制液壓板料剪切生產(chǎn)線-成品接料架設(shè)計【液壓剪板機】

自動控制液壓板料剪切生產(chǎn)線-成品接料架設(shè)計【液壓剪板機】,液壓剪板機,自動控制液壓板料剪切生產(chǎn)線-成品接料架設(shè)計【液壓剪板機】,自動控制,液壓,板料,剪切,生產(chǎn)線,成品,架設(shè),板機 柔性臂架自行式起重機傾翻載荷 柔性臂架自行式起重機傾翻載荷 S. KILICASLAN, T. BALKAN AND S. K. IDER 土耳其安卡拉 中東科技大學(xué)機械工程學(xué)院 16531 摘 要 在這一項研究中，自行式起重機的特性是利用基于柔性多體動力學(xué)理論建模分析得到的，用以確定起重機不產(chǎn)生傾翻危險的起重量。只有起重機的臂架被假設(shè)成柔性的，因為臂架是唯一一個在吊載過程中撓度較大的部件。用數(shù)學(xué)方程描述起重機的剛性部件與柔性部件之間的相互作用與耦合，并且開發(fā)了用于進行動力學(xué)分析的應(yīng)用軟件。變幅油缸的推力因臂架的仰角不同而不同，臂架運動的同時，變幅油缸的推力被計算出來，同時繪制出載荷曲線，并且將結(jié)果與用一臺10噸的自行式工程起重機實驗得到的數(shù)據(jù)進行比較。 1.緒論 作為機械系統(tǒng)的起重機一般被看成包含柔性部件的閉環(huán)的機械裝置。在以臂架的角度位置作為自變量來確定起重機額定起重量的時候，就需要用到動力學(xué)的解決方法。 目前只有極少數(shù)關(guān)于自行式起重機的動力學(xué)分析控制的研究，并且這些研究中的大部分都沒有考慮部件的撓性。Sato和Sakawa建立了一個動力學(xué)模型，用于控制一個同時進行三種運動（回轉(zhuǎn)運動、起升運動、變幅運動）的柔性回轉(zhuǎn)起重機[1]。只有主臂與副臂間的連接被假設(shè)成為柔性的，這樣假設(shè)的目的是能夠使載荷得到恰當(dāng)?shù)霓D(zhuǎn)化，這樣在轉(zhuǎn)化結(jié)束時，載荷幅值的擺動能夠以最快的速度衰減。這個能夠?qū)ζ鹬貦C進行實際載荷和額定載荷對比顯示并加以限制的系統(tǒng)在巴爾干半島上得到應(yīng)用，這個基于微型計算機的控制系統(tǒng)是通過油壓和臂架仰角來確定當(dāng)前吊鉤的實際載荷的。 在這篇文章中，起重機的特性是通過柔性多體動力學(xué)來分析得到的。給出了柔性多體動力學(xué)的動力學(xué)和運動學(xué)方程，同時開發(fā)了對起重機進行動力學(xué)分析的應(yīng)用軟件。在多體動力學(xué)分析時，系統(tǒng)的剛性聯(lián)接和柔性運動通過運用絕對的耦合和形參變量用公式來表達[3, 4]。然后，部件間的連接和指定的運動用約束方程來描述。柔性體是通過有限元方法來模擬的，形參變量是通過模型轉(zhuǎn)化得到的，用以代替彈性變量。 變幅油缸的推力因臂架的角度位置不同而不同，臂架運動的同時，變幅油缸的推力被計算出來，通過開發(fā)的應(yīng)用軟件來圖示彈性效果。同時繪制出臂架不同仰角的載荷曲線，并且將結(jié)果與起重機制造廠家提供的數(shù)據(jù)進行比較。 2.起重機建模 我們用一臺930型自行式起重機的傾翻載荷控制的實驗數(shù)據(jù)用于開發(fā)的軟件的計算，并且用了這臺起重機的結(jié)構(gòu)和參數(shù)[5]。但這個分析方法只要經(jīng)過簡單的修改就可以用于類似的起重機的分析。 重物吊在吊鉤上，臂架起升，因為重物起升過高非常危險，重物的高度是通過測量起升繩的長度來控制的。在重物起升、下降和運輸時，起重機是不能回轉(zhuǎn)的，這是由于一些限制條件，比如重物非常大或非常重、空間問題等。 在臂架的每一個角度位置都對應(yīng)著一個最大載荷，超過這個載荷后起重機很有可能傾翻。因為臂架的仰角只有在重物起升或下降時才會發(fā)生改變，而且恰巧這又是個平面運動，所以建模和分析計算都是在二維空間進行的。 圖1.試驗用起重機示意圖 圖2.試驗用起重機動力學(xué)模型 單位：mm；A0G 2000;GC 17500;A0A 5823;A0D 5850;AD 565;BD 3455;OB0 2350;OA0 805 圖1表示實驗用起重機的示意圖，圖2表示它的動力學(xué)模型，模型共分成五部分。 橫截面、材料特性及各部分的尺寸是通過技術(shù)數(shù)據(jù)文件和直接對實驗起重機進行測量得到的。部件1（臂架）的橫截面是個空心的多邊形，壁厚為t，如圖3所示。這個橫截面的尺寸線性地從A0增大到G，又線性地從G減小到C，截面A0、G和C的尺寸如圖3所示。部件2是個圓柱形的桿，直徑25mm。部件3是個液壓活塞桿，直徑180mm壁厚20mm。部件4是變幅油缸筒，內(nèi)徑230mm，外徑246mm， 長3440mm。部件1的彈性模量和密度分別為200GPa和5750kg/m3。其它部件的密度均認為是7850kg/m3。 圖3.部件1（臂架）橫截面.（單位：mm） 當(dāng)考慮了各部件的尺寸（長度和橫截面）和彈性模量，就可以只認為部件1（臂架）是彈性的，這樣的話其他部件均被假設(shè)成是剛性的。 在對起重機進行分析時要考慮以下假設(shè)： 1.液壓油的質(zhì)量包含在液壓缸（部件4）的質(zhì)量中，液壓缸質(zhì)量的改變是由于考慮了缸內(nèi)的液壓油質(zhì)量的改變； 2.液壓油被認為是不可壓縮的； 3.吊重被看成是集中質(zhì)量，并且通過一根被看成是剛性桿的繩子連接到臂架頭部。這根繩子在平面內(nèi)可以繞C點自由旋轉(zhuǎn)。只要這個桿相對于垂直位置的擺動很小并且這個桿仍然處于張緊狀態(tài)時這個假設(shè)就是正確的。在正常操作速度和吊重下這些條件都符合； 4.臂架結(jié)構(gòu)的阻尼比用Rayleigh衰減法來計算確定； 5.吊重與地面間的距離假設(shè)保持不變，這是通過在起重機臂架升降過程中改變繩長來實現(xiàn)的。 3.動力學(xué)方程 令表示一個部件的結(jié)構(gòu)，部件k的變形通過這個系數(shù)來定義，表示一個確定的結(jié)構(gòu)。令表示結(jié)構(gòu)的原點的位置，表示部件k的角速度。 利用有限單元法，部件k的i單元上的任一點P的變形位移向量為： （1） 其中是變形的單元形函數(shù)矩陣，是單元間聯(lián)系的坐標(biāo)變換矩陣，單元節(jié)點位移向量。 點P的速度可表示為： （2） 其中是中點Q到P的變形后的位置向量，是的變形協(xié)調(diào)矩陣，Tk是從到的坐標(biāo)變換矩陣，，是用于減小彈性變形的模型轉(zhuǎn)換變量，是模型形變向量。方程（2）可寫成： （3） 其中是影響系數(shù)矩陣，是部件k的速度向量。 連接各個部件的系統(tǒng)N的連接處和角度指示在速度水平上用運動學(xué)約束方程表示為： Cy=g （4） 其中C是雅客比約束矩陣，y是系統(tǒng)的速度向量，由下式確定： yT=[y(1)T…y(N)T] （5） Kane方程用于確定系統(tǒng)的運動方程： My+CTλ=Q+Fs+Fd+F （6） 其中λ是約束反力向量，M是質(zhì)量矩陣，Q、Fs、Fd及F分別是Coriolis力向量、彈性力向量、阻尼力向量和實際力向量，分別為： ，,, , （7） 質(zhì)量矩陣Mk和部件k的Coriolis向量Qk為 （8） （9） 其中Ek是部件k中的有限單元的數(shù)量，Vki是單元的體積，是它的密度。 Fsk和Fdk可按下式給出： ， （10） 其中Kk是部件k的結(jié)構(gòu)剛度矩陣，Dk是阻尼矩陣。在仿真時，結(jié)構(gòu)的質(zhì)量和用于組成Dk的結(jié)構(gòu)的剛度復(fù)數(shù)有2%的衰減。 在平面系統(tǒng)中，減小到，其中是個標(biāo)量。變成，其中是根據(jù)變換而來的。 當(dāng)公式（8）和（9）中由空間決定的變量分離后，就能獲得[3,4]不隨時間改變的矩陣。 （11） ，，，，， （12） （13） （14） （15） ， （16，17） （18） （19） 臂架是靠駕駛員控制的液壓缸來驅(qū)動的，一般來說，運動的整個過程中液壓缸以恒定的速度運動，所以臂架和活塞的振動都能控制在一個很小的水平上。為了避免沖擊載荷的產(chǎn)生，活塞的開始運動時速度從0增大到以及最終停止時從減小到0，速度的變化假設(shè)成隨時間呈擺線形變化。這個理想的速度曲線如圖4所示，并可用下列方程描述： （20） 如果部件1和2的絞點在不同的位置，這個系統(tǒng)將變成一個不能動的結(jié)構(gòu)。系統(tǒng)之所以能夠運動是因為二者絞點位置相同。所以，部件1和2的約束方程是線性相關(guān)的。由于這個原因，其中的一個約束方程可以分解以減少線性損耗。 圖4.擺線形加速度的速度曲線 臂架仰角（度） 活塞反力（kN） 圖5.不同臂架仰角的活塞反力（吊重32.4kN，起升時間30s） 臂架仰角（度） 活塞反力（kN） 圖6.不同臂架仰角的活塞反力（吊重32.4kN，起升時間10s） 臂架仰角（度） 橫向位移（m） 圖7.不同臂架仰角的節(jié)點3、8、13的橫向位移（載荷32.4kN，起升時間30s） 橫向位移（m） 臂架仰角（度） 圖8.不同臂架仰角的節(jié)點3、8、13的橫向位移（載荷32.4kN，起升時間30s） 圖9.（a）節(jié)點13的橫向位移的時間響應(yīng) （b）節(jié)點13的橫向位移的快速傅氏變換算法 4.起重機特性的計算機仿真及與實驗數(shù)據(jù)的對比 用于分析實驗用起重機的應(yīng)用軟件已經(jīng)開發(fā)完畢。在這個軟件中，部件1（臂架）的其中任一個有限單元的形函數(shù)可以代換其他任何一個。 Balkan已經(jīng)對臂架的工作范圍內(nèi)，臂架運動的30s過程做了實驗[2],選取這個速度是為了減小彈性變形的作用。在研究中測量了液壓系統(tǒng)的壓力和臂架的仰角位置。由于臂架的振動而引起的油壓的振動已經(jīng)通過控制系統(tǒng)過濾掉，所以在測量的數(shù)據(jù)中是看不到的。實驗起重機在起升過程中的吊重為32.4kN，液壓系統(tǒng)油壓的變化是在臂架運動的30s過程中測定的。所以在臂架起升的30s過程中液壓活塞的反力變化與臂架仰角的變化有關(guān)，所以可以根據(jù)吊重32.4kN來計算臂架處于不同仰角位置時的變幅油缸活塞反力。 臂架在起升的30s過程中不同仰角位置的變幅油缸活塞反力通過利用計算機代碼模擬出來，并且在圖5中給出。 臂架運動的30s過程的實驗結(jié)果也同時在圖5中給出。這些數(shù)據(jù)不包括活塞加速及減速過程。并且，由于臂架振動帶來的影響已經(jīng)被濾除，所以在圖中是不能看到的。從圖中可以看出，臂架運動過程中仿真的結(jié)果與實驗數(shù)據(jù)非常接近。 吊重32.4kN，起升時間10s時臂架不同位置的變幅活塞反力也通過計算機程序計算出來了，為的是模擬更有意義的彈性效果，如圖6所示。 在仿真時，臂架被離散成12個單元。其中兩個在A0G之間，臂架橫截面積從A0到G線性地增加。另外的十個在GC之間，臂架的橫截面積從G到C線性地減小。部件1的阻尼是通過在最初兩種模式基礎(chǔ)上依次減小2%來近似計算的。在模擬起升時間為30s時，假設(shè)臂架起升最初1.5s為加速過程，最后1.5s為減速過程。在臂架起升時間為10s的情況下，加速和減速的時間分別假設(shè)為1s。 吊重為32.4kN，臂架起升時間為30s和10s的兩個工況時，節(jié)點3（節(jié)點3在臂架節(jié)點A0和A之間），節(jié)點8（位于臂架節(jié)點A和C之間），及節(jié)點13（對應(yīng)于臂架尖端節(jié)點C）的橫向位移都根據(jù)臂架的不同仰角位置計算出來了，分別表達在圖7和圖8中。因為節(jié)點3的橫向位移的數(shù)量級為10-5m，所以這個節(jié)點的橫向位移在圖中是看不到的。 從圖5-8可以看出活塞反力的幅度和平均值及節(jié)點的橫向位移在起升時間為10s工況時的值要大于起升時間為30s的工況。因此，臂架的彈性效果可以清晰地看出來。 在所有的仿真過程中，當(dāng)臂架在起升過程中，變幅活塞的反力像期望的一樣隨之減小。在變幅活塞加速運動過程時，活塞反力、橫向位移量及它們的振動幅度都比勻速運動期間的要大。在臂架勻速起升的過程中，活塞反力的大小、橫向位移及各自的振動幅度都平穩(wěn)地減小。在臂架減速起升的過程中，活塞反力和橫向位移減小，但振動的幅度卻增加了。但在減速的過程中，活塞反力數(shù)值及幅值、橫向位移的變化卻比加速過程的小。 從仿真中還可以看出存在兩種類型的振動。一種是由臂架的振動產(chǎn)生的，周期較小，另一種由于載荷的振動產(chǎn)生的，周期較大。吊重32.4kN，起升時間為10s工況時，臂架尖端（節(jié)點13）的橫向振動的時間響應(yīng)如圖9（a）所示。圖9（a）中數(shù)據(jù)快速傅氏變換算法數(shù)值曲線在圖9（b）中給出。 頻率小的振動是由于系統(tǒng)的激勵產(chǎn)生的，它的大小變化符合這個規(guī)律。當(dāng)臂架向上運動時，它朝垂直的位置變化，從而引起臂架的橫向偏轉(zhuǎn)量減小。由于載荷振動而產(chǎn)生的振動的頻率也降到這個頻率范圍。頻率大于1.5Hz的振動是由于臂架以它的自然頻率振動。自然頻率隨時間變化是多體系統(tǒng)的一個特性，并且這引起了我們在圖9（a）中看到的尖銳的信號。 5.起升能力的仿真 起重機傾翻的模擬是在閉環(huán)的條件下來完成的，其目的是為了能夠看到在起重機臂架向上及向下運動的過程中究竟何時發(fā)生傾翻。當(dāng)起重機的一個支腿受到地面的支反力變成零時，傾翻便發(fā)生了。利用起重機底盤的自由體受力圖，如圖10所示，傾翻時可用方程（21）表示： （21） 式中，、和、是起重機底盤對臂架及變幅油缸支反作用力的分力，、是地面對起重機兩個支腿的反作用力，是起重機底盤的重力。 圖10. 起重機底盤的自由體受力圖.單位：m.A5.50；C3.37；A12.57；A24.60；B11.77；B22.25 圖11.起升能力曲線. 30s仿真；10s仿真；廠家的數(shù)據(jù) 當(dāng)小于或等于零時，就滿足了傾翻的條件。對于臂架起升時間分別為10s和30s的兩種工況，對于不同吊重，對應(yīng)的力變成零的臂架位置是通過開發(fā)的應(yīng)用軟件計算出來的。實驗用起重機的仿真結(jié)果在圖11中給出了。試驗用起重機生產(chǎn)廠家提供的額定起重量也同時列于圖11中。其中的幅度R是這樣定義的，起重機在傾翻方向上由回轉(zhuǎn)中心到臂架端點，即載荷所在位置的水平距離，通過式計算得出。 從圖11中可以看出，當(dāng)臂架運動時間減少時，相同幅度時對應(yīng)的許用起重量隨之減少了。雖然在算得許用起重量的時候沒有包含任何類似臂架運動時間的信息，但與臂架運動時間為30s時的許用起重量的曲線非常相似。除此之外，制造廠家還注釋說這些許用載荷數(shù)據(jù)應(yīng)該以一個安全系數(shù)來用，安全系數(shù)取為1.5。 6.結(jié)論 在研究中，起重機的起重能力是通過柔性多體動力學(xué)分析得到的。為了達到這個分析目的，我們開發(fā)了一個能夠?qū)ζ鹬貦C進行動力學(xué)分析的應(yīng)用軟件。系統(tǒng)的剛性連接和彈性運動通過運用絕對的耦合和形參變量用公式來表達[3, 4]。然后，部件間的鉸接和特定的運動是通過約束方程來表達的。柔性體是通過有限元方法來模擬的，形參變量是通過模型轉(zhuǎn)化得到的，用以代替彈性變量。并且利用計算機程序?qū)Φ踔貫?2.4kN，臂架起升運動時間分別為30s和10s兩種工況的變幅油缸活塞反力進行了仿真，運動速度的加速和減速時的加速度都是圓滑的擺線形。并用臂架運動時間為30s的仿真結(jié)果與實驗的結(jié)果進行了對比。臂架運動時間分別為30s和10s兩種工況下，節(jié)點3、節(jié)點8和節(jié)點13在臂架處于不同仰角位置的橫向位移都計算出來了。最后，形成了臂架起升時間分別為10s和30s的兩種工況的載荷曲線，并將其與起重機生產(chǎn)廠家提供的數(shù)據(jù)進行對比。 從分析中可以看出臂架運動時間將很顯著地影響起重機的動力學(xué)特性。在活塞運動速度較低時（比如臂架起升時間為30s），臂架的柔性起的作用很小，所以臂架在這時可以看成是剛性體。然而，當(dāng)活塞運動速度提高時（比如臂架起升時間為10s），臂架的柔性起顯著的作用。在臂架運動時間為10s時，從活塞反力和起升能力的仿真結(jié)果中可以看出臂架的柔性起的作用。 另外需要注意的是，這里計算出來的載荷曲線是在活塞運動速度曲線為擺線形的情況下得出的，這是為了近似地模擬有經(jīng)驗的駕駛員操作時的速度特性。其他形狀的速度曲線對應(yīng)的載荷曲線也同樣可以通過開發(fā)的計算機程序計算生成。變幅活塞突然加速和減速時將會顯著的降低起升能力，這是由于這將會產(chǎn)生較大的橫向位移和慣性沖擊載荷。 參考文獻：略 - 13 -