均衡亞克力通透試驗和參數摹擬
實驗方法實驗裝置如示，用高速轉鏡掃描相機和單狹縫技術來測定爆轟波透出端面的波形。平面波透鏡炸藥選用JO9159，內部采用惰性材料有機玻璃調整波形，底端覆蓋一層高能炸藥JO9159，起爆方式為單點起爆。
相機的掃描轉速為15萬轉/分鐘，狹縫像在底片上的掃描速度為7.5mms-1，中小方框示底片處理后得到的Rt實驗數據。
數值模擬狀態方程在計算模型中，對JO9159炸藥采用三項式的點平面波透鏡示意微秒級高壓電雷管，2JO9159炸藥，3有機玻璃，4反光鏡，5高速轉鏡掃描相機火增長模型，對有機玻璃采用Gruneisen方程。
三項式點火增長模型<7，8>：ddt=I（1-）b0-1-ax+G1（1-）cdpy+G2（1-）egpz（1）式中，是炸藥反應度，t是時間，是密度，0是初始密度，單位為gcm-3，p是壓力，單位為GPa，I，G1，G2，a，b，x，c，d，y，e，g和z是常數。
反應產物的狀態方程采用JWL形式的狀態方程：p=A1-wR1Ve-R1V+B1-wR2Ve-R2V+wE0V（2）式中的系數均為待定參數。
Gruneisen方程<9，10>：p=0C21+1-02-a21-（S1-1）+S2+1-S3（+1）2+（0+a）E（3）式中，C為usup曲線的截距，S1、S2、S3是usup曲線斜率的系數，0是Gruneisen系數，a是對0的一階體積修正。
JO9159炸藥化學反應速率方程的參數及爆轟產物狀態方程的參數如2所示，有機玻璃狀態方程參數如所示，2、3中參數的下標0示材料初始值。JO9159炸藥點火增長模型反應速率方程參數計算模型與結果分析在笛卡兒坐標系中建立三維計算模型如所示。
由于測試實驗裝置為軸對稱，故將其剖為原來四分之一的形狀，在其兩個剖面分別加上關于x，y軸的對稱約束。
ABCD區為JO9159炸藥部分，CDEF區為有機玻璃，EFGH區為JO9159炸藥。采用單點起爆方式，起爆點為A點。
為模型的初分網格，網格劃分為六面體形，邊長取為1mm.為1.42s時模型的等壓分布，由于軸向稀疏波的作用，爆轟產物向后飛散，起爆位置成小山坡狀，這時平面波透鏡內傳播的爆轟波形近似為球面波形。為4s時模型的等壓分布，這時有機玻璃已被壓縮，由于有機玻璃中的沖擊波速度小于高速層炸藥的爆速值，造成在同一軸向距離的位置，轟波到達有機玻璃與高速層炸藥界面的時間與沖擊波沿有機玻璃軸線傳播的時間幾乎相同，因而使有機玻璃中傳播的沖擊波波形近似為一平面。為7.32s時模型的等壓分布，這時底層JO9159已被起爆，由于側向稀疏波的作用，平面波透鏡成為凸出的腰鼓狀，軸向稀疏波的作用使起爆部分的山坡形狀進一步擴大。等壓分布中靠近軸線位置的爆轟波近似為平面，而邊側的等壓分布則由于徑向稀疏波的作用向后計算模型的主視計算模型初分網格的主視1.42s時JO9159炸藥及有機玻璃中的等壓分布4s時JO9159炸藥及有機玻璃中的等壓分布2s時有機玻璃及底層JO9159炸藥中的等壓分布爆轟波透出端面波形的計算值彎曲。為爆轟波透出端面波形的計算值與實驗值的比較，計算結果與實驗值基本相符合，只是在靠近邊側位置的計算結果與實驗結果相差較大，這是由于在建模過程中為了計算穩定，對靠近平面波透鏡側面的有機玻璃進行光滑處理后而造成的波形計算失真。
結論應用三維非線性有限元流體動力學程序（ANSYS/LSDYNA）對平面波透鏡實驗進行了數值模擬，數值模擬中對JO9159炸藥采用點火增長模型，炸藥爆轟產物采用JWL狀態方程，有機玻璃采用Gruneisen方程，計算結果與實驗基本相符。所得炸藥及有機玻璃參數可用于平面波透鏡的設計。
