劉建軍

摘要：根據(jù)仿生學(xué)原理，對(duì)長(zhǎng)根莖類中藥材收獲機(jī)械的關(guān)鍵部件挖掘犁刀進(jìn)行了創(chuàng)新設(shè)計(jì)，同時(shí)建立了深耕挖掘犁刀的三維實(shí)體模型，運(yùn)用有限元分析軟件ＡＮＳＹＳ對(duì)其進(jìn)行了仿真分析，確定了深耕挖掘犁刀在挖掘時(shí)的應(yīng)變分布圖。結(jié)果表明，深耕挖掘犁刀產(chǎn)生的最大整體變形量為９．６２0 4 ｍｍ，發(fā)生在深耕挖掘犁刀的刀頭；最大主彈性應(yīng)變強(qiáng)度為２．５０５ ７ ｅ－２，最大彈性應(yīng)變強(qiáng)度為３．６４６ ３ ｅ－２，最大等效彈性應(yīng)變強(qiáng)度為２．４４１ ８ ｅ－２，發(fā)生在深耕挖掘犁刀螺母1上。

關(guān)鍵詞：仿生；挖掘犁刀；建模；仿真

中圖分類號(hào)：Ｓ２２５．７文獻(xiàn)標(biāo)識(shí)碼：Ａ文章編號(hào)：0439－８114（２０12）14－3084-04

Ｄｙｎａｍｉｃ Ｍｏｄｅｌｉｎｇ ａｎｄ Ｓｉｍｕｌａｔｉｎｇ Ａｎａｌｙｓｉｓ ｏｆ Ｂｉｏｎｉｃｓ Ｐｌｏｕｇｈｉｎｇ Ｃｏｌｔｅｒ

ＬＩＵ Ｊｉａｎ－ｊｕｎ

（Ｄｅｐａｒｔｍｅｎｔ ｏｆ Ｐｈｙｓｉｃｓ， Ｍｅｃｈａｎｉｃａｌ ａｎｄ Ｅｌｅｃｔｒｉｃ Ｅｎｇｉｎｅｅｉｎｇ， Ｓａｎｍｉｎｇ Ｕｎｉｖｅｒｓｉｔｙ， Ｓａｎｍｉｎｇ ３６５００４， Ｆｕｊｉａｎ，Ｃｈｉｎａ）

Ａｂｓｔｒａｃｔ： Ａｃｃｏｒｄｉｎｇ ｔｏ ｔｈｅ ｐｒｉｎｃｉｐｌｅ ｏｆ ｂｉｏｎｉｃｓ， ａ ｎｅｗ ｔｙｐｅ ｏｆ ｄｅｅｐ ｐｌｏｕｇｈｉｎｇ ｃｏｌｔｅｒ， ｔｈｅ ｋｅｙ ｏｆ ｈａｒｖｅｓｔｉｎｇ ｍａｃｈｉｎｅｒｙ ｆｏｒ ｍｅｄｉｃｉｎａｌ ｈｅｒｂｓ ｗｉｔｈ ｌｏｎｇ ｒｏｏｔ ｏｒ ｓｔａｌｋ was designed． Ｔｈｅ ３Ｄ ｓｏｌｉｄ ｍｏｄｅｌ ｏｆ ｄｅｅｐ ｅｘｃａｖａｔｉｏｎ ｐｌｏｕｇｈｉｎｇ ｃｏｌｔｅｒ ｗａｓ ｅｓｔａｂｌｉｓｈｅｄ； ａｎｄ ｔｈｅ ｖｉｓｕａｌ ｓｉｍｕｌａｔｉｏｎ ｗａｓ ｃｏｎｄｕｃｔｅｄ ｂｙ ｕｓｉｎｇ ｆｉｎｉｔｅ ｅｌｅｍｅｎｔ ａｎａｌｙｓｉｓ ｓｏｆｔｗａｒｅ ＡＮＳＹＳ ｔｏ ｃｏｎｓｔｒｕｃｔ ｔｈｅ ｓｔｒａｉｎ－ｒｅｌｅａｓｅ ｍａｐ ｏｆ ｄｅｅｐ ｅｘｃａｖａｔｉｏｎ ｐｌｏｕｇｈｉｎｇ ｃｏｌｔｅｒ ｗｈｅｎ ｄｉｇｇｉｎｇ． Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ ｔｈｅ ｇｒｅａｔｅｓｔ ｔｏｔａｌ ｄｅｆｏｒｍａｔｉｏｎ ｆｏｒ ｔｈｅ ｄｅｅｐ ｐｌｏｕｇｈｉｎｇ ｃｏｕｌｔｅｒ ｗａｓ ９．６２0 4 ｍｍ， ｏｃｃｕｒｒｉｎｇ ｉｎ ｔｈｅ ｃｕｔｔｅｒ ｈｅａｄ； ｔｈｅ ｍａｘｉｍｕｍ ｐｒｉｎｃｉｐａｌ ｅｌａｓｔｉｃ ｓｔｒａｉｎ ｗａｓ ２．５０５ ７ｅ－２； ｍａｘｉｍｕｍ ｅｌａｓｔｉｃ ｓｔｒａｉｎ ｓｔｒｅｎｇｔｈ ｗａｓ ３．６４６ ３ｅ－２； ｔｈｅ ｍａｘｉｍｕｍ ｅｑｕｉｖａｌｅｎｔ ｅｌａｓｔｉｃ ｓｔｒａｉｎ ｗａｓ ２．４４１ ８ｅ－２， ｏｃｃｕｒｒｉｎｇ ａｔ ｔｈｅ ｆｒｏｎｔ ｎｕｔ1．

Ｋｅｙ ｗｏｒｄｓ： ｂｉｏｎｉｃｓ； ｐｌｏｕｇｈｉｎｇ ｃｏｌｔｅｒ； ｅｓｔａｂｌｉｓｈｍｅｎｔ ｏｆ ｓｏｌｉｄ ｍｏｄｅｌ； ｓｉｍｕｌａｔｉｏｎ ａｎａｌｙｓｉｓ

中藥材是中國(guó)幾千年來治療、保健的主要資源，目前中國(guó)中藥資源物種數(shù)高達(dá)１３ １１２種。２０１１年中國(guó)藥材種植面積已達(dá)２０７．６萬ｈｍ２，中成藥年產(chǎn)值為８２０多億元，為促進(jìn)農(nóng)村經(jīng)濟(jì)發(fā)展、增加農(nóng)民收入發(fā)揮了重要作用。但中藥材收獲卻具有嚴(yán)格的時(shí)節(jié)性，如果不能及時(shí)、按時(shí)收獲，將大大影響藥材的品質(zhì)。目前收獲的方式主要是人工挖收或用簡(jiǎn)易挖掘機(jī)械挖掘，然后進(jìn)行人工揀拾，此項(xiàng)作業(yè)存在勞動(dòng)強(qiáng)度大、勞動(dòng)效率低、收獲損失高、藥材的品質(zhì)差等缺陷，嚴(yán)重影響經(jīng)濟(jì)效益，極大地限制了中藥材種植業(yè)的發(fā)展。

以長(zhǎng)根莖類中藥材收獲機(jī)的關(guān)鍵部件挖掘犁刀為研究對(duì)象，以仿生學(xué)的原理設(shè)計(jì)一種新型的減阻降耗犁刀，并對(duì)此進(jìn)行相關(guān)的仿真分析。

１仿生深耕挖掘犁刀的設(shè)計(jì)

１．１仿生減阻機(jī)理研究

從生命科學(xué)的角度來看，生物體經(jīng)過漫長(zhǎng)的歷史進(jìn)化進(jìn)程，大大優(yōu)化了自身的結(jié)構(gòu)，使其適應(yīng)環(huán)境的能力大大提高，這樣就得到了生存和繁衍。許多困擾人類的科技難題都把焦點(diǎn)聚集在了生物和自然身上，仿生學(xué)就是以生物與自然中的基本特征或形態(tài)為藍(lán)本，借助于高科技的技術(shù)手段，對(duì)其進(jìn)行“形”或“神”的模仿，實(shí)現(xiàn)特定的技術(shù)功能或目標(biāo)［１，２］。這種仿生減阻技術(shù)為耕作部件的減阻技術(shù)研究提供了一種新的思路。

仿生減阻的研究大體分為兩個(gè)步驟。第一步是對(duì)現(xiàn)有生物系統(tǒng)的結(jié)構(gòu)、功能、過程或行為機(jī)理進(jìn)行內(nèi)涵研究，獲取面向仿生應(yīng)用的生物信息，這是仿生減阻研發(fā)的基礎(chǔ)和前提［３］。第二步是基于生物結(jié)構(gòu)、功能、過程或行為機(jī)理解決科學(xué)技術(shù)難題，研制具有類似于生物系統(tǒng)減阻功能的技術(shù)或裝備。

郭志軍等［４］使用ＸＴＬ３０－ＣＴＶ體視顯微鏡分析了公雞、家鼠、螻蛄以及田鼠等幾種動(dòng)物爪趾的幾何特征。對(duì)于每一個(gè)爪趾而言，其縱剖面上都具有２條輪廓線，位于觸土面一側(cè)的輪廓線稱為內(nèi)輪廓線，而把位于背部的輪廓線稱作外輪廓線。內(nèi)、外輪廓線的相對(duì)位置數(shù)據(jù)可以通過最小二乘法來分析與擬合，把每條輪廓線的四次擬合公式換入標(biāo)準(zhǔn)的曲率方程，可計(jì)算出每條曲線在不同位置坐標(biāo)下的曲率。擬合結(jié)果表明，各爪趾的輪廓線特別是內(nèi)輪廓線表現(xiàn)出明顯的變曲率特征，以田鼠爪趾為例，圖１給出了其中的一只田鼠左前爪中趾的體視放大照片，測(cè)量過程中保持爪趾輪廓線所在平面與鏡頭垂直，則四次多項(xiàng)式擬合方程為＝０．０００ ２ｘ４－０．００１ ４ｘ３＋０．００７ ５ｘ２－０．３４４ ２ｘ－１．６８９ ２，Ｒ２＝０．９９０ ６，圖２給出了與圖１所示爪趾內(nèi)輪廓線相應(yīng)的曲率變化趨勢(shì)線，橫坐標(biāo)為測(cè)量過程中爪趾的相對(duì)位置，縱坐標(biāo)為曲率，原點(diǎn)則表示爪趾與肢體相連接的位置。由圖２可知，爪趾內(nèi)輪廓線的曲率趨勢(shì)線具有明顯的波動(dòng)變化的特征，其變化規(guī)律可看作兩條拋物線的近似耦合。

深耕挖掘犁刀是長(zhǎng)根莖類中藥材收獲機(jī)最重要的工作部件，主要由刀頭和刀柄兩部分組成，工作過程為：刀頭切開并松動(dòng)根莖底層土壤，刀柄刃切開并松動(dòng)行間土壤，實(shí)現(xiàn)降低后續(xù)收獲阻力和減少挖掘刀面磨損的功能。挖掘犁刀在工作時(shí)利用拖拉機(jī)的牽引力對(duì)土壤進(jìn)行切削，挖掘犁刀與土壤之間的相互作用不僅與各自的材料特性有關(guān)，而且還與挖掘犁刀的結(jié)構(gòu)參數(shù)和運(yùn)動(dòng)參數(shù)有關(guān)。

１．２挖掘刀頭的設(shè)計(jì)

１．２．１刀頭的受力分析刀面的受力模型見圖３，工作時(shí)，機(jī)具以速度ｖｍ前進(jìn)，與水平面成α角的斜面ＡＢ相當(dāng)于一個(gè)單面楔。土壤對(duì)其作用力Ｎ垂直于ＡＢ面，近似作用于ＡＢ段的中點(diǎn)，行進(jìn)中除阻力Ｎ外，尚有一個(gè)沿ＡＢ面指向后方的摩擦力Ｆ，Ｎ與Ｆ合成為阻力Ｒ。阻力Ｒ可分解成一個(gè)垂直分力Ｒｚ和一個(gè)水平分力Ｒｘ。垂直分力Ｒｚ和水平分力Ｒｘ為工作時(shí)必須要克服的牽引阻力［５］。

Ｒｘ＝Ｒｚ·ｃｏｔβ

而作用力Ｒｚ的大小隨土壤的土質(zhì)、含水量及前茬作物的種類不同而不同，在保證Ｒｚ一定的情況下，β角越大，牽引阻力Ｒｘ越小，且有關(guān)系式：

α＋β＋φ＝９０°

式中，α為起土角，β為合力作用角，φ為摩擦角。

１．２．２刀頭的形式及參數(shù)確定深耕挖掘犁刀通常具有３種形式：鑿形、箭形和雙翼形，鑿形有利于作物行間深松土壤，不翻轉(zhuǎn)表層土，故設(shè)計(jì)選用鑿形（圖４）。深耕挖掘犁刀工作參數(shù)包括結(jié)構(gòu)參數(shù)和作業(yè)參數(shù)，結(jié)構(gòu)參數(shù)主要有刀面長(zhǎng)度Ｌ、刀面寬度Ｂ、刀板厚度δ和刀板前角θ等，大小分別為１６５、４０、１０ ｍｍ和２０°；作業(yè)參數(shù)主要有挖掘犁刀工作速度ｖ０、起土角α、挖掘深度等。

深耕挖掘犁刀刀頭的起土角α一般為２０°左右，此角過大將增大阻力。試驗(yàn)結(jié)果表明，起土角α從４５°左右減小到２０°時(shí)，牽引阻力隨著α減小而降低；當(dāng)α＞５０°時(shí)，刀頭將粘結(jié)上由壓實(shí)土壤構(gòu)成的土楔，設(shè)計(jì)中實(shí)際采用的α＝２３°。挖掘深度按照中藥材的成長(zhǎng)深度選定為４００ ｍｍ。挖掘犁刀采用ＧＢ／Ｔ ７１１中規(guī)定的６５Ｍｎ鋼制造，刃部應(yīng)進(jìn)行熱處理。淬火區(qū)寬度為２０～３０ ｍｍ，硬度應(yīng)為ＨＲＣ４８～５６。

１．２．３刀頭最大橫向影響范圍土壤受到窄型工作部件擠壓時(shí)，工件兩側(cè)的土壤破裂界面在水平投影面內(nèi)呈扇形，深耕挖掘犁刀的最大橫向影響范圍計(jì)算公式為：

b≤d

式中，ｄ為深耕挖掘犁刀的最大橫向?qū)挾?；ａ為深耕挖掘犁刀的作業(yè)深度；θ為土壤的前切滑移角，一般為５０°左右；φ為土壤的剪切滑移破裂面與法線方向的夾角，一般為土壤與金屬的摩擦角，大約為２５°；α為深耕挖掘犁刀刀頭的起土角。將上述數(shù)值帶入公式得到挖掘刀頭的最大橫向影響范圍為ｂ≤１１１．１ ｃｍ。由理論計(jì)算可以得出深耕挖掘犁刀橫向影響范圍比較寬，有利于行間土壤的松動(dòng)，降低根莖收獲時(shí)的阻力。

１．３挖掘刀柄的設(shè)計(jì)

１．３．１刀柄的受力分析刀柄的截面形狀可看作雙面楔，其受力模型如圖５所示。

式中，Ｐ１為楔面上的正壓力產(chǎn)生的阻力分力；Ｐ２為楔面上的剪力產(chǎn)生的阻力分力；Ｐ３為側(cè)面上的剪力產(chǎn)生的阻力分力；μ為摩擦系數(shù)。

由于正壓力是土壤抵抗變形產(chǎn)生的力，其大小難確定，令Ｎ＝Ｋ１Ｆ１，Ｎ１＝Ｋ２Ｆ２。式中，Ｋ１為土壤抵抗變形的單位阻力；Ｋ２為土壤對(duì)于楔子側(cè)面的單位壓力；Ｆ１為楔面面積；Ｆ２為楔子側(cè)面面積；Ｎ為楔面上的正壓力；Ｎ１為側(cè)面上的剪力。Ｋ２與Ｋ１有所不同，其不同之處在于犁柄不運(yùn)動(dòng)時(shí)Ｋ１為零而Ｋ２不變。Ｋ２反映出土壤彈性恢復(fù)特性。所以：

１．３．２刀柄的參數(shù)確定依據(jù)中藥材根莖的分布寬度、生長(zhǎng)深度、土壤土質(zhì)、根系的抓土程度來綜合確定挖掘刀柄的參數(shù)。利用三面楔原理碎土，當(dāng)土壤受到刀面的力大于土壤的內(nèi)聚力時(shí)會(huì)產(chǎn)生剪切破壞［６］。如圖６所示，刀柄的入土部分采用弧形最為省力，還可以減少雜草纏結(jié)。挖掘犁刀觸土面導(dǎo)線形式為復(fù)合型，中部為圓弧，圓弧半徑設(shè)計(jì)為３２０ ｍｍ，弧形角為２３°。為更好地減少阻力，刀柄與刀頭連接處采用與水平面成２３°角的切線過度。為配合刀頭入土４００ ｍｍ，刀柄上設(shè)計(jì)有６個(gè)Φ２５ ｍｍ等距孔，這和機(jī)架連接使用，同時(shí)還可以作為入土深度調(diào)整孔。考慮到刀尖處受到不平衡外力，應(yīng)設(shè)計(jì)足夠的抗彎、抗扭強(qiáng)度。因此，刀柄垂直部分上部截面為２５ ｍｍ×６０ ｍｍ，耕層內(nèi)圓弧呈棱角形，棱角為６０°。刀柄采用不低于ＧＢ／Ｔ ７００規(guī)定的Ｑ２７５鋼制造，且應(yīng)保證其剩余變形量不超過３ ｍｍ。深耕挖掘犁刀和挖掘刀柄的其他技術(shù)條件應(yīng)符合標(biāo)準(zhǔn)ＪＢ／Ｔ ９７８８—１９９９中的相應(yīng)規(guī)定。

２仿生深耕挖掘犁刀的仿真分析

２．１仿生深耕挖掘犁刀的三維實(shí)體模型

基于Ｐｒｏ／Ｅ建立深耕挖掘犁刀的三維模型，如圖７所示，仿生深耕挖掘犁刀主要由刀柄、刀頭、沉頭螺栓和螺母構(gòu)成。犁刀材料為優(yōu)質(zhì)碳素結(jié)構(gòu)鋼６５Ｍｎ，其彈性模量為２．１×１０８ ｋＰａ，剪切模量為８．０１６×１０７ ｋＰａ，泊松比為０．３，密度為７ ８５０ ｋｇ／ｍ３。為了研究深耕挖掘犁刀的強(qiáng)度特性，計(jì)算工作載荷下的應(yīng)變強(qiáng)度和變形量，把模型進(jìn)行簡(jiǎn)化和抽象，忽略一些圓角。

２．２仿生深耕挖掘犁刀的有限元分析

采用理想彈塑性材料本構(gòu)關(guān)系來模擬切削工具—土壤接觸系統(tǒng)，不考慮切削速度對(duì)切削阻力的影響。土壤彈性模量為１ ２５０ ｋＰａ，剪切模量為５４３．４７８ ｋＰａ，泊松比為０．１５［７］。土壤模型的網(wǎng)格劃分采用自動(dòng)劃分方式，單元大小采用默認(rèn)值。整個(gè)土壤模型單元有２２ ９１９個(gè)，節(jié)點(diǎn)有２５ ８４４個(gè)［８］。挖掘犁刀模型網(wǎng)格也采用自動(dòng)劃分的方式，單元有２４９ ７７７個(gè)，節(jié)點(diǎn)有６１ ７５２個(gè)（圖８）。在保證計(jì)算精度的前提下，為了降低計(jì)算量，將靠近接觸區(qū)域部分的有限元網(wǎng)格劃分得比較密集，而將遠(yuǎn)離接觸區(qū)域的網(wǎng)格劃分得比較稀疏。

根據(jù)深耕挖掘犁刀作業(yè)的實(shí)際測(cè)試和土槽試驗(yàn)的結(jié)果，在刀柄的上部垂直于挖掘刀面方向施加１７８ ＭＰａ的壓強(qiáng)，接觸部分用懲罰函數(shù)法進(jìn)行計(jì)算。通過ＡＮＳＹＳ軟件靜態(tài)分析得到刀頭整體的變形量如圖９所示，分析結(jié)果見表１。

分析結(jié)果表明，在試驗(yàn)條件下深耕挖掘犁刀產(chǎn)生的最大整體變形量為９．６２０ ４ ｍｍ，發(fā)生在深耕挖掘犁刀的刀頭；最大主彈性應(yīng)變強(qiáng)度為２．５０５ ７ｅ－２，發(fā)生在深耕挖掘犁刀螺母1上。最大彈性應(yīng)變強(qiáng)度為３．６４６ ３ｅ－２，最大等效彈性應(yīng)變強(qiáng)度為２．４４１ ８ｅ－２，發(fā)生在深耕挖掘犁刀螺母1上。通過以上仿真分析可知最大應(yīng)變值小于材料許用應(yīng)變值，說明該結(jié)構(gòu)在此靜載荷作用下是安全的。

３結(jié)論

通過對(duì)仿生深耕挖掘犁刀的設(shè)計(jì)與仿真分析，可以得出如下結(jié)論。

１）利用仿生學(xué)原理，模擬生態(tài)環(huán)境中逐步形成的具有特定曲率的土壤動(dòng)物爪趾，以刀頭、刀柄的受力分析為基礎(chǔ)，確定了挖掘犁刀的各項(xiàng)基本參數(shù)，設(shè)計(jì)了新型的仿生深耕挖掘犁刀，該部件挖掘過程中所受切削阻力較小，降低了磨損。

２）建立了仿生深耕挖掘犁刀的三維實(shí)體模型，并進(jìn)行了有限元的仿真分析。結(jié)果表明深耕挖掘犁刀產(chǎn)生的最大整體變形量為９．６２０ ４ ｍｍ，發(fā)生在深耕挖掘犁刀的刀頭；最大主彈性應(yīng)變強(qiáng)度為２．５０５ ７ ｅ－２，發(fā)生在深耕挖掘犁刀螺母1上。最大彈性應(yīng)變強(qiáng)度為３．６４６ ３ ｅ－２，最大等效彈性應(yīng)變強(qiáng)度為２．４４１ ８ ｅ－２，發(fā)生在深耕挖掘犁刀螺母1上。仿真結(jié)果表明得到的最大應(yīng)變值小于材料許用應(yīng)變值，說明該結(jié)構(gòu)在此靜載荷作用下是合理并且安全的。

３）新型的深耕挖掘犁刀設(shè)計(jì)以仿生學(xué)為基礎(chǔ)，可以顯著減少挖掘阻力，降低能源消耗，為中藥材收獲機(jī)的研制提供了科學(xué)的理論依據(jù)，具有較大的工程實(shí)用意義。

參考文獻(xiàn)：

［１］ 李守仁，林金天． 驅(qū)動(dòng)型土壤耕作機(jī)械的理論與計(jì)算［Ｍ］．北京：機(jī)械工業(yè)出版社，１９９７．

［２］ 朱鳳武，佟金． 土壤深松技術(shù)及高效節(jié)能仿生研究的發(fā)展［Ｊ］．吉林大學(xué)學(xué)報(bào)（工學(xué)版），２００３，３３（２）：９５－９９．

［３］ 郭志軍，周志立，任露泉． 達(dá)烏爾黃鼠爪趾幾何特征分析［Ｊ］．河南科技大學(xué)學(xué)報(bào)（自然科學(xué)版），２００３，２４（１）：１－４．

［４］ 郭志軍，周志立，佟金，等．拋物線型削面刀具切削性能二維有限元分析［Ｊ］．河南科技大學(xué)學(xué)報(bào)（自然科學(xué)版），２００２，２３（４）：１－４．

［５］ 郭志軍，周志立，徐東，等．高效節(jié)能仿生深松部件的試驗(yàn)［Ｊ］． 河南科技大學(xué)學(xué)報(bào)（自然科學(xué)版），２００３，２４（３）：１－３．

［６］ 郭志軍，鄧志強(qiáng)，陳海燕． 仿生挖掘斗－土壤接觸二維有限元分析［Ｊ］．建筑機(jī)械：上半月，２００７（3）：７４－７７．

［７］ 陸懷民． 切土部件與土壤相互作用的粘彈塑性有限元分析［Ｊ］．土木工程學(xué)報(bào)，２００２，３５（６）：７９－８１．

［８］ 郭志軍，佟金，周志立，等． 耕作部件－土壤接觸問題研究方法分析［Ｊ］． 農(nóng)業(yè)機(jī)械學(xué)報(bào)，２００１，３２（４）：１０２－１０４．