干了這么多年機(jī)械結(jié)構(gòu)仿真曲柄滑塊機(jī)構(gòu)可以說(shuō)是我接觸最多的機(jī)構(gòu)之一。從內(nèi)燃機(jī)的活塞-連桿-曲軸系統(tǒng)到?jīng)_壓機(jī)的傳動(dòng)機(jī)構(gòu)再到自動(dòng)生產(chǎn)線上的送料裝置曲柄滑塊幾乎無(wú)處不在。它結(jié)構(gòu)簡(jiǎn)單——一根曲柄、一根連桿、一個(gè)滑塊再加上機(jī)架就能把旋轉(zhuǎn)運(yùn)動(dòng)轉(zhuǎn)成直線往復(fù)運(yùn)動(dòng)或者反過(guò)來(lái)把直線運(yùn)動(dòng)轉(zhuǎn)成旋轉(zhuǎn)運(yùn)動(dòng)。但就是這么一個(gè)“看似簡(jiǎn)單”的四桿機(jī)構(gòu)真要在設(shè)計(jì)階段把它吃透光靠手算公式和經(jīng)驗(yàn)系數(shù)遠(yuǎn)遠(yuǎn)不夠。這也是為什么現(xiàn)在做機(jī)械產(chǎn)品研發(fā)幾乎離不開(kāi)多體動(dòng)力學(xué)仿真系統(tǒng)。把曲柄滑塊機(jī)構(gòu)放進(jìn)仿真環(huán)境里建立剛體模型、定義運(yùn)動(dòng)副和約束、施加驅(qū)動(dòng)力矩和工作載荷然后求解整個(gè)運(yùn)動(dòng)周期內(nèi)的位移、速度、加速度和受力情況——這套流程聽(tīng)起來(lái)很常規(guī)但真正用好了它能回答很多在圖紙階段回答不了的問(wèn)題連桿銷(xiāo)軸處的受力峰值到底多大軸承壽命能不能滿足要求機(jī)構(gòu)的急回特性對(duì)生產(chǎn)節(jié)拍有什么影響在轉(zhuǎn)速變化和負(fù)載波動(dòng)下機(jī)構(gòu)會(huì)不會(huì)出現(xiàn)振動(dòng)或沖擊這些問(wèn)題都直接決定設(shè)計(jì)能不能落地。這篇文章我想從“仿真系統(tǒng)功能分析”的角度切入把曲柄滑塊機(jī)構(gòu)多體動(dòng)力學(xué)仿真系統(tǒng)的核心功能、使用邏輯和實(shí)操細(xì)節(jié)系統(tǒng)地拆一遍。文章不僅適合剛開(kāi)始接觸多體仿真的機(jī)械工程師也適合需要引入或選型仿真工具的技術(shù)管理者。無(wú)論你用的是哪款商業(yè)軟件還是開(kāi)源求解器下面講的功能邏輯和分析思路都是通用的。1. 先說(shuō)清楚多體動(dòng)力學(xué)仿真系統(tǒng)在這個(gè)機(jī)構(gòu)上到底分析什么很多新手拿到曲柄滑塊機(jī)構(gòu)第一反應(yīng)就是直接建模、直接跑仿真、直接看動(dòng)畫(huà)。這樣做不是不行但很容易跑完只得到一堆五顏六色的結(jié)果云圖卻說(shuō)不清自己到底在分析什么。在做功能拆解之前我們必須先把“仿真系統(tǒng)在這個(gè)機(jī)構(gòu)上到底用來(lái)解決什么問(wèn)題”這個(gè)底層邏輯理清。1.1 曲柄滑塊機(jī)構(gòu)的運(yùn)動(dòng)本質(zhì)旋轉(zhuǎn)-往復(fù)的轉(zhuǎn)換與慣性沖擊曲柄滑塊機(jī)構(gòu)的核心運(yùn)動(dòng)功能是把曲柄的整周旋轉(zhuǎn)通過(guò)連桿的平面復(fù)合運(yùn)動(dòng)轉(zhuǎn)換為滑塊在導(dǎo)軌上的往復(fù)直線運(yùn)動(dòng)。這個(gè)轉(zhuǎn)換過(guò)程看起來(lái)順暢但隱含著一個(gè)工程上必須重視的現(xiàn)象滑塊的速度和加速度在整個(gè)行程內(nèi)劇烈變化尤其在曲柄轉(zhuǎn)角接近0度和180度附近時(shí)速度方向反轉(zhuǎn)加速度達(dá)到峰值。舉個(gè)例子一個(gè)典型的內(nèi)燃機(jī)曲柄連桿機(jī)構(gòu)曲柄轉(zhuǎn)速6000轉(zhuǎn)/分曲柄半徑50mm連桿長(zhǎng)度150mm時(shí)滑塊活塞的峰值加速度可以達(dá)到數(shù)千m/s2遠(yuǎn)遠(yuǎn)超過(guò)重力加速度。這意味著什么意味著滑塊和連桿在換向點(diǎn)附近承受著巨大的慣性力。慣性力不解決機(jī)構(gòu)振動(dòng)、噪聲、螺栓松動(dòng)、軸承磨損都會(huì)提前找上門(mén)。多體動(dòng)力學(xué)仿真系統(tǒng)在這里的第一個(gè)核心功能就是在給定機(jī)構(gòu)幾何尺寸、質(zhì)量和慣量參數(shù)之后精確計(jì)算出整個(gè)運(yùn)動(dòng)周期內(nèi)滑塊的位移、速度、加速度曲線以及各個(gè)構(gòu)件的角速度、角加速度。這些運(yùn)動(dòng)學(xué)結(jié)果是一切后續(xù)強(qiáng)度分析和動(dòng)力學(xué)分析的基礎(chǔ)。1.2 從運(yùn)動(dòng)學(xué)到動(dòng)力學(xué)系統(tǒng)功能分析的兩大層次我在給團(tuán)隊(duì)做內(nèi)訓(xùn)時(shí)經(jīng)常強(qiáng)調(diào)一個(gè)觀點(diǎn)多體動(dòng)力學(xué)仿真系統(tǒng)的功能至少要分成兩個(gè)層次來(lái)理解。第一層是運(yùn)動(dòng)學(xué)分析。運(yùn)動(dòng)學(xué)分析不涉及質(zhì)量也不涉及力只關(guān)心幾何約束下機(jī)構(gòu)的運(yùn)動(dòng)狀態(tài)。分析目標(biāo)是計(jì)算各個(gè)構(gòu)件的位置、速度和加速度。運(yùn)動(dòng)學(xué)分析的輸入是曲柄轉(zhuǎn)速或驅(qū)動(dòng)函數(shù)輸出是滑塊運(yùn)動(dòng)曲線和連桿姿態(tài)變化。這個(gè)層次主要用于驗(yàn)證機(jī)構(gòu)能否實(shí)現(xiàn)預(yù)期的運(yùn)動(dòng)規(guī)律比如行程是否滿足要求、速度是否均勻、有沒(méi)有運(yùn)動(dòng)干涉。第二層是動(dòng)力學(xué)分析。動(dòng)力學(xué)分析要計(jì)入構(gòu)件的質(zhì)量、轉(zhuǎn)動(dòng)慣量、重力、驅(qū)動(dòng)力矩、工作載荷、摩擦力、接觸力等一系列因素求解系統(tǒng)的運(yùn)動(dòng)方程得到各個(gè)運(yùn)動(dòng)副約束反力、構(gòu)件內(nèi)力以及運(yùn)動(dòng)響應(yīng)。這個(gè)層次的分析結(jié)果直接為后續(xù)的有限元強(qiáng)度分析提供邊界條件也為選電機(jī)、定軸承、校核螺栓提供依據(jù)。一套完整的曲柄滑塊機(jī)構(gòu)仿真系統(tǒng)功能分析必須以這兩個(gè)層次為主線。只跑得動(dòng)運(yùn)動(dòng)學(xué)而算不了動(dòng)力學(xué)的系統(tǒng)只能算是一個(gè)幾何分析工具談不上真正的多體動(dòng)力學(xué)仿真系統(tǒng)。1.3 功能分析的服務(wù)對(duì)象設(shè)計(jì)驗(yàn)證與故障復(fù)現(xiàn)除了常規(guī)的設(shè)計(jì)階段仿真驗(yàn)證仿真系統(tǒng)還有一類(lèi)非常實(shí)用的功能場(chǎng)景——故障復(fù)現(xiàn)與問(wèn)題排查。我在實(shí)際項(xiàng)目中遇到過(guò)這樣的情況一臺(tái)沖壓設(shè)備在高速運(yùn)轉(zhuǎn)時(shí)滑塊的導(dǎo)向?qū)к壆惓０l(fā)熱同時(shí)伴有明顯的沖擊噪聲。拆開(kāi)檢查導(dǎo)軌磨損嚴(yán)重。在設(shè)計(jì)時(shí)所有人都覺(jué)得導(dǎo)軌承受的載荷不大沒(méi)有仔細(xì)校核。后來(lái)在仿真系統(tǒng)里精確建模加上了運(yùn)動(dòng)副間隙和導(dǎo)軌接觸模型結(jié)果發(fā)現(xiàn)滑塊在換向瞬間對(duì)導(dǎo)軌的側(cè)向沖擊力比理論計(jì)算值大了近三倍而且沖擊力方向隨曲柄轉(zhuǎn)角周期性變化。這個(gè)發(fā)現(xiàn)直接改變了導(dǎo)軌的設(shè)計(jì)方案。所以說(shuō)多體動(dòng)力學(xué)仿真系統(tǒng)的功能分析不能只停留在“能跑出曲線”這個(gè)層面更要理解它作為工程驗(yàn)證平臺(tái)在設(shè)計(jì)優(yōu)化、故障復(fù)現(xiàn)、壽命預(yù)測(cè)這些場(chǎng)景中的真實(shí)價(jià)值。2. 系統(tǒng)功能模塊全拆解從幾何建模到結(jié)果后處理的完整鏈路一套合格的曲柄滑塊機(jī)構(gòu)多體動(dòng)力學(xué)仿真系統(tǒng)在功能架構(gòu)上通常由五個(gè)核心模塊構(gòu)成幾何建模模塊、約束與連接定義模塊、驅(qū)動(dòng)與載荷施加模塊、求解器模塊、后處理模塊。這五個(gè)模塊形成一個(gè)完整的分析鏈路任何一個(gè)環(huán)節(jié)的功能缺失或邏輯錯(cuò)誤都會(huì)直接污染最終仿真結(jié)論。2.1 幾何建模與剛體定義仿真精度的第一道門(mén)檻幾何建模是多體動(dòng)力學(xué)仿真的起點(diǎn)。在當(dāng)前主流軟件中幾何建模通常有三種途徑第一種是直接在仿真系統(tǒng)內(nèi)置的建模界面中繪制簡(jiǎn)化幾何體。這種方法適合機(jī)構(gòu)概念設(shè)計(jì)階段因?yàn)榇藭r(shí)不需要很精確的零件外形只需要把曲柄、連桿、滑塊、機(jī)架這些關(guān)鍵構(gòu)件的拓?fù)潢P(guān)系和各鉸點(diǎn)位置確定下來(lái)就可以進(jìn)行運(yùn)動(dòng)學(xué)驗(yàn)證。比如先畫(huà)一條線代表曲柄一條線代表連桿一個(gè)方塊代表滑塊快速驗(yàn)證機(jī)構(gòu)的運(yùn)動(dòng)學(xué)邏輯。第二種是從CAD系統(tǒng)導(dǎo)入三維幾何模型。這種方式精度較高適合詳細(xì)設(shè)計(jì)階段可以直接使用裝配體中已經(jīng)定好的尺寸和配合。但要注意一個(gè)問(wèn)題導(dǎo)入的幾何模型往往包含大量對(duì)多體動(dòng)力學(xué)分析無(wú)關(guān)緊要的特征比如圓角、倒角、螺紋孔、裝飾性凸臺(tái)等。這些特征會(huì)增加接觸檢測(cè)的計(jì)算量卻不會(huì)提高求解精度。第三種是參數(shù)化建模。這是目前比較推薦的做法尤其在做機(jī)構(gòu)系列化設(shè)計(jì)時(shí)通過(guò)參數(shù)化尺寸驅(qū)動(dòng)能夠在幾秒鐘內(nèi)生成不同桿長(zhǎng)、不同偏心距的機(jī)構(gòu)模型為后續(xù)優(yōu)化掃清障礙。無(wú)論用哪種方式建模都必須保證一個(gè)核心要素正確鉸點(diǎn)位置曲柄與機(jī)架的旋轉(zhuǎn)中心、曲柄與連桿的連接點(diǎn)、連桿與滑塊的連接點(diǎn)以及滑塊與機(jī)架導(dǎo)軌的相對(duì)位置。在曲柄滑塊機(jī)構(gòu)中這些鉸點(diǎn)坐標(biāo)直接決定了機(jī)構(gòu)的運(yùn)動(dòng)學(xué)特性哪怕只差0.1mm對(duì)高速機(jī)構(gòu)來(lái)說(shuō)都會(huì)產(chǎn)生明顯的加速度誤差。剛體定義模塊還需要為每個(gè)構(gòu)件指定質(zhì)量、質(zhì)心位置和轉(zhuǎn)動(dòng)慣量。很多人在這個(gè)環(huán)節(jié)習(xí)慣采用軟件默認(rèn)的均勻密度材料或者粗略估一個(gè)質(zhì)量填進(jìn)去。初期驗(yàn)證可以接受但到了動(dòng)力學(xué)分析階段尤其是要做受力峰值評(píng)估時(shí)質(zhì)量特性的誤差會(huì)直接傳導(dǎo)到慣性力和支座反力的計(jì)算結(jié)果上造成設(shè)計(jì)偏危險(xiǎn)或偏保守。2.2 運(yùn)動(dòng)副與約束定義把零件“粘”成機(jī)構(gòu)的關(guān)鍵一步幾何模型建好之后零件之間還是相互獨(dú)立的需要靠運(yùn)動(dòng)副約束把它們組建成一個(gè)能夠運(yùn)動(dòng)的機(jī)構(gòu)。曲柄滑塊機(jī)構(gòu)中最常用到的運(yùn)動(dòng)副類(lèi)型有旋轉(zhuǎn)副、移動(dòng)副和圓柱副。旋轉(zhuǎn)副Revolute Joint約束兩個(gè)構(gòu)件之間只能繞公共軸線相對(duì)轉(zhuǎn)動(dòng)用在曲柄與機(jī)架之間、曲柄與連桿之間以及連桿與滑塊之間。移動(dòng)副Translational Joint約束兩個(gè)構(gòu)件之間只能沿指定方向相對(duì)滑動(dòng)用在滑塊與導(dǎo)軌/機(jī)架之間。在某些設(shè)計(jì)場(chǎng)景下連桿與滑塊之間也可以使用圓柱副Cylindrical Joint允許沿軸向的一定滑動(dòng)以模擬實(shí)際關(guān)節(jié)間隙或油膜效果。這個(gè)模塊的功能設(shè)計(jì)直接決定機(jī)構(gòu)的自由度是否正確。曲柄滑塊機(jī)構(gòu)在平面內(nèi)的自由度是1——曲柄轉(zhuǎn)角決定一切位置。如果在定義約束時(shí)多加了一個(gè)約束方向機(jī)構(gòu)就變成了過(guò)約束系統(tǒng)少加一個(gè)約束機(jī)構(gòu)就會(huì)多出不該有的自由度導(dǎo)致仿真時(shí)零件飛出去或出現(xiàn)異常漂移。我在給初學(xué)者指導(dǎo)時(shí)通常建議每創(chuàng)建一個(gè)運(yùn)動(dòng)副就問(wèn)自己三個(gè)問(wèn)題——它約束了哪個(gè)構(gòu)件相對(duì)于哪個(gè)構(gòu)件的哪幾個(gè)自由度這符合實(shí)際機(jī)械結(jié)構(gòu)中該處連接件的真實(shí)約束條件嗎整個(gè)機(jī)構(gòu)加起來(lái)的自由度等于1嗎這三個(gè)問(wèn)題想清楚約束定義就不會(huì)犯低級(jí)錯(cuò)誤。另外現(xiàn)代仿真系統(tǒng)還有一個(gè)重要功能——運(yùn)動(dòng)副摩擦定義。旋轉(zhuǎn)副中的摩擦力矩、移動(dòng)副中的滑動(dòng)摩擦力都可以通過(guò)設(shè)定摩擦系數(shù)和預(yù)緊力來(lái)模擬。別小看摩擦的作用在很多機(jī)構(gòu)中摩擦是系統(tǒng)阻尼的主要來(lái)源特別是在低速重載的沖壓設(shè)備中摩擦對(duì)機(jī)構(gòu)的停止位置和動(dòng)態(tài)響應(yīng)有決定性影響。2.3 驅(qū)動(dòng)與載荷施加工況真實(shí)性決定仿真價(jià)值約束定義完成后機(jī)構(gòu)本身還是靜止的。要讓機(jī)構(gòu)動(dòng)起來(lái)并模擬真實(shí)工作狀態(tài)需要在運(yùn)動(dòng)副或構(gòu)件上施加驅(qū)動(dòng)和載荷。驅(qū)動(dòng)Motion通常加在曲柄的旋轉(zhuǎn)副上有兩種常見(jiàn)方式一種是給定運(yùn)動(dòng)規(guī)律例如曲柄角速度恒定或按某條角速度曲線變化。這種方式適合運(yùn)動(dòng)學(xué)分析與剛體動(dòng)力學(xué)分析能夠精確計(jì)算給定運(yùn)動(dòng)規(guī)律下的機(jī)構(gòu)受載。另一種是給定驅(qū)動(dòng)力矩讓系統(tǒng)自己求解機(jī)構(gòu)在這種驅(qū)動(dòng)力矩下的實(shí)際運(yùn)動(dòng)狀態(tài)。這種方式更接近真實(shí)電機(jī)驅(qū)動(dòng)行為尤其適合分析啟停過(guò)程、負(fù)載突加情況下的動(dòng)態(tài)響應(yīng)。比如內(nèi)燃機(jī)曲柄機(jī)構(gòu)是由燃?xì)鈮毫νㄟ^(guò)活塞推動(dòng)的這時(shí)就不能直接給曲柄轉(zhuǎn)速而應(yīng)該給活塞頂部一個(gè)隨時(shí)間變化的壓力函數(shù)讓仿真系統(tǒng)自己求解曲軸轉(zhuǎn)速波動(dòng)。載荷施加需要區(qū)分力和力矩。對(duì)于曲柄滑塊機(jī)構(gòu)來(lái)說(shuō)最常見(jiàn)的工作載荷是滑塊端的工作阻力比如沖壓時(shí)的沖裁力、壓縮機(jī)中的氣體壓力、注塑機(jī)中的合模力。這些載荷通常與滑塊位置和速度相關(guān)需要通過(guò)定義函數(shù)或查表的方式施加。這里特別強(qiáng)調(diào)一個(gè)功能設(shè)計(jì)上的關(guān)鍵點(diǎn)具備函數(shù)編輯器和高階插值能力的仿真系統(tǒng)在工程實(shí)踐中的價(jià)值遠(yuǎn)遠(yuǎn)超過(guò)只支持常值載荷的系統(tǒng)。真實(shí)的設(shè)備工作載荷幾乎都是隨位置和時(shí)間變化的如果系統(tǒng)只能施加恒定載荷很多工況就模擬不了分析結(jié)果很容易失真。2.4 求解器配置數(shù)值穩(wěn)定性與精度的平衡藝術(shù)很多剛接觸多體仿真的人會(huì)忽略求解器配置直接采用默認(rèn)設(shè)置殊不知這就像一臺(tái)精密儀器用了不合適的分辨率去測(cè)量數(shù)據(jù)量和計(jì)算量都上去了但結(jié)果的可靠性反而不高。曲柄滑塊機(jī)構(gòu)仿真中求解器配置最核心的兩個(gè)參數(shù)是積分步長(zhǎng)和誤差容限。積分步長(zhǎng)直接影響求解精度和計(jì)算耗時(shí)。步長(zhǎng)過(guò)大會(huì)導(dǎo)致數(shù)值積分誤差增大甚至出現(xiàn)能量不守恒的可笑結(jié)果——比如一個(gè)不含阻尼的理想機(jī)構(gòu)仿真后系統(tǒng)的總能量竟然不斷增長(zhǎng)。步長(zhǎng)過(guò)小則會(huì)導(dǎo)致計(jì)算時(shí)間大幅增加尤其在仿真多個(gè)運(yùn)動(dòng)周期時(shí)耗時(shí)非常可觀。合理的做法是先采用較大的步長(zhǎng)做預(yù)分析觀察主要曲線形態(tài)然后在速度變化劇烈的區(qū)間滑塊換向點(diǎn)附近加密步長(zhǎng)或者在軟件中啟用自適應(yīng)步長(zhǎng)功能。誤差容限Tolerance決定求解器在每個(gè)時(shí)間步上是否認(rèn)為計(jì)算已經(jīng)收斂。容差設(shè)置得太寬松結(jié)果曲線會(huì)出現(xiàn)鋸齒狀波動(dòng)太嚴(yán)格則計(jì)算時(shí)間顯著增加且可能不收斂。工程經(jīng)驗(yàn)值一般設(shè)置在10的負(fù)4次方到負(fù)6次方之間具體需要結(jié)合模型規(guī)模來(lái)調(diào)整。求解器的選擇也很重要。對(duì)于剛體多體系統(tǒng)當(dāng)前主流求解器多為隱式積分算法如Gear算法或BDF算法。這類(lèi)算法數(shù)值穩(wěn)定性好適合強(qiáng)耦合的非線性系統(tǒng)求解。曲柄滑塊機(jī)構(gòu)的運(yùn)動(dòng)方程屬于帶有強(qiáng)約束的微分代數(shù)方程DAE不是一般的常微分方程O(píng)DE因此仿真系統(tǒng)必須支持DAE求解功能。如果在選型時(shí)發(fā)現(xiàn)某款軟件只能處理ODE系統(tǒng)而無(wú)法直接處理約束方程那它處理曲柄滑塊這種機(jī)構(gòu)問(wèn)題時(shí)就會(huì)很吃力甚至?xí)l繁報(bào)錯(cuò)。2.5 后處理與結(jié)果輸出從曲線到設(shè)計(jì)決策的最后一步后處理模塊的功能豐富程度直接影響仿真成果的表達(dá)和數(shù)據(jù)再利用率?；A(chǔ)的后處理必須包含位移、速度、加速度曲線繪制運(yùn)動(dòng)副反力曲線繪制驅(qū)動(dòng)力矩曲線繪制以及機(jī)構(gòu)運(yùn)動(dòng)動(dòng)畫(huà)回放。進(jìn)階的后處理還需要支持頻譜分析、數(shù)據(jù)導(dǎo)出到通用格式、曲線疊加對(duì)比、云圖生成等。在設(shè)計(jì)評(píng)審時(shí)動(dòng)畫(huà)回放配合特征曲線是說(shuō)服力最強(qiáng)的交付形式。評(píng)審人員不關(guān)心你求解器配置了什么算法、誤差容差設(shè)了多少他們關(guān)心的是機(jī)構(gòu)在一個(gè)運(yùn)動(dòng)周期內(nèi)最高加速度出現(xiàn)在哪個(gè)位置最大受力出現(xiàn)在哪個(gè)角度滑塊速度是否滿足節(jié)拍要求如果后處理模塊能在幾分鐘內(nèi)把這些問(wèn)題用清晰的圖表回答出來(lái)仿真的價(jià)值就得到了體現(xiàn)。另外數(shù)據(jù)導(dǎo)出功能千萬(wàn)不要忽略。很多企業(yè)做聯(lián)合仿真時(shí)需要把曲柄滑塊機(jī)構(gòu)的動(dòng)力學(xué)分析結(jié)果輸出為有限元軟件能識(shí)別的載荷邊界條件或者把速度、位置數(shù)據(jù)輸出給控制系統(tǒng)仿真工具。如果后處理模塊不支持格式靈活的數(shù)據(jù)導(dǎo)出整個(gè)仿真鏈路就會(huì)在這個(gè)節(jié)點(diǎn)卡住嚴(yán)重影響研發(fā)效率。3. 影響分析結(jié)果可信度的關(guān)鍵功能點(diǎn)約束間隙、接觸與求解細(xì)節(jié)前面章節(jié)講的是系統(tǒng)有哪些功能模塊這一節(jié)我要講的是就算功能模塊都用了結(jié)果也未必可信。曲柄滑塊機(jī)構(gòu)的多體動(dòng)力學(xué)仿真真正的難點(diǎn)和風(fēng)險(xiǎn)點(diǎn)往往集中在下面幾個(gè)細(xì)節(jié)里。這些細(xì)節(jié)能不能處理好正是區(qū)分“功能分析”和“模擬動(dòng)畫(huà)”的分水嶺。3.1 運(yùn)動(dòng)副間隙理想約束與真實(shí)機(jī)構(gòu)之間的鴻溝到目前為止絕大多數(shù)仿真分析使用理想運(yùn)動(dòng)副來(lái)建模。也就是說(shuō)仿真系統(tǒng)假設(shè)旋轉(zhuǎn)副的軸與軸承孔完美同軸沒(méi)有間隙移動(dòng)副的滑塊與導(dǎo)軌之間沒(méi)有側(cè)向空間。這種理想化假設(shè)給仿真帶來(lái)了極大的便利但也掩蓋了一個(gè)實(shí)際問(wèn)題真實(shí)機(jī)構(gòu)中運(yùn)動(dòng)副必然存在間隙。間隙的影響在低速機(jī)構(gòu)中可能不顯著但在高速工況下會(huì)被急劇放大。我在分析一個(gè)高速送料機(jī)構(gòu)時(shí)發(fā)現(xiàn)僅僅在連桿兩端的旋轉(zhuǎn)副中各加入0.05mm的間隙滑塊在換向點(diǎn)的加速度峰值就比理想模型高出了百分之四十多同時(shí)還在時(shí)域曲線中引入了明顯的高頻振蕩成分。這個(gè)現(xiàn)象在機(jī)構(gòu)學(xué)中被稱(chēng)為間隙引起的沖擊效應(yīng)。如果系統(tǒng)不具備運(yùn)動(dòng)副間隙定義功能或者只在特定模塊中支持那這類(lèi)研究就基本無(wú)法開(kāi)展。值得說(shuō)明的是如果只是為了常規(guī)強(qiáng)度校核采用理想運(yùn)動(dòng)副建模的結(jié)果偏于保守還是偏于危險(xiǎn)并沒(méi)有一個(gè)絕對(duì)答案。間隙的存在會(huì)改變載荷在各運(yùn)動(dòng)副上的分配關(guān)系有時(shí)減少峰值有時(shí)則產(chǎn)生碰撞沖擊造成更高的峰值。因此在做關(guān)鍵結(jié)構(gòu)的載荷邊界提取時(shí)建議至少在其中一個(gè)方案中考慮間隙的影響。3.2 接觸模型滑塊與導(dǎo)軌、銷(xiāo)軸與襯套的相互作用曲柄滑塊機(jī)構(gòu)中涉及接觸的位置主要有兩類(lèi)一類(lèi)是銷(xiāo)軸與襯套/軸承孔的圓柱面接觸另一類(lèi)是滑塊底面與導(dǎo)軌平面的面接觸。接觸功能在仿真系統(tǒng)中是否健全直接影響機(jī)構(gòu)受力分析的準(zhǔn)確性。在多體動(dòng)力學(xué)仿真系統(tǒng)中接觸力的計(jì)算通常采用罰函數(shù)法或拉格朗日乘子法。罰函數(shù)法通過(guò)定義接觸剛度和穿透量來(lái)計(jì)算接觸力優(yōu)點(diǎn)是計(jì)算效率高、實(shí)現(xiàn)簡(jiǎn)單缺點(diǎn)是接觸剛度選取不當(dāng)會(huì)造成結(jié)果失真或數(shù)值振蕩。拉格朗日乘子法精度更高但計(jì)算消耗更大且對(duì)接觸狀態(tài)的突變比較敏感。在工程仿真中我一般建議分段處理分析機(jī)構(gòu)整體受力時(shí)可以不啟用接觸模型直接采用運(yùn)動(dòng)副約束反力作為參考但分析局部磨損、沖擊載荷和振動(dòng)響應(yīng)時(shí)必須開(kāi)啟接觸功能并且在正式分析前做一個(gè)接觸剛度敏感性驗(yàn)證。簡(jiǎn)單的做法是把接觸剛度提高一倍和降低一半各計(jì)算一次如果結(jié)果峰值變化在5%以?xún)?nèi)說(shuō)明接觸剛度選取基本合理如果變化超過(guò)20%就要回頭檢查接觸參數(shù)設(shè)置。接觸阻尼的設(shè)置也同樣重要。阻尼過(guò)小接觸過(guò)程中的高頻振蕩無(wú)法衰減結(jié)果曲線像噪聲一樣阻尼過(guò)大則會(huì)人為地吸收能量導(dǎo)致接觸力偏低。實(shí)際處理中通常先根據(jù)材料恢復(fù)系數(shù)推算接觸阻尼再根據(jù)仿真曲線微調(diào)。3.3 仿真步長(zhǎng)與收斂性那些“跑飛”的模型到底是誰(shuí)的鍋在曲柄滑塊機(jī)構(gòu)仿真中新手遇到最多的報(bào)錯(cuò)就是“求解失敗”或“計(jì)算發(fā)散”。模型在某一時(shí)刻突然飛出去或者速度曲線瞬間飆升到離譜的數(shù)值通常不是系統(tǒng)功能缺失而是模型或求解參數(shù)出了問(wèn)題。常見(jiàn)的發(fā)散原因有以下幾類(lèi)一是初始狀態(tài)不合理。機(jī)構(gòu)初始位置與驅(qū)動(dòng)函數(shù)不匹配比如給定了一個(gè)在初始時(shí)刻加速度無(wú)窮大的驅(qū)動(dòng)函數(shù)系統(tǒng)在第一步就無(wú)法收斂。最典型的錯(cuò)誤是在曲柄初始角速度不為零的前提下又讓角速度從零加速到目標(biāo)值這兩組條件沖突求解器根本無(wú)法滿足。二是約束冗余沖突。在一個(gè)旋轉(zhuǎn)副上同時(shí)定義驅(qū)動(dòng)和額外約束或者滑塊同時(shí)被兩個(gè)移動(dòng)副限位導(dǎo)致系統(tǒng)方程過(guò)約束。三是接觸剛度與時(shí)間步長(zhǎng)不匹配。接觸力對(duì)穿透非常敏感如果時(shí)間步長(zhǎng)太大求解器可能在一步之內(nèi)就出現(xiàn)過(guò)度穿透造成接觸力爆炸式增長(zhǎng)進(jìn)而把機(jī)構(gòu)彈飛。排查發(fā)散問(wèn)題的方法其實(shí)很簡(jiǎn)單先降低驅(qū)動(dòng)速度跑通運(yùn)動(dòng)學(xué)模型確認(rèn)約束和幾何沒(méi)有問(wèn)題再加接觸逐步增加轉(zhuǎn)速和載荷每一步確認(rèn)曲線形態(tài)正常后再進(jìn)入下一步。這種由簡(jiǎn)到繁的建模策略能幫你快速定位導(dǎo)致發(fā)散的具體功能環(huán)節(jié)而不是在一個(gè)復(fù)雜的失敗模型里大海撈針。4. 仿真系統(tǒng)在參數(shù)化分析與優(yōu)化設(shè)計(jì)中的功能支撐曲柄滑塊機(jī)構(gòu)的研發(fā)通常不是一錘子買(mǎi)賣(mài)。設(shè)計(jì)前期要比較不同桿長(zhǎng)比的方案中期要做尺寸公差敏感性分析后期還要結(jié)合樣機(jī)測(cè)試結(jié)果修正模型。仿真系統(tǒng)的參數(shù)化分析能力和優(yōu)化功能在這些環(huán)節(jié)中發(fā)揮著不可替代的作用。4.1 參數(shù)化建模如何讓機(jī)構(gòu)方案對(duì)比從一周縮短到半天想象一個(gè)典型的方案設(shè)計(jì)場(chǎng)景客戶要求滑塊行程100mm機(jī)構(gòu)急回特性系數(shù)K1.4同時(shí)要求機(jī)構(gòu)的外廓尺寸盡量小。這三個(gè)需求交織在一起曲柄半徑、連桿長(zhǎng)度、偏心距之間存在多組可行解。如果每個(gè)方案都在三維CAD中重新畫(huà)圖一個(gè)方案畫(huà)一天一周時(shí)間只能比五個(gè)方案。有了參數(shù)化建模功能就不一樣了。把曲柄半徑、連桿長(zhǎng)度和偏距設(shè)置為設(shè)計(jì)變量把行程和急回特性系數(shù)設(shè)置為優(yōu)化目標(biāo)仿真系統(tǒng)可以在很短時(shí)間內(nèi)批量計(jì)算幾十組參數(shù)組合并自動(dòng)輸出滿足約束條件的可行域。這個(gè)環(huán)節(jié)真正把工程師從重復(fù)建模中解放出來(lái)把精力集中在理解參數(shù)與性能之間的耦合關(guān)系上。在參數(shù)化過(guò)程中有一點(diǎn)值得特別提醒參數(shù)化參數(shù)不僅包括幾何尺寸還應(yīng)該包括質(zhì)量屬性。因?yàn)閹缀纬叽绺淖兒蟾鳂?gòu)件的質(zhì)量和轉(zhuǎn)動(dòng)慣量也會(huì)改變進(jìn)而影響系統(tǒng)的動(dòng)力學(xué)響應(yīng)。手工同步更新質(zhì)量和慣量非常繁瑣且容易遺漏只有系統(tǒng)支持完整參數(shù)化幾何和質(zhì)量屬性同步聯(lián)動(dòng)動(dòng)力學(xué)的優(yōu)化結(jié)果才是可信的。4.2 目標(biāo)函數(shù)與約束條件的設(shè)定邏輯做機(jī)構(gòu)優(yōu)化不能只看單一指標(biāo)。比如一味降低連桿質(zhì)量可以減少慣性力但連桿剛度也會(huì)下降實(shí)際結(jié)構(gòu)中可能引起失穩(wěn)或變形過(guò)大。所以在優(yōu)化設(shè)計(jì)中目標(biāo)函數(shù)要盡量貼近工程綜合需求。以曲柄滑塊機(jī)構(gòu)為例一個(gè)合理的優(yōu)化問(wèn)題可以設(shè)定為在滿足滑塊行程、最小傳動(dòng)角和連桿最大應(yīng)力約束的前提下最小化曲柄的峰值驅(qū)動(dòng)力矩。這個(gè)目標(biāo)函數(shù)把運(yùn)動(dòng)學(xué)指標(biāo)行程、傳力特性傳動(dòng)角和結(jié)構(gòu)強(qiáng)度應(yīng)力都納入了約束體系優(yōu)化結(jié)果才具有工程實(shí)用意義。不過(guò)這里也要潑一盆冷水目前的通用多體動(dòng)力學(xué)仿真系統(tǒng)自帶優(yōu)化算法在機(jī)械工程領(lǐng)域的實(shí)用性參差不齊。對(duì)于簡(jiǎn)單的參數(shù)掃描和響應(yīng)面分析內(nèi)置功能基本夠用但涉及復(fù)雜的多目標(biāo)優(yōu)化比如同時(shí)最小化沖擊載荷和機(jī)構(gòu)體積更推薦的做法是將多體仿真系統(tǒng)與專(zhuān)業(yè)優(yōu)化工具聯(lián)合使用通過(guò)數(shù)據(jù)接口傳遞分析結(jié)果由外部?jī)?yōu)化算法迭代尋優(yōu)。4.3 模型標(biāo)定把仿真結(jié)果拉回真實(shí)世界的功能接口任何仿真系統(tǒng)的功能強(qiáng)弱最終都要通過(guò)模型與實(shí)際測(cè)量數(shù)據(jù)的一致性來(lái)驗(yàn)證。工程上把這個(gè)過(guò)程稱(chēng)為模型標(biāo)定。標(biāo)定通常需要實(shí)測(cè)數(shù)據(jù)包括滑塊位移曲線、關(guān)鍵點(diǎn)加速度曲線、驅(qū)動(dòng)力矩或運(yùn)動(dòng)副反力信號(hào)。仿真系統(tǒng)提供的曲線導(dǎo)出功能要與常見(jiàn)數(shù)據(jù)采集格式兼容并能方便地進(jìn)行時(shí)域?qū)R。我在標(biāo)定一臺(tái)高速?zèng)_壓機(jī)構(gòu)的模型時(shí)發(fā)現(xiàn)仿真與測(cè)試在換向點(diǎn)附近的加速度曲線差異很大。排查到最后問(wèn)題出在兩個(gè)環(huán)節(jié)一是連桿在變形狀態(tài)下質(zhì)心位置與剛體模型不一致二是旋轉(zhuǎn)副的摩擦在低速換向區(qū)間的非線性行為沒(méi)有在仿真中準(zhǔn)確模擬。調(diào)整了連桿的等效質(zhì)量和當(dāng)量摩擦系數(shù)后仿真曲線與實(shí)測(cè)曲線的吻合度有了顯著提升。模型標(biāo)定是一項(xiàng)非?？简?yàn)工程師功力與系統(tǒng)靈活性的工作。一個(gè)優(yōu)秀的仿真系統(tǒng)應(yīng)該在建模模塊中允許用戶自定義本構(gòu)關(guān)系和摩擦特性在后處理中提供多曲線對(duì)比和誤差統(tǒng)計(jì)功能同時(shí)支持將標(biāo)定參數(shù)保存為模板方便以后同類(lèi)機(jī)構(gòu)直接套用。5. 從實(shí)際項(xiàng)目出發(fā)曲柄滑塊仿真系統(tǒng)的常見(jiàn)“踩坑”復(fù)盤(pán)這一節(jié)我梳理幾個(gè)真實(shí)項(xiàng)目中的典型問(wèn)題每一個(gè)都對(duì)應(yīng)仿真系統(tǒng)某一特定功能的使用誤區(qū)。寫(xiě)出來(lái)是想讓大家少走彎路遇到類(lèi)似現(xiàn)象時(shí)能快速定位問(wèn)題根源。5.1 初始位形“卡死”運(yùn)動(dòng)學(xué)仿真跑不動(dòng)的排查過(guò)程某次做一個(gè)偏心曲柄滑塊機(jī)構(gòu)的運(yùn)動(dòng)學(xué)驗(yàn)證。模型建完后一運(yùn)行就報(bào)錯(cuò)提示“機(jī)構(gòu)奇異”或者“位置初始化失敗”。我當(dāng)時(shí)的第一反應(yīng)是檢查運(yùn)動(dòng)副方向是否定義正確。隨后逐一檢查旋轉(zhuǎn)副軸線方向和移動(dòng)副滑動(dòng)方向并未發(fā)現(xiàn)明顯錯(cuò)誤。接著把目標(biāo)轉(zhuǎn)速?gòu)暮愣ǜ臑閺牧憔徛龁?dòng)問(wèn)題依舊。后來(lái)我把模型簡(jiǎn)化到只保留幾何線和鉸點(diǎn)通過(guò)手算機(jī)構(gòu)在初始轉(zhuǎn)角下的裝配條件發(fā)現(xiàn)初始位形處于機(jī)構(gòu)的極限位置附近——連桿與曲柄近乎共線。在這個(gè)奇異位形下機(jī)構(gòu)從靜止開(kāi)始運(yùn)動(dòng)的速度方向極不穩(wěn)定數(shù)值求解困難。解決辦法很簡(jiǎn)單把初始曲柄轉(zhuǎn)角調(diào)整到一個(gè)非奇異位置比如與連桿呈75度夾角仿真器很快跑通然后再通過(guò)瞬態(tài)分析過(guò)渡到目標(biāo)工況。這個(gè)問(wèn)題的本質(zhì)是啟動(dòng)位置選擇不當(dāng)而不是系統(tǒng)功能缺陷。此后我在所有曲柄滑塊模型中都固定了一個(gè)好習(xí)慣——初始位形避開(kāi)機(jī)構(gòu)的死點(diǎn)區(qū)和近似死點(diǎn)區(qū)。5.2 質(zhì)量屬性缺失引起的“幽靈慣性力”還有一個(gè)我印象很深的案例。一個(gè)同事做機(jī)構(gòu)受力分析跑出來(lái)的運(yùn)動(dòng)副反力曲線峰值比理論估算大了將近一倍而且反力方向非常異常。他反復(fù)檢查載荷設(shè)置甚至懷疑是軟件計(jì)算出錯(cuò)。我把模型調(diào)出來(lái)看了一下發(fā)現(xiàn)他把所有零件的材料密度都保留了默認(rèn)值而默認(rèn)值對(duì)應(yīng)的是一種幾乎不存在的超輕材料。這意味著實(shí)際零件質(zhì)量遠(yuǎn)高于模型賦值慣性力被嚴(yán)重低估。但奇怪的是他的結(jié)果偏大這不符合邏輯。進(jìn)一步檢查后才發(fā)現(xiàn)他在CAD導(dǎo)入時(shí)沒(méi)有正確設(shè)置單位零件把毫米模型的體積按米為單位計(jì)算導(dǎo)致質(zhì)量被放大了十的九次方倍。這是一個(gè)非常低級(jí)但又極其常見(jiàn)的錯(cuò)誤。多體動(dòng)力學(xué)軟件對(duì)單位系統(tǒng)的處理非常講究很多系統(tǒng)內(nèi)部采用SI單位制但允許用戶在建模時(shí)使用其他單位。如果你從CAD系統(tǒng)導(dǎo)入的是毫米單位模型而仿真系統(tǒng)默認(rèn)使用米單位導(dǎo)入時(shí)不統(tǒng)一幾何尺寸乘了一千倍質(zhì)量乘了十的九次方倍整個(gè)仿真結(jié)果自然面目全非。進(jìn)入正式分析前先粗略核對(duì)一下模型的總體質(zhì)量和各構(gòu)件的質(zhì)心位置是一種成本極低但收益極高的良好習(xí)慣。5.3 后處理曲線“毛刺”與采樣頻率陷阱在查看動(dòng)力學(xué)分析結(jié)果時(shí)很多工程師會(huì)遇到一個(gè)困擾速度或加速度曲線在平滑的主趨勢(shì)上疊加了大量高頻毛刺。有人會(huì)懷疑模型出了問(wèn)題也有人嘗試通過(guò)增大積分誤差來(lái)解決。毛刺的來(lái)源主要有三個(gè)第一是接觸碰撞引起的真實(shí)物理響應(yīng)尤其在有間隙的機(jī)構(gòu)中這種高頻振蕩是真實(shí)存在的不能為了“好看”而強(qiáng)行濾波第二是求解器容差過(guò)大導(dǎo)致數(shù)值誤差累積這種毛刺幅度通常不大但無(wú)規(guī)律可循適當(dāng)收緊容差可以減少第三是后處理輸出采用了過(guò)高的采樣頻率把求解器數(shù)值迭代過(guò)程中的微小振蕩也一并繪制了出來(lái)這種情況下可以降低輸出采樣頻率或進(jìn)行合理的低通濾波。判斷毛刺是否有效的簡(jiǎn)單方法是減小步長(zhǎng)重新計(jì)算觀察毛刺是否隨步長(zhǎng)變化。如果毛刺幅度和形態(tài)基本不隨步長(zhǎng)變化可以判定為物理現(xiàn)象如果毛刺在減小步長(zhǎng)后大幅減弱則屬于數(shù)值噪聲。這個(gè)判斷方法在多個(gè)項(xiàng)目中幫我避免了對(duì)真實(shí)故障的誤判和對(duì)數(shù)值噪聲的過(guò)度解讀。我在處理沖壓機(jī)構(gòu)的振動(dòng)問(wèn)題時(shí)正是靠這個(gè)判斷方法確認(rèn)了高頻毛刺屬于運(yùn)動(dòng)副間隙引發(fā)的物理沖擊順藤摸瓜找到了連桿襯套磨損異常的根本原因最終通過(guò)更換自潤(rùn)滑襯套和調(diào)整配合間隙解決了問(wèn)題。6. 最后再分享三點(diǎn)關(guān)于仿真系統(tǒng)的選型與使用體會(huì)第一點(diǎn)體會(huì)關(guān)于仿真系統(tǒng)的定位。曲柄滑塊機(jī)構(gòu)雖然結(jié)構(gòu)簡(jiǎn)單但它涵蓋了多體動(dòng)力學(xué)仿真的大部分核心功能需求旋轉(zhuǎn)副、移動(dòng)副、接觸、摩擦、參數(shù)化、后處理。一套仿真系統(tǒng)如果能把曲柄滑塊機(jī)構(gòu)分析得準(zhǔn)確、穩(wěn)定、高效那它應(yīng)對(duì)更復(fù)雜的平面或空間機(jī)構(gòu)問(wèn)題通常也不會(huì)太吃力。因此把曲柄滑塊機(jī)構(gòu)作為仿真系統(tǒng)的能力和精度測(cè)試基準(zhǔn)是一個(gè)很務(wù)實(shí)的做法。第二點(diǎn)體會(huì)關(guān)于過(guò)程管理。很多工程師過(guò)分關(guān)注最終結(jié)果曲線卻忽略了仿真過(guò)程的規(guī)范管理。建議在仿真項(xiàng)目中建立模板化的建模規(guī)范統(tǒng)一的單位制、統(tǒng)一的構(gòu)件命名規(guī)則、統(tǒng)一的鉸點(diǎn)坐標(biāo)系方向、統(tǒng)一的輸出文件命名格式。這些規(guī)范看起來(lái)花時(shí)間但在項(xiàng)目協(xié)作和模型復(fù)用階段能節(jié)省大量返工時(shí)間。我自己就吃過(guò)“模型文件名帶_2_3_final”這種坑后來(lái)強(qiáng)制規(guī)范之后項(xiàng)目流轉(zhuǎn)順暢了很多。第三點(diǎn)體會(huì)關(guān)于工程師判斷力的不可替代性。再完善的仿真系統(tǒng)功能分析能力也不能替代工程師對(duì)機(jī)構(gòu)的物理直覺(jué)。每次仿真前先用靜力學(xué)或運(yùn)動(dòng)學(xué)手算估算關(guān)鍵值每次仿真后再反問(wèn)自己“這個(gè)結(jié)果合不合理”“為什么會(huì)出現(xiàn)這個(gè)峰值”“如果結(jié)構(gòu)參數(shù)變化結(jié)果趨勢(shì)是否符合物理判斷”。這種習(xí)慣能幫你攔住大量因?yàn)閰?shù)誤設(shè)、單位錯(cuò)誤、約束矛盾而導(dǎo)致的失敗仿真而不是把時(shí)間浪費(fèi)在為一個(gè)錯(cuò)誤模型反復(fù)調(diào)參上。多體動(dòng)力學(xué)仿真的意義不在于讓結(jié)果看起來(lái)高深而在于讓機(jī)構(gòu)設(shè)計(jì)從“憑經(jīng)驗(yàn)”過(guò)渡到“有依據(jù)”工具越強(qiáng)大使用者越要保持清醒。