基于虚单元法(VEM)的固体推进剂药柱结构完整性分析  

Structural integrity analysis of solid propellant grain based on virtual element method(VEM)

作　　者：崔辉如 CUI Huiru(College of Defense Engineering,Army Engineering University of PLA,Nanjing210007,China;Postdoctoral Station of Civil Engineering,Army Engineering University of PLA,Nanjing210007,China)

机构地区：[1]陆军工程大学国防工程学院,南京210007 [2]陆军工程大学土木工程博士后科研流动站,南京210007

出　　处：《固体火箭技术》2023年第4期528-534,共7页Journal of Solid Rocket Technology

基　　金：国家自然科学基金(12202496);中国博士后科学基金资助项目(2022M723866);江苏省青年基金(BK20210435)。

摘　　要：为利用多边形单元在网格划分时的灵活特性,采用虚单元法(VEM)开展固体推进剂药柱结构完整性分析。在虚单元空间的基础上,介绍了虚单元法中刚度矩阵的主要计算流程。针对固体推进剂典型的线粘弹性本构,开展了一系列相关的算例验证工作。仿真结果表明,VEM可以精确地模拟出不同阶跃应变条件下粘弹性薄板的应力松弛响应;在粘弹性薄板双侧蠕变测试中,VEM计算出的薄板双侧应变响应与解析解一致;在点火阶段,利用VEM计算出的圆管药柱环向、径向应力-应变与解析解、ABAQUS计算出的数值解相互吻合。结合VEM在处理多边形单元领域的优势,可以进一步利用VEM进行药柱结构完整性分析领域的深入研究。To make use of the flexible characteristics of polygon elements in mesh generation,the virtual element method(VEM)was utilized to analyze the structural integrity of the solid propellant grain.Based on the virtual element space,the main cal⁃culation flow of the stiffness matrix in the virtual element method was introduced.A series of related examples were performed to veri⁃fy the method for the typical linear viscoelastic constitutive model of the solid propellant.The simulation results show that VEM can accurately simulate the stress relaxation response of viscoelastic thin plate under different step strain conditions.In the bilateral creep test of viscoelastic thin plate,the bilaterod strain responses calculated by VEM are consistent with the analytical solution.In the igni⁃tion stage,the hoop/radial stress and strain of the cylindrical grain calculated by VEM are in agree with the analytical and numerical solutions calculated by ABAQUS.Combined with the advantages of VEM in the field of processing polygonal elements,VEM can be further utilized for in⁃depth research of propellant grain structural integrity analysis.

关 键 词：虚单元法 固体推进剂药柱 结构完整性 粘弹性 

分 类 号：V435[航空宇航科学与技术—航空宇航推进理论与工程]