飛機(jī)維修專業(yè)畢業(yè)設(shè)計(jì)

1、飛機(jī)維修（發(fā)動(dòng)機(jī)維修）專業(yè)畢業(yè)設(shè)計(jì) 燃?xì)廨啓C(jī)動(dòng)力渦輪組件的強(qiáng)度分析 摘要：燃?xì)廨啓C(jī)動(dòng)力渦輪組件是燃?xì)廨啓C(jī)的主要組件之 一, 由于其不僅工作溫度高，而且還要承受燃?xì)廨啓C(jī)在起動(dòng) 和停機(jī)時(shí)，因溫度劇烈變化引起的熱沖擊，工作條件惡劣， 故渦輪組件是決定燃?xì)廨啓C(jī)壽命的關(guān)鍵部件，為確保有足夠 的壽命。本文重點(diǎn)對(duì)某型燃?xì)廨啓C(jī)低壓渦輪壓氣機(jī)轉(zhuǎn)子葉片 的強(qiáng)度特性進(jìn)行了分析，并驗(yàn)證了其穩(wěn)定性及可靠性。本文 使用 SAmcEF/Field 軟件的轉(zhuǎn)子動(dòng)力學(xué)分析模塊對(duì)該轉(zhuǎn)子葉 片進(jìn)行了分析計(jì)算。根據(jù)機(jī)組實(shí)際運(yùn)行的條件，計(jì)算了該機(jī) 組轉(zhuǎn)子具有較小的臨界轉(zhuǎn)速、穩(wěn)態(tài)不平衡響應(yīng)、葉片丟失瞬 態(tài)響應(yīng)等。計(jì)算結(jié)果表明臨界轉(zhuǎn)速安全

2、系數(shù)合理；轉(zhuǎn)子系 統(tǒng)選取的平衡量具有較小的振動(dòng)幅值；轉(zhuǎn)子的瞬態(tài)響應(yīng)結(jié)果 驗(yàn)證了結(jié)構(gòu)方案的合理性，轉(zhuǎn)子系統(tǒng)具有較好的穩(wěn)定性。得 出了此轉(zhuǎn)子結(jié)構(gòu)方案能保證低壓渦輪壓氣機(jī)穩(wěn)定運(yùn)行的結(jié) 論，從而為燃?xì)廨啓C(jī)的可靠性、 維修性、保障性提供了參考。關(guān)鍵詞：低壓渦輪；強(qiáng)度；轉(zhuǎn)速；振動(dòng)；可靠性 、八 、-前言 上世紀(jì)二十年代，德國(guó)人霍爾茨瓦特制成第一臺(tái)實(shí)用的 燃?xì)廨啓C(jī)，其效率為 13%、功率為 370 千瓦，自此之后，燃 氣輪機(jī)就逐漸進(jìn)入了人們的生活領(lǐng)域。燃?xì)廨啓C(jī)、汽輪機(jī)、 發(fā)電機(jī)及電動(dòng)機(jī)等都是典型的旋轉(zhuǎn)機(jī)械，都以轉(zhuǎn)子作為工作主體。轉(zhuǎn)子連同軸承、支座等稱為轉(zhuǎn)子 支承系統(tǒng)。機(jī)組運(yùn)轉(zhuǎn)時(shí)轉(zhuǎn)子系統(tǒng)常常發(fā)生振動(dòng)，而振動(dòng)產(chǎn)

3、生噪聲。降低了 機(jī)組的工作效率；嚴(yán)重的振動(dòng)會(huì)導(dǎo)致轉(zhuǎn)子斷裂，造成重大危 害。轉(zhuǎn)子 支承系統(tǒng)的振動(dòng)是多樣的，包括轉(zhuǎn)軸的扭曲振動(dòng)、彎曲振動(dòng)和輪盤(pán)葉片的振動(dòng)等。轉(zhuǎn)軸的振動(dòng)較為復(fù)雜， 牽涉的因素較多。轉(zhuǎn)子動(dòng)力學(xué)就是以轉(zhuǎn)軸的彎曲振動(dòng)作為主 要研究對(duì)象的。轉(zhuǎn)子的不平衡量所引起的振動(dòng)屬于“強(qiáng)迫振動(dòng)” ，它的 角頻率和轉(zhuǎn)動(dòng)角速度相等。對(duì)于高速轉(zhuǎn)子，除了不平衡質(zhì)量 引起的振動(dòng)以外，還有頻率與轉(zhuǎn)動(dòng)角速度不相等的振動(dòng)，稱 為“渦動(dòng)”。轉(zhuǎn)子連接件配合面的摩擦、轉(zhuǎn)軸的內(nèi)阻、軸承 油膜力或密封引起的氣動(dòng)力等都是產(chǎn)生渦動(dòng)的因素。理論上， 轉(zhuǎn)軸的這種渦動(dòng)屬于“自激振動(dòng)” 1 。渦動(dòng)甚至可以破壞 轉(zhuǎn)軸或軸承。本文研究的對(duì)象是該型

4、燃?xì)廨啓C(jī)低壓渦輪壓氣機(jī)轉(zhuǎn)子葉片，通過(guò)對(duì)轉(zhuǎn)子-支承建模，使用 SAmcEF專業(yè)轉(zhuǎn)子動(dòng) 力學(xué)分析軟件，采用有限元素法分析了其轉(zhuǎn)子動(dòng)力學(xué)特性， 包括轉(zhuǎn)子的臨界轉(zhuǎn)速計(jì)算、穩(wěn)態(tài)不平衡響應(yīng)分析、轉(zhuǎn)子穩(wěn)定 性分析等。驗(yàn)證了其在工程應(yīng)用方面的可用性及可靠性。同 時(shí)得出了分析其動(dòng)力學(xué)特性的基本方法及結(jié)論。2 結(jié)構(gòu)簡(jiǎn)介該型燃?xì)廨啓C(jī)低壓渦輪壓氣機(jī)轉(zhuǎn)子呈軸流輪轂式整體結(jié)構(gòu)，它由 9 級(jí)低壓壓氣機(jī)和 l 級(jí)低壓渦輪組成，低壓壓氣 機(jī)與低壓渦輪之問(wèn)通過(guò)低壓渦輪軸連接并以花鍵傳遞扭矩（見(jiàn)圖1）。08級(jí)低壓壓氣機(jī)輪盤(pán)、葉片材料為鈦合金； 低壓渦輪盤(pán)、葉片材料為高溫合金；低壓渦輪軸材料為馬氏 體不銹鋼。低壓渦輪壓氣機(jī)轉(zhuǎn)子呈 3

5、 點(diǎn)支承結(jié)構(gòu)。前支承采 用徑向止推滾珠軸承，支承點(diǎn)位于壓氣機(jī) o 級(jí)輪盤(pán)的前段， 由彈性支承、擠壓油膜阻尼器、滾珠軸承組成；中間支承采 用滾柱軸承由擠壓油膜阻尼器、滾珠軸承組成，支承點(diǎn)位 于低壓壓氣機(jī)后軸徑后段；后支承采用滾柱軸承，由彈性支 承、擠壓油膜阻尼器、滾珠軸承組成，支承點(diǎn)位于低壓渦輪 軸后段。3 計(jì)算模型該型燃機(jī)低壓渦輪壓氣機(jī)轉(zhuǎn)子的 3D 模型較大，這將導(dǎo) 致在網(wǎng)格劃分以及計(jì)算過(guò)程中花費(fèi)大量時(shí)間，因此我們對(duì)本 機(jī)組的計(jì)算采用 2D 軸對(duì)稱模型。在總體直角坐標(biāo)系下建立 二維軸對(duì)稱單元，其種類有 3 節(jié)點(diǎn)或高階 6 節(jié)點(diǎn)的三角形單 元、 4 節(jié)點(diǎn)或高階 8 節(jié)點(diǎn)的四邊形單元。每個(gè)節(jié)點(diǎn)有

6、9 個(gè)自 由度，前 6 個(gè)自由度與梁?jiǎn)卧粯樱謩e為沿旋轉(zhuǎn)軸線方向 的拉伸和扭轉(zhuǎn)，以及由彎曲而引起的其他 2 個(gè)方向的線位移 和角位移。另外， 3 個(gè)自由度與旋轉(zhuǎn)軸的橫截面變形有關(guān)， 分別為拉伸引起的徑向位移和彎曲引起的 2 個(gè)切向位移。同 時(shí)使用這類單元可以很好地模擬轉(zhuǎn)子的“渦動(dòng)效應(yīng)”建立總體直角坐標(biāo)系（ X,y,Z ）和局部圓柱坐標(biāo)系。設(shè) 定為局部坐標(biāo)系下的位移分量 , 將它們沿 0 向按傅立葉級(jí)數(shù) 展開(kāi)如下：式中： n 為波數(shù) ;m 為相數(shù)。每個(gè)諧函數(shù)對(duì)應(yīng) 1 個(gè)結(jié)構(gòu)的 基本變形 2 。低壓渦輪壓氣機(jī)轉(zhuǎn)子的壓氣機(jī)葉盤(pán)間的連接方式為焊 接或螺栓緊固渦輪葉盤(pán)與渦輪軸之間也是以螺栓緊固的方 式

7、連接，低壓壓氣機(jī)與低壓渦輪軸之間的連接為花鍵。以此 為依據(jù)進(jìn)行整體建模。低壓渦輪壓氣機(jī)轉(zhuǎn)子葉片較長(zhǎng)，其質(zhì) 量及轉(zhuǎn)動(dòng)慣量都很大，會(huì)產(chǎn)生較大的離心力和回轉(zhuǎn)力矩，這 種回轉(zhuǎn)效應(yīng)會(huì)對(duì)轉(zhuǎn)子的臨界轉(zhuǎn)速產(chǎn)生比較大的影響葉片等 效簡(jiǎn)化為集中質(zhì)量加到計(jì)算模型中。4 計(jì)算內(nèi)容4.1 支承剛度和阻尼計(jì)算 現(xiàn)代航空發(fā)動(dòng)機(jī)高速轉(zhuǎn)子系統(tǒng)大都采用“柔軸”設(shè)計(jì)， 使發(fā)動(dòng)機(jī)工作轉(zhuǎn)速高于轉(zhuǎn)子系統(tǒng)的臨界轉(zhuǎn)速，該型燃機(jī)低壓 渦輪壓氣機(jī)轉(zhuǎn)子就采用了這種設(shè)計(jì)理念。由轉(zhuǎn)子在臨界轉(zhuǎn)速 時(shí)的撓度和傳遞率可知，此類轉(zhuǎn)子工作在臨界轉(zhuǎn)速或者通過(guò) 臨界轉(zhuǎn)速時(shí)，其振幅正比于不平衡量和轉(zhuǎn)子的臨界角速度， 反比于阻尼系數(shù) 3 。顯然，對(duì)于“柔軸”類轉(zhuǎn)子系統(tǒng)

8、采取 適當(dāng)措施減小轉(zhuǎn)子的不平衡量、降低臨界轉(zhuǎn)速和增大阻尼都 是可以達(dá)到減小轉(zhuǎn)子振動(dòng)幅值的目的。而在發(fā)動(dòng)機(jī)總體結(jié)構(gòu)設(shè)計(jì)已經(jīng)確定的前提下，采用彈性支承或者阻尼器，或者兩 者兼而有之，就成為了最為有效和可行的方法。該型燃機(jī)低 壓渦輪壓氣機(jī)轉(zhuǎn)子就采用了這種減振方式。機(jī)組所采用的彈 性支承就是所謂低剛性彈性支承，可以使轉(zhuǎn)子在通過(guò)臨界轉(zhuǎn) 速時(shí)的振幅和通過(guò)軸承座的外傳載荷均很小，這就可使轉(zhuǎn)子 工作在低階臨界轉(zhuǎn)速以上。因此確定該型燃機(jī)低壓壓氣機(jī)轉(zhuǎn) 子支承剛度是計(jì)算轉(zhuǎn)子 支承系統(tǒng)臨界轉(zhuǎn)速所必需的前提條件。4.1.1 前支承剛度和阻尼計(jì)算 低壓壓氣機(jī)轉(zhuǎn)子前支承位于 0 級(jí)輪盤(pán)前段，由鼠籠式彈 性支撐、擠壓油膜阻尼

9、器、滾珠軸承組成。）滾珠軸承的剛度 4式中：klxx為徑向剛度；Fr為徑向外力;n為滾珠數(shù)量; d為滾珠直徑；B為接觸角。2）鼠籠川度 以有限元法計(jì)算了鼠籠的靜剛度。3）擠壓油膜阻尼器的剛度及阻尼計(jì)算轉(zhuǎn)子 支承系統(tǒng)采用擠壓油膜阻尼器后，便會(huì)出現(xiàn)油膜剛度和油膜阻尼，在分析其動(dòng)力特性時(shí)，除轉(zhuǎn)子一支承系 統(tǒng)自身的剛度外，還要計(jì)及油膜剛度。油膜阻尼遠(yuǎn)遠(yuǎn)大于轉(zhuǎn) 子一支承系統(tǒng)其他阻尼，存在油膜阻尼時(shí)，其他阻尼甚至可 以忽略不計(jì)。因此對(duì)于前支承來(lái)說(shuō)，油膜阻尼就是其支承系統(tǒng)的阻尼。依據(jù)其結(jié)構(gòu)分析，前支承的阻尼器應(yīng)按短軸承 半油膜計(jì)算。4）前支承總剛度、阻尼N/m；N*s/m 。4.1.2 中支承剛度和阻尼計(jì)算

10、中間支承位于壓氣機(jī)后軸頸后段，由擠壓油膜阻尼器、 滾珠軸承組成。）滾珠軸承的剛度計(jì)算。式中：klxx為徑向剛度；Fr為徑向外力;n為滾珠數(shù)量; L為滾子有效長(zhǎng)度；B為接觸角。2）中間軸承總剛度及總阻尼計(jì)算N/m；N*s/m 。4.1.3 后支承剛度和阻尼計(jì)算 后支承位于低壓渦輪軸后段，由彈性支承、擠壓油膜阻 尼器、滾珠軸承組成。后支承總剛度及總阻尼計(jì)算如下：N/m；N*s/m 。4.2 低壓渦輪壓氣機(jī)轉(zhuǎn)子臨界轉(zhuǎn)速計(jì)算低壓渦輪壓氣機(jī)轉(zhuǎn)子是由低壓壓氣機(jī)轉(zhuǎn)子、低壓渦輪轉(zhuǎn) 子組成的軸系掌握軸系廠卜單獨(dú)轉(zhuǎn)子的動(dòng)力學(xué)特性對(duì)分析 整體軸系的動(dòng)力學(xué)特性是十分必要的。由于低壓渦輪前端以 花鍵與低壓壓氣機(jī)后軸徑相

11、連接，并依靠低壓后軸徑上的中 問(wèn)軸承作為其前支承點(diǎn)，即低壓渦輪并沒(méi)有獨(dú)立的前支承， 也就沒(méi)有單獨(dú)計(jì)算低壓渦輪轉(zhuǎn)子臨界轉(zhuǎn)速的必要。因此在計(jì) 算整體軸系的臨界轉(zhuǎn)速之前首先對(duì)低壓壓氣機(jī)轉(zhuǎn)子進(jìn)行了 臨界轉(zhuǎn)速計(jì)算分析。4.2.1 材料屬性低壓渦輪壓氣機(jī)轉(zhuǎn)子是由幾種不同的主體材料結(jié)構(gòu)組 成，對(duì) 2 種不同材料的屬性分別進(jìn)行定義。低壓壓氣機(jī)輪盤(pán) 材料為鈦合金、低壓壓氣機(jī)后軸徑為馬氏體不銹鋼、低壓渦 輪軸為馬氏體不銹鋼、低壓渦輪為高溫合金，分別計(jì)算低壓 壓氣機(jī)葉片與低壓渦輪轉(zhuǎn)子葉片的等效質(zhì)量和等效轉(zhuǎn)動(dòng)慣 量后以集中質(zhì)量的方式對(duì)其加載。圖2施加各項(xiàng)邊界條件后的轉(zhuǎn)子4.2.2 計(jì)算結(jié)果及分析）計(jì)算結(jié)果低壓渦輪壓氣

12、機(jī)轉(zhuǎn)子臨界轉(zhuǎn)速計(jì)算結(jié)果如表1 所示。1-3階臨界轉(zhuǎn)速的振型如（圖 3-圖 5）所示。2）結(jié)果分析由于低壓渦輪壓氣機(jī)的慢車轉(zhuǎn)速N1=2400r/min，工作轉(zhuǎn)速為N。=7700r / min，根據(jù)本文的計(jì)算結(jié)果，可知 臨界轉(zhuǎn)速的安全系數(shù)：k1=Np1/N1=36802400=l.53 k2=Np2/N0=123307700=1.60 超速時(shí)的臨界轉(zhuǎn)速的安全系數(shù)： k2'=Np2/=12330/=1.394.3 穩(wěn)態(tài)諧波響應(yīng)分析穩(wěn)態(tài)不平衡響應(yīng)的計(jì)算是轉(zhuǎn)子動(dòng)力學(xué)分析中與臨界轉(zhuǎn) 速計(jì)算同等重要的基木任務(wù)。穩(wěn)態(tài)不平衡響應(yīng)分析也可以用 來(lái)確定系統(tǒng)的臨界轉(zhuǎn)速，但它更重要的任務(wù)是不是用來(lái)求解 在轉(zhuǎn)子系

13、統(tǒng)中可能存在的不平衡量作用下，轉(zhuǎn)子一支承系統(tǒng) 的穩(wěn)態(tài)不平衡響應(yīng)分析研究如何采取措施。限制最大不平 衡響應(yīng)及減小不平衡響應(yīng)。在任一轉(zhuǎn)速下的不平衡響應(yīng)，可能是若干個(gè)主振型的線 性疊加，隨著轉(zhuǎn)速的改變可由某一階振型向另一階振型過(guò)渡。 小平衡量的分布規(guī)律的改變可引起穩(wěn)態(tài)不平衡響應(yīng)的變化， 不僅僅是量的變化。也會(huì)引起相應(yīng)“振型”的變化。因此， 穩(wěn)態(tài)不平衡響應(yīng)計(jì)算比臨界轉(zhuǎn)速計(jì)算能提供更多的有用信 息。轉(zhuǎn)子的激振力就是轉(zhuǎn)子的轉(zhuǎn)速，振動(dòng)幅值與轉(zhuǎn)速的平方 成比例。分析該轉(zhuǎn)子的不平衡響應(yīng)對(duì)于掌握其在不平衡量作 用下的動(dòng)態(tài)響應(yīng)特性是非常重要的。4.3.1 計(jì)算模型按 G1 級(jí)精度確定低壓渦輪轉(zhuǎn)子的單位質(zhì)量允許的剩余

14、 不平衡量。加重位置位于第 o 級(jí)壓氣機(jī)輪盤(pán)及渦輪盤(pán)上第 0 級(jí)壓氣機(jī)輪盤(pán)上的不平衡量為 2.54g, 渦輪盤(pán)上的不平衡量 為 l.37g 。4.3.2 計(jì)算結(jié)果 在不平衡響應(yīng)計(jì)算結(jié)果中，可以得到轉(zhuǎn)子在計(jì)算頻域范 圍內(nèi)的最大響應(yīng)值，同時(shí)也可以得到工作轉(zhuǎn)速范圍內(nèi)的最大 響應(yīng)值。最大響應(yīng)值包括最大相對(duì)變形以及最大相對(duì)旋轉(zhuǎn)角 度。4.4瞬態(tài)響應(yīng)分析轉(zhuǎn)子系統(tǒng)的瞬態(tài)響應(yīng)分析主要是轉(zhuǎn)子系統(tǒng)不平衡突然 變化，作用在轉(zhuǎn)子系統(tǒng)上的外載荷突然變化或轉(zhuǎn)子系統(tǒng)在變 轉(zhuǎn)速下工作等情況下轉(zhuǎn)子系統(tǒng)的響應(yīng)分析包括轉(zhuǎn)子系統(tǒng) 的位移、變形以及支承結(jié)構(gòu)的傳遞載荷分布。對(duì)于帶擠壓油 膜阻尼器的轉(zhuǎn)子系統(tǒng)，在機(jī)動(dòng)載荷作用下的響應(yīng)需進(jìn)行瞬

15、態(tài) 分析。飛機(jī)機(jī)動(dòng)飛行時(shí)引起的機(jī)動(dòng)載荷為單位階躍函數(shù)。轉(zhuǎn)子變速過(guò)程通常近似為等加速或等減速過(guò)程，此時(shí)有：式中：& amp;為轉(zhuǎn)子轉(zhuǎn)角；w0為初始角速度。對(duì)運(yùn)動(dòng)方程式進(jìn)行直接積分可以求得系統(tǒng)在變速時(shí)的 不平衡響應(yīng)。這時(shí)。式中油膜耦合力應(yīng)采用瞬態(tài)力表達(dá)式，不平衡離心力的表達(dá)式改為：（對(duì)于平面）（對(duì)于平面）對(duì)于離心力F,應(yīng)根據(jù)不同的不平衡量的相對(duì)位置來(lái)確 定5 。在瞬態(tài)響應(yīng)計(jì)算結(jié)果中，可以得到轉(zhuǎn)子在計(jì)算時(shí)域范圍 內(nèi)支承點(diǎn)以及葉片丟失點(diǎn)的最大響應(yīng)值，其中包括位移、速 度、加速度等。也可以得到軸心軌跡等數(shù)據(jù)；還可以從計(jì)算 結(jié)果中得到轉(zhuǎn)子穩(wěn)定性的結(jié)論。同時(shí)也可以得到計(jì)算時(shí)域內(nèi) 的響應(yīng)變化曲線等。5

16、 結(jié)論） 臨界轉(zhuǎn)速計(jì)算結(jié)論 本文計(jì)算的該型燃機(jī)低壓渦輪壓氣機(jī)轉(zhuǎn)子的慢車轉(zhuǎn)速 的臨界轉(zhuǎn)速安全系數(shù)為 1.53 ，全工況工作轉(zhuǎn)速的臨界轉(zhuǎn)速安 全系數(shù)為 1.60 ，超速時(shí)的臨界轉(zhuǎn)速安全系數(shù)為 1.39 ?？梢?jiàn)， 此結(jié)構(gòu)方案能滿足低壓渦輪壓氣機(jī)轉(zhuǎn)子穩(wěn)定運(yùn)行的要求且 有較高裕度。2） 不平衡響應(yīng)計(jì)算結(jié)論 在工作頻率范圍內(nèi)，轉(zhuǎn)子的動(dòng)態(tài)響應(yīng)較小，最大變形為 o 055mm（ 62.5Hz） 。證明改進(jìn)技術(shù)研究達(dá)到預(yù)期效果，精 度得到很大程度的提高。因此，低壓渦輪壓氣機(jī)轉(zhuǎn)子選取的 平衡精度是合理的，同時(shí)具有較小的振動(dòng)幅值。3） 瞬態(tài)響應(yīng)計(jì)算結(jié)論分析了轉(zhuǎn)子由 0 轉(zhuǎn)速升至工作轉(zhuǎn)速過(guò)程中，第 0 級(jí)葉片 丟失狀

17、態(tài)下轉(zhuǎn)子的瞬態(tài)響應(yīng)。由轉(zhuǎn)了的軸心軌跡判斷。轉(zhuǎn)子 系統(tǒng)是穩(wěn)定的，己形成了穩(wěn)態(tài)的圓軌跡；轉(zhuǎn)子的振動(dòng)幅值響 應(yīng)是在系統(tǒng)阻尼間隙范圍內(nèi)。經(jīng)以上分析證明，該轉(zhuǎn)子系統(tǒng)的設(shè)計(jì)是合理的，此轉(zhuǎn)子 結(jié)構(gòu)方案能保證低壓渦輪壓氣機(jī)的穩(wěn)定運(yùn)行。參考文獻(xiàn)1 顧XX,等.轉(zhuǎn)子動(dòng)力學(xué)m.北京，國(guó)防工業(yè)出版社， 1985.2 梅 XX .二維有限元模型在燃?xì)廨啓C(jī)轉(zhuǎn)子 支承系統(tǒng)臨界轉(zhuǎn)速計(jì)算中的應(yīng)用j. 燃?xì)鉁u輪實(shí)驗(yàn)與研究，XX,16（ 3）：42-44.3 鐘XX.轉(zhuǎn)子動(dòng)力學(xué)m.北京：清華大學(xué)出版社，1987.4 航空發(fā)動(dòng)機(jī)設(shè)計(jì)手冊(cè)總編委會(huì)，航空發(fā)動(dòng)機(jī)設(shè)計(jì)手 冊(cè)，（第 19 冊(cè)） m. 北京：航空工業(yè)出版社， XX.5 聞XX.高等轉(zhuǎn)子動(dòng)力學(xué)m.北京：機(jī)械工業(yè)出版社， 1999.致謝本論文是在丁 XX 老師的悉心指導(dǎo)下完成的。丁老師和藹可親