半圓形拉深筋圓角半徑及金屬流動方向?qū)罱钭枇Φ挠绊?/h1>
2015-06-14 07:38:06那景新劉海鵬閆亞坤
吉林大學(xué)學(xué)報(工學(xué)版) 2015年5期
那景新，白 霜，劉海鵬，閆亞坤
（吉林大學(xué) 汽車仿真與控制國家重點實驗室，長春130022）
0 引 言
近年來隨著沖壓工藝的發(fā)展，沖壓制件日趨復(fù)雜，作為拉深工藝重要控制手段之一的拉深筋逐步成為拉深模具設(shè)計中的一個關(guān)鍵環(huán)節(jié)［1］。在汽車覆蓋件成形過程中，零件各部分成形條件不同，因此各部分成形所需的成形力不同，通常需要采用拉深筋來進(jìn)行控制，拉深筋參數(shù)的合理取值是控制金屬流動、防止出現(xiàn)起皺和破裂的重要手段［2－3］。
拉深筋最直接的作用是增加板料成形時的阻力［4］。拉深筋阻力有許多影響因素，國內(nèi)外都進(jìn)行了相關(guān)研究。例如，Demeri［5］通過試驗研究表明，無論何種材料，拉深筋高度和壓邊力的增加都能引起拉深筋阻力的增加；邢忠文等［6］系統(tǒng)地研究了拉深筋形式、參數(shù)、材料特性、潤滑條件、變形速度以及壓邊力等對拉深筋阻力的影響。然而，在沖壓過程中，而關(guān)于金屬流動方向并不都與拉深筋法向平行，金屬流動方向與拉深筋法向的夾角這一因素對拉深筋阻力的影響的研究卻很少。本文結(jié)合試驗重點研究金屬流動方向?qū)罱钭枇Φ挠绊?，并且同時考慮拉深筋圓角半徑這一影響因素。
本文專門設(shè)計了一套測量拉深筋阻力的試驗裝置，分別測取了以圓角半徑和金屬流動方向為變量的試驗數(shù)據(jù)，分析得出圓角半徑及金屬流動方向?qū)罱钭枇Φ挠绊懸?guī)律：圓角半徑對拉深筋阻力的影響呈現(xiàn)近似線性關(guān)系；金屬流動方向?qū)罱钭枇Φ挠绊懗尸F(xiàn)近似的拋物線關(guān)系。綜合兩個參數(shù)的影響規(guī)律并且基于回歸分析方法建立了拉深筋阻力的近似解析模型，并且驗證了回歸分析結(jié)果與試驗數(shù)據(jù)的一致性。最后對回歸系數(shù)進(jìn)行T 檢驗得出，相比于圓角半徑，金屬流動方向?qū)罱钭枇Φ挠绊懜@著，并且金屬流動方向和圓角半徑的耦合效應(yīng)較小，為拉深筋的進(jìn)一步研究奠定了基礎(chǔ)。
1 試驗裝置及方案
本文設(shè)計的試驗裝置及拉伸試驗機(jī)（測力裝置）的實物圖如圖1所示，試驗裝置的CATIA 效果圖如圖2所示。
圖1 拉深筋試驗裝置及拉伸試驗機(jī)實物Fig.1 Entity of experiment device of drawbead and tensile testing machine
圖2 拉深筋試驗裝置效果圖Fig.2 CATIA picture of experiment device of drawbead
根據(jù)實際沖壓需要，限定圓角半徑R 和金屬流動方向與拉深筋法向的夾角θ 的取值范圍分別為R＝4～16mm，θ＝0°～30°，基于差值原理，每個變量選 取3、4 個數(shù)值，其中R 取6、8、10、12 mm，θ 取0°、15°、30°。試驗所用的板料尺寸為190mm×40mm×1mm，材料為st14鋼板，性能如下：彈性模量為2.07×105MPa，屈服強(qiáng)度為210MPa，抗拉強(qiáng)度為270 MPa，泊松比為0.28，伸長率為0.36。
為了減少試驗過程中的裝卸操作工作量，試驗按照如下方式分組進(jìn)行：根據(jù)拉深筋圓角半徑R 將試驗分成四大組，每大組再按照金屬流動方向與拉深筋法向的夾角θ分成三小組，每小組安排3個試件，即每組參數(shù)進(jìn)行3次試驗，通過采用多次試驗取平均值的方法來減小試驗誤差所帶來的影響。
在本試驗裝置中，將水平拉深筋的中點作為圓心，通過將拉深筋扭轉(zhuǎn)到不同角度來模擬不同的金屬流動方向，從而可在不改變試件運動方向情況下形成金屬流動方向與拉深筋法向不同的夾角，如圖3所示，節(jié)約了試驗成本。
圖3 拉深筋角度布置Fig.3 Layout of drawbead angle
2 圓角半徑及金屬流動方向?qū)罱钭枇Φ挠绊?/h2>
由于試驗得到的拉深筋阻力是沿著板料寬度方向（40mm）的總阻力，將其除以板料寬度，即可得單位長度的拉深筋阻力。結(jié)合試驗數(shù)據(jù)，如表1所示，分析拉深筋的圓角半徑及金屬流動方向?qū)罱钭枇Φ挠绊懸?guī)律。
表1 拉深筋阻力試驗數(shù)據(jù)Table 1 Experimental data of drawbead restraining force
2.1 拉深筋圓角半徑對拉深筋阻力的影響
根據(jù)三組不同金屬流動方向下圓角半徑分別為6、8、10、12mm 時所測得的拉深筋阻力的試驗數(shù)據(jù)，分別繪制了金屬流動方向分別為0°、15°、30°時，以圓角半徑為橫坐標(biāo)，單位拉深筋阻力為縱坐標(biāo)的關(guān)系曲線，如圖4所示。
拉深筋阻力主要是由板料通過拉深筋時，發(fā)生彎曲和反彎曲及其接觸面的摩擦而形成的［7－8］。由圖4可以看出：當(dāng)金屬流動方向與拉深筋法向的夾角一定時，隨著半徑的增大，拉深筋阻力減小，這主要是由于拉深筋在金屬流動方向的截面的圓弧曲率減小，從而導(dǎo)致彎曲和反彎曲的阻力減小所引起的。圓角半徑對拉深筋阻力的影響呈現(xiàn)近似線性關(guān)系，并且影響程度相差不大。
圖4 拉深筋阻力與拉深筋圓角半徑的關(guān)系Fig.4 Relationships between drawbead restraining force and corner radius
2.2 金屬流動方向?qū)罱钭枇Φ挠绊?/h3>
為進(jìn)一步研究金屬流動方向?qū)罱钭枇Φ挠绊懸?guī)律，本文又分別依據(jù)四組不同圓角半徑下金屬流動方向與拉深筋法向夾角分別為0°、15°、30°時所測得的拉深筋阻力的試驗數(shù)據(jù)，分別繪制了圓角半徑分別為6、8、10、12 mm 時，以夾角為橫坐標(biāo)，單位拉深筋阻力為縱坐標(biāo)的關(guān)系曲線，如圖5所示。
圖5 拉深筋阻力與金屬流動方向的關(guān)系Fig.5 Relationships between drawbead restraining force and metal flow direction
圖5表明，金屬流動方向?qū)罱钭枇Φ挠绊懗尸F(xiàn)近似的拋物線關(guān)系，且在夾角為15°時達(dá)到最大值。這主要是因為在金屬流動方向與拉深筋法向的夾角不同時，與金屬流動方向平行的拉深筋的截面的圓弧曲率的不同引起的彎曲和反彎曲的阻力變化幅度不同，而且板料與拉深筋的接觸面積的不同引起的摩擦力變化幅度也不同，拉深筋阻力大小的改變是這兩個主要因素的綜合影響結(jié)果，并且不同夾角時的影響程度不同。
例如，當(dāng)夾角為0°～15°時，夾角的增大導(dǎo)致板料與拉深筋的接觸面積不斷增大，從而導(dǎo)致摩擦力增大；金屬流經(jīng)拉深筋時，與金屬流動方向平行的拉深筋的截面為橢圓形，夾角的增大導(dǎo)致橢圓的曲率減?。ㄇ拾霃皆龃螅?，那么由彎曲反彎曲所產(chǎn)生的阻力也會減小。
綜上并結(jié)合圖5分析可知，當(dāng)夾角為0°～15°時，金屬流經(jīng)拉深筋和筋槽產(chǎn)生的摩擦力增加的幅度大于由彎曲反彎曲產(chǎn)生的阻力減小的幅度，導(dǎo)致拉深筋阻力提高，并且提高幅度較小，表現(xiàn)為曲線上升趨勢緩慢。當(dāng)圓角半徑為12mm 時，上升幅度最小，為1.25N；當(dāng)圓角半徑為6 mm 和10mm 時，上升幅度最大，也僅為7.5N。當(dāng)夾角為15°～30°時，彎曲反彎曲產(chǎn)生阻力減小的幅度大于摩擦力增加的幅度，導(dǎo)致拉深筋阻力減小，并且減小幅度較大，表現(xiàn)為曲線下降趨勢明顯，尤其是當(dāng)拉深筋圓角半徑為10 mm 和12 mm 時，下降幅度分別為21.666N 和16.667N。
3 拉深筋阻力近似解析模型的建立
由圖4可以看出：圓角半徑對拉深筋阻力的影響呈現(xiàn)近似線性關(guān)系，采用二次函數(shù)便可以表示拉深筋阻力與圓角半徑間的關(guān)系。由圖5可以看出：金屬流動方向?qū)罱钭枇Φ挠绊懗尸F(xiàn)近似拋物線關(guān)系，可以采用開口向下的二次拋物線函數(shù)表示拉深筋阻力與金屬流動方向間的關(guān)系。
綜上分析，為了進(jìn)一步研究圓角半徑與金屬流動方向的耦合效應(yīng)對拉深筋阻力的影響，可用一個包含交叉項的二元二次數(shù)學(xué)模型來描述：
式中：Ff為單位拉深筋阻力；a0～a5為系數(shù)。
本項目組張萬才［9］曾經(jīng)采用獨立編寫的遺傳算法進(jìn)行曲線擬合，得到拉深筋阻力的近似解析模型，但是遺傳算法不僅編程繁瑣而且不是很穩(wěn)定，本文通過Excel進(jìn)行二元非線性回歸分析得到拉深筋阻力的近似解析模型，此方法不僅分析過程簡單而且分析精度較高且穩(wěn)定。利用Excel進(jìn)行多元非線性回歸分析需要將非線性問題轉(zhuǎn)化為線性問題才能進(jìn)行［10］，即將二元二次非線性回歸問題轉(zhuǎn)化為五元一次線性回歸問題，轉(zhuǎn)化后的回歸方程如下：
表2為線性化后的數(shù)據(jù)格式，表3為線性回歸分析的方差分析結(jié)果，F(xiàn) 值為36.8963，其顯著性水平為0.0002，小于0.05，表明回歸方程的線性關(guān)系顯著。當(dāng)檢驗整個回歸效果顯著時，還需檢驗每個變量xi對y 有無顯著的線性影響［11］。表4為回歸系數(shù)的檢驗結(jié)果，其中回歸系數(shù)a2和a4的顯著水平大于0.05，回歸系數(shù)不顯著，優(yōu)先剔除回歸方程的高次項x4（R2），可得四元一次線性回歸方程如式（3）所示：
改進(jìn)后進(jìn)行回歸分析，分析結(jié)果如表5和表6所示。
表2 線性化后的數(shù)據(jù)格式Table 2 Data format after linearization
表3 線性回歸分析的方差分析表Table 3 Analysis of variance of the linear regression analysis
表4 回歸方程回歸系數(shù)的檢驗Table 4 Testing of regression coefficient of regression equation
表5 改進(jìn)后線性回歸分析的方差分析表Table 5 Analysis of variance of the linear regression analysis of the improved equation
表6 改進(jìn)后回歸方程回歸系數(shù)的檢驗Table 6 Testing of regression coefficient of the improved regression equation
通過表6可得到四元一次線性回歸方程為：
將其還原回去即得二元二次非線性回歸方程為：
公式（5）即為考慮圓角半徑及金屬流動方向兩個因素的拉深筋阻力的近似解析模型。為了說明基于回歸分析方法得到的解析模型的準(zhǔn)確程度，將試驗數(shù)據(jù)與分析結(jié)果進(jìn)行對比，如表7 所示?？梢钥闯?，各組殘差很小，只有兩個百分比大于2%，但是也不超過5%，并且絕大多數(shù)都小于1%，可見回歸分析結(jié)果與試驗數(shù)據(jù)基本一致，即該近似解析模型是可以接受的。
結(jié)合解析模型進(jìn)一步分析圓角半徑及金屬流動方向?qū)罱钭枇Φ挠绊?。其中，式?）中的一次項θ（t＝5.2156）的系數(shù)為正，二次項θ2（t＝－6.6121）系數(shù)為負(fù)，說明單獨增大金屬流動方向與拉深筋法向的夾角可以提高拉深筋阻力，但是夾角過大時拉深筋阻力反而會減??；一次項R（t＝－3.7793）的系數(shù)為負(fù)，說明單獨增大拉深筋圓角半徑會減小拉深筋阻力；一次項θ系數(shù)的｜t｜＝5.2156大于一次項R系數(shù)的｜t｜＝3.7793，說明夾角對拉深筋阻力的影響要比圓角半徑顯著；交叉項θR（t＝－2.6258）的系數(shù)為－0.0701，說明夾角和圓角半徑的耦合效應(yīng)較小。
表7 試驗數(shù)據(jù)與分析結(jié)果對比Table 7 Contrast of experimental data and regression analysis result
圖6為利用Matlab繪制的拉深筋阻力近似解析模型的三維曲面圖，可以清晰看出圓角半徑與金屬流動方向?qū)罱钭枇Φ挠绊戁厔荨?/p>
圖6 回歸分析結(jié)果三維圖Fig.6 3Dgraph of result of regression analysis
4 結(jié) 論
（1）拉深筋圓角半徑對拉深筋阻力的影響呈近似的線性關(guān)系，單獨增大拉深筋圓角半徑會減小拉深筋阻力。
（2）金屬流動方向?qū)罱钭枇Φ挠绊懗式频膾佄锞€關(guān)系，單獨增大金屬流動方向與拉深筋法向的夾角可以提高拉深筋阻力，但是夾角過大時拉深筋阻力反而會減小。
（3）圓角半徑和金屬流動方向與拉深筋法向的夾角兩因素相比，夾角對拉深筋阻力的影響更顯著。
（4）金屬流動方向與拉深筋法向的夾角和圓角半徑的耦合效應(yīng)較小。
［1］邢宇.拉深筋的力能特性的實驗研究［D］.秦皇島：燕山大學(xué)機(jī)械工程學(xué)院，2006.Xing Yu.Experimental research on the mechanics characters of drawbead［D］.Qinhuangdao：College of Mechanical Engineering，Yanshan University，2006.
［2］Sun Guang－yong，Li Guang－yao.Multiobjective robust optimization method for drawbead design in sheet metal forming［J］.Materials and Design，2010，31（4）：1917－1929.
［3］余海燕，陳關(guān)龍，張衛(wèi)剛，等.板料成型技術(shù)中拉深筋的研究進(jìn)展［J］.塑性工程學(xué)報，2004，11（3）：77－81.Yu Hai－yan，Chen Guan－long，Zhang Wei－gang，et al.Development of drawbead investigation in sheet metal forming technology［J］.Journal of Plasticity Engineering，2004，11（3）：77－81.
［4］Halkaci H S，Turkoz M，Dilmec M，et al.Enhancing formability in hydromechanical deep drawing process adding a shallow drawbead to the blank holder［J］.Journal of Materials Processing Technology，2014，214（8）：1638－1646.
［5］Demeri M Y.Drawbeads in sheet metal forming［J］.Journal of Materials Engineering and Performance，1993，2（6）：863－866.
［6］邢忠文，楊玉英，郭剛.薄板沖壓成形中的拉深筋阻力及其影響因素研究［J］.模具工業(yè)，1994，158（4）：30－33.Xing Zhong－wen，Yang Yu－ying，Guo Gang.Research on drawbead restraining force and its influence factors in sheet metal forming［J］.Die ＆Mould Industry，1994，158（4）：30－33.
［7］Nine H D.The applicability of coulomb＇s friction law to drawbeads in sheet metal forming［J］.Journal of Applied Metal Working，1982，2（3）：200－210.
［8］Yang J H，Chen J，He D N，et al.Experimental testing of drawbead restraining force in sheet metal forming［J］.Acta Metallurgica Sinica，2003，16（1）：46－50.
［9］張萬才.基于實驗擬合的拉深筋的等效模型研究及其在一步逆算中的應(yīng)用［D］.長春：吉林大學(xué)汽車工程學(xué)院，2012.Zhang Wan－cai.The research of the equivalent drawbead model based on the fitting of experim－ental data and its application in the ore step inverse algorithm［D］.Changchun：College of Automobile Engineering，Jilin University，2012.
［10］龔江，石培春，李春燕.巧用Excel解決多元非線性回歸分析［J］.農(nóng)業(yè)網(wǎng)絡(luò)信息，2011（1）：46－48.Gong Jiang，Shi Pei－chun，Li Chun－yan.Application of excel software in multi－nonlinear regress analysis［J］.Agriculture Network Information，2011（1）：46－48.
［11］劉建波，魏宗康.石英加速度計誤差系數(shù)顯著性分析［J］.中國慣性技術(shù)學(xué)報，2011，19（5）：615－620.Liu Jian－bo，Wei Zong－kang.Significance analysis of QFPA＇s error model＇s coefficients［J］.Journal of Chinese Inertial Technology，2011，19（5）：615－620.
猜你喜歡
CAD直接建模中圓角特征識別技術(shù)的研究與應(yīng)用電子技術(shù)與軟件工程(2021年8期)2021-06-16 06:30:58 采用直線回歸方程預(yù)測桑癭蚊防治適期科學(xué)與財富(2021年36期)2021-05-10 04:54:37 線性回歸方程的求解與應(yīng)用中學(xué)生數(shù)理化·高一版(2021年2期)2021-03-19 08:32:06 線性回歸方程要點導(dǎo)學(xué)中學(xué)生數(shù)理化·高一版(2021年2期)2021-03-19 08:32:02 探究鐘表上的夾角中學(xué)生數(shù)理化·七年級數(shù)學(xué)人教版(2020年12期)2021-01-18 06:57:42 小圓角棱線車門外板工藝設(shè)計裝備制造技術(shù)(2020年1期)2020-12-25 05:18:50 求解異面直線夾角問題的兩個路徑語數(shù)外學(xué)習(xí)·高中版上旬(2020年8期)2020-09-10 07:22:44 連續(xù)展成磨削小半徑齒頂圓角的多刀逼近法制造技術(shù)與機(jī)床(2019年6期)2019-06-25 10:17:18 走進(jìn)回歸分析，讓回歸方程不再是你高考的絆腳石中學(xué)生數(shù)理化(高中版.高二數(shù)學(xué))(2019年6期)2019-06-24 03:37:52 任意夾角交叉封閉邊界內(nèi)平面流線計算及應(yīng)用西南石油大學(xué)學(xué)報(自然科學(xué)版)(2018年4期)2018-08-02 05:42:38
那景新，白 霜，劉海鵬，閆亞坤
（吉林大學(xué) 汽車仿真與控制國家重點實驗室，長春130022）
0 引 言
近年來隨著沖壓工藝的發(fā)展，沖壓制件日趨復(fù)雜，作為拉深工藝重要控制手段之一的拉深筋逐步成為拉深模具設(shè)計中的一個關(guān)鍵環(huán)節(jié)［1］。在汽車覆蓋件成形過程中，零件各部分成形條件不同，因此各部分成形所需的成形力不同，通常需要采用拉深筋來進(jìn)行控制，拉深筋參數(shù)的合理取值是控制金屬流動、防止出現(xiàn)起皺和破裂的重要手段［2－3］。
拉深筋最直接的作用是增加板料成形時的阻力［4］。拉深筋阻力有許多影響因素，國內(nèi)外都進(jìn)行了相關(guān)研究。例如，Demeri［5］通過試驗研究表明，無論何種材料，拉深筋高度和壓邊力的增加都能引起拉深筋阻力的增加；邢忠文等［6］系統(tǒng)地研究了拉深筋形式、參數(shù)、材料特性、潤滑條件、變形速度以及壓邊力等對拉深筋阻力的影響。然而，在沖壓過程中，而關(guān)于金屬流動方向并不都與拉深筋法向平行，金屬流動方向與拉深筋法向的夾角這一因素對拉深筋阻力的影響的研究卻很少。本文結(jié)合試驗重點研究金屬流動方向?qū)罱钭枇Φ挠绊?，并且同時考慮拉深筋圓角半徑這一影響因素。
本文專門設(shè)計了一套測量拉深筋阻力的試驗裝置，分別測取了以圓角半徑和金屬流動方向為變量的試驗數(shù)據(jù)，分析得出圓角半徑及金屬流動方向?qū)罱钭枇Φ挠绊懸?guī)律：圓角半徑對拉深筋阻力的影響呈現(xiàn)近似線性關(guān)系；金屬流動方向?qū)罱钭枇Φ挠绊懗尸F(xiàn)近似的拋物線關(guān)系。綜合兩個參數(shù)的影響規(guī)律并且基于回歸分析方法建立了拉深筋阻力的近似解析模型，并且驗證了回歸分析結(jié)果與試驗數(shù)據(jù)的一致性。最后對回歸系數(shù)進(jìn)行T 檢驗得出，相比于圓角半徑，金屬流動方向?qū)罱钭枇Φ挠绊懜@著，并且金屬流動方向和圓角半徑的耦合效應(yīng)較小，為拉深筋的進(jìn)一步研究奠定了基礎(chǔ)。
1 試驗裝置及方案
本文設(shè)計的試驗裝置及拉伸試驗機(jī)（測力裝置）的實物圖如圖1所示，試驗裝置的CATIA 效果圖如圖2所示。
圖1 拉深筋試驗裝置及拉伸試驗機(jī)實物Fig.1 Entity of experiment device of drawbead and tensile testing machine
圖2 拉深筋試驗裝置效果圖Fig.2 CATIA picture of experiment device of drawbead
根據(jù)實際沖壓需要，限定圓角半徑R 和金屬流動方向與拉深筋法向的夾角θ 的取值范圍分別為R＝4～16mm，θ＝0°～30°，基于差值原理，每個變量選 取3、4 個數(shù)值，其中R 取6、8、10、12 mm，θ 取0°、15°、30°。試驗所用的板料尺寸為190mm×40mm×1mm，材料為st14鋼板，性能如下：彈性模量為2.07×105MPa，屈服強(qiáng)度為210MPa，抗拉強(qiáng)度為270 MPa，泊松比為0.28，伸長率為0.36。
為了減少試驗過程中的裝卸操作工作量，試驗按照如下方式分組進(jìn)行：根據(jù)拉深筋圓角半徑R 將試驗分成四大組，每大組再按照金屬流動方向與拉深筋法向的夾角θ分成三小組，每小組安排3個試件，即每組參數(shù)進(jìn)行3次試驗，通過采用多次試驗取平均值的方法來減小試驗誤差所帶來的影響。
在本試驗裝置中，將水平拉深筋的中點作為圓心，通過將拉深筋扭轉(zhuǎn)到不同角度來模擬不同的金屬流動方向，從而可在不改變試件運動方向情況下形成金屬流動方向與拉深筋法向不同的夾角，如圖3所示，節(jié)約了試驗成本。
圖3 拉深筋角度布置Fig.3 Layout of drawbead angle
2 圓角半徑及金屬流動方向?qū)罱钭枇Φ挠绊?/h2>
由于試驗得到的拉深筋阻力是沿著板料寬度方向（40mm）的總阻力，將其除以板料寬度，即可得單位長度的拉深筋阻力。結(jié)合試驗數(shù)據(jù)，如表1所示，分析拉深筋的圓角半徑及金屬流動方向?qū)罱钭枇Φ挠绊懸?guī)律。
表1 拉深筋阻力試驗數(shù)據(jù)Table 1 Experimental data of drawbead restraining force
2.1 拉深筋圓角半徑對拉深筋阻力的影響
根據(jù)三組不同金屬流動方向下圓角半徑分別為6、8、10、12mm 時所測得的拉深筋阻力的試驗數(shù)據(jù)，分別繪制了金屬流動方向分別為0°、15°、30°時，以圓角半徑為橫坐標(biāo)，單位拉深筋阻力為縱坐標(biāo)的關(guān)系曲線，如圖4所示。
拉深筋阻力主要是由板料通過拉深筋時，發(fā)生彎曲和反彎曲及其接觸面的摩擦而形成的［7－8］。由圖4可以看出：當(dāng)金屬流動方向與拉深筋法向的夾角一定時，隨著半徑的增大，拉深筋阻力減小，這主要是由于拉深筋在金屬流動方向的截面的圓弧曲率減小，從而導(dǎo)致彎曲和反彎曲的阻力減小所引起的。圓角半徑對拉深筋阻力的影響呈現(xiàn)近似線性關(guān)系，并且影響程度相差不大。
圖4 拉深筋阻力與拉深筋圓角半徑的關(guān)系Fig.4 Relationships between drawbead restraining force and corner radius
2.2 金屬流動方向?qū)罱钭枇Φ挠绊?/h3>
為進(jìn)一步研究金屬流動方向?qū)罱钭枇Φ挠绊懸?guī)律，本文又分別依據(jù)四組不同圓角半徑下金屬流動方向與拉深筋法向夾角分別為0°、15°、30°時所測得的拉深筋阻力的試驗數(shù)據(jù)，分別繪制了圓角半徑分別為6、8、10、12 mm 時，以夾角為橫坐標(biāo)，單位拉深筋阻力為縱坐標(biāo)的關(guān)系曲線，如圖5所示。
圖5 拉深筋阻力與金屬流動方向的關(guān)系Fig.5 Relationships between drawbead restraining force and metal flow direction
圖5表明，金屬流動方向?qū)罱钭枇Φ挠绊懗尸F(xiàn)近似的拋物線關(guān)系，且在夾角為15°時達(dá)到最大值。這主要是因為在金屬流動方向與拉深筋法向的夾角不同時，與金屬流動方向平行的拉深筋的截面的圓弧曲率的不同引起的彎曲和反彎曲的阻力變化幅度不同，而且板料與拉深筋的接觸面積的不同引起的摩擦力變化幅度也不同，拉深筋阻力大小的改變是這兩個主要因素的綜合影響結(jié)果，并且不同夾角時的影響程度不同。
例如，當(dāng)夾角為0°～15°時，夾角的增大導(dǎo)致板料與拉深筋的接觸面積不斷增大，從而導(dǎo)致摩擦力增大；金屬流經(jīng)拉深筋時，與金屬流動方向平行的拉深筋的截面為橢圓形，夾角的增大導(dǎo)致橢圓的曲率減?。ㄇ拾霃皆龃螅?，那么由彎曲反彎曲所產(chǎn)生的阻力也會減小。
綜上并結(jié)合圖5分析可知，當(dāng)夾角為0°～15°時，金屬流經(jīng)拉深筋和筋槽產(chǎn)生的摩擦力增加的幅度大于由彎曲反彎曲產(chǎn)生的阻力減小的幅度，導(dǎo)致拉深筋阻力提高，并且提高幅度較小，表現(xiàn)為曲線上升趨勢緩慢。當(dāng)圓角半徑為12mm 時，上升幅度最小，為1.25N；當(dāng)圓角半徑為6 mm 和10mm 時，上升幅度最大，也僅為7.5N。當(dāng)夾角為15°～30°時，彎曲反彎曲產(chǎn)生阻力減小的幅度大于摩擦力增加的幅度，導(dǎo)致拉深筋阻力減小，并且減小幅度較大，表現(xiàn)為曲線下降趨勢明顯，尤其是當(dāng)拉深筋圓角半徑為10 mm 和12 mm 時，下降幅度分別為21.666N 和16.667N。
3 拉深筋阻力近似解析模型的建立
由圖4可以看出：圓角半徑對拉深筋阻力的影響呈現(xiàn)近似線性關(guān)系，采用二次函數(shù)便可以表示拉深筋阻力與圓角半徑間的關(guān)系。由圖5可以看出：金屬流動方向?qū)罱钭枇Φ挠绊懗尸F(xiàn)近似拋物線關(guān)系，可以采用開口向下的二次拋物線函數(shù)表示拉深筋阻力與金屬流動方向間的關(guān)系。
綜上分析，為了進(jìn)一步研究圓角半徑與金屬流動方向的耦合效應(yīng)對拉深筋阻力的影響，可用一個包含交叉項的二元二次數(shù)學(xué)模型來描述：
式中：Ff為單位拉深筋阻力；a0～a5為系數(shù)。
本項目組張萬才［9］曾經(jīng)采用獨立編寫的遺傳算法進(jìn)行曲線擬合，得到拉深筋阻力的近似解析模型，但是遺傳算法不僅編程繁瑣而且不是很穩(wěn)定，本文通過Excel進(jìn)行二元非線性回歸分析得到拉深筋阻力的近似解析模型，此方法不僅分析過程簡單而且分析精度較高且穩(wěn)定。利用Excel進(jìn)行多元非線性回歸分析需要將非線性問題轉(zhuǎn)化為線性問題才能進(jìn)行［10］，即將二元二次非線性回歸問題轉(zhuǎn)化為五元一次線性回歸問題，轉(zhuǎn)化后的回歸方程如下：
表2為線性化后的數(shù)據(jù)格式，表3為線性回歸分析的方差分析結(jié)果，F(xiàn) 值為36.8963，其顯著性水平為0.0002，小于0.05，表明回歸方程的線性關(guān)系顯著。當(dāng)檢驗整個回歸效果顯著時，還需檢驗每個變量xi對y 有無顯著的線性影響［11］。表4為回歸系數(shù)的檢驗結(jié)果，其中回歸系數(shù)a2和a4的顯著水平大于0.05，回歸系數(shù)不顯著，優(yōu)先剔除回歸方程的高次項x4（R2），可得四元一次線性回歸方程如式（3）所示：
改進(jìn)后進(jìn)行回歸分析，分析結(jié)果如表5和表6所示。
表2 線性化后的數(shù)據(jù)格式Table 2 Data format after linearization
表3 線性回歸分析的方差分析表Table 3 Analysis of variance of the linear regression analysis
表4 回歸方程回歸系數(shù)的檢驗Table 4 Testing of regression coefficient of regression equation
表5 改進(jìn)后線性回歸分析的方差分析表Table 5 Analysis of variance of the linear regression analysis of the improved equation
表6 改進(jìn)后回歸方程回歸系數(shù)的檢驗Table 6 Testing of regression coefficient of the improved regression equation
通過表6可得到四元一次線性回歸方程為：
將其還原回去即得二元二次非線性回歸方程為：
公式（5）即為考慮圓角半徑及金屬流動方向兩個因素的拉深筋阻力的近似解析模型。為了說明基于回歸分析方法得到的解析模型的準(zhǔn)確程度，將試驗數(shù)據(jù)與分析結(jié)果進(jìn)行對比，如表7 所示?？梢钥闯?，各組殘差很小，只有兩個百分比大于2%，但是也不超過5%，并且絕大多數(shù)都小于1%，可見回歸分析結(jié)果與試驗數(shù)據(jù)基本一致，即該近似解析模型是可以接受的。
結(jié)合解析模型進(jìn)一步分析圓角半徑及金屬流動方向?qū)罱钭枇Φ挠绊?。其中，式?）中的一次項θ（t＝5.2156）的系數(shù)為正，二次項θ2（t＝－6.6121）系數(shù)為負(fù)，說明單獨增大金屬流動方向與拉深筋法向的夾角可以提高拉深筋阻力，但是夾角過大時拉深筋阻力反而會減??；一次項R（t＝－3.7793）的系數(shù)為負(fù)，說明單獨增大拉深筋圓角半徑會減小拉深筋阻力；一次項θ系數(shù)的｜t｜＝5.2156大于一次項R系數(shù)的｜t｜＝3.7793，說明夾角對拉深筋阻力的影響要比圓角半徑顯著；交叉項θR（t＝－2.6258）的系數(shù)為－0.0701，說明夾角和圓角半徑的耦合效應(yīng)較小。
表7 試驗數(shù)據(jù)與分析結(jié)果對比Table 7 Contrast of experimental data and regression analysis result
圖6為利用Matlab繪制的拉深筋阻力近似解析模型的三維曲面圖，可以清晰看出圓角半徑與金屬流動方向?qū)罱钭枇Φ挠绊戁厔荨?/p>
圖6 回歸分析結(jié)果三維圖Fig.6 3Dgraph of result of regression analysis
4 結(jié) 論
（1）拉深筋圓角半徑對拉深筋阻力的影響呈近似的線性關(guān)系，單獨增大拉深筋圓角半徑會減小拉深筋阻力。
（2）金屬流動方向?qū)罱钭枇Φ挠绊懗式频膾佄锞€關(guān)系，單獨增大金屬流動方向與拉深筋法向的夾角可以提高拉深筋阻力，但是夾角過大時拉深筋阻力反而會減小。
（3）圓角半徑和金屬流動方向與拉深筋法向的夾角兩因素相比，夾角對拉深筋阻力的影響更顯著。
（4）金屬流動方向與拉深筋法向的夾角和圓角半徑的耦合效應(yīng)較小。
［1］邢宇.拉深筋的力能特性的實驗研究［D］.秦皇島：燕山大學(xué)機(jī)械工程學(xué)院，2006.Xing Yu.Experimental research on the mechanics characters of drawbead［D］.Qinhuangdao：College of Mechanical Engineering，Yanshan University，2006.
［2］Sun Guang－yong，Li Guang－yao.Multiobjective robust optimization method for drawbead design in sheet metal forming［J］.Materials and Design，2010，31（4）：1917－1929.
［3］余海燕，陳關(guān)龍，張衛(wèi)剛，等.板料成型技術(shù)中拉深筋的研究進(jìn)展［J］.塑性工程學(xué)報，2004，11（3）：77－81.Yu Hai－yan，Chen Guan－long，Zhang Wei－gang，et al.Development of drawbead investigation in sheet metal forming technology［J］.Journal of Plasticity Engineering，2004，11（3）：77－81.
［4］Halkaci H S，Turkoz M，Dilmec M，et al.Enhancing formability in hydromechanical deep drawing process adding a shallow drawbead to the blank holder［J］.Journal of Materials Processing Technology，2014，214（8）：1638－1646.
［5］Demeri M Y.Drawbeads in sheet metal forming［J］.Journal of Materials Engineering and Performance，1993，2（6）：863－866.
［6］邢忠文，楊玉英，郭剛.薄板沖壓成形中的拉深筋阻力及其影響因素研究［J］.模具工業(yè)，1994，158（4）：30－33.Xing Zhong－wen，Yang Yu－ying，Guo Gang.Research on drawbead restraining force and its influence factors in sheet metal forming［J］.Die ＆Mould Industry，1994，158（4）：30－33.
［7］Nine H D.The applicability of coulomb＇s friction law to drawbeads in sheet metal forming［J］.Journal of Applied Metal Working，1982，2（3）：200－210.
［8］Yang J H，Chen J，He D N，et al.Experimental testing of drawbead restraining force in sheet metal forming［J］.Acta Metallurgica Sinica，2003，16（1）：46－50.
［9］張萬才.基于實驗擬合的拉深筋的等效模型研究及其在一步逆算中的應(yīng)用［D］.長春：吉林大學(xué)汽車工程學(xué)院，2012.Zhang Wan－cai.The research of the equivalent drawbead model based on the fitting of experim－ental data and its application in the ore step inverse algorithm［D］.Changchun：College of Automobile Engineering，Jilin University，2012.
［10］龔江，石培春，李春燕.巧用Excel解決多元非線性回歸分析［J］.農(nóng)業(yè)網(wǎng)絡(luò)信息，2011（1）：46－48.Gong Jiang，Shi Pei－chun，Li Chun－yan.Application of excel software in multi－nonlinear regress analysis［J］.Agriculture Network Information，2011（1）：46－48.
［11］劉建波，魏宗康.石英加速度計誤差系數(shù)顯著性分析［J］.中國慣性技術(shù)學(xué)報，2011，19（5）：615－620.Liu Jian－bo，Wei Zong－kang.Significance analysis of QFPA＇s error model＇s coefficients［J］.Journal of Chinese Inertial Technology，2011，19（5）：615－620.
猜你喜歡
電子技術(shù)與軟件工程(2021年8期)2021-06-16 06:30:58
科學(xué)與財富(2021年36期)2021-05-10 04:54:37
中學(xué)生數(shù)理化·高一版(2021年2期)2021-03-19 08:32:06
中學(xué)生數(shù)理化·高一版(2021年2期)2021-03-19 08:32:02
中學(xué)生數(shù)理化·七年級數(shù)學(xué)人教版(2020年12期)2021-01-18 06:57:42
裝備制造技術(shù)(2020年1期)2020-12-25 05:18:50
語數(shù)外學(xué)習(xí)·高中版上旬(2020年8期)2020-09-10 07:22:44
制造技術(shù)與機(jī)床(2019年6期)2019-06-25 10:17:18
中學(xué)生數(shù)理化(高中版.高二數(shù)學(xué))(2019年6期)2019-06-24 03:37:52
西南石油大學(xué)學(xué)報(自然科學(xué)版)(2018年4期)2018-08-02 05:42:38
