李美春，鄧運(yùn)來(lái)，唐建國(guó)，萬(wàn) 里，張新明

（中南大學(xué)1.材料科學(xué)與工程學(xué)院，長(zhǎng)沙 410083；2.有色金屬材料科學(xué)與工程教育部重點(diǎn)實(shí)驗(yàn)室，長(zhǎng)沙 410012）

0 引 言

6061鋁合金是可熱處理強(qiáng)化的中高強(qiáng)度鋁合金，具有良好的機(jī)械加工性能和抗腐蝕能力，廣泛應(yīng)用于要求具有一定強(qiáng)度和高抗蝕能力的各種結(jié)構(gòu)件中，如汽車(chē)、船舶、軌道車(chē)輛等的結(jié)構(gòu)件［1］。改進(jìn)鑄造、軋制及熱處理等工藝是提高該合金綜合性能的主要途徑。固溶和時(shí)效熱處理是提高力學(xué)性能的主要手段，但是合金在固溶時(shí)的淬火過(guò)程中會(huì)產(chǎn)生較大的殘余應(yīng)力［2－3］，殘余應(yīng)力的存在會(huì)極大地影響板材的加工性能、疲勞壽命以及耐蝕性能等［4］。目前，消減殘余應(yīng)力的主要方法有機(jī)械法和熱處理方法［5－7］。機(jī)械法在工程實(shí)際應(yīng)用中需要較高的成本。因此，對(duì)于6061鋁合金，研究合適的熱處理工藝在提高其力學(xué)性能的同時(shí)降低殘余應(yīng)力具有重要的理論和實(shí)際意義。

目前國(guó)內(nèi)外對(duì)鋁合金板材的變形過(guò)程、熱處理過(guò)程、顯微組織及性能之間的關(guān)系做了大量的研究［8－10］，對(duì)鋁合金厚板淬火殘余應(yīng)力的相關(guān)研究報(bào)道也很多［5，11］，但大部分文獻(xiàn)對(duì)中厚板殘余應(yīng)力的研究都從宏觀(guān)角度出發(fā)，很少有人研究材料顯微組織與淬火殘余應(yīng)力之間的關(guān)系。研究表明鋁鎂硅合金在固溶熱處理過(guò)程中會(huì)發(fā)生顯著的再結(jié)晶［12］，再結(jié)晶過(guò)程中織構(gòu)有可能發(fā)生改變，而織構(gòu)的改變意味著材料的性能將發(fā)生改變。為此，作者以22mm厚的6061鋁合金熱軋中厚板為研究對(duì)象，研究了不同時(shí)間固溶處理后板材的組織、性能與表面殘余應(yīng)力，并計(jì)算Schmidt因子，從晶體學(xué)織構(gòu)的角度定性解釋了淬火殘余應(yīng)力差異的原因。

1 試樣制備與試驗(yàn)方法

試驗(yàn)材料為中國(guó)精美鋁業(yè)公司提供的厚度為22mm的6061鋁合金熱軋板。合金的化學(xué)成分如表1所示，熱軋總變形量為95%。為了減小固溶后淬火工藝對(duì)殘余應(yīng)力及性能的影響，取3塊尺寸為500mm×400mm×22mm的試樣在鹽浴爐中加熱到530℃后固溶保溫2h時(shí)，再放入另外3塊相同尺寸的試樣，繼續(xù)加熱保溫2h后再淬入室溫水中，淬火轉(zhuǎn)移時(shí)間小于5s，淬火后水溫不高于50℃。熱軋態(tài)試樣標(biāo)為HR，固溶保溫2h和4h的試樣分別記為S1、S2，將固溶試樣在180℃下進(jìn)行16h時(shí)效處理，分別記為A1、A2。

表1 6061鋁合金熱軋板的化學(xué)成分（質(zhì)量分?jǐn)?shù)）Tab.1 Chemical composition of 6061 Al alloy hot rolled plate（mass） %

采用盲孔法測(cè)試試樣的表面殘余應(yīng)力，選取中厚板對(duì)角線(xiàn)上1/4，1/2，3/4處的三個(gè)點(diǎn)為殘余應(yīng)力測(cè)試點(diǎn)，取三組試樣的平均值；將試樣制成半徑為5mm的圓柱試樣，在Instrong 8032型萬(wàn)能材料力學(xué)試驗(yàn)機(jī)上以2mm·min－1的速度測(cè)試?yán)煨阅埽挥肏V-5型小負(fù)荷維氏硬度計(jì)進(jìn)行硬度測(cè)試，載荷為29.4N，保載時(shí)間15s；在熱軋、固溶和時(shí)效態(tài)板材的表層和心部（1/2厚度處）取金相試樣，尺寸為10mm×10mm×2mm，經(jīng)粗磨、細(xì)磨、拋光后用體積分?jǐn)?shù)為40%HF溶液浸蝕，在XJP-6A型光學(xué)顯微鏡下觀(guān)察顯微組織；在試樣S1、S2、A1、A2的表層和中部分別切取薄片制備透射電鏡試樣，先預(yù)磨成厚約0.08mm的薄片，再?zèng)_成直徑3mm圓片后進(jìn)行雙噴減薄，電解液為25%HNO3＋75%CH3OH，采用液氮冷卻，溫度控制在－20℃以下，在Tecnai G220型透射電鏡上進(jìn)行微細(xì)結(jié)構(gòu)分析，加速電壓為200kV，并對(duì)組織中的粒子進(jìn)行能譜分析；織構(gòu)分析在金相試樣上采用Brucker D8Discovery型X射線(xiàn)儀測(cè)｛111｝、｛200｝、｛220｝、｛113｝4個(gè)不完全極圖，運(yùn)用球諧函數(shù)級(jí)數(shù)展開(kāi)法求取向分布函數(shù)（ODF）［13］，采用PSO 織構(gòu)組分法計(jì)算織構(gòu)體積分?jǐn)?shù)及真 ODF［14］。

2 試驗(yàn)結(jié)果與討論

2.1 熱軋板的組織與織構(gòu)

從圖1可以看出，熱軋板中的再結(jié)晶現(xiàn)象不明顯，表層和心部的組織有較大差異，表層組織模糊混亂，晶界不清晰；而心部為典型的變形組織，晶界較清晰。由圖2可知，熱軋板中存在粗大的板條狀粒子，且表層中粗大粒子尺寸和數(shù)目明顯大于心部的。經(jīng)能譜分析，該粒子為含鉻、錳和鐵的β化合物相。

圖1 熱軋?jiān)嚇覱M形貌Fig.1 The OMmorphology of hot-rolled specimen（a）surface and（b）core

從圖3中可見(jiàn)，表層的織構(gòu)與心部的差異較大。表層的比較散亂，而心部從 Goss（｛110｝〈001〉）至B（｛011｝〈112〉）到反 Goss（｛011｝〈011〉）的α取向線(xiàn)織構(gòu)組分明顯，為典型的軋制織構(gòu)；二者的β取向線(xiàn)組分都不明顯；除此以外還有少量的旋轉(zhuǎn)立方織構(gòu)（｛001｝〈110〉）（一種典型的剪切織構(gòu)），且試樣表層的旋轉(zhuǎn)立方織構(gòu)比心部的略強(qiáng)。

由以上分析可知，該板材在軋制過(guò)程中變形不均勻。由于軋輥與軋件之間的摩擦使得板材表面發(fā)生了沿軋向的剪切變形，其應(yīng)變模式為在平錘壓縮應(yīng)變的基礎(chǔ)上疊加一個(gè)沿軋向的剪切應(yīng)變，因而表層織構(gòu)中出現(xiàn)了剪切織構(gòu)，文獻(xiàn)［15－16］也得到相類(lèi)似的結(jié)果。對(duì)于中厚板，由于軋輥接觸弧長(zhǎng)ld與板材平均厚度h的比值比較小，心部幾乎不受摩擦力的影響，為典型的平錘壓縮變形模式，因此織構(gòu)表現(xiàn)為典型的軋制織構(gòu)。從表層到心部剪切應(yīng)變逐漸減小，因而表層的應(yīng)變模式復(fù)雜，其相應(yīng)的變形組織也較復(fù)雜。

圖2 熱軋?jiān)嚇拥腡EM形貌Fig.2 TEMmorphology of hot-rolled specimen：（a）surface and（b）core

2.2 不同時(shí)間固溶板的組織與織構(gòu)

從圖4中可以看出，板材在固溶時(shí)發(fā)生了不同程度的再結(jié)晶，與圖1熱軋時(shí)相比再結(jié)晶顯著。試樣S1的表層晶粒大小不均勻，在大量粗大晶粒中間有一些尺寸較小的晶粒，而其心部為類(lèi)似于變形組織的長(zhǎng)條形晶粒。試樣S2表層和心部不再呈細(xì)長(zhǎng)的纖維組織，表層出現(xiàn)了許多邊界為曲線(xiàn)的細(xì)小晶粒，心部晶粒向等軸晶發(fā)展，表明當(dāng)固溶時(shí)間從2h延長(zhǎng)到4h，板材再結(jié)晶程度增大。

圖3 熱軋?jiān)嚇拥腛DF圖Fig.3 The ODF of hot-rolled specimen：（a）surface and（b）core

由圖5可知，試樣S1中存在兩種不同形狀的粗大粒子，經(jīng)能譜分析知長(zhǎng)寬比較大的為含鉻、錳和鐵的β化合物相，而近似等軸狀粒子為含鐵、錳、鉻和硅的α相；試樣S1表層和心部的化合物相形態(tài)和尺寸差異較大，表層β相粗大粒子較多，心部幾乎全部為等軸的α相。與圖2比較可知，試樣S1、S2厚度方向組織中粗大粒子的差異是由于熱軋板厚度方向組織的不均勻性導(dǎo)致的。固溶4h后，粗大含鉻、錳和鐵的β化合物相基本都溶入了基體，α相變化不明顯。由此可見(jiàn)，固溶時(shí)間的延長(zhǎng)促進(jìn)了厚度方向組織的均勻化。

圖4 固溶試樣S1和S2的OM形貌Fig.4 OMmorphology of solution specimens（S1，S2）：（a）S1，surface；（b）S1，core；（c）S2，surface and（d）S2，core

圖5 固溶試樣S1和S2的TEM形貌Fig.5 TEMmorphology of solution specimens（S1，S2）：（a）S1，surface；（b）S1，core；（c）S2，surface and（d）S2，core

圖6是試樣S1、S2表層和心部的取向分布函數(shù)（ODF）的 Cube-RD 和 Cube-ND 取向線(xiàn)、β－取向線(xiàn)（β-fiber）、α－取向線(xiàn)（α-fiber）分析結(jié)果。Cube-ND和Cube-RD取向線(xiàn)分別表示了Cube－織構(gòu)沿軋向轉(zhuǎn)動(dòng)和沿板法向轉(zhuǎn)動(dòng)的取向密度f(wàn)（g）分布情況。從圖中可以看出，試樣S1中典型再結(jié)晶織構(gòu)｛001｝〈100〉Cube-織構(gòu)和｛001｝〈110〉旋轉(zhuǎn)立方織構(gòu)的表層取向密度分別為2.48，2.58，心部取向密度分別為1.88，0.14，差別明顯。可知表層再結(jié)晶織構(gòu)取向密度大于心部的；而β－取向線(xiàn)和α－取向線(xiàn)上的｛011｝〈211〉（B-），｛123｝〈634〉（S-），｛112｝〈111〉（C-）和｛011｝〈100〉（G-）等典型變形織構(gòu)組分的取向密度，心部均強(qiáng)于表層的，即表面再結(jié)晶程度大于心部的。試樣S2表層和心部｛001｝〈100〉Cube-織構(gòu)的f（g）分別為0.98，1.02，｛001｝〈110〉旋轉(zhuǎn)立方織構(gòu)的f（g）均為－0.82，即固溶4h后，材料中的｛001｝〈100〉Cube-織構(gòu)、｛001｝〈110〉旋轉(zhuǎn)立方織構(gòu)接近隨機(jī)織構(gòu)分布取向密度為1。試樣S2的β－取向線(xiàn)和α－取向線(xiàn)上的B-、S-、C-和 G-等典型織構(gòu)組分的取向密度均強(qiáng)于試樣S1的，且表層和心部各織構(gòu)組分的取向密度差均遠(yuǎn)小于后者的，即厚度方向上織構(gòu)分布均勻。

固溶熱處理過(guò)程中板材中發(fā)生了一定程度的再結(jié)晶，因此織構(gòu)有可能發(fā)生改變。研究表明［17－18］，無(wú)第二相粒子的再結(jié)晶主要由立方織構(gòu)構(gòu)成，但對(duì)含有大量第二相粒子的鋁合金，由于第二相粒子促進(jìn)了再結(jié)晶晶核的形成（即PSN效應(yīng)），再結(jié)晶形核機(jī)制有所改變，故再結(jié)晶織構(gòu)組分也會(huì)發(fā)生改變。關(guān)于再結(jié)晶過(guò)程由PSN所控制的鋁合金的再結(jié)晶織構(gòu)，Humphreys［19］、Qrsund等［17］采用 X射線(xiàn)極圖方法研究了第二相粒子對(duì)再結(jié)晶織構(gòu)的影響，結(jié)果表明PSN導(dǎo)致隨機(jī)織構(gòu)。而Engler［18］等認(rèn)為PSN機(jī)制導(dǎo)致取向密度較高且繞法向旋轉(zhuǎn)一定角度（20°～25°）的Cube-ND立方織構(gòu)。試樣S1和S2表層和心部都有Cube-ND立方織構(gòu)的存在，試樣S1的Cube-ND立方織構(gòu)的取向密度大于試樣S2的，且表層取向密度大于心部的。Cube-ND立方織構(gòu)的出現(xiàn)表明固溶過(guò)程中再結(jié)晶的形核機(jī)制主要是PSN。

圖6 固溶試樣S1和S2的ODF取向線(xiàn)分析結(jié)果Fig.6 The analysis results of orientation line of the ODF of solution specimens（S1，S2）

2.3 不同時(shí)間固溶后中厚板的殘余應(yīng)力與力學(xué)性能

由表2可見(jiàn)，試樣S2沿軋向的表面殘余應(yīng)力σ1與試樣S1相比，平均下降了42.5%，最大下降率為57%；沿橫向表面殘余應(yīng)力σ2平均下降22.7%，最大降低率為49.4%。

中厚板表面殘余應(yīng)力處于平面應(yīng)力狀態(tài)，應(yīng)用第四強(qiáng)度理論計(jì)算的綜合殘余應(yīng)力σr進(jìn)行評(píng)估較為恰當(dāng)，其計(jì)算公式為

計(jì)算得出試樣S1、S2的表面綜合殘余應(yīng)力如表2所示，試樣S2表面綜合殘余應(yīng)力明顯低于S1的。

表2 固溶試樣S1和S2的表面殘余應(yīng)力Tab.2 The surface residual stress of solution specimens S1，S2

由表3可看出，固溶時(shí)間對(duì)中厚板力學(xué)性能的影響并不明顯，板材軋向力學(xué)性能略?xún)?yōu)于橫向的。試樣A1、A2的軋向和橫向屈服強(qiáng)度均大于300MPa，可知試樣S2的殘余應(yīng)力約等于0.2σ0.2。龔海等［20］用裂紋柔度法對(duì)鋁合金厚板的淬火應(yīng)力進(jìn)行了研究，結(jié)果表明鋁合金厚板殘余應(yīng)力的分布規(guī)律與Hossain［21］所研究的結(jié)果相符，殘余應(yīng)力表現(xiàn)出外壓應(yīng)力內(nèi)拉應(yīng)力的分布特征，殘余應(yīng)力在厚度方向呈M型分布，且壓應(yīng)力絕對(duì)值大于拉應(yīng)力絕對(duì)值，即表面殘余應(yīng)力絕對(duì)值最大。故內(nèi)部殘余應(yīng)力小于0.2σ0.2，滿(mǎn)足產(chǎn)品的工藝要求，不需要再進(jìn)行特定的殘余應(yīng)力消減處理。

表3 時(shí)效試樣A1和A2的力學(xué)性能Tab.3 The mechanical properties of aging specimens A1，A2

由圖7可知，試樣A1、A2硬度沿厚度方向存在一定的差異，心部硬度大于表層硬度。相對(duì)于試樣A2，試樣A1的表層硬度略高，心部硬度略低。比較硬度隨厚度變化的方差可知，試樣A2厚度方向硬度較均勻。

圖7 時(shí)效試樣A1和A2沿厚度方向的硬度分布及方差Fig.7 Hardness distribution（a）along the thickness of aging specimens（A1，A2）and their variances（b）

2.4 織構(gòu)對(duì)殘余應(yīng)力的影響

淬火后試樣的殘余應(yīng)力主要是由快速冷卻過(guò)程中熱應(yīng)力產(chǎn)生的不均勻塑性應(yīng)變?cè)斐?，由前面分析可知固溶熱處理時(shí)間主要影響合金的織構(gòu)與含鐵相，其中對(duì)織構(gòu)的影響較大，并且含鐵相對(duì)塑性變形幾乎沒(méi)有影響，且淬火工藝相同，因此，可認(rèn)為殘余應(yīng)力的差異主要是由于織構(gòu)的差異而導(dǎo)致的。

熱應(yīng)力是淬火殘余應(yīng)力產(chǎn)生的原因，因此從熱應(yīng)力入手分析。表層由于冷卻速率快，因此以拉應(yīng)力為主，為了分析問(wèn)題的方便，假設(shè)其為平面應(yīng)力狀態(tài)且兩主軸方向的應(yīng)力大小相等，即

式中：σij為淬火過(guò)程中所產(chǎn)生的熱應(yīng)力。

同樣心部由于冷卻速率較慢，以壓應(yīng)力為主，假設(shè)其為平面應(yīng)力狀態(tài)且兩主軸方向的應(yīng)力大小相等，即

由于在變形的起始階段，晶粒間的交互作用較小，可采用Sachs模型的假設(shè)，即假設(shè)各晶粒所承受的應(yīng)力狀態(tài)相同，表層所有晶粒的應(yīng)力狀態(tài)見(jiàn)式（2），心部所有晶粒的應(yīng)力狀態(tài)見(jiàn)式（3）。同時(shí)由于淬火熱應(yīng)力所引起的塑性應(yīng)變非常小，對(duì)應(yīng)塑性變形的起始階段，為塑性變形的易滑移階段，單個(gè)滑移系參與塑性變形［22］。根據(jù)Schmidt定律［23］，分切應(yīng)力最大的滑移系其滑移面為n＝（n1，n2，n3）、滑移方向?yàn)閎＝（b1，b2，b3），即為參與塑性變形的滑移系，其Schmidt因子為

多晶體材料的F值為其各織構(gòu)組分的Fi值按體積分?jǐn)?shù)φi加權(quán)平均，即：

分別對(duì)試樣S1和S2表層和心部織構(gòu)進(jìn)行織構(gòu)體積分?jǐn)?shù)的定量分析和F值的定量分析。試樣S1、S2表層和心部的 F 值分別約為0.13，0.21，0.16，0.19。試樣S1心部F值最大，表層的F值最小。由式（5）可推知材料屈服強(qiáng)度與F值成反比，故在淬火過(guò)程中試樣S1心部較容易變形，表層最難發(fā)生變形，這會(huì)導(dǎo)致淬火快速冷卻過(guò)程中由熱應(yīng)力產(chǎn)生的塑性應(yīng)變的不均勻性加劇，而淬火殘余應(yīng)力的產(chǎn)生主要是由快速冷卻過(guò)程中熱應(yīng)力產(chǎn)生的不均勻塑性應(yīng)變?cè)斐傻?；故淬火過(guò)程中表面會(huì)形成較大的殘余應(yīng)力，與試樣S1的表層殘余應(yīng)力較大相符合。試樣S2表層與心部的F值差異較小，與其厚度方向力學(xué)性能較均勻符合。

顯然，中厚板淬火殘余應(yīng)力的產(chǎn)生還受到其他眾多因素的影響，實(shí)際情況遠(yuǎn)比計(jì)算假設(shè)復(fù)雜得多，但從晶體學(xué)織構(gòu)的角度出發(fā)，為淬火殘余應(yīng)力的研究提供了一個(gè)新的方向。

3 結(jié) 論

（1）固溶保溫時(shí)間從2h延長(zhǎng)到4h，熱軋中厚板厚度方向組織和性能更均勻，同時(shí)促進(jìn)了再結(jié)晶過(guò)程的進(jìn)行，使得板材織構(gòu)明顯強(qiáng)化。

（2）在平面應(yīng)力狀態(tài)下只經(jīng)過(guò)2h固溶后淬火板材心部的Schmidt因子最大，為0.21，因而在淬火過(guò)程中最容易發(fā)生變形，同時(shí)其表層和心部的Schmidt因子差異比經(jīng)過(guò)4h固溶保溫處理試樣的大；這與試樣S1的表層殘余應(yīng)力較大及試樣A2厚度方向力學(xué)性能更均勻相符合。

（3）采用合適的熱處理方法可以在不影響材料性能的前提下，實(shí)現(xiàn)降低表面殘余應(yīng)力的目的。

［1］王祝堂.國(guó)際協(xié)定牌號(hào)的加工純鋁及加工鋁合金［J］.輕合金加工技術(shù)，1996，24（4）：24-27.

［2］MICHAEL B P，MICHAEL R H.Residual stress，stress relief，and Inhomogeneity in aluminum plate［J］.Scripta Materialia，2002，46（4）：77-82.

［3］NOYAN I C ，COHEN J B .Residual stress［M］.United States：Springer-Verlag，New York Inc，1986：220-230.

［4］WALKER D M，HOMR Y.Residual stress analysis of aircraft aluminum forgings［J］.Advanced Materials ＆ Processes，2002，160（6）：57-60.

［5］ROBINSON J S，TANNER D A.Reducing residual stress in 7050aluminum alloy die forgings by heat treatment［J］.Journal of Engineering Materials and Technology，2008，130：1-8.

［6］KOC M，CULP J，ALTAN T.Prediction of residual stresses in quenched aluminum blocks and their reduction through cold working processes［J］.Journal of Materials Processing Technology，2006，174（1/3）：342-354.

［7］LADOS D A，APELIAN D，WANG L.Minimization of residual stress in heat-treated Al-Si-Mg cast alloysusing uphill quenching：mechanisms and effects on static and dynamic properties［J］.Materials Science and Engineering：A，2010，527：3159-3165.

［8］MAISONNETTE D，SUERY M，NELIAS D，et al.Effects of heat treatments on the microstructure and mechanical properties of a 6061aluminium alloy［J］.Materials Science and Engi-neering：A，2011，528：2718-2724.

［9］LEE S H，SAITO Y，SAKAI T，et al.Microstructures and mechanical properties of 6061aluminum alloy processed by accumulative roll-bonding［J］.Materials Science and Engineering：A，2002，325：228-235.

［10］李衛(wèi)超，鄒勇志，曾建民.固溶溫度和時(shí)間對(duì)ZL114A合金組織的影響［J］.機(jī)械工程材料，2008，32（11）：25-27.

［11］LIAO Kai，JING Xin-xing ，ZHOU Ya，et al.Research of micromechanism on residual stress in aluminum alloy thick plate［J］.Advanced Materials Research，2010，146/147：233-237.

［12］GJESTLAND H，LOHNE O，RYUMN.The effect of quenching on recrystallization after hot working［J］.Metallurgical Transactions：A，1978，9（11）：1982-1983.

［13］BUNGE H J.Mathematische methoden der texturanalyse［M］.Berlin：Akademie-Verlag，1969.

［14］TANG J G，ZHANG X M，DENG Y L，et al.Texture decomposition with particle swarm optimization method［J］.Computational Materials Science，2006，38（2）：395-399.

［15］唐建國(guó)，張新明，徐敏，等.3104鋁合金熱粗軋板的織構(gòu)梯度及其熱壓縮變形后退火織構(gòu)演變的影響［J］.中國(guó)有色金屬學(xué)報(bào)，2012，32（10）：1932-1940.

［16］鄧運(yùn)來(lái)，張新明，劉瑛，等.幾何因素與摩擦耦合對(duì)高純鋁箔剪切織構(gòu)的影響［J］.中國(guó)有色金屬學(xué)報(bào)，2002，12（4）：634-638.

［17］ORSUND R，NES E.Effect of particles on recrystallization textures in aluminium-manganese alloys［J］.Scripta Metallurgic，1988，22：665-669.

［18］ENGLER O，KONG X W，LUCKE K.Recrystallisation textures of particle containing Al-Cu and Al-Mn single crystals［J］.Acta Materialia，2001，49（8）：1701-1715.

［19］HUMPHREYS F J，HATHERLY M.Recrystallization and related annnealing phenomena［M］.Oxford，New York：Pergamon，1996：235-266.

［20］龔海，吳運(yùn)新，廖凱.不同淬火工藝對(duì)7075鋁合金厚板殘余應(yīng)力的影響［J］.中南大學(xué)學(xué)報(bào)：自然科學(xué)版，2010，41（4）：1354-1359.

［21］HOSSAIN S，DAYMOND MR，TRUMAN C E，et al.Prediction and measurement of residual stresses in quenched stainless-steel Spheres［J］.Materials Science and Engineering：A，2004，373：339-349.

［22］王自強(qiáng)，段祝平.塑性細(xì)觀(guān)力學(xué)［M］.北京：北京大學(xué)出版社，1995.

［23］REID C N.Deformation geometry for materials scientists［M］.Oxford，New York：Pergamon Press，1973.