胡 耿，張 定，譚武中

（1.中國(guó)航發(fā)湖南動(dòng)力機(jī)械研究所，湖南 株洲 412002；2.直升機(jī)傳動(dòng)技術(shù)國(guó)防科技重點(diǎn)實(shí)驗(yàn)室，湖南 株洲 412002）

1 引言

膜盤聯(lián)軸器是一種金屬?gòu)椥月?lián)軸器，具有可靠性高、無(wú)需潤(rùn)滑、補(bǔ)償量大等優(yōu)點(diǎn)，特別適用于高速重載場(chǎng)合，在航空、船舶等領(lǐng)域得到了廣泛應(yīng)用。膜盤型面是聯(lián)軸器的關(guān)鍵彈性部位，需要具備足夠的彈性以補(bǔ)償兩端轉(zhuǎn)子的不對(duì)中［1-2］。由于工作時(shí)型面應(yīng)力狀態(tài)復(fù)雜，在實(shí)際使用中以高周疲勞失效為主［3］，因此在設(shè)計(jì)時(shí)需要重點(diǎn)考慮膜盤型面的疲勞強(qiáng)度。

目前對(duì)膜盤型面的研究主要以靜態(tài)應(yīng)力為主，對(duì)膜盤疲勞的研究較為少見：部分學(xué)者分別應(yīng)用Timoshenko 模型和有限元方法研究了在不同載荷工況下膜盤型面的靜態(tài)應(yīng)力分布［4-6］，并在靜強(qiáng)度分析的基礎(chǔ)上，開展了膜盤優(yōu)化研究［7-8］；Calistrat、楊憲峰等通過(guò)試驗(yàn)手段研究了膜盤的疲勞特性［3，9］。

為提高膜盤抗疲勞性能，通過(guò)分析膜盤在復(fù)合工況下的應(yīng)力，確定了膜盤高周疲勞應(yīng)力計(jì)算方法，并在此基礎(chǔ)上研究了型面參數(shù)對(duì)膜盤高周疲勞的影響，然后以高周疲勞安全系數(shù)最大為目標(biāo)，建立了基于有限元分析的膜盤型面優(yōu)化設(shè)計(jì)方法，最后以某膜盤為算例，完成了優(yōu)化設(shè)計(jì)和試驗(yàn)驗(yàn)證，證實(shí)了該方法的可行性和有效性。

2 膜盤型面結(jié)構(gòu)及參數(shù)

膜盤型面通常為變厚薄板，其側(cè)面一般為特殊曲線，結(jié)構(gòu)，如圖1所示。膜盤型面曲線方程［2］為：

圖1 膜盤結(jié)構(gòu)Fig.1 Cross-Section of A Diaphragm

式中：t—任意半徑上型面厚度值；ta—型面上最薄處厚度值；r—型面任意一點(diǎn)半徑值；a—型面最薄處半徑值；n—型面指數(shù)，用于表征膜盤型面厚度變化梯度。

3 膜盤高周疲勞應(yīng)力分析

膜盤聯(lián)軸器在工作過(guò)程中需承受扭矩、離心載荷、軸向偏移及角向偏斜，各載荷作用下膜盤的應(yīng)力分布，如圖2所示。在扭矩、離心載荷和軸向偏移作用下，膜盤型面應(yīng)力周向是均勻分布的，而在角向偏斜作用下，膜盤型面應(yīng)力為軸對(duì)稱分布（非周向均勻），原因在于角向偏斜使膜盤型面產(chǎn)生了彎曲內(nèi)應(yīng)力。

圖2 膜盤應(yīng)力分布Fig.2 Stress Distribution of A Diaphragm

利用ANSYS瞬態(tài)分析方法，將角向偏斜作為瞬態(tài)分析的初始條件，得到在各載荷共同作用下膜盤型面頂部區(qū)域應(yīng)力隨時(shí)間的變化曲線結(jié)果，如圖3所示。根據(jù)圖3，膜盤在扭矩、離心載荷及軸向偏移作用下，應(yīng)力并不隨時(shí)間變化（穩(wěn)態(tài)應(yīng)力），而施加角向偏斜載荷后，應(yīng)力隨時(shí)間周期波動(dòng)（動(dòng)態(tài)應(yīng)力），該應(yīng)力即為膜盤的高周疲勞應(yīng)力。由于聯(lián)軸器的轉(zhuǎn)速大，角向偏斜導(dǎo)致的循環(huán)應(yīng)力變化頻率相對(duì)較高。

圖3 膜盤型面應(yīng)力變化曲線Fig.3 Stress Fluctuation of Diaphragm Profile

根據(jù)S-N曲線，其高周疲勞安全系數(shù)Ks可以表示為：

式中：σa、σm、σb、σN—膜盤型面的循環(huán)應(yīng)力幅值、平均應(yīng)力、膜盤型面材料強(qiáng)度極限和疲勞極限。

根據(jù)圖3，平均應(yīng)力σm為扭矩、離心載荷及周向偏移共同作用下的應(yīng)力值，而應(yīng)力幅值σa為疊加角向偏斜作用后產(chǎn)生的應(yīng)力波動(dòng)值。

4 型面參數(shù)對(duì)疲勞應(yīng)力的影響分析

最小厚度ta和型面指數(shù)n是型面設(shè)計(jì)的兩個(gè)關(guān)鍵參數(shù)，對(duì)型面整體的穩(wěn)態(tài)和動(dòng)態(tài)應(yīng)力分布有重要影響。以某膜盤型面為算例，利用ANSYS軟件計(jì)算不同型膜盤面的高周疲勞安全系數(shù)結(jié)果，如圖4、圖5所示。

圖4 型面指數(shù)n對(duì)高周疲勞安全系數(shù)的影響Fig.4 Relation of Fatigue Safety Factor to Profile Index n

圖5 型面最小厚度ta對(duì)高周疲勞安全系數(shù)的影響Fig.5 Relation of Fatigue Safety Factor to Minimum Thickness ta

根據(jù)圖4，當(dāng)最小厚度ta相同時(shí)，隨型面指數(shù)增大，安全系數(shù)并非單調(diào)遞增或遞減，且在全局范圍內(nèi)存在多處極值。當(dāng)取不同ta值時(shí)，各極值點(diǎn)對(duì)應(yīng)的型面指數(shù)值也不相同。當(dāng)ta=0.4mm時(shí)，最大值安全系數(shù)對(duì)應(yīng)的型面指數(shù)為n=0.55；當(dāng)ta=0.6mm，最大安全系數(shù)對(duì)應(yīng)的型面指數(shù)n=1.1。

根據(jù)圖5，當(dāng)型面指數(shù)n相同時(shí)，疲勞安全系數(shù)KS不隨ta值單調(diào)變化，且不同型面指數(shù)對(duì)應(yīng)的變化趨勢(shì)也不相同，說(shuō)明型面指數(shù)會(huì)影響到ta和KS之間的變化關(guān)系。

因此，最小厚度ta和型面指數(shù)n對(duì)膜盤型面高周疲勞影響不是獨(dú)立的，兩者存在一定的相關(guān)性，設(shè)計(jì)時(shí)需要相互匹配才能實(shí)現(xiàn)膜盤型面高周疲勞強(qiáng)度最優(yōu)化。

5 膜盤型面優(yōu)化設(shè)計(jì)思路

由于型面疲勞安全系數(shù)并不隨型面最小厚度和型面指數(shù)單調(diào)變化，且ta和n兩參數(shù)存在相互關(guān)聯(lián)性，因此需要進(jìn)行優(yōu)化來(lái)實(shí)現(xiàn)參數(shù)的匹配設(shè)計(jì)。膜盤型面優(yōu)化設(shè)計(jì)基本思路為：初步確定膜盤型面結(jié)構(gòu)參數(shù)，利用ANSYS軟件進(jìn)行膜盤聯(lián)軸器參數(shù)化建模，再開展有限元計(jì)算分析，根據(jù)計(jì)算結(jié)果對(duì)型面參數(shù)進(jìn)行優(yōu)化，并重新有限元建模計(jì)算，直至優(yōu)化收斂，最終得到最優(yōu)型面參數(shù)。

5.1 參數(shù)化建模

利用ANSYS 自帶的APDL 語(yǔ)言對(duì)膜盤進(jìn)行參數(shù)化建模。在確定膜盤型面參數(shù)后，建立膜盤型面方程，并在徑向均勻選取多節(jié)點(diǎn)，采用B樣條進(jìn)行擬合，得到型面曲線。對(duì)于膜盤聯(lián)軸器的空心軸，采用beam 單元模擬。膜盤參數(shù)化有限元模型，如圖6所示。

圖6 膜盤和空心軸有限元模型Fig.6 FEM Model of a Diaphragm and a Hollow Shaft

5.2 有限元分析

根據(jù)膜盤聯(lián)軸器工況特點(diǎn)，利用ANSYS多載荷步計(jì)算，得到膜盤型面在不同載荷作用下的應(yīng)力結(jié)果以及聯(lián)軸器整體的臨界轉(zhuǎn)速。對(duì)高周疲勞應(yīng)力，需根據(jù)應(yīng)力特點(diǎn)，對(duì)各單元在穩(wěn)態(tài)工況和角向偏斜工況下的應(yīng)力進(jìn)行合成，確定型面各單元的應(yīng)力均值和幅值，并篩選出其中的“危險(xiǎn)單元”，得到其疲勞安全系數(shù)。

5.3 型面參數(shù)優(yōu)化

5.3.1 優(yōu)化設(shè)計(jì)變量

根據(jù)式1和圖1，膜盤型面結(jié)構(gòu)的主要參數(shù)包括：型面最小厚度ta，膜盤型面最小厚度半徑值a，型面指數(shù)n，型面根部半徑b。其中，a和b可根據(jù)接口尺寸和安裝空間確定，因此型面優(yōu)化的主要設(shè)計(jì)變量為：

（1）型面最小厚度ta，可行域一般?。?.2，1.4］［10］；

（2）型面指數(shù)n，可行域?yàn)椋?，2］。

5.3.2 目標(biāo)函數(shù)

為提高膜盤抗疲勞性能，將型面高周疲勞強(qiáng)度安全系數(shù)KS最大化作為優(yōu)化設(shè)計(jì)目標(biāo)。當(dāng)膜盤材料及載荷條件確定的情況下，平均應(yīng)力和應(yīng)力幅值為型面參數(shù)的函數(shù)，即：

5.3.3 約束條件

優(yōu)化設(shè)計(jì)約束條件包括：靜強(qiáng)度條件和臨界轉(zhuǎn)速條件。對(duì)于靜強(qiáng)度，膜盤在最大瞬態(tài)載荷狀態(tài)下的安全系數(shù)K1應(yīng)不小于1.15，即：

式中：σmax—膜盤型面最大應(yīng)力；

σs—膜盤材料的屈服極限。

為避免聯(lián)軸器在工作過(guò)程中產(chǎn)生共振破壞，膜盤型面須具備一定的剛度確保聯(lián)軸器在工作時(shí)避開臨界轉(zhuǎn)速，通常第一階臨界轉(zhuǎn)速安全裕度K2通常應(yīng)不低于20%，即：

式中：n1—聯(lián)軸器的第一階臨界轉(zhuǎn)速；

nS—聯(lián)軸器工作轉(zhuǎn)速。

5.3.4 優(yōu)化工具

利用ANSYS軟件中的優(yōu)化方法進(jìn)行優(yōu)化。該軟件提供了兩種優(yōu)化設(shè)計(jì)方法：零階方法和一階方法，其中零階方法是一個(gè)較均衡的優(yōu)化方法，迭代速度塊，適用于大多數(shù)工程優(yōu)化問(wèn)題，一階方法是基于目標(biāo)函數(shù)對(duì)設(shè)計(jì)變量的偏導(dǎo)數(shù)進(jìn)行優(yōu)化，優(yōu)化精度高，但是收斂速度慢，且可能收斂于局部最優(yōu)解。為提高優(yōu)化設(shè)計(jì)的速度和適應(yīng)性，采用零階方法進(jìn)行優(yōu)化。

6 優(yōu)化算例及試驗(yàn)驗(yàn)證

6.1 優(yōu)化算例

某傳動(dòng)軸為帶膜盤結(jié)構(gòu)，其一側(cè)為兩個(gè)薄壁膜盤結(jié)構(gòu)，如圖7所示。傳動(dòng)軸傳遞的穩(wěn)態(tài)扭矩為460N·m，瞬態(tài)扭矩為550N·m，工作轉(zhuǎn)速為19800rpm，穩(wěn)態(tài)軸向偏角1°，瞬態(tài)軸向偏角1.5°，穩(wěn)態(tài)軸向偏移0.5mm，瞬態(tài)軸向偏移1mm。

臨床分期主要與淋巴瘤侵及的范圍有關(guān)，本組患者中Ⅰ～Ⅱ期淋巴瘤共52例（50.5%），SUVmax為8.32±4.58；Ⅲ～Ⅳ期51例（49.5%），SUVmax為10.45±6.43，二者間差異無(wú)統(tǒng)計(jì)學(xué)意義（P＞0.05）。

圖7 動(dòng)力傳動(dòng)軸結(jié)構(gòu)簡(jiǎn)圖Fig.7 Cross-Section of a Transmission Shaft

6.2 優(yōu)化結(jié)果及討論

借助ANSYS優(yōu)化工具，開展膜盤型面參數(shù)優(yōu)化，共得到15個(gè)滿足設(shè)計(jì)要求的膜盤型面（可行解15組），其中一個(gè)組為最優(yōu)解（ta=0.46，n=1.2）。優(yōu)化設(shè)計(jì)得到的膜盤和原膜盤（優(yōu)化前）的型面參數(shù)和計(jì)算結(jié)果，如表1所示。

表1 優(yōu)化前后膜盤型面參數(shù)及結(jié)果Tab.1 Parameters and Results before and after Optimization

相比原膜盤，優(yōu)化設(shè)計(jì)膜盤的靜強(qiáng)度安全系數(shù)提高了7%，臨界轉(zhuǎn)速安全裕度基本保持不變，高周疲勞強(qiáng)度安全系數(shù)提高了27%。

根據(jù)優(yōu)化設(shè)計(jì)中的有限元結(jié)果文件，提取膜盤型面危險(xiǎn)區(qū)域（動(dòng)靜復(fù)合應(yīng)力最大單元）的應(yīng)力數(shù)據(jù)，和原膜盤型面有限元計(jì)算結(jié)果進(jìn)行對(duì)比結(jié)果，如表2所示。

表2 優(yōu)化前后膜盤型面疲勞應(yīng)力Tab.2 Fatigue Stresses Before and After Optimization

應(yīng)力結(jié)果顯示，原膜盤的疲勞危險(xiǎn)區(qū)域位于型面頂部（應(yīng)力均值205MPa，應(yīng)力幅值179MPa），優(yōu)化設(shè)計(jì)膜盤的危險(xiǎn)區(qū)域位于型面根部，其應(yīng)力均值為196MPa，和原膜盤型面頂部應(yīng)力基本相當(dāng)，但是應(yīng)力幅值減小到115MPa，相比原膜盤型面減小了35.7%，有較大改善。

6.3 試驗(yàn)驗(yàn)證

為了驗(yàn)證優(yōu)化設(shè)計(jì)方法的有效性，加工了四件膜盤試驗(yàn)件（原膜盤和優(yōu)化設(shè)計(jì)膜盤試驗(yàn)件各兩件），開展膜盤聯(lián)軸器疲勞運(yùn)轉(zhuǎn)試驗(yàn)。試驗(yàn)時(shí)，除了施加穩(wěn)態(tài)扭矩外，通過(guò)設(shè)置拉桿來(lái)施加軸向偏移，并通過(guò)調(diào)整兩端安裝座的徑向偏移來(lái)施加角向偏斜。試驗(yàn)方案，如圖8所示。

圖8 膜盤試驗(yàn)件運(yùn)轉(zhuǎn)試驗(yàn)方案Fig.8 Schematic of Coupling’s Fatigue Testing Bench

優(yōu)化前兩件膜盤分別在運(yùn)轉(zhuǎn)613h 和659h 后，型面表面出現(xiàn)了裂紋，而優(yōu)化后的兩件膜盤在運(yùn)行750h 后，型面表面經(jīng)檢查無(wú)異常。試驗(yàn)結(jié)果表明：優(yōu)化設(shè)計(jì)膜盤的疲勞性能優(yōu)于原膜盤。

7 結(jié)論

（1）膜盤型面最小厚度ta和型面指數(shù)n對(duì)膜盤高周疲勞應(yīng)力有重要影響，且ta和n存在相互關(guān)聯(lián)，型面設(shè)計(jì)時(shí)需要進(jìn)行匹配優(yōu)化以取得最佳高周疲勞安全系數(shù)；

（2）建立以高周疲勞安全系數(shù)為目標(biāo)的型面優(yōu)化設(shè)計(jì)方法，并以某膜盤為例，通過(guò)試驗(yàn)驗(yàn)證了該方法的可行性和有效性。