固体火箭技术2026，Vol.49Issue(3)：409-417,9.DOI:10.7673/j.issn.1006-2793.2026.03.010
基于真实空间分布的固体推进剂三维有限元建模方法研究
Three-dimensional finite element reconstruction method for solid propellant based on real spatial distribution
摘要
Abstract
To address the key challenges in three-dimensional(3D)mesostructure modeling of novel composite solid propellants,such as insufficient component segmentation accuracy,truncation and interlayer defects easily caused by the traditional method of planar segmentation followed by 3D reconstruction,and the difficulty in effectively distinguishing components with similar grayscale values,a mesoscopic finite element reconstruction method integrating μ-CT experiments and spatial digital image processing technology was proposed.This method adopts an optimized technical route of 3D reconstruction prior to spatial segmentation,which preserves the spatial topological relationships of mesoscopic components.Initially,two-dimensional slice images of the propellant specimen at different heights are acquired via μ-CT scanning,and preprocessed to eliminate scanning noise and invalid background regions.Subsequently,all preprocessed slices are stacked along the height direction to construct a three-dimensional grayscale array containing complete spatial position and grayscale information.On this basis,targeted segmentation strategies are implemented for different components respectively,in combination with the grayscale ranges,spatial morphological characteristics and size distributions of each component obtained from SEM/EDS tests.For initial defects,background interference is eliminated by combining the lower grayscale limit with connected domain volume screening.For lead(Pb)particles with independent,non-overlapping grayscale ranges,segmentation is completed directly via grayscale thresholding.For ammonium perchlorate(AP)and aluminum(Al)particles,sequential spatial morphological operations are performed on the basis of preliminary grayscale screening,coupled with screening by connected domain equivalent diameter and voxel count,to remove interference from burrs,noise and agglomerates.For HMX and the viscoelastic matrix with overlapping grayscale values,the mixture region is first obtained through preliminary grayscale range screening.The optimal segmentation threshold is then determined using the Otsu adaptive threshold method.Finally,based on the morphological characteristics and spatial dimensions of HMX,spatial connected domain analysis is applied to filter out undersized matrix noise and oversized agglomerates,realizing effective separation between the two components.The segmentation results show that the volume fractions of AP,Al,Pb,initial defects,HMX,and the viscoelastic matrix are 8.28%,4.37%,0.254%,0.16%,39.74%,and 47.20%,respectively,which are in good agreement with experimental measurements.Compared with conventional methods,the method proposed effectively avoids truncation and interlayer defects caused by planar segmentation,and more realistically reflects the 3D spatial distribution characteristics of each component.Furthermore,a 3D finite element model of the representative volume element(RVE)is established based on the segmentation results,which provides more reasonable structural input and reliable technical support for multiscale mechanical analysis,damage evolution simulation of solid propellants,and structural integrity assessment of propellant grains.关键词
固体推进剂/药柱结构完整性/细观建模/μ-CT试验/空间数字图像处理/组分分割Key words
solid propellant/grain structural integrity/mesoscale modeling/μ-CT test/spatial digital image processing/component segmentation分类
航空航天引用本文复制引用
王沿朝,石峰,赵志鹏,白建方,段磊光,王学仁..基于真实空间分布的固体推进剂三维有限元建模方法研究[J].固体火箭技术,2026,49(3):409-417,9.基金项目
国家自然科学基金(52108387) （52108387）
陕西省自然科学基金(2025JC-YBMS-508) （2025JC-YBMS-508）
固体推进全国重点实验室研究基金资助项目(SPL2025JJ002). （SPL2025JJ002）
