亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        流固耦合作用下制動氣室動態(tài)特性建模與試驗

        2021-11-17 12:04:26朱茂桃
        計算機(jī)仿真 2021年9期
        關(guān)鍵詞:有限元模型

        邵 瑜,朱茂桃

        (江蘇大學(xué)汽車與交通工程學(xué)院,江蘇 鎮(zhèn)江 212013)

        1 引言

        由于氣壓制動的制動力矩大、操縱輕便,絕大部分載貨汽車及客車安裝氣壓制動系統(tǒng)[1-2]。制動氣室作為氣壓制動系統(tǒng)中的主要裝置,用于連接汽車制動回路與制動執(zhí)行機(jī)構(gòu)。制動氣室的壓力特性影響制動器的制動力輸出和制動性能[3],而目前國內(nèi)在制動氣室設(shè)計開發(fā)方面的研究較少,因此,研究制動氣室的動態(tài)特性對高性能氣壓制動系統(tǒng)的開發(fā)具有重要意義。

        氣壓制動系統(tǒng)已經(jīng)得到了廣泛的研究,許多學(xué)者建立了相應(yīng)的模型預(yù)測制動腔內(nèi)的瞬態(tài)壓力。周佳瑋等[4]根據(jù)制動氣室工作原理,建立制動氣室數(shù)學(xué)模型,在此基礎(chǔ)上采用AMESim軟件分析其在不同壓力源下的動態(tài)相應(yīng)特性及推桿行程變化。朱薔等[5]依據(jù)氣體流動的狀態(tài)方程、連續(xù)性方程及動力學(xué)方程對制動氣室輸出壓力特性進(jìn)行數(shù)學(xué)描述,并進(jìn)行試驗驗證。李興麗等[6]通過對制動氣室的流量特性方程、狀態(tài)方程及活塞盤的運動方程無因次化,得到制動氣室無因次解析模型,并通過試驗進(jìn)行驗證。

        雖然上述學(xué)者們在制動氣室動態(tài)特性方面做了很多研究,但他們的研究主要是對其進(jìn)行數(shù)學(xué)建模,缺乏對流體和固體之間相互耦合作用的考慮,誤差較大。本文基于流固耦合有限元分析的理論與方法,運用ADINA仿真軟件建立了制動氣室的流固耦合有限元模型,分析制動氣室在壓縮空氣作用下的動態(tài)特性,并與試驗結(jié)果對比,驗證了建模方法的正確性。文中對制動氣室的建模方法和試驗方法,可為實際產(chǎn)品的開發(fā)提供理論支持。

        2 制動氣室結(jié)構(gòu)及試驗方法

        制動氣室的基本結(jié)構(gòu)如圖1所示,主要由進(jìn)氣口、橡膠膜片、殼體、回位彈簧、推桿、夾箍和螺栓等組成。橡膠膜片將氣室分成兩個腔室,當(dāng)汽車制動時,壓縮空氣從進(jìn)氣口進(jìn)入制動氣室壓縮腔,在壓縮空氣的作用下使膜片發(fā)生變形,推動推桿,并帶動制動調(diào)整臂,將制動蹄摩擦片壓向制動鼓從而產(chǎn)生制動;當(dāng)推桿達(dá)到極限位移時,制動氣室仍將持續(xù)沖入壓縮氣體,直至壓力達(dá)到最大;制動結(jié)束,制動氣室壓縮腔內(nèi)的氣體通過快放閥釋放,在回位彈簧的作用下氣室壓縮腔的體積不斷減小,直至回歸原位。

        圖1 制動氣室結(jié)構(gòu)圖

        根據(jù)QCT790-2007《制動氣室性能要求及臺架試驗方法》,建立制動氣室實驗裝置,如圖2和圖3所示。壓縮機(jī)為制動系統(tǒng)提供壓縮空氣,可最大提供0.9 Mpa的壓縮空氣,與壓縮機(jī)集成一個容量為90公升的儲氣罐,并使用調(diào)壓閥調(diào)節(jié)儲氣罐輸出的壓力,壓縮空氣通過軟管輸送至制動氣室。壓力傳感器安裝在制動氣室前端蓋打孔處,用于測量壓縮腔內(nèi)的空氣壓力變化情況。采用激光位移傳感器測得推桿的位移變化情況。所得到的電磁閥信號、制動氣室壓強信號、推桿位移信號通過數(shù)據(jù)采集儀,最終送達(dá)信號處理系統(tǒng)。

        圖2 試驗現(xiàn)場圖

        圖3 試驗系統(tǒng)布置圖

        3 流固耦合理論

        流固耦合運動過程中,結(jié)構(gòu)在流體載荷的作用下發(fā)生運動,其變形進(jìn)而改變流體載荷的大小和分布[7]。為分析壓縮空氣作用下制動氣室的動態(tài)響應(yīng),在流體-結(jié)構(gòu)界面的運動學(xué)和動力學(xué)條件為[8]

        df=ds

        (1)

        nτf=nτs

        (2)

        其中df和ds分別表示流體和結(jié)構(gòu)的邊界位移,τf和τs分別表示流體和結(jié)構(gòu)的邊界應(yīng)力,n為耦合邊界上的外法線矢量。

        根據(jù)動力條件,流體牽引力沿流體-結(jié)構(gòu)界面整合成流體力,施加于結(jié)構(gòu)節(jié)點上,其大小為

        (3)

        其中,hd表示固體位移的虛量。

        耦合系統(tǒng)的解向量為X=(Xf,Xs),其中Xf和Xs分別為分別定義在流體節(jié)點和固體節(jié)點上的流體和結(jié)構(gòu)域的解向量。此時ds=ds(Xs),τf=τf(Xf)。耦合流固耦合體系的有限元方程可以表示為

        (4)

        其中,F(xiàn)f和Fs分別為流體和結(jié)構(gòu)的流固耦合邊界上關(guān)于時間積分的有限元代數(shù)方程組。

        運用ADINA軟件進(jìn)行流固耦合計算,分別在ANDINA結(jié)構(gòu)模塊和流體模塊中建立結(jié)構(gòu)模型和流體模型,將兩個模型放入ADINA流固耦合求解器中進(jìn)行求解。采用Newton-Raphson迭代算法對兩個模型進(jìn)行耦合求解即對式(4)進(jìn)行求解,得到流場和結(jié)構(gòu)場的位移情況。

        4 制動氣室流固耦合模型建立

        4.1 有限元模型

        由于結(jié)構(gòu)和流體幾何模型均具有軸對稱特點,且載荷和邊界條件也具有軸對稱特點,因此將三維計算分析簡化成二維計算分析。同時建立結(jié)構(gòu)幾何模型和流體幾何模型,在建立流固耦合計算模型時,在軟件結(jié)構(gòu)模塊中劃分結(jié)構(gòu)部分的網(wǎng)格,在軟件流體模塊中劃分流體部分的網(wǎng)格。

        圖4 制動氣室流固耦合有限元模型

        4.2 材料參數(shù)

        流體區(qū)域的材料為空氣。流體的邊界是完全封閉,當(dāng)施加進(jìn)氣壓力時,流固耦合邊界發(fā)生變形,因此需要考慮空氣的可壓縮性。對于可壓縮流體模型,材料參數(shù)包括:粘度、等壓比熱、等容比熱、熱傳導(dǎo)系數(shù)、參考溫度。流體的密度是由與溫度和壓力相關(guān)的狀態(tài)方程計算,因此不考慮體積模量和密度[9]。

        表1 流體材料參數(shù)

        橡膠材料為超彈性材料,一般通過某種簡單的變形情況下的應(yīng)力應(yīng)變屬性描述橡膠力學(xué)性能。采用回歸分析,以一個適當(dāng)?shù)膽?yīng)變能函數(shù)實驗得到的應(yīng)力應(yīng)變數(shù)據(jù)進(jìn)行擬合,并將擬合得到的有關(guān)數(shù)據(jù)作為有限元分析的輸入[10]。本文采用Ogden函數(shù)描述橡膠的力學(xué)性能,其多項式應(yīng)變能函數(shù)為:

        (5)

        其中λ1、λ2、λ3是三個主拉伸比,μi、αi為由實驗數(shù)據(jù)決定的材料常數(shù),Di為表征可壓縮性的體積模量。

        采用GOTECH AI-7000 M拉伸試驗機(jī)對橡膠材料進(jìn)行單軸拉伸、壓縮試驗,用于拉伸實驗的啞鈴形試片,其工作標(biāo)線間距25±0.5 mm,工作部分寬度為6.0±0.4 mm,厚度為2.00±0.03 mm;用于壓縮實驗的圓柱形試塊,其直徑為29 mm,厚度為12.5 mm。表2為橡膠材料的應(yīng)力應(yīng)變數(shù)據(jù),軟件可自動擬合應(yīng)力應(yīng)變數(shù)據(jù),生成Ogden材料本構(gòu)模型系數(shù)及體積模量,并繪制擬合曲線。

        表2 橡膠材料應(yīng)力應(yīng)變數(shù)據(jù)

        采用微機(jī)控制彈簧拉壓試驗機(jī)對芯軸回位彈簧進(jìn)行壓縮試驗,每壓縮10 mm記錄一次負(fù)載值,彈簧的計算剛度偏差在5 %左右,假設(shè)芯軸回位彈簧為線性彈簧,取剛度的平均值為2.882 N/mm。

        4.3 邊界條件

        圖5為固體模型的邊界條件,對固體模型中的推桿施加由安裝過程中產(chǎn)生的彈簧預(yù)應(yīng)力,并約束推盤除運動方向外所有的自由度,固體模型與流體模型中的流固耦合接觸面一一對應(yīng)。對于軸對稱模型,不論是結(jié)構(gòu)模型還是流體模型,對稱中心線可以不施加任何邊界,軟件可以識別對稱中心線并施加相應(yīng)的約束邊界。

        圖5 結(jié)構(gòu)模型的邊界條件

        對于可壓縮流體的流固耦合計算,流體的壓力、溫度和密度的初始條件必須滿足狀態(tài)方程,如下式

        p=(cp-cv)ρθ

        (6)

        其中,p為流體初始壓力,pa;θ為流體溫度,K;ρ為流體初始密度,kg/m3。

        圖6為流體模型的邊界條件,設(shè)置流體的初始壓強為0.1 Mpa,溫度為338 K。流體壁面定義為非滑移壁面,對進(jìn)氣口施加壓力載荷。

        圖6 流體模型的邊界條件

        流體模型采用Turbulent K-Epsilon 湍流模型,橡膠膜片上下表面設(shè)置為流固耦合接觸面, 結(jié)構(gòu)仿真模型采用隱式計算方法,流體仿真模型采用瞬態(tài)仿真計算方法,耦合時間步設(shè)為1×10-5s。

        5 制動氣室動態(tài)特性分析

        5.1 制動氣室動態(tài)特性仿真與試驗結(jié)果對比

        圖7和圖8為制動氣室在0.4 Mpa的進(jìn)氣壓力下,制動氣室壓力曲線和推桿位移曲線。由圖7可知。制動氣室的運動可分為三個部分。運動初期,制動氣室內(nèi)的壓強快速增加;達(dá)到0.1368 MPa時在壓強的作用下橡膠膜片發(fā)生變形,推動推桿進(jìn)行運動,隨著推桿的運動,制動氣室內(nèi)的體積快速增大,制動氣室內(nèi)的壓強在一定范圍內(nèi)發(fā)生波動;當(dāng)推桿到達(dá)極限位置時,制動氣室的壓強迅速提高,直至最大。制動氣室在充氣過程中的計算結(jié)果與實驗結(jié)果的曲線總體趨勢一致,證明本文建模方法與計算方法的正確性。

        圖7 制動氣室壓強變化曲線

        圖8 制動氣室推桿位移曲線

        壓力響應(yīng)時間為制動氣室的壓力達(dá)到穩(wěn)態(tài)最大壓力值的75%時所經(jīng)歷的時間,表3為制動氣室壓力響應(yīng)時間情況,由表3可知,不同供氣壓力下制動氣室壓力響應(yīng)時間的試驗值和計算值誤差均小于15 %,所見的計算模型精度較高。

        表3 制動氣室壓力響應(yīng)時間計算值與試驗結(jié)果

        5.2 制動氣室流場特性仿真結(jié)果

        制動氣室內(nèi)部流場試驗測試與理論分析較為困難,通過流固耦合有限元分析計算可以直觀地獲得制動氣室內(nèi)部流場的分布情況,同時能掌握制動氣室結(jié)構(gòu)變化對其性能的影響情況

        圖10為不同推桿位移情況下制動氣室流場情況。從圖中可以看出:推桿位移從20 mm增加到60 mm,制動氣室的壓強變化不大,這是由于在壓縮空氣的作用下,制動氣室密封腔內(nèi)的體積變大,橡膠膜片變形較復(fù)雜,制動氣室的壓強緩慢地增加。入口處流體壓力場急劇變化,壓力梯度大,在入口側(cè)存在小范圍地負(fù)壓區(qū)。

        圖10 不同推桿位移情況下流體壓力場

        5.3 入口直徑對制動氣室動態(tài)響應(yīng)的影響

        圖9和圖10為0.4 Mpa供氣壓力下不同入口直徑時制動氣室壓力響應(yīng)曲線及推桿位移曲線。由圖可知,入口直徑對制動氣室壓力響應(yīng)的影響較大,入口直徑越大,制動氣室的壓力響應(yīng)時間越短。入口直徑為12 mm時,制動氣室的壓力響應(yīng)時間最短。入口直徑為8 mm時,制動氣室的壓力響應(yīng)時間最長。然而隨著入口直徑的增大,入口直徑的變大對壓力響應(yīng)的影響越不明顯。

        圖11 不同入口直徑下制動氣室壓強曲線

        圖12 不同入口直徑下制動氣室推桿位移曲線

        6 結(jié)論

        1)本文建立了膜片制動氣室的流固耦合有限元模型,計算分析其動態(tài)特性,并與試驗結(jié)果進(jìn)行對比分析,驗證本文建模方法和計算方法的正確性。

        2)分析了入口直徑對制動氣室壓力響應(yīng)的影響。結(jié)果表明,提高入口直徑可以有效增加制動氣室的壓力響應(yīng),但入口直徑越大,提高效果越不明顯。

        3)本文在諸多假設(shè)與簡化條件的基礎(chǔ)上展開建模,在往后的研究中應(yīng)進(jìn)一步完善制動氣室的計算模型,更加精準(zhǔn)地預(yù)測制動氣室地動態(tài)特性。

        猜你喜歡
        有限元模型
        一半模型
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
        新型有機(jī)玻璃在站臺門的應(yīng)用及有限元分析
        基于有限元的深孔鏜削仿真及分析
        基于有限元模型對踝模擬扭傷機(jī)制的探討
        3D打印中的模型分割與打包
        FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
        磨削淬硬殘余應(yīng)力的有限元分析
        基于SolidWorks的吸嘴支撐臂有限元分析
        女人张开腿让男桶喷水高潮 | 少妇aaa级久久久无码精品片| 曰本无码人妻丰满熟妇5g影院| 国产精品国产午夜免费福利看| 91人妻人人做人人爽九色| 国产精品一区二区性色| 欧美黑吊大战白妞| 在线精品无码一区二区三区| 给我播放的视频在线观看| 干日本少妇一区二区三区| 中字幕人妻一区二区三区| 国产自精品| 久久精品网站免费观看| 日本高清一级二级三级| 久久国内精品自在自线图片| 亚洲精品第一国产麻豆 | 亚洲中文字幕一区二区在线| 好紧好爽免费午夜视频| 国产精品久久婷婷六月丁香| 亚洲色欲色欲大片WWW无码| 亚洲国产精品日韩av专区| 2019最新中文字幕在线观看| 国产av无码专区亚洲av手机麻豆 | 日本伊人精品一区二区三区 | 欧美一级三级在线观看| 日本精品人妻一区二区三区| 日本免费大片一区二区| 久久露脸国产精品| 中文字幕精品一二三区| 日本一区二区三区清视频| 色费女人18毛片a级毛片视频 | 人人妻人人澡人人爽欧美一区九九| 久久青草伊人精品| 日本高清无卡一区二区三区| 三级全黄裸体| 国产精品久久久久国产a级| 99久久久久久亚洲精品| 户外精品一区二区三区| 中文字幕天天躁日日躁狠狠躁免费 | 中日韩精品视频在线观看| 加勒比在线一区二区三区|