張 寧，王耀奇，侯紅亮，張艷苓，董曉萌，李志強

(1 北京航空制造工程研究所，北京 100024；2 塑性成形技術航空科技重點實驗室，北京 100024；3 數(shù)字化塑性成形技術及裝備北京市重點實驗室，北京 100024；4 北京科技大學，北京 100083)

7B04鋁合金超塑性變形行為

張 寧1,2,3，王耀奇1,2,3，侯紅亮1,2,3，張艷苓1,2,3，董曉萌4，李志強1

(1 北京航空制造工程研究所，北京 100024；2 塑性成形技術航空科技重點實驗室，北京 100024；3 數(shù)字化塑性成形技術及裝備北京市重點實驗室，北京 100024；4 北京科技大學，北京 100083)

針對7B04鋁合金開展了變形溫度為470～530℃，應變速率為0.0003～0.01s-1的高溫超塑性拉伸實驗，研究了材料的超塑性變形行為和變形機制。結果表明，7B04鋁合金的流動應力隨著變形溫度的升高和應變速率的降低而逐漸減小，伸長率隨之增加；在變形溫度為530℃，應變速率為0.0003s-1時，7B04鋁合金的伸長率達到最大1105%，超塑性能最佳；應變速率敏感性指數(shù)m值均大于0.3，且隨變形溫度的升高而增加；在500～530℃的變形溫度范圍內，m值大于0.5，表明7B04鋁合金超塑性變形以晶界滑動為主要變形機制；變形激活能Q為190kJ/mol，表明7B04鋁合金的超塑性變形主要受晶內擴散控制；7B04鋁合金超塑性變形中在晶界附近有液相產生，且適量的液相有利于提高材料的超塑性能。

7B04鋁合金；超塑性；變形行為；變形機制

7XXX系(Al-Zn-Mg-Cu系)合金屬于超高強變形鋁合金，具有比重小、強度高、熱加工性能好以及良好的斷裂韌度和耐腐蝕性等優(yōu)點，被廣泛應用于航空航天等工業(yè)領域，尤其在航空工業(yè)中占有非常重要的地位，是飛機的主體結構材料之一[1-3]。目前，A380飛機、空客A320、波音B777客機等世界先進飛機上的許多機體結構(如機翼蒙皮、機翼梁、機身隔框、主框架附件等)已普遍采用7050，7010，7475，7055等高純高強高韌鋁合金材料[3-5]。

7B04鋁合金(美國牌號7475鋁合金)是我國在7A04鋁合金基礎上經過改進而成的一種高純高強鋁合金，屬于Al-Zn-Mg-Cu系，主要應用于航空航天領域[6,7]。相關研究表明，7475鋁合金具有比其他7XXX系合金更加優(yōu)良的超塑性能[8]，各國學者對7475鋁合金的超塑性進行了大量的研究[8-13]。在美國利用7475合金的超塑性已成功制成T-39前機身隔柜、B-lB轟炸機襟翼翼肋和機身構件；洛克韋爾國際公司則用7475合金成形Apache直升機艙板；諾斯洛普用7475改型合金超塑成形F-5E 型戰(zhàn)斗機前部控制系統(tǒng)的下艙板[13]。然而目前針對國產7B04，我國學者主要開展了材料的制備、熱處理和蠕變時效工藝研究[14-19]，對于該材料的超塑成形工藝研究較少，還不能采用超塑成形工藝制造7B04鋁合金復雜構件，嚴重制約了該合金在我國飛機制造業(yè)中的應用。

本工作通過高溫超塑性拉伸實驗對7B04鋁合金開展超塑性變形行為研究，分析變形溫度和應變速率對材料超塑性能的影響規(guī)律，獲得材料的最佳超塑性成形工藝參數(shù)；根據(jù)實驗結果，利用理論計算求解應變速率敏感性指數(shù)m和變形激活能Q值，同時結合斷口形貌，研究7B04鋁合金的超塑性變形機制。

1 實驗材料與方法

實驗所采用的材料是由中南大學提供的7B04鋁合金細晶板材，板材厚度為2.0mm，合金主要化學成分如表1所示。7B04鋁合金細晶板材的制備工藝：首先將25mm厚的熱軋板經460～485℃固溶4h后水淬，然后在300～400℃，8～48h條件下進行時效處理，經300℃預熱1～2h后軋制，且經多道次軋制得到2.0mm厚板，最后在480℃，30min條件下再結晶退火，從而得到平均晶粒尺寸為10μm的細晶板材，如圖1所示。

表1 7B04鋁合金主要化學成分(質量分數(shù)/%)

7B04鋁合金的高溫超塑性拉伸實驗在LETRYDL-20T型電子萬能拉伸試驗機上進行，實驗變形溫度為470，485，500，515，530℃，應變速率為0.0003，0.001，0.003，0.01s-1，拉伸過程中保持夾頭的速度恒定不變。拉伸試樣尺寸示意圖如圖2所示，拉斷后快速水淬，以保留其高溫變形組織，并采用掃描電鏡(SEM)對拉斷后試樣的斷口形貌進行觀察。

2 結果與分析

2.1 超塑性變形行為

2.1.1 真應力-真應變曲線

圖1 7B04鋁合金板材微觀組織Fig.1 Microstructure of 7B04 Al alloy sheet

圖2 超塑性拉伸試樣尺寸示意圖Fig.2 Schematic diagram of dimensions of tensile specimen

圖3所示為 7B04 鋁合金試樣超塑性拉伸變形前后的宏觀形貌?？梢姡?B04鋁合金超塑性拉伸之后基本沒有縮頸現(xiàn)象產生，變形相對比較均勻，這說明7B04鋁合金具有較好的超塑性，超塑性的均勻變形有效抑制了頸縮的出現(xiàn)。對比圖3中7B04鋁合金在不同變形條件下拉伸后的試樣尺寸，不難發(fā)現(xiàn)當變形溫度為530℃，應變速率為0.0003s-1時，材料的伸長率達到最大，為1105%，說明在此變形條件下材料具有最佳超塑性。

圖3 不同變形條件下超塑性拉伸后試樣宏觀形貌；；；Fig.3 Elongated specimens under different deformation =0.0003s-1；=0.001s-1；=0.003s-1；=0.01s-1

圖4 7B04鋁合金在不同變形條件下超塑性拉伸后的真應力-真應變曲線 (a)470℃；(b)485℃；(c)500℃；(d)515℃；(e)530℃Fig.4 True stress-true strain curves of 7B04 Al alloy after superplastic tensile deformation under different deformation parameters(a)470℃；(b)485℃；(c)500℃；(d)515℃；(e)530℃

7B04鋁合金在不同變形條件下超塑性拉伸后的真應力-真應變曲線，如圖4所示，從圖中可以看出，在不同變形條件下，7B04鋁合金超塑性變形的應力-應變曲線形狀基本相同：變形初期，應力隨應變的增加而迅速升高，表現(xiàn)出明顯的硬化效應；隨著變形的進行，由于應變軟化效應，應力緩慢增加，進入穩(wěn)態(tài)流動階段，且應變速率越低，穩(wěn)態(tài)流動階段越長；最后，當斷裂發(fā)生時應力急劇下降。

2.1.2 變形參數(shù)對超塑性的影響

7B04鋁合金超塑性變形過程中峰值應力隨變形溫度的變化如圖5所示，從圖中可以看出，在相同應變速率下，材料的峰值應力隨著變形溫度的升高而減小，這是因為升高變形溫度，提高了原子的平均動能，促進了位錯的運動，使晶界滑移易于發(fā)生，降低了變形過程中的流動應力；同時高溫下動態(tài)回復產生一定的軟化作用使材料的流動應力減小。

圖5 7B04鋁合金超塑性變形過程中峰值應力與變形溫度關系曲線Fig.5 Relationship curves between peak stress and temperature of 7B04 Al alloy during superplastic deformation

圖6是7B04鋁合金伸長率隨變形溫度的變化曲線，當應變速率在0.0003s-1和0.001s-1之間時，相同應變速率下，隨著變形溫度逐漸上升，伸長率呈陡峭上升趨勢，且當變形溫度為530℃時，伸長率達到最大值1105%，說明在此變形條件下，升高變形溫度能夠顯著提高材料的超塑性能；當應變速率為0.003s-1時，伸長率隨變形溫度的升高而緩慢增加，但變化范圍較小，在210%～370%之間，說明在此應變速率下，變形溫度對材料的超塑性能影響較?。划攽兯俾蕿?.01s-1時，伸長率隨著變形溫度的升高而先增加后減小，但波動不大，在130%～240%之間，說明此時材料的超塑性能較差，且對變形溫度的變化不敏感。

圖6 7B04鋁合金伸長率-變形溫度曲線Fig.6 Elongation-deformation temperature curves of 7B04 Al alloy

7B04鋁合金超塑性變形過程中峰值應力隨應變速率的變化如圖7所示，可見，在相同溫度下，材料的峰值應力隨著應變速率的提高而逐漸增加，這是由于在高應變速率下，位錯密度迅速增加，導致位錯塞積，而變形過程中需要擴散蠕變和位錯滑移來調節(jié)晶界的滑移，所以此時晶界滑移受阻，易于在三角晶界處、晶內的析出相周圍產生應力集中，從而使拉伸過程中的流動應力增加。圖8是7B04鋁合金伸長率隨應變速率變化的曲線圖?？梢?，相同變形溫度下，隨著應變速率的增加，材料的伸長率逐漸減??；當應變速率在0.0003s-1和0.001s-1之間時，相同溫度下， 材料的伸長率變化緩慢，且均具有較高的數(shù)值，說明在此變形條件下，材料具有較好的超塑性；當應變速率為0.003s-1時，伸長率急速下降，但其最小值為210%，仍有超塑性；當應變速率為0.01s-1時，不同變形溫度下的伸長率均較小，且相差不大，說明應變速率升高到一定程度后，材料的超塑性變差。因此，7B04鋁合金的超塑性變形應變速率應小于0.01s-1。

圖7 7B04鋁合金峰值應力-應變速率曲線Fig.7 Peak stress-strain rate curves of 7B04 Al alloy

圖8 7B04鋁合金伸長率-應變速率曲線Fig.8 Elongation-strain rate curves of 7B04 Al alloy

2.2 應變速率敏感性指數(shù)m值和變形激活能Q值

鋁合金在高溫超塑性變形過程中流動應力一般較低，穩(wěn)態(tài)流動階段應變速率與流動應力之間的關系可采用Arrhenius關系式[20]來描述：

(1)

式中：A為材料常數(shù)；n為應力指數(shù)，n=1/m，m為應變速率敏感性指數(shù)；Q為變形激活能；R為氣體常數(shù)，R=8.314J/(mol·K)；T為熱力學溫度，即絕對溫度。

對式(1)兩邊取自然對數(shù)，得出：

(2)

在一定的變形溫度和應變下，式(2)中l(wèi)nA-Q/RT可視為一常數(shù)Km，則m值可通過下式求解：

(3)

圖9 7B04鋁合金在不同變形溫度下的lnσ 曲線Fig.9 lnσ curves of 7B04 Al alloy at different deformation temperatures

相關文獻指出[21]，晶界滑動對變形總量的貢獻越大，應變速率敏感指數(shù)m值也越大，而且當m值達到0.5時，超塑性變形的主要機制便是晶界滑動。因此，由圖9可知，當變形溫度為470℃時，m值僅為0.3，說明此時晶界滑動對變形總量的貢獻較小，可能存在其他變形協(xié)調機制；當變形溫度從470℃上升到500℃時，m值增加，晶界滑動對變形總量的貢獻逐漸增大，直至變形溫度大于500℃時，晶界滑動成了主要的變形機制，對變形總量的貢獻最大。

(4)

(5)

繪制lnσ-1/T曲線，并進行線性擬合，其結果如圖10所示。根據(jù)式(5)求取各應變速率下的變形激活能，m值取不同變形溫度下的平均值，從而獲得材料的超塑性變形激活能為190kJ/mol，高于純鋁的晶格自擴散激活能142kJ/mol，這說明7B04鋁合金在超塑變形中的晶界滑動主要受晶內擴散控制[19]。

圖10 7B04鋁合金lnσ -1/T擬合曲線Fig.10 lnσ -1/Tfitting curves of 7B04 Al alloy

2.3 斷口形貌

圖11是7B04鋁合金在同一應變速率0.0003s-1，不同變形溫度470℃和485℃下超塑性拉伸變形后斷口的掃描電鏡照片，從圖中可以看出斷口處存在一些須狀物，說明材料在低應變速率下變形時晶界附近有液相產生；對比圖11(a)，(b)不難發(fā)現(xiàn)，當變形溫度為485℃時，須狀物的數(shù)量有所增加，且較密集。

圖11 7B04鋁合金在同一應變速率0.0003s-1不同變形溫度下拉伸斷口的SEM照片 (a)470℃；(b)485℃Fig.11 SEM images of fracture of 7B04 Al alloy at the strain rate of 0.0003s-1 at different deformation temperatures (a)470℃；(b)485℃

圖12 7B04鋁合金在同一變形溫度530℃不同應變速率下拉伸斷口的SEM照片 (a)0.0003s-1；(b)0.001s-1Fig.12 SEM images of fracture of 7B04 Al alloy at the deformation temperature of 530℃ at different strain rates (a)0.0003s-1；(b)0.001s-1

7B04鋁合金在同一變形溫度530℃，不同應變速率0.0003s-1和0.001s-1下超塑性變形后斷口的掃描電鏡照片如圖12所示，可以看出，在高應變速率下，只有部分區(qū)域存在須狀物，且少而短。

上述現(xiàn)象說明7B04鋁合金在超塑性變形過程中，晶界處的低熔點相會熔化形成須狀物[9,22]；且隨著變形溫度的升高和應變速率的降低，須狀物的數(shù)量和長度都有所增加，說明須狀物的數(shù)量和形貌與變形溫度和應變速率有關。由于須狀物的存在能夠在一定程度上提高晶界的結合強度、促進應力集中的釋放[12,23,24]，從而使晶界易于滑動，說明在高溫和低速下，7B04鋁合金超塑變形過程中存在液相協(xié)調機制，提高了材料的超塑性。但相關研究表明[23,24]，液相過多時，晶界、界面由于失去了結合力，材料的超塑性能會有所降低。因此，只有適量的液相才能起到協(xié)調晶界滑動的作用，從而在一定程度上提高7B04鋁合金的超塑性能，過多過少都不利于超塑變形。

3 結論

(1)7B04鋁合金在超塑性變形中，隨著變形溫度的升高和應變速率的降低，峰值應力逐漸減小，而伸長率逐漸增加；當應變速率在0.0003s-1和0.001s-1之間時，升高變形溫度能夠顯著提高材料的超塑性能；在高應變速率下，材料的超塑性能較差；7B04鋁合金的超塑性變形溫度范圍為470～530℃，應變速率要小于0.01s-1；當變形溫度為530℃，應變速率為0.0003s-1時，材料具有最佳的超塑性能。

(2)7B04鋁合金的應變速率敏感性指數(shù)m值均大于0.3，且隨著變形溫度的升高而逐漸增加，說明低溫下晶界滑動對超塑變形總量的貢獻較少，存在其他超塑性變形機制，而當變形溫度高于500℃時，晶界滑動為主要的變形機制。變形激活能Q=190kJ/mol，高于純鋁的晶格自擴散激活能，表明7B04鋁合金超塑性變形主要受晶內擴散控制；7B04鋁合金超塑性拉伸后的斷口處有須狀物存在，說明變形過程中晶界附近有液相產生，且一定的液相有助于調節(jié)晶界的滑動，從而提高材料的超塑性能。

[1] 黃敏，陳軼，李超，等. 7A12-T7352鋁合金高溫力學性能及斷裂行為研究[J]. 航空材料學報，2014，34(1)：82-85.

HUANG M，CHEN Y，LI C，et al. High temperature mechanical properties and fracture characteristics of 7A12-T7352[J]. Journal of Aeronautical Materials，2014，34(1)：82-85.

[2] 佘歡，疏達，儲威，等. Fe和Si雜質元素對7XXX系高強航空鋁合金組織及性能的影響[J]. 材料工程，2013，(6)：92-98.

SHE H，SHU D，CHU W，et al. Effects of Fe and Si impurities on the microstructure and properties of 7XXX high strength aircraft aluminum alloys[J]. Journal of Materials Engineering，2013，(6)：92-98.

[3] 姬浩. 7000系高強鋁合金的發(fā)展及其在飛機上的應用[J]. 航空科學技術，2015，26(6)：75-78.

JI H. Development and application of 7000 high strength aluminum alloys on airplane [J]. Aeronautical Science & Technology，2015，26(6)：75-78.

[4] 劉昌斌，夏長清，戴曉元. 高強高韌鋁合金的研究現(xiàn)狀及發(fā)展趨勢[J]. 礦冶工程，2003，23(6)：74-82.

LIU C B，XIA C Q，DAI X Y. Present states of research and developing trends of high strength and high toughness aluminum alloy[J]. Mining and Metallurgical Engineering，2003，23(6)：74-82.

[5] 謝優(yōu)華，陸 政，戴圣龍. Al-6Mg-0.2Sc鋁合金高溫力學行為研究[J]. 材料工程，2007，(3)：49-52.

XIE Y H，LU Z，DAI S L. Mechanical behavior of Al-6Mg-0.2Sc aluminum alloys at elevated temperature[J]. Journal of Materials Engineering，2007，(3)：49-52.

[6] 寧愛林，曾蘇民. 時效制度對7B04鋁合金組織和性能的影響[J]. 中國有色金屬學報，2004，14(6)：922-927.

NING A L，ZENG S M. Effects of ageing system on microstructure and mechanical properties of 7B04 aluminium alloy[J]. The Chinese Journal of Nonferrous Metals，2004，14(6)：922-927.

[7] 張智慧，左玉婷，劉淑鳳. 7B04鋁合金板沿厚度方向顯微組織、織構及力學性能的研究[J]. 電子顯微學報，2011，30(4-5)：322-327.

ZHANG Z H，ZUO Y T，LIU S F. Study on microstructure, mixture and tensile property of 7B04 aluminium alloy plate along thickness direction [J].Journal of Chinese Electron Microscopy Society，2011，30(4-5)：322-327.

[8] 喬志偉. 變形條件對7475鋁合金組織和超塑性性能的影響[D]. 湖南：湖南大學，2007.

QIAO Z W. The effect of deformation processing on the microstructure and superplastic properties of 7475 Al alloy[D]. Hunan: Hunan university，2007.

[9] TAKAYAMA Y, TOZAMA T, KATO H. Superplasticity and thickness of liquid phase in the vicinity of solidus temperature in a 7475 aluminum alloy[J]. Acta Materialia, 1999, 47(4): 1263-1270.

[10] SMOLEJ A, GNAMUS M, SLACK E. The influence of the thermomechanical processing and forming parameters on superplastic behavior of the 7475 aluminium alloy[J]. Journal of Materials Processing Technology, 2001, 118(1-3): 397-402.

[11] CHEN C L, TAN M J. Effect of grain boundary character distribution (GBCD) on the cavitation behavior during superplastic deformation of Al 7475 [J]. Materials Science & Engineering: A, 2002, 338(1-2): 243-252.

[12] IWASAKI H, MABUCHI M, HIGASHI K. Plastic cavity growth during superplastic flow in AA 7475Al alloy containing a small amount of liquid[J]. Acta Materialia, 2001, 49(12): 2269-2275.

[13] 張林. 噴射成形7475鋁合金熱處理及超塑性研究[D]. 上海：上海交通大學，2010.

ZHANG L. Study on heat treatment process and superplasticity of 7475 alloy by spray forming[D]. Shanghai: Shanghai Jiao Tong University，2010.

[14] WANG T, YIN Z M, SUN Q. Effect of homogenization treatment on microstructure and hot workability of high strength 7B04 aluminium alloy[J]. Transactions of Nonferrous Metals Society of China, 2007, 17(2): 335-339.

[15] LI Z H, XIONG B Q, ZHANG Y A, et al. Effects of the two-step ageing treatment on the microstructure and properties of 7B04 alloy pre-stretched thick plates[J]. Rare Metals, 2007, 26(3): 193-199.

[16] LI Z H, XIONG B Q, ZHANG Y A, et al. Microstructural evolution of aluminum alloy 7B04 thick plate by various thermal treatments[J]. Transactions of Nonferrous Metals Society of China, 2008, 18(1): 40-45.

[17] GAO F H, LI N K, TIAN N, et al. Overheating temperature of 7B04 highstrength aluminum alloy[J]. Transactions of Nonferrous Metals Society of China, 2008, 18(2): 321-326.

[18] 景武，趙懌，許廣興，等. 7B04鋁合金薄板成形過程中的退火工藝優(yōu)化[J]. 沈陽理工大學學報，2012，31(2)：91-94.

JING W，ZHAO Y，XU G X，et al. Optimum annealing process scheme in aluminum alloy 7B04 sheet formation[J]. Journal of Shenyang Ligong University，2012，31(2)：91-94.

[19] 李超，萬敏，黃霖. 7B04鋁合金蠕變過程中析出相的影響因素[J]. 航空材料學報，2009，29(2)：13-17.

LI C，WAN M，HUANG L. Effect on precipitation phases of 7B04 aluminum alloy during creep process[J]. Journal of Aeronautical Materials，2009，29(2)：13-17.

[20] LIN D, SUN F. Superplasticity in a large-grained TiAl alloy[J]. Intermetallics, 2004, 12(7): 875-883.

[21] 張凱鋒，尹德良，王國峰，等. 熱軋AZ31鎂合金超塑變形中的微觀組織演變及斷裂行為[J]. 航空材料學報，2005，25(1)：5-10.

ZHANG K F，YIN D L，WANG G F，et al. Microstructure evolution and fraction behavior in superplastic deformation of hot-rolled AZ31 Mg alloy[J]. Journal of Aeronautical Materials，2005，25(1)：5-10.

[22] CHEN K C, TAN M J. Cavity growth and filament formation of superplastically deformed Al 7475 alloy[J]. Materials Science & Engineering: A, 2001, 298(1): 235-244.

[23] 張福全，喬志偉，陳振華，等. 微量液相與7475鋁合金的超塑性變形[J].礦冶工程，2007，27(4)：71-73.

ZHANG F Q，QIAO Z W，CHEN Z H，et al. Effects of small amount of liquid on the superplasicity deformation of 7475 aluminium alloy[J]. Mining and Metallurgical Engineering，2007，27(4)：71-73.

[24] 許曉靜. 微量液相與鋁基復合材料高應變速率超速性[J]. 宇航材料工藝，1999，29(5)：41-43.

XU X J. Role of liquid in high strain rate superplasticity of aluminum matrix composites[J]. Aerospace Materials and Technology，1999，29(5)：41-43.

(本文責編：齊書涵)

Superplastic Deformation Behavior of 7B04 Al Alloy

ZHANG Ning1,2,3，WANG Yao-qi1,2,3，HOU Hong-liang1,2,3，ZHANG Yan-ling1,2,3，DONG Xiao-meng4，LI Zhi-qiang1

(1 Beijing Aeronautical Manufacturing Technology Research Institute，Beijing 100024，China；2 Aeronautical Key Laboratory for Plastic Forming Technology，Beijing 100024，China；3 Beijing Key Laboratory of Digital Plasticity Forming Technology and Equipment，Beijing 100024，China；4 University of Science and Technology Beijing，Beijing 100083，China)

The superplastic tensile tests were carried out for 7B04 Al alloy under a temperature range of 470-530℃ at a strain rate range of 0.0003-0.01s-1. Then,the superplastic deformation behavior and deformation mechanism of the material were investigated. The results show that with the increase of temperature and decrease of strain rate, the flow stress of 7B04 Al alloy is decreased and the elongation is increased.When the deformation temperature is at 530℃ and the strain rate is 0.0003s-1, the elongation reaches the maximum value of 1105%, which indicates that at this process parameter the optimal superplasticity is obtained; the strain rate sensitivity indexmvalue is above 0.3 and increases as the temperature going up; when the deformation temperature is at the range of 500-530℃, themvalue is above 0.5,the superplastic deformation mechanism of 7B04 Al alloy predominated by grain boundary sliding; the superplastic deformation activation energyQvalue is 190kJ/mol which indicates that the superplastic deformation is significantly controlled by the intercrystalline diffusion; the liquid phase will be produced along the grain boundaries of 7B04 Al alloy during superplastic deformation,which can enhance the superplasticity with a moderate amount.

7B04 Al alloy；superplasticity；deformation behavior；deformation mechanism

10.11868/j.issn.1001-4381.2016.000686

TG135.3；TG146.2

A

1001-4381(2017)04-0027-07

國家自然科學基金項目(51334006,51405457); 國家重點基礎研究發(fā)展計劃項目(2011CB012803)

2016-06-05;

2017-01-03

侯紅亮(1963-)，男，研究員，主要從事鈦合金、鋁合金等先進航空材料鈑金成形技術研究工作，聯(lián)系地址：北京市340信箱106分箱(100024)，E-mail: hou_ll@163.com