張子瀟，孔德仁，王　芳，楊麗俠

(1. 南京理工大學(xué) 機(jī)械工程學(xué)院，江蘇 南京 210094； 2. 中國(guó)兵器工業(yè)204所，陜西 西安 710065)

?

基于準(zhǔn)動(dòng)態(tài)校準(zhǔn)的Φ3.5×8.75 mm柱形銅柱力學(xué)性能仿真分析*

張子瀟1，孔德仁1，王芳1，楊麗俠2

(1. 南京理工大學(xué) 機(jī)械工程學(xué)院，江蘇 南京 210094； 2. 中國(guó)兵器工業(yè)204所，陜西 西安 710065)

摘要:為了改善銅柱測(cè)壓中靜標(biāo)體系的缺陷，我國(guó)曾設(shè)計(jì)了采用準(zhǔn)動(dòng)態(tài)校準(zhǔn)的Φ3.5×8.75 mm銅柱測(cè)壓器材，而在實(shí)際使用過(guò)程中該型銅柱測(cè)壓系統(tǒng)在高壓測(cè)試中出現(xiàn)了一些問(wèn)題，其結(jié)構(gòu)設(shè)計(jì)需要進(jìn)行優(yōu)化，這必然涉及到對(duì)該型銅柱測(cè)壓范圍的界定與預(yù)壓值的選擇. 應(yīng)用非線性動(dòng)力學(xué)軟件LS-DYNA對(duì)銅柱測(cè)壓系統(tǒng)在不同預(yù)壓與溫度下的力學(xué)性能行了仿真分析，得到了預(yù)壓與溫度對(duì)銅柱變形量的影響，研究成果為Φ3.5×8.75 mm柱形銅柱測(cè)壓器的優(yōu)化提供了依據(jù)以及參考.

關(guān)鍵詞:塑性測(cè)壓； 準(zhǔn)動(dòng)態(tài)校準(zhǔn)； 力學(xué)性能分析； 壓力測(cè)量； 有限元分析

塑性測(cè)壓法與電測(cè)法是測(cè)試火炮等常規(guī)武器膛內(nèi)火藥燃?xì)鈮毫?下稱膛壓)的主要手段，其中塑性測(cè)壓法因?yàn)槠洳僮骱?jiǎn)單、測(cè)試結(jié)果穩(wěn)定、用于大口徑武器時(shí)無(wú)需開測(cè)壓孔等優(yōu)點(diǎn)，目前仍然廣泛應(yīng)用于武器驗(yàn)收等工作中. 塑性法根據(jù)塑性敏感元件不同可分為銅柱測(cè)壓法與銅球測(cè)壓法，目前我國(guó)的銅柱測(cè)壓器材主要還在使用引進(jìn)自前蘇聯(lián)的靜態(tài)標(biāo)定體系，由于銅材存在應(yīng)變率效應(yīng)，在同樣壓力情況下靜壓變形量大于動(dòng)壓變形，這就會(huì)引入相應(yīng)的系統(tǒng)誤差，即所謂“靜動(dòng)差”[1]. 此外，靜標(biāo)銅柱因覆蓋壓力范圍窄所以器材規(guī)格繁多，且不同規(guī)格銅柱對(duì)同一壓力源的測(cè)壓結(jié)果無(wú)可比性，這對(duì)于塑性測(cè)壓器材的生產(chǎn)、管理與使用都是極為不利的[2].

為了消除或減少測(cè)壓銅柱的靜動(dòng)差，國(guó)外從20世紀(jì)70年代開始，對(duì)塑性測(cè)壓的標(biāo)定體系進(jìn)行重大改革，即用具有與膛壓上升沿類似動(dòng)態(tài)特性的半正弦壓力脈沖對(duì)塑性測(cè)壓器材進(jìn)行準(zhǔn)動(dòng)態(tài)標(biāo)定，這也是塑性測(cè)壓法的發(fā)展趨勢(shì)[3-4].

針對(duì)我國(guó)靜態(tài)標(biāo)定體系中存在的靜動(dòng)差與器材規(guī)格過(guò)多等問(wèn)題，我國(guó)曾經(jīng)設(shè)計(jì)了Φ3.5×8.75 mm系列準(zhǔn)動(dòng)態(tài)校準(zhǔn)銅柱測(cè)壓器材. 該系列銅柱測(cè)壓器件采用的是一次統(tǒng)一預(yù)壓的準(zhǔn)動(dòng)態(tài)校準(zhǔn)，較大程度上減少了銅柱測(cè)壓中靜動(dòng)差，其測(cè)量值與電測(cè)值更為一致，大大提高了銅柱測(cè)壓的精度. 但是這一規(guī)格器材在高溫高壓測(cè)試中會(huì)出現(xiàn)測(cè)壓油滲漏導(dǎo)致銅柱壓偏的情況，而且低溫高壓下也存在著活塞組件卡滯的問(wèn)題； 因此需要對(duì)Φ3.5×8.75 mm系列銅柱測(cè)壓器進(jìn)行優(yōu)化設(shè)計(jì)，這必然涉及到對(duì)Φ3.5×8.75 mm測(cè)壓銅柱的測(cè)壓范圍進(jìn)行合理界定與統(tǒng)一預(yù)壓值的選擇問(wèn)題. 目前初步確定Φ3.5×8.75mm系列動(dòng)標(biāo)銅柱的測(cè)壓范圍為20～700 MPa. 因?yàn)殂~柱測(cè)壓法中通常用錐形銅柱進(jìn)行低壓測(cè)試，柱形銅柱則用于高壓測(cè)試，為了對(duì)測(cè)壓器設(shè)計(jì)進(jìn)行優(yōu)化需要界定錐形銅柱與柱形銅柱測(cè)壓范圍的界限，并據(jù)此為柱形銅柱選擇合理的預(yù)壓值.

1Φ3.5×8.75 mm銅柱測(cè)壓的工作原理

圖1　銅柱測(cè)壓器結(jié)構(gòu)示意圖Fig.1　The diagram of copper crusher gauge

銅柱測(cè)壓器基本工作原理如圖1 所示，測(cè)壓器在使用時(shí)置于火炮炮膛中，被測(cè)武器膛壓p(t)經(jīng)測(cè)壓油作用于活塞桿工作面上，通過(guò)活塞桿將壓力施加于銅柱上，引起銅柱軸向壓縮變形. 由于銅柱材料屈服極限較低，在壓力作用下發(fā)生永久塑性應(yīng)變，在膛壓作用完全結(jié)束后銅柱依然存在有一定殘余變形. 由于銅柱在測(cè)試中的變形過(guò)程無(wú)法測(cè)量，所以采殘余變形量(后稱變形量)的大小Y用來(lái)度量膛壓p(t)的峰值大小pm[5].

為了消除銅柱個(gè)體間差異并且縮小銅柱在工作中的變形量以減少測(cè)量中的靜動(dòng)差，測(cè)壓銅柱在使用前需要進(jìn)行預(yù)壓. 對(duì)于Φ3.5×8.75 mm 柱形銅柱采用的是一次統(tǒng)一預(yù)壓，預(yù)壓后進(jìn)行準(zhǔn)動(dòng)態(tài)標(biāo)定并編制壓力表. 進(jìn)行過(guò)預(yù)壓后銅柱的高度稱為壓前高，完成膛壓測(cè)試后銅柱高則稱為壓后高，用壓前高減去壓后高即為銅柱變形量.

2有限元模型的建立與計(jì)算

本文通過(guò)有限元仿真對(duì)預(yù)壓值變化對(duì)銅柱變形規(guī)律的影響以及高低溫情況下銅柱的變形情況進(jìn)行了分析，為預(yù)壓值的選擇與測(cè)壓器的設(shè)計(jì)提供了參考. 由于塑性測(cè)壓法涉及到瞬態(tài)以及非線性問(wèn)題，所以需要使用非線性動(dòng)力學(xué)分析軟件，本文中使用的軟件為L(zhǎng)S-DYNA.

2.1仿真邊界條件

在實(shí)際測(cè)試中，銅柱一端與測(cè)壓器活塞接觸，另一端與測(cè)壓器砧體接觸. 由于測(cè)壓器所用40CrNiMoA鋼材經(jīng)過(guò)加工后剛度遠(yuǎn)遠(yuǎn)大于銅柱，所以在計(jì)算中可視為剛體. 而且測(cè)壓器砧體面積大于銅柱，故銅柱與砧體的接觸可認(rèn)看作對(duì)銅柱一端施加一個(gè)軸向約束.

因?yàn)樘艍航?jīng)過(guò)測(cè)壓油傳播至活塞桿表面上時(shí)可近似看作均勻分布的壓力，所以在加載時(shí)視作一個(gè)均布?jí)毫d荷.

2.2有限元模型建立

由于外殼等其他部件在正常工作時(shí)對(duì)銅柱變形影響很小，所以在本文的分析中僅關(guān)注銅柱與活塞桿. 由于這兩個(gè)部件均為旋轉(zhuǎn)體，所以從節(jié)省計(jì)算資源與方便網(wǎng)格劃分的角度上，將模型簡(jiǎn)化為1/4后視為對(duì)稱結(jié)構(gòu)進(jìn)行計(jì)算. 為了減少計(jì)算中的沙漏效應(yīng)與局部應(yīng)力異常，用幾何模型劃分為大小合適且形狀規(guī)則的體單元，得到有限元模型. 模型中活塞桿有1 385個(gè)單元、1 968個(gè)節(jié)點(diǎn)，銅柱有1 188個(gè)單元、1 628個(gè)節(jié)點(diǎn).

圖2　柱形銅柱及活塞桿的幾何模型以與有限元模型Fig.2　The geometric and finite element model of copper cylinder and piston rod

2.3銅柱本構(gòu)關(guān)系的選擇

測(cè)壓銅柱在使用前先要經(jīng)過(guò)一次靜態(tài)預(yù)壓，在受到膛壓變形時(shí)則是一個(gè)動(dòng)態(tài)過(guò)程； 在這兩個(gè)變形階段，銅柱應(yīng)當(dāng)產(chǎn)生彈性變形與塑性變形. 所以對(duì)銅柱需要使用能夠反映材料在靜態(tài)與動(dòng)態(tài)兩種情況下發(fā)生非線性變形的本構(gòu)關(guān)系，所以本文中本構(gòu)關(guān)系采用的為冪法則模型，該模型為含有Cowper-Symonds應(yīng)變率模型的Ramburgh-Osgood模型，在LS-DYNA中稱為Power-Law模型，可用于靜態(tài)與動(dòng)態(tài)的金屬塑性變形仿真[6]，不同溫度下該模型參數(shù)有所不同. 其塑性段應(yīng)力應(yīng)變關(guān)系可通過(guò)式(1)表示.

(1)

仿真中單位制為kg-mm-ms，所用材料參數(shù)如表2 所示.

表2　銅柱Power-Law本構(gòu)模型參數(shù)

2.4模擬加載設(shè)置

由于半正弦脈沖與膛壓上升沿有著類似的頻率特性，所以國(guó)內(nèi)外通常使用脈寬為6±10% ms的半正弦壓力脈沖作為準(zhǔn)動(dòng)態(tài)校準(zhǔn)的壓力源[7]，本文仿真中使用的模擬壓力載荷為脈寬6 ms的半正弦壓力脈沖[8].

由于銅柱在使用中需要進(jìn)行預(yù)壓，所以在施加膛壓載荷之前需要先施加預(yù)壓載荷，在完成預(yù)壓基礎(chǔ)上再施加膛壓載荷. 目前暫定Φ3.5×8.75 mm測(cè)壓范圍為20～700 MPa，由于較高預(yù)壓不宜測(cè)量低膛壓，所以本文中對(duì)預(yù)壓值175～185 Mpa的柱形銅柱加載峰值200～750 Mpa的半正弦壓力脈沖，對(duì)預(yù)壓值較高的銅柱則加載峰值為250～750 MPa的壓力脈沖.

圖4 中為部分半正弦壓力載荷的曲線，壓力載荷于33 ms處達(dá)到峰值.

圖3　半正弦壓力脈沖與膛壓曲線示意圖Fig.3　The diagram of half sine pressure pulse and chamber pressure curve

圖4　仿真動(dòng)態(tài)壓力載荷曲線Fig.4　The loading curves in simulation

3仿真結(jié)果與分析

觀察銅柱發(fā)生變形的時(shí)程圖，記錄變形量大小，分析不同預(yù)壓值下銅柱的變形情況.

圖5　預(yù)壓185 MPa銅柱變形曲線Fig.5　The deformation curves of copper cylinder with preload of 185 MPa

圖5 中是常溫下預(yù)壓值為185 MPa柱形銅柱在不同峰值壓力脈沖下的變形時(shí)程圖，可觀察到在壓力脈沖下降沿存在可見的彈性回復(fù)，且銅柱的變形量峰值較壓力脈沖峰值在時(shí)間上存在一定滯后，這與銅柱的實(shí)際變形情況相符.

表3 中為不同預(yù)壓值下不同壓力峰值銅柱的變形量Y值，由于高預(yù)壓不適用于低膛壓，所以不同預(yù)壓值所采用的最小壓力峰值有所區(qū)別

為了方便進(jìn)行分析，對(duì)不同預(yù)壓值的變形量進(jìn)行二次回歸，即ε=f(P)=A0+B1P+B2P2，回歸精度通過(guò)相關(guān)系數(shù)r和剩余標(biāo)準(zhǔn)偏差RMSE考核，回歸結(jié)果見表4.

從表4 中可見，銅柱變形量回歸質(zhì)量較好，其變形規(guī)律較為符合二次多項(xiàng)式.

表4　銅柱變預(yù)壓變形量回歸分析

在銅柱測(cè)壓中，測(cè)壓范圍下限一般取銅柱預(yù)壓值加15 MPa，將不同預(yù)壓值銅柱預(yù)壓后高度(壓前高)、測(cè)壓范圍與最大變形量即壓力達(dá)到測(cè)壓上限時(shí)的變形量匯總?cè)氡?，作為測(cè)壓器設(shè)計(jì)參考.

表5　銅柱預(yù)壓高、測(cè)壓范圍與最大變形量

根據(jù)表4 中回歸方程可求得壓力靈敏度K，壓力靈敏度定義為單位變形量表征的壓力，單位為MPa/mm[9]將表4中回歸方程參數(shù)帶入式(2)可求得銅柱變形的壓力靈敏度方程以及壓力靈敏度范圍. 目前測(cè)壓范圍上限暫定為700 MPa，根據(jù)所用壓力峰值范圍，對(duì)于175 MPa, 180 MPa, 185 MPa預(yù)壓值統(tǒng)計(jì)200～700 MPa靈敏度，其余預(yù)壓值則統(tǒng)計(jì)250～700 MPa范圍.

(2)

圖6　銅柱靈敏度曲線Fig.6　Sensitivity curves of copper cylinder

將柱形銅柱壓力靈敏度繪制成曲線圖可更加直觀地觀察預(yù)壓值變化對(duì)靈敏度曲線的影響，見圖6.

從圖6 線中可以看出，采用較高的預(yù)壓值能夠得到線性更好的靈敏度曲線，但是從表5中可見采用高預(yù)壓值會(huì)縮小銅柱的測(cè)壓范圍. 為此，應(yīng)當(dāng)選擇一個(gè)折衷方案，從靈敏度曲線來(lái)看，采用較低預(yù)壓值銅柱的靈敏度差異較小，且靈敏度變化在可接受范圍內(nèi)，而采用較高預(yù)壓值可得到線性更好的靈敏度曲線，綜合分析后暫定柱形銅柱的測(cè)壓范圍為200～700 MPa統(tǒng)一預(yù)壓值為185 MPa.

暫時(shí)確定測(cè)壓范圍及預(yù)壓值后，需要為測(cè)壓器設(shè)計(jì)提供尺寸參考，而銅柱在不同溫度下變形量存在差異，為此需要分析在高低溫工況下銅柱的變形情況. 表7中為不同溫度條件下，預(yù)壓值為185 MPa銅柱的變形量Y仿真數(shù)據(jù).

表7　銅柱變溫度變形量

從表7 中可以看出，高低溫下銅柱的變形量發(fā)生了一定變化，而且符合材料在溫度升高時(shí)軟化、低溫硬化的情況[10]； 對(duì)不同溫度下銅柱變形量進(jìn)行二次回歸，數(shù)據(jù)填入表8 中.

表8　柱形銅柱變溫度回歸

將變形量與預(yù)壓值的關(guān)系繪制為曲線圖，見圖7，以更為直觀地分析溫度變化的影響.

圖7　柱形銅柱高低溫變預(yù)壓變形量回歸曲線Fig.7　The regression curve of deformation amount of copper cylinder in different tempreture

從圖7 中可見，溫度變化在一定程度上影響了銅柱的變形情況： 高溫情況下變形量增大，相應(yīng)的低溫情況下變形量減少，且不同溫度曲線互不相交，可見不同溫度下銅柱有類似的變形規(guī)律. 若采用186 MPa 作為統(tǒng)一預(yù)壓值，則對(duì)測(cè)壓器進(jìn)行設(shè)計(jì)，活塞行程應(yīng)大于仿真得到的50 ℃時(shí)700 MPa變形量值3.209 mm，并且保留一定余量.

4結(jié)論

1) 通過(guò)LS-DYNA對(duì)柱形測(cè)壓銅柱采用不同預(yù)壓值時(shí)的變形規(guī)律與靈敏度進(jìn)行了仿真分析，為界定柱形銅柱測(cè)壓范圍與統(tǒng)一預(yù)壓值提供了參考.

2) 通過(guò)仿真分析暫定柱形銅柱測(cè)壓范圍為200～700 MPa，采用186 MPa的預(yù)壓值并且分析了采用該預(yù)壓值的柱形銅柱在高溫與低溫環(huán)境下的變形情況，為之后的相關(guān)實(shí)驗(yàn)設(shè)計(jì)提供了邊界條件.

參考文獻(xiàn)：

[1]伍鷹，朱明武. 測(cè)壓圓柱靜動(dòng)態(tài)特性計(jì)算方法[J]. 兵工學(xué)報(bào)，1993(4)： 7-69-74.

Wu Ying, Zhu Mingwu. Calculation of the static and dynamic characteristics of pressure-measuring cylinders [J]. ActaArmamentarii, 1993(4)： 7-69-74. (in Chinese)

[2]楊麗俠. 我國(guó)火炮膛壓塑性測(cè)壓的計(jì)量發(fā)展研究[C]. 中國(guó)航空工業(yè)集團(tuán)公司第三○四研究所. 2013年力學(xué)計(jì)量測(cè)試技術(shù)學(xué)術(shù)交流會(huì)論文專集. 南平： 北京長(zhǎng)城計(jì)量測(cè)試技術(shù)研究所，2013： 1-4.

[3]孔德仁, 朱明武. 銅柱、銅球準(zhǔn)動(dòng)態(tài)校準(zhǔn)的機(jī)理[J]. 南京理工大學(xué)學(xué)報(bào)，2004，28(4)： 375-379.

Kong Deren, Zhu Mingwu. Mechanism of quasi-dynamic calibration for copper cylinder and copper ball[J]. Journal of Nanjing University of Science and Technology，2004，28(4)： 375-379. (in Chinese)

[4]張淑寧. 準(zhǔn)動(dòng)態(tài)校準(zhǔn)技術(shù)實(shí)驗(yàn)研究[D]. 南京： 南京理工大學(xué)，2002.

[5]孔德仁，狄長(zhǎng)安，范啟勝. 塑性測(cè)壓技術(shù)[M]. 北京： 兵器工業(yè)出版社，2006.

[6]王澤鵬，胡仁喜，康士廷，等. ANSYS 13. 0/LS-DYNA非線性有限元分析[M]. 北京： 機(jī)械工業(yè)出版社，2011.

[7]朱明武. 半正弦壓力源的模型及設(shè)計(jì)[C]. 中國(guó)計(jì)量測(cè)試學(xué)會(huì). 2001全國(guó)壓力計(jì)量測(cè)試技術(shù)年會(huì)論文集. 太原： 中國(guó)計(jì)量測(cè)試學(xué)會(huì)，2001： 112-118.

[8]李波.φ8×13 mm測(cè)壓銅柱靜動(dòng)差研究[D]. 南京： 南京理工大學(xué)，2013.

[9]劉素潔. 測(cè)壓銅柱的一致性及靜動(dòng)差研究[D]. 南京： 南京理工大學(xué)，2008.

[10]孔德仁，朱明武，張淑寧，等. 測(cè)壓銅柱溫度統(tǒng)一修正公式研究[J]. 兵工學(xué)報(bào)，2003，24(2)： 145-148.

Kong Deren, Zhu Mingwu, Zhang Shuning, et al. A study on a temperature correction equation for pressure— measuring copper cylinders [J]. ActaArmamentarii, 2003, 24(2)： 145-148. (in Chinese)

Mechanical Property Simulation Analysis ofΦ3.5×8.75 mm Copper Cylinder Based on Quasi-Dynamic Calibration

ZHANG Zixiao1, KONG Deren1, WANG Fang1, YANG Lixia2

(1. Colloge of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China;2.Xi’an Modern Chemistry Institute, Xi’an 710065, China)

Abstract:The Φ3.5×8.75 mm copper cylinder pressure-measuring elements have been designed which is based on quasi-Dynamic calibration technique in order to improve the defects of the static calibration system in the pressure-measuring using copper cylinder. But there are problems with this model of copper crusher gauge in high pressure measurement. So the structure design of the copper crusher gauge should be optimized that involves the defining of the measuring range and the preloading of the copper cylinder. By the nonlinear finite element analysis software ANSYS/LS-DYNA, the mechanical property of the copper cylinder pressure measuring system in different preloading and temperature situations is simulated to analyze the influence of preloading and temperature to the deformation amount. The research results provide reference to the optimal design of Φ3.5×8.75 mm copper cylinder crusher gauges.

Key words:plastic pressure-measurement; quasi-dynamic calibration; mechanical property analysis; pressure-measurement; finite element analysis

文章編號(hào):1671-7449(2016)04-0299-07

收稿日期：2015-10-21

基金項(xiàng)目：國(guó)防計(jì)量基礎(chǔ)預(yù)研資助項(xiàng)目

作者簡(jiǎn)介：張子瀟(1992-)，男，碩士生，主要從事塑性測(cè)壓技術(shù)研究.

中圖分類號(hào):O344.3

文獻(xiàn)標(biāo)識(shí)碼:A

doi:10.3969/j.issn.1671-7449.2016.04.004