• Synthesis and Characterization of 4，4，8，8-Tetranitro-2，6-dioxaadamantane

  LI Huan, ZHOU Qi, HOU Tian-jiao, WANG Gui-xiang, LUO Jun

  Online:March 13, 2025  DOI: 10.11943/CJEM2025012

  Abstract(10) HTML(94) PDF(1.38 M)( 37)

  Abstract:A novel cage-like energetic compound， 4，4，8，8-tetranitro-2，6-dioxaadamantane， was synthesized via four steps involving oxidative cyclization， oxidation， oximation and gem-dinitration by using 9-oxabicyclo［3.3.1］nona-2，6-diene as raw material. Its structure was characterized by nuclear magnetic resonance（NMR）， Fourier transform infrared spectroscopy （FT‐IR） and elemental analysis （EA）， and single crystal X-ray diffraction （SC-XRD）was adopted to further confirm its crystal structure. The thermal stability was investigated by differential scanning calorimetry‐thermogravimetry（DSC‐TG）analysis， and detonation parameters were predicted by EXPLO5. Results show that 4，4，8，8-tetranitro-2，6-dioxaadamantane has a crystal density of 1.75 g·cm^-3 and belongs to monoclinic system， space groups P2₁/c. Its thermal decomposition temperature is 190.6 ℃， theoretical detonation velocity is 7705 m·s^-1， and detonation pressure is 25.75 GPa.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
• Preparation and properties of surface-modified boron and its composite powder

  WU Jun-ying, WANG Jian-yu, LIU Xin-hang, HU Liang, SHANG Yi-ping, LIU Dan-yang, Chen Lang

  Online:March 13, 2025  DOI: 10.11943/CJEM2024290

  Abstract(20) HTML(106) PDF(1.77 M)( 64)

  Abstract:Boron powder is often used as a combustible agent in energetic materials due to its high calorific value， volume calorific value and clean combustion products. However， the surface oxide layer of boron powder makes it difficult to ignite and brings low combustion efficiency. In order to improve the ignition and combustion performances， boron powder was wetly milled in hot acetonitrile to remove the surface oxide layer for obtaining pre-treated boron powder with high activity， according to the good solubility of boron oxide in acetonitrile solvent. Acetonitrile and n-hexane were used as a dual control agent， and then the pretreated boron powder and highly active aluminum were performed a secondary ball milling to finally prepare the boron-aluminum composite powder with surface-activated boron. The morphological characteristics， thermogravimetric， ignition and combustion characteristics of boron and composite powder were studied. The results showed that the content of surface boron oxide of boron powder was reduced after pretreatment with acetonitrile， pretreatment boron powder was easier to react with oxygen when heated in air， and the percentage of mass increase was 25.6% more than that of untreated boron powder. After pretreatment with hot acetonitrile， the surface boron oxide content decreased， the active boron content increased， and the ignition and combustion performances were significantly improved. The mass of composite powder with boron-aluminum mass ratio of 60/40 increased by 93% when heated in air， the ignition temperature was 738.1 ℃ at low heating rate， and the particle combustion time was 11.2 ms.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
• e-mail：luming@njust.edu.cn Developing Energetic Materials with High Hydrogen Content， Nitrogen Rich and High Enthalpy of Formation

  LU Ming, YANG Feng, WANG Xian-feng

  Online:January 20, 2025  DOI: 10.11943/CJEM2024272

  Abstract(135) HTML(190) PDF(531.46 K)( 295)

  Abstract:

  • 1
  • 2
  • 3
  • 4
• Process Parameters of Single-Layer Stacking Process through 3D Printing of Solid Propellants

  WANG Jia-wei, SHI Hong-bin, LIU Xuan-jie, SONG Shi-xiong, SHI Jia-wei, WANG Qi-hu

  Online:January 03, 2025  DOI: 10.11943/CJEM2024210

  Abstract(128) HTML(54) PDF(2.22 M)( 434)

  Abstract:The process parameters have a direct impact on the 3D printing quality of solid propellant grains. To more reasonably adjust the 3D printing process parameters and improve printing quality， based on the single-layer stacking process， a numerical simulation method was employed to conduct an orthogonal experimental study on three influencing factors： extrusion speed， printing height， and printing temperature. The degree of influence of each factor was calculated through variance and range analysis. The grey relational analysis method was adopted for comparison， and the optimal combination of process parameters was selected after comprehensively considering the printing accuracy of special points. A method for calculating the printing line spacing based on single-line cross-sectional data was proposed for the first time， and simulation and experimental verification were conducted. The results indicated that the extrusion speed had the greatest impact on printing quality. When the extrusion speed was set to 12 mm·s^-1， the nozzle height was 1.2 mm， and the printing temperature was 55 ℃， the printed part exhibited optimal quality. After parameter adjustment， the tensile strength of the specimen increased from 0.21 MPa to 0.43 MPa， and the density rose from 1.43×103 kg·m^-3 to 1.65×103 kg·m^-3. Single-layer printing simulations and experiments demonstrated a significant improvement in molding quality after parameter adjustment.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
• Reaction Characteristics of Perfluoropolyether-Functionalized Micro/Nano Aluminum in Hexanitrohexaazaisowurtzitane

  WU Cheng-cheng, SUN Sen, LI Shi-wei, GUO Xue-yong

  Online:December 27, 2024  DOI: 10.11943/CJEM2024230

  Abstract(123) HTML(56) PDF(3.43 M)( 320)

  Abstract:In order to investigate the reaction characteristics of oxidizer-fuel composite materials with different particle sizes in constant-volume combustion， laser ignition and detonation environment of high energy explosive hexanitrohexaazaisowurtzitane （CL-20）， three particle sizes of perfluoropolyether-functionalized micro/nano aluminum （nAl_150@xPEPE， μAl_1@xPEPE and μAl_5@xPEPE， where x=2.5%， 5.0%， 7.5%） was constructed by particle suspension method， and CL-20 based aluminized explosive was prepared by kneading granulation method. The pressure-time curve， laser-induced ignition process， energy release rate and efficiency of samples in CL-20 were studied by means of closed constant-volume explosive device， laser ignition， detonation velocity and detonation heat test equipment， respectively. The results showed that with the increase of PFPE mass fraction， the peak pressure and pressurization rate of nAl_150@xPEPE samples and μAl_1@xPEPE samples increased gradually， while the peak pressure of μAl_1@7.5%PEPE sample reached 4138.4 kPa and its pressurization rate reached 0.216 MPa·ms^-1. However， when the PFPE mass fraction exceeded 5.0%， the pressurization rate seemed to slow down. At the same time， with the increase of PFPE mass fraction， the burning rate of PFPE-functionalized micro/nano aluminum in CL-20 increased gradually. When x=7.5%， the burning rate of all the three samples with different particle sizes in CL-20 increased by 2.1 cm·s^-1， 1.8 cm·s^-1 and 2.3 cm·s^-1， respectively. In addition， four kinds of fuel-rich CL-20 based aluminized explosives were designed. Among them， the measured detonation velocity of JWL-3 explosive （62% CL-20/32% μAl_1@5.0%PEPE/6% binder） was 8125 m·s^-1， the measured detonation heat was 8049.8 kJ·kg^-1， and the energy release efficiency reached 86.10% （measured by detonation heat）.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
• Synthesis and Properties of 6-Amino-4-（trinitromethyl）-2-carbonyl-1H-1，3，5-triazine

  ZHANG Tong-wei, XU Yuan-gang, LU Ming

  Online:December 12, 2024  DOI: 10.11943/CJEM2024239

  Abstract(155) HTML(192) PDF(1.41 M)( 395)

  Abstract:A novel high-energy compound， 6-Amino-4-（trinitromethyl）-2-carbonyl-1H-1，3，5-triazine， was synthesized in one step. The crystal structure of this compound was characterized by X-ray single crystal diffraction. Its structure and properties were characterized by ¹H and ¹³C NMR， FT-IR and DSC. The detonation performance was calculated by EXPLO5. The sensitivity testing was performed according to the BAM standard method.The compound crystallizes in orthorhombic space group C 2/c， a=10.183（4） ?， b=9.388（3） ?， c=21.324（8） ?， V=2005.9（13） ?³， α=90°， β=100.246（10）°， γ=90°， Z=8. The calculated detonation velocity and pressure for compound 1 are 8167 m·s^-1 and 27.6 GPa， respectively， with measured impact sensitivity of 6 J and friction sensitivity of 210 N.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
• Thermal Behavior and Decomposition Mechanism of 2，2-Azobi［4，5-bis （tetrazole-5-yl）］-1，2，3-triazole

  LIU Shu-liang, CAI Tao, ZHANG Li-nan, QI Yuan, MA Hui-chao, LIN Qiu-han

  Online:November 27, 2024  DOI: 10.11943/CJEM2024229

  Abstract(150) HTML(200) PDF(1.93 M)( 366)

  Abstract:To study the thermal decomposition behavior of 2，2-azobi［4，5-bis（tetrazole-5-yl）］-1，2，3-triazole （NL24）， the structure， morphology and thermal decomposition characteristics of NL24 were studied by means of scanning electron microscopy， thermogravimetric analyzer， differential scanning calorimeter and thermogravimetric infrared mass spectrometry. The kinetic parameters such as apparent activation energy and pre-exponential factor were calculated by Kissinger， Ozawa and ?atava-?esták method， and the thermal decomposition mechanism of NL24 was speculated. Results show that NL24 has two main weight loss stages at the heating rate of 10 ℃·min-1 . The first weight loss stage occurs at about 180 ℃， which belongs to the volatile endothermal process of dimethyl sulfoxide. The violent thermal decomposition of NL24 occurs at the second weight loss stage between 270 ℃ and 300 ℃， which has not only rapid gas generation rate， but also belongs to autocatalytic reaction. The main gaseous products are N2， HCN， HN3， etc. The apparent activation energy and pre-exponential factor of the decomposition process are 174.69 kJ·mol-1 and 1016.60 s-1 ， respectively. The reaction model of thermal decomposition stage of NL24 is random nucleation and subsequent growth.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
• Synthesis and Properties of N²，N⁶-Dimethyl-N²，N⁴，N⁶，3，5-Pentanitro-2，4，6-Pyridinetriamine

  ZHANG Rong-zheng, LU Ming, XU Yuan-gang

  Online:November 21, 2024  DOI: 10.11943/CJEM2024211

  Abstract(129) HTML(142) PDF(1.56 M)( 398)

  Abstract:A pyridine energetic molecule， N²，N⁶-dimethyl-N²，N⁴，N⁶，3，5-pentanitro-2，4，6-pyridinetriamine （NNDP）， has been synthesized in two steps from 4-amino-2，6-dichloropyridine. The process was found to be effective and simple. The structure of this compound is characterized by ¹H and ¹³C NMR， FT-IR and DSC. The crystal structure of this compound is characterized by X-ray single crystal diffraction. Results shows that compound NNDP belongs to the monoclinic space group P 2₁/c， a=16.3215（17） ?， b=7.9819（8） ?， c=13.1954（13） ?， V=1712.3（3） ?³， α=90（6）^o， β=95.093（3）^o， γ=90（7）^o， Z=4. The presence of multiple nitro and nitramine groups contributes to a low decomposition temperature. Its detonation performance was predicted using EXPLO5， and sensitivity testing was conducted using the BAM standard method.It was found that the detonation performance and impact sensitivity of NNDP（D=8762 m·s^-1， p=34.5 GPa， IS=7.7 J） are comparable to those of RDX.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
• Theoretical Study of the Electronic Structure， Bonding properties and Aromaticity of Hexazine Anion ［N₆］^4-

  JIANG Shuai-jie, ZHANG Guang-yuan, XU Yuan-gang, WANG Peng-cheng, LU Ming

  Online:November 13, 2024  DOI: 10.11943/CJEM2024205

  Abstract(97) HTML(140) PDF(1.86 M)( 355)

  Abstract:To understand the properties of the novel polynitrogen compound hexazine anion ［N₆］^4-， computational chemical methods were used to study the electronic structure， bonding properties and aromaticity of N₆， ［N₆］^2- and ［N₆］^4-. The M06-2X method combined with the def2-TZVP basis set was used to optimized the structures and calculated the electronic structure features， such as bond length， bond angle， dihedral angle， molecular size and so on. Subsequently， multiple bond orders were calculated， using the atoms-in-molecules （AIM） theory to calculate multiple bond properties， and drawing the electron deformation density map to directly show the bond behavior. Finally， various aromatic indices were calculated to show the aromatic characteristics of three hexazine rings. The calculation results show that by comparing with the electronic structure optimized by CCSD， the M06-2X method in the common DFT method is suitable for studying the current system. Mayer bond order shows that the N—N bond has a certain degree of σ bond characteristics. The aromaticity study shows that the［N₆］^4- is aromatic， with the aromatic harmonic oscillator model （HOMA） value at 0.96 and the nuclear independent chemical shift （NICS_ZZ（1）） at -18.97 ppm. The IR， Raman and UV-Visible spectra of ［N₆］^4- were simulated to provide reference for experimental detection.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
• Synthesis and Performance of 5-Nitro-3-（trinitromethyl）-1H-1，2，4- triazole Nitrogen-Rich Energetic Ionic Salts

  WANG Xian-feng, YANG Feng, XU Yuan-gang, LU Ming

  Online:September 18, 2024  DOI: 10.11943/CJEM2024224

  Abstract(197) HTML(102) PDF(2.31 M)( 545)

  Abstract:To further balance the energy and safety of 5-nitro-3-（trinitromethyl）-1H-1，2，4-triazole， four nitrogen-rich energetic ionic salts were synthesized using 2-（5-amino-1H-1，2，4-triazole-3-yl） acetic acid as a starting material through a silver salt substitution reaction. The structures of all new compounds were characterized using nuclear magnetic resonance， Fourier transform infrared spectroscopy， differential scanning calorimetry， thermogravimetric analysis， and single crystal X-ray diffraction. The results indicate that the ammonium salt， hydrazine salt， and guanidine salt of 5-nitro-3-（trinitromethyl）-1H-1，2，4-triazole exhibit higher initial decomposition temperature than that of the precursor. Moreover， the hydrazine salt and guanidine salt belong to the different crystal systems with distinct crystal packing arrangements and densities. However， they share consistent characteristics in terms of intermolecular weak interactions， with the H…O interaction being the predominant contributor. With the decreasing of the ratios of N…O and O…O interactions， the sensitivity of the nitrogen-rich energetic ionic salts to mechanical stimuli decreases. Finally， the analysis of the distribution of molecular electrostatic potential supplements the explanation for the change in impact sensitivity of 5-nitro-3-（trinitromethyl）-1H-1，2，4-triazole after salt formation. Among the four ionic compounds， the hydrazine salt exhibits outstanding detonation performance （D=8634 m·s^-1， p=30.2 GPa， I_sp=263.5 s） with relative high sensitivity. In contrast， the triaminoguanidine salt demonstrates excellent overall performance. It has a detonation velocity comparable to that of the hydrazine salt （D=8627 m·s^-1）， a heat of formation nearly 1.4 times greater than that of the precursor （ΔH_f=0.644 kJ·g^-1）， and a low mechanical sensitivity （IS=10.3 J， FS=150 N）.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25

• Combustion Characteristics of Ternary Reactive Metal Fuels Al/B/Mg and Al/B/MgH₂

  CHENG Zhi-peng, XIA Yu, LUO Yi-min, MA Teng, XU Fei-yang, ZHANG Yu, WU Xing-liang, XU Sen

  Online:August 16, 2024  DOI: 10.11943/CJEM2024100

  Abstract(370) HTML(280) PDF( 1678)

  Abstract:To investigate the combustion characteristics of ternary active metal fuels Al/B/Mg （ABM） and Al/B/MgH₂ （ABM-H）， the heat of combustion and minimum ignition energy were studied by using an oxygen bomb calorimeter and a Hartmann tube， respectively. The sub-transient process of flame propagation and the spatiotemporal distribution characteristics of temperature fields were determined by using a high-speed camera system and a high-speed infrared camera system. The results indicate that the calorific values of ABM and ABM-H are 34.1 and 32.2 MJ·kg^-1， respectively， exhibiting increases of 14.4% and 8.1% over pure Al （29.8 MJ·kg^-1）. The minimum ignition energies of ABM， ABM-H， and Al are 160-170， 100-110， and 20-30 mJ， respectively. Compared to pure Al， the combustion duration of ABM and ABM-H increase by 65.5%， 34.5% and the peak flame propagation velocities increase by 12.6%， 23.0%， respectively， at a mass concentration of 625 g·m^-3. At a mass concentration of 500 g·m^-3， ABM-H and ABM exhibit the largest peak flame propagation velocities by 45.05， 38.7 m·s^-1， and the maximum temperatures peak of flame surface by 1856， 1717 ℃， respectively， where ABM-H shows a 7.6% improvement on temperatures peak of flame surface and a faster heating-rate compared to ABM. It suggests that the ABM and ABM-H formulations significantly reduce the explosion risk of the dust/air mixture， and significantly improving the combustion performance. ABM demonstrates superior thermal effects in calorific value and duration of combustion， whereas ABM-H exhibits higher reactivity in terms of minimum ignition energy， flame propagation speed， and temperature rise rate.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
  • 109
  • 110
  • 111
  • 112
  • 113
  • 114
• Research on Combustion Characteristics and Injury Effects of Methane Vapor Clouds in Tunnels

  GU Lin-lin, XU Yong-hang, ZHU Huang-hao, WANG Zhen

  Online:August 20, 2024  DOI: 10.11943/CJEM2024155

  Abstract(272) HTML(140) PDF( 1357)

  Abstract:In order to explore the propagation patterns and characteristics of methane vapor cloud combustion waves in tunnels， the CE/SE （space-time conservation element and solution element） method in LS-DYNA software was employed to establish a pre-mixed explosion model of methane and air in the tunnel， which was validated through experimental data. In this paper， typical combustion waveforms of methane vapor cloud with a concentration 9.5% in different test positions were demonstrated by numerical simulation. The propagation and evolution law of overpressure and temperature was analyzed. The injury effects of overpressure and thermal radiation on human in tunnel were investigated. It was revealed that the combustion pressure wave along the tunnel can be divided into four stages： free expansion， reflection dissipation， wall acceleration， and Mach propagation. The pressure variation presented three characteristics： wall impact rise， reflective decay， and stable propagation. The pressure wave presented a sort of periodical reflection propagation mode radially， while the intensity was declining according to the consumption of methane. The temperature field evolved symmetrically from the ignition point to the tunnel entrance and the peak temperature decayed rapidly along the path. The temperature field radiated from the ignition point to the bottom of the tunnel， leading to a gradual convergence of in a certain section and decreased slowly over time. For the injury effects caused by a combination of combustion overpressure and thermal radiation， the fatal distance was 13.51m， the severe injury distance was 13.51~23.51m， the moderate injury distance was 23.51-160 m while the concentration of methane vapor cloud was 5%. For the methane vapor cloud with a concentration 6.5%， those distances were 16.46 m， 16.46~45.36 m and 45.36~160 m respectively. As for a concentration 9.5%， the fetal distance was 20.58m and the severe injury distance was 20.58~160m.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
• Experimental and Numerical sdudy on the Explosion Characteristics of Hydrogen-Methane-Air in a Closed Vessel

  MEI Liang, GUO Jin, HUANG Shi-kai, WANG Jin-gui

  Online:October 25, 2024  DOI: 10.11943/CJEM2024186

  Abstract(101) HTML(124) PDF( 340)

  Abstract:In order to guide the application of hydrogen/methane mixture as fuel， explosion experiments were carried out using a cylindrical closed vessel with an inner diameter and length of 300 mm. The effects of hydrogen fraction （X_H2） from 0 to 100% and equivalence ratio （Φ） from 0.6 to 1.4 on the flame evolution and explosion pressure were investigated. Meanwhile， CHEMKIN software was introduced to analyze the laminar burning velocity and sensitivity coefficient of the H₂-CH₄-air premixed gas. The results showed that， for a certain Φ， the maximum explosion pressure （p_max）， the maximum pressure rise rate （（dp/dt）_max）， the explosion index （K_G）， and the laminar burning velocity increased monotonically with the increase of X_H2. The duration to reach p_max and （dp/dt）_max ， named t_A and t_B， respectively， decreased gradually. After ignition， the flame surface gradually transformed from a smooth structure to a honeycomb flame lattice structure. With a constant Φ and an increasing X_H2， the duration from ignition to the termination of the explosion decreased dramatically. Meanwhile， at the same moment the flame radius increased but the fold on the flame surface increased. The simulation results showed that the elementary reactions R35 and R52 had the most significant influence on the laminar burning velocity. The maximum molar fractions of the key radicals （H， O， and OH） had a positive correlation with the laminar burning velocity， and the increase of X_H2 lead to a significant increase in the maximum molar fractions of the key radicals. The primitive reactions R38 and R84 were the dominant reactions affecting the rate of production （ROP） of key radicals.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
  • 109
  • 110
  • 111
  • 112
  • 113
  • 114
• Design and Optimization of Multi-layer Composite Structure under Combined Loading of Shock Wave and Fragments

  ZHOU Meng, LIANG Min-zu, CHEN Rong, LIN Yu-liang, ZHANG Yu-wu

  Online:February 28, 2025  DOI: 10.11943/CJEM2024175

  Abstract(63) HTML(40) PDF( 144)

  Abstract:To effectively mitigate the combined impulsive and fragmentary loads exerted by near-field explosion， a multi-layer composite protective structure has been developed， incorporating an anti-penetration layer， a coordinated support layer， and a cushioning energy-absorption layer. A finite element analysis model was established and validated through near-field explosion experiment. Utilizing the outcomes of finite element simulations， a response surface approximation model for the composite protective structure was constructed. With the areal density and overall thickness of the composite protective structure as the optimization targets， a multi-objective optimization of the thickness distribution was conducted under the individual and combined effects of fragments and shock waves using the Non-dominated Sorting Genetic Algorithm II （NSGA-II）. And the Pareto optimal solution set was obtained. The findings demonstrate that， in comparison to the initial design， the areal density of the optimized composite protective structure subjected to the individual action of fragments can be decreased by up to 19.2%， with a maximum thickness reduction of 10.0%. Under the individual action of shock waves， the areal density can be reduced by up to 34.9%， and the thickness by up to 27.5%. Under the combined action of shock waves and fragments， the areal density can be reduced by up to 19.2%， and the thickness by up to 10.0%. For application scenarios where the thickness is constrained to no more than 40 mm， the optimized composite protective structure exhibits an approximately 17.5% reduction in areal density and a 9.1% reduction in total thickness compared to the initial design. It was noted that the Pareto optimal solution sets obtained from the individual fragment action and the combined action are nearly identical， indicating that the composite protective structure significantly diminishes the influence of shock waves on the subsequent fragment impact， thereby effectively mitigating the combined effects of shock waves and fragments.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
• Plasma Characteristics of Typical Nitrogen-containing Compounds in Nitrogen Atmosphere

  LIU Hai-qing, XIANG Shu-jie, FANG Pu-yixing, LI Chun-tian, SHEN Rui-qi, ZHANG Wei

  Online:October 15, 2024  DOI: 10.11943/CJEM2024101

  Abstract(139) HTML(112) PDF( 420)

  Abstract:Nitrogen-containing compounds， acting as nitrogen donors， directly influence the types of high nitrogen compounds formed under laser irradiation. To understand the impact of various nitrogen-containing compounds on the formation of high-nitrogen compounds， three representative compounds NaN₃， Si₃N₄ and P₃N₅ were ablated using a pulsed Nd： YAG laser in a nitrogen atmosphere. The plasma characteristics and the evolution of the transient intermediates generated by laser sputtering were investigated by transient spectrometer. The findings indicate that the laser ablation of NaN₃ yields the highest number of nitrogen atoms （N Ⅰ）， monovalent nitrogen ions （N Ⅱ）， and trivalent nitrogen ions （N Ⅲ）， with the longest duration of nitrogen plasma. The lifetimes of N Ⅰ， N Ⅱ， and N Ⅲ reached 39，400 ns， 39，400 ns， and 19，450 ns， respectively. Among the three nitrogen donors， the laser ablation of NaN₃ in a nitrogen atmosphere is most likely to result in the formation of high-nitrogen or all-nitrogen compounds.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
  • 109
  • 110
  • 111
  • 112
  • 113
  • 114
  • 115
  • 116
  • 117
  • 118
  • 119
  • 120
  • 121
  • 122
  • 123
  • 124
  • 125
  • 126
  • 127
  • 128
  • 129
  • 130
  • 131
  • 132
  • 133
  • 134
  • 135
  • 136
  • 137
  • 138
  • 139
  • 140
  • 141
  • 142
  • 143
  • 144
  • 145
  • 146
  • 147
  • 148
  • 149
  • 150
  • 151
  • 152
  • 153
  • 154
  • 155
  • 156
  • 157
  • 158
  • 159
  • 160
  • 161
  • 162
  • 163
  • 164
  • 165
  • 166
  • 167
  • 168
  • 169
  • 170
  • 171
  • 172
  • 173
  • 174
  • 175
  • 176
  • 177
  • 178
  • 179
  • 180
  • 181
  • 182
  • 183
  • 184
  • 185
• Using BDD Electrode for Electrochemical Treatment of Nitrogen in Real Energetic Material Wastewater

  YOU Jia-jun, WANG Bing, XIONG Ying, WANG Jian

  Online:June 04, 2024  DOI: 10.11943/CJEM2024092

  Abstract(476) HTML(132) PDF( 1834)

  Abstract:The pilot-scale energetic material wastewater is a kind of wastewater which is extremely difficult to degrade， containing high concentrations of various nitrogen-containing compounds， such as ammonia nitrogen （NH₃-N）， nitrite （NO₂^-）， nitrate （NO₃^-）， and other organic pollutants . To realize the efficient and directional removal of these nitrogen-containing compounds， boron-doped diamond （BDD） electrodes were prepared by the hot-filament chemical vapor deposition （HFCVD） method and utilized to degrade the wastewater. The effects of electrolyte composition and concentration， modified electrode type， and electrolysis device structure on the degradation efficiency were investigated. It demonstrated that adding 0.1 M sodium chloride （NaCl） electrolyte to energetic material wastewater could improve the selectivity of NH₃-N direct conversion to nitrogen （N₂）. Using Cu/BDD and Ni/BDD cathodes accelerate the conversion process of high-valent nitrogen to NH₃-N. Under the dual electrolysis cell structure system， employing Cu/BDD and Ni/BDD electrodes as anodes improve the degradation rate of NH₃-N conversion to N₂. Therefore， the approach utilizes metal-modified BDD electrodes as anodes， is expected to be a highly effective method in the rapid and selective degradation of energy material wastewater， especially when using 0.1 M NaCl as electrolyte.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
• Initial Thermal Aging Behavior of PNIMMO Studied by Isothermal Gasometric Method

  GONG Xue-ling, GUAN Jian, LIU Hong-ni, MO Hong-chang, ZHANG Qing-yuan, PENG Ru-fang, JIN Bo

  Online:November 19, 2024  DOI: 10.11943/CJEM2024157

  Abstract(102) HTML(96) PDF( 448)

  Abstract:The initial isothermal aging behavior of poly （3-nitratomethyl-3-methyloxetane） （PNIMMO） was studied. The aging kinetic parameters and thermal aging mechanism of PNIMMO at 100-120 ℃ were investigated using an isothermal gas measuring device. The storage life of PNIMMO was determined by the Berthelot equation. The results indicate that the activation energy （E_a） is 156.42 kJ·mol^-1 and the logarithm of the pre-exponential factor （lgA） is 16.86 s^-1 when the aging depth of PNIMMO reaches 0.1%. Conversely， at an aging depth of 0.5%， E_a is measured at 156.05 kJ·mol^-1 and lgA at 16.03 s^-1， as derived from the Arrhenius equation. According to the mode matching method， the thermal aging of PNIMMO at 100-120 ℃ conforms to the mechanism function No.28， that is， the reaction order n=1/4， E_a=154.33 kJ·mol^-1. Using an aging depth of 0.1% as the evaluation criterion， PNIMMO can be stored at room temperature for 51.6 years. During the initial phase of thermal decomposition， the side chain ─O─NO₂ bond undergoes cleavage followed by hydrogenation， subsequently leading to gradual degradation of the main chain into stable polyaromatic compounds.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
• Simulation Study on the Anisotropy Rule of Impact Sensitivity of ε-CL-20

  LIANG Lin, WANG Ya-jun, GAN Qiang, ZHANG Wen-bo, REN Shu, LI Gen, FENG Chang-gen

  Online:November 19, 2024  DOI: 10.11943/CJEM2024121

  Abstract(133) HTML(56) PDF( 414)

  Abstract:To investigate the anisotropy of impact sensitivity of the cage-like energetic material hexanitrohexaazaisowurtzitane （ε-CL-20），this work used the ReaxFF-lg reactive force field and molecular dynamics method， multiscale impact loading simulations were performed on six typical crystallographic planes： （010）， （110）， （20）， （011）， （11）， and （001）. The correlation between stress， temperature， chemical reactions， and the direction of impact was analyzed. Results indicate a pronounced anisotropy in the impact sensitivity of ε-CL-20， with the sensitivity ranking of the planes as （010）>（110）>（20）≈（011）>（11）> （001）. The system exhibits the strongest thermo-mechanical and chemical responses when impacted perpendicular to the （010） plane， implying the highest sensitivity. In contrast， the weakest responses and lowest sensitivity occurs when impacted perpendicular to the （001） plane. Based on these findings， for planar layered energetic materials， impacts parallel to the molecular layers yield the highest sensitivity， while the impacts perpendicular to the molecular layer opposite result in low sensitivity.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
  • 109
  • 110
  • 111
  • 112
• Solubility and Thermodynamic Modeling of 3-Nitro-l，2，4-triazole-5-one （NTO） in Different Binary Solvents

  GUO Hao-qi, YANG Yu-lin

  Online:December 25, 2024  DOI: 10.11943/CJEM2024184

  Abstract(84) HTML(140) PDF( 297)

  Abstract:Using a dynamic laser monitoring technique， the solubility of 3-nitro-l，2，4-triazole-5-one （NTO） was investigated in two different binary systems， namely hydroxylamine nitrate （HAN）-water and boric acid （HB）-water ranging from 278.15 K to 318.15 K. The solubility in each system was found to be positively correlated with temperature. Furthermore， solubility data were analyzed using four equations： the modified Apelblat equation， Van’t Hoff equation， λh equation and CNIBS/R-K equations， and they provided satisfactory results for both two systems. The average root-mean-square deviation （10⁵RMSD） values for these models were less than 13.93. Calculations utilizing the Van’t Hoff equation and Gibbs equations facilitated the derivation of apparent thermodynamic properties of NTO dissolution in the two systems， including values for Gibbs free energy， enthalpy and entropy. The %ζ_H is larger than %ζ_TS， and all the %ζ_H data are ≥ 58.63%， indicating that the enthalpy make a greater contribution than entropy to the ΔG_soln^o.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
• 3D Printing Technology in the Manufacturing of Solid Propellant Grain： A review of Applications and Prospects

  SONG Shi-xiong, REN Quan-bin, WANG Jia-wei, PANG Ai-min, TANG Min

  Online:September 14, 2024  DOI: 10.11943/CJEM2024129

  Abstract(298) HTML(300) PDF( 937)

  Abstract:3D printing technology has the characteristics of customization， mold free， and flexibility， which can provide an effective approach for the shaping of special structure solid propellant grains in multi-thrust or multi -pulse solid rocket motors. At present， research on 3D printing of solid propellant grain has been conducted both domestically and internationally. This article focuses on the application of typical 3D printing processes such as binder jetting， photopolymerization curing， and material extrusion in the formation of heterogeneous solid propellant grains with contained complex structures， gradients， and multi material integration. It summarizes the key issues that exist in the 3D printing of these three types of solid propellant grains. The future research directions were prospected， and it was emphasized that the future manufacturing of heterogeneous solid propellant grains should focus on low sensitivity specialized solid propellant slurries， printing equipment for large grain forming， and insulation coating printing technology.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66

Vol, 33, No.2, 2025    

>Preparation and Property

• Synthesis， Structure， and Performance of Potassium 1-（Dinitromethyl）-3-nitro-5-azidopyrazole

  KANG Xue-meng, DING Ning, XU Xu-dong, SUN Qi, LI Sheng-hua, PANG Si-ping

  2025,33(2):103-109, DOI: 10.11943/CJEM2024296

  Abstract(230) HTML(19) PDF(1007.16 K)( 183)

  Abstract:1-（Trinitromethyl）-3-nitro-5-azidopyrazole （compound 1） was synthesized from 1-acetonyl-3，5-dinitropyrazole via nitration. The 1-（dinitromethyl）-3-nitro-5-azidopyrazole potassium salt （compound 2） was further prepared through a salification reaction with KI. The structures of both two compounds were characterized by Fourier-transform infrared spectroscopy （FT-IR）， nuclear magnetic resonance （NMR）， and single crystal X-ray diffraction. The results show that compound 2 contains both azido and dinitromethylene groups， exhibiting superior overall properties of high decomposition temperature （151.41 ℃） and detonation velocity （8369 m∙s^-1）. The relatively high sensitivity towards external mechanical stimuli make it suitable to be a primary explosive. The detonation performance of compound 2 was further analyzed using Gaussian and Explo5 software. This research provides valuable insights for the design of novel green explosives.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
• Preparation and Characterization of CL-20/MNP Co-crystal

  MA Fan, LIU Ze-ning, NING Ting-yao, QIAO Shen, YANG Zong-wei, LIANG Jin-hu

  2025,33(2):110-116, DOI: 10.11943/CJEM2025004

  Abstract(199) HTML(31) PDF(2.12 M)( 136)

  Abstract:In order to reduce the mechanical sensitivity of hexanitrohexaazaisowurtzitane（CL-20）， the co-crystal explosives CL-20/3-MNP （co-crystal 1） and CL-20/4-MNP （cocrystal 2） were prepared by cocrystal technology by combining the isomeric molecules of CL-20 and methyl-nitropyrazole（MNP）：1-methyl-3-nitropyrazole（3-MNP） and 1-methyl-4-nitropyrazole （4-MNP）. The structure of the single crystal was determined by single crystal X-ray diffractometer （SXRD）， and the thermal properties and impact sensitivity of the cocrystal were measured by differential scanning calorimetry （DSC） and BAM drop weight impact sensitivity instrument， and the detonation performance was predicted by EXPLO5. The results show that cocrystal 1 is a triclinic crystal system with a P space group. Cocrystal 2 is an orthorhombic crystal system with a Pbca space group. The melting points of cocrystal 1 and cocrystal 2 are 113.2 ℃ and 147.2 ℃， respectively， which are 25.2 ℃ and 53.7 ℃ higher than those of 3-MNP and 4-MNP， respectively， and the decomposition temperature is higher than that of CL-20， which has better thermal stability， and the impact sensitivity of cocrystal 1 and cocrystal 2 is 20 J and 18 J， respectively， which effectively reduces the sensitivity of CL-20 （2.5 J）. The theoretical detonation velocity is 8060 m·s^-1 and 8643 m·s^-1， and the theoretical detonation pressure is 26.16 GPa and 32.61 GPa， respectively， which are lower than those of CL-20， and the cocrystal 2 sensitivity and detonation performance are comparable to those of LLM-105（17 J， 8560 m·s^-1， 34.99 GPa）， which is expected to be used as a new high-energy and low-inductance cocrystal explosive.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
• Preparation and Properties of TKX-50/AP Composite Microspheres

  WANG Xin-fei, YANG Jun-qing, ZHAO Han-xiao, WANG Ting-hui, HAO Ga-zi, XIAO Lei, JIANG Wei

  2025,33(2):117-126, DOI: 10.11943/CJEM2024267

  Abstract(179) HTML(19) PDF(1.81 M)( 138)

  Abstract:To obtain new energetic materials with high energy， low sensitivity and liable storage， TKX-50/AP composite microspheres were prepared by spray drying method. The effects of different feed rates on the morphology of composite microspheres were investigated， and the optimum process conditions were obtained. The morphology and crystal form of the samples were characterized by scanning electron microscope， X-ray diffractometer and Fourier infrared spectrometer. The thermal decomposition performance， impact sensitivity and hygroscopicity of the microspheres were also studied. Results show that the samples are spherical particles with sphericity of Φ=96.40% and d₅₀=1.33 μm at the feed rate of 3.6 mL·min^-1. The crystal structure is different from that of mechanical mixing of TKX-50 and AP. The thermal decomposition activation energy is 184.43 kJ·mol^-1， the characteristic drop is 50.1 cm， and the hygroscopicity is 37% lower than that of the physical mixing， showing excellent safety performance.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
• Preparation and Properties of Modified Spherical ADN with Amine Compounds

  QU Wei-chen, YAN Jia-wei, XIA Wen-tao, LI Lei, DU Fang, WU Peng, LIN Yu-hui, HE Jian-ming, TAO Bo-wen

  2025,33(2):127-135, DOI: 10.11943/CJEM2024222

  Abstract(150) HTML(15) PDF(3.19 M)( 102)

  Abstract:Ammonium dinitramide （ADN）， as a high-energy green oxidizer， faces significant challenges in engineering applications within solid propellants due to its high surface polarity and strong hygroscopicity. To improve the hygroscopicity of ADN， both emulsion and liquid-phase methods were employed to graft hydrophobic amine compounds， oleylamine and 4-fluorobenzylamine， onto the surface of spherical ADN （PADN）， thereby preparing the modified spherical ADN（PMADN）. The structure and properties of the modified ADN were characterized using scanning electron microscopy-energy dispersive spectroscopy （SEM-EDS）， thermogravimetric-differential thermal analysis （TG-DTA）， and X-ray photoelectron spectroscopy （XPS）. Additionally， hygroscopicity studies were conducted using the desiccator equilibrium method. The results confirmed both amine compounds were grafted successfully and the spherical morphology of the modified ADN was intact. The thermal decomposition temperature of the modified ADN increased from 196.3 ℃ （raw ADN） to 198.3 ℃ and 200.7 ℃， and the impact sensitivity improved from 17.95 J to 24.35 J and 28.80 J （average value）， respectively. After 144 h under 25 ℃ and 57% relative humidity， the moisture rates of the two modified spherical ADN were 1.37% and 1.07%， decreasing 74.95% and 80.44% compared to that of the raw ADN. Additionally， no caking or liquid water was observed on the surface， indicating the modified ADN had excellent anti-hygroscopic properties.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38

>Analysis and Testing

• Machine Learning Recognition of Impact Sensitivity of Energetic Materials Based on Acoustic Signals

  ZHANG Bing-ru, LI Ming, WEN Yu-shi, GUO Feng, YU Shao-jun, JI Chun-liang, LIN Cong-mei, HAO Li-xiao, HAN Yong, DENG Chuan, DAI Xiao-gan

  2025,33(2):136-147, DOI: 10.11943/CJEM2024300

  Abstract(191) HTML(30) PDF(3.20 M)( 145)

  Abstract:To improve the accuracy and objectivity of explosives’ impact sensitivity testing， machine learning methods were applied in the intelligent recognition of explosives’ impact response acoustic signals. Experiments on mixed explosives were conducted using a drop-weight impact sensitivity test device， with an audio acquisition system synchronously capturing acoustic signals during the impact process. One-dimensional time domain and frequency domain features， such as maximum value and bandwidth， were extracted. Short-Time Fourier Transform （STFT） was employed to convert audio data into frequency spectrograms. Data augmentation for one-dimensional data was performed using a Conditional Generative Adversarial Network（cGAN）， while a Deep Convolutional Generative Adversarial Network（DCGAN） was applied to enhance spectrogram data. Multiple machine learning models were employed for explosion classification， including Random Forest （RF）， eXtreme Gradient Boosting （XGBoost）， Back-propagation Neural Network （BPNN）， Support Vector Machine （SVM）， k-means clustering， Convolutional Neural Network （CNN）， and Vision Transformer （ViT）. Results demonstrate that RF， XGBoost， BPNN， and SVM achieve accuracy rates exceeding 99.5% on the real dataset and achieve 100% on the cGAN-augmented dataset. In contrast， k-means clustering initially reaches an accuracy of 98.5% on the real dataset， but accuracy shows a trend of increase followed by decline on augmented data. CNN and ViT achieve accuracies of 98.1% and 98.4% on the real dataset， respectively， and improved to 98.4% and 98.9% on augmented data. However， their performance remaine slightly lower than traditional models due to the constraints of small sample sizes and minor overfitting issues. The proposed deep learning-based intelligent recognition method for explosives’ impact sensitivity in this study achieved a high level of accuracy， demonstrating its reliability and practicality in the task of detecting explosive sound signals. At the same time， it effectively mitigates the subjectivity and inefficiency associated with traditional manual recognition methods， providing a reliable technical solution for the safety use of explosives.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
  • 109
  • 110
  • 111
  • 112
  • 113
  • 114
• Detection of Hydrazoic Acid Gas by Carboxylated Multi-walled Carbon Nanotube Modified Screen Printed Electrode

  LIN Jun-chi, BA Shu-hong, HAN Ji-min, YANG Li

  2025,33(2):148-155, DOI: 10.11943/CJEM2024257

  Abstract(144) HTML(7) PDF(1.42 M)( 95)

  Abstract:For the fast and easy detection of hydrazoic acid gas， an electrosensor for in-situ detection of hydrazoic acid was prepared based on the principle of electrochemical analysis using a screen-printed electrode as a substrate and modified with carboxylated multi-walled carbon nanotubes. The electrochemical detection of hydrazoic acid was constructed by optimizing the solvent of the modification solution， the pH value of the detection solution， and the scan rate. The morphology and performance were characterized. The results show that the in-situ detection electrosensor for hydrazoic acid was prepared as a microelectrode with carboxylated multi-walled carbon nanotubes/glacial acetic acid modification solution， and its response current is about 121% higher than that of the unmodified electrode. The detection sensitivity is high at the solution pH 7.5. The square root of the scan rate is linearly related to the oxidation peak current， and the electrochemical oxidation of N₃^- is a diffusion-controlled process with good selectivity， stability， and reproducibility. The detection limit of N₃^- was 10.4 µmol·L^-1 in the concentration range of 5×10^-5-1×10^-3 mol·L^-1 N₃^- using the differential pulse voltammetry method. The prediction equations for the concentration of HN₃ gas produced by different NaN₃ feedstock contents at different times were derived from the online detection of the actual synthesis of HN₃ gas， and the recoveries of HN₃ are 96.8%-99.5%. In addition， the relationship between the concentration of synthesis HN₃ gas and the response current has been established.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64

>Propulsion and Projection

• Molecular Dynamics Simulations of the Condensation Behavior of Nitroglycerine-containing Volatiles on Solid Surfaces

  HUANG Chao-ran, PAN Bao, LI Guo, WANG Yu, XIE Lin-sheng, LIU Suo-en, ZHAO Jun-bo, LIU Yun-zhang, TAN Kai-xin

  2025,33(2):156-164, DOI: 10.11943/CJEM2024178

  Abstract(154) HTML(8) PDF(2.40 M)( 77)

  Abstract:The condensation and accumulation of Nitroglycerin （NG）-containing volatiles on various solid surfaces during the propellant rolling process， which pose safety hazards， were investigated using molecular dynamics simulation methods. The study was conducted by constructing a hybrid system model consisting of NG volatiles and solid surfaces， examining the effects of solid surface material， surface roughness， and NG content on molecular dynamics characteristic parameters such as radial distribution function， mean square displacement， diffusion coefficient， and relative density distribution of NG volatiles in the hybrid system. The findings demonstrate that as the mass fraction of NG increases， the size of volatile condensate clusters on the solid surface progressively diminishes. Conversely， the condensation ratio of volatiles exhibits a trend of initial increase followed by a decrease， with the maximum condensation ratio occurring at 70% NG， corresponding to a diffusion coefficient of 0.0364. The diffusion coefficient for the condensation of volatiles containing NG on a silica （SiO₂） surface is 2.1228， which is substantially greater than that on surfaces composed of copper （Cu）， calcium oxide （CaO）， and ferrum （Fe）. However， the uniformity of the SiO₂ surface condensate cluster is poor. The introduction of surface roughness factors has opposite effects on the condensation amount of volatiles on the SiO₂ and Fe surfaces. When the SiO₂ surface goes from smooth to roughness of 0.4 nm， the diffusion coefficient increases from 2.1228 to 10.7156， and the condensation amount of volatiles on the surface increases； however， when the Fe surface goes from smooth to roughness of 0.4 nm， the diffusion coefficient decreases from 17.5673 to 1.8462， and the condensation amount of the surface volatiles decreases.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
  • 109
  • 110
  • 111
  • 112
  • 113
  • 114
  • 115
  • 116
  • 117
  • 118
  • 119
  • 120
  • 121
  • 122
  • 123
  • 124
  • 125
  • 126
  • 127
  • 128
  • 129
  • 130
  • 131
  • 132
  • 133
  • 134
  • 135
  • 136
• Response Characteristics Analysis and Parameters Calibration of Ignition and Growth Model of NEPE Propellants under Impact Load

  GUO Zong-tao, XU Jin-sheng, CHEN Xiong, CAO Xin-yu, PANG Song-lin, WANG Jin-dong

  2025,33(2):165-177, DOI: 10.11943/CJEM2024246

  Abstract(154) HTML(12) PDF(2.05 M)( 83)

  Abstract:In order to ensure the safe application of solid rocket motors in complex battlefield environments， it is necessary to conduct in-depth research on the response characteristics of NEPE high-energy solid propellants under impact loads. A one-dimensional Lagrange test system was established to measure the pressure wave at different Lagrange location. The unreacted shock adiabatic curve of the measured propellant was obtained through the momentum conservation relationship before and after the shock wavefront， and the JWL EOS of the unreacted NEPE propellant were obtained by using a genetic algorithm. A Φ50 mm NEPE propellant probe-type cylinder test platform was constructed and a 12-channel probe with radial displacement difference was used to record the time of copper cylinder expansion to different probe positions and the time curve of cylinder expansion velocity was obtained. Based on the results of the Φ50 mm cylinder test， the JWL-Miller EOS parameters of the detonation product of NEPE propellant was calibrated by Gurney model and genetic algorithm. Finally， the pressure curves at different Lagrange position were fitted with the ignition and growth model. The results show that the fitting correlation coefficients of JWL EOS parameter curves for unreacted propellant and detonation products are high enough and the obtained ignition and growth model parameters well simulate shock initiation experimental results and the obtained parameters can provide reference for the safety evaluation of solid rocket motors.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
• Inhibition of Hydrogen Release of AlH₃ under HTPB System

  WANG Bang-qing, HU Da-shuang, JIN Guo-rui, XIAO Xu, LUO Cong, YANG Gen, CAO Rong, ZHAO Cheng-yuan, CHI Xu-hui

  2025,33(2):178-187, DOI: 10.11943/CJEM2024245

  Abstract(142) HTML(15) PDF(4.27 M)( 116)

  Abstract:In order to achieve the inhibition of hydrogen release from aluminum trihydride （AlH₃）， the hydroxyl polybutadiene （HTPB） was chosen as hydrogen absorber to undergo the hydrogenation reaction with unsaturated C═C double bonds under the catalytic action of CA catalyst， eliminating the hydrogen released by the decomposition of AlH₃ in time. The catalytic absorption characteristics of hydrogen released by AlH₃ under HTPB-CA system were studied by micron three-dimensional image microscopy （micron CT）， elemental analysis， thermogravimetry-differential thermal analysis （TG-DTA）， scanning electron microscopy （SEM）， X-ray diffraction （XRD） and other characterization methods. And the effect of hydrogen absorption on the decomposition of AlH₃ was further analyzed. It indicates that CA could catalyze the absorption of hydrogen released from AlH₃ in the HTPB system， inhibiting the generation of voids in AlH₃/HTPB grain which significantly prolongs its volume cracking time. The addition of CA extends the induction period AlH₃ decomposition in the AlH₃/HTPB at 70 ℃ from 6 days to 15 days. And the H content of AlH₃ of the AlH₃/HTPB/CA mixture is 8.94% after 203 days of aging at 60 ℃， while that H content of AlH₃ of the AlH₃/HTPB mixture is only 3.12%. It suggests the inhibition of the decomposition of AlH₃ after the hydrogen reacting with C═C double bond of HTPB to form C─H bond under the catalytic action of CA， demonstrating a hydrogen autocatalytic decomposition mechanism of AlH₃.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90
  • 91
  • 92
  • 93
  • 94
  • 95
  • 96
  • 97
  • 98
  • 99
  • 100
  • 101
  • 102
  • 103
  • 104
  • 105
  • 106
  • 107
  • 108
  • 109
  • 110
  • 111
  • 112
  • 113
  • 114
  • 115
  • 116
  • 117
  • 118
  • 119
  • 120
  • 121
  • 122
  • 123
  • 124
  • 125
  • 126
  • 127
  • 128
  • 129
  • 130
  • 131
  • 132
  • 133
  • 134
  • 135
  • 136
  • 137
  • 138
  • 139
  • 140
• Mechanical Properties and Combustion Start-Stop Performance of Semi-Interpenetrating-Network HAN Based Electrically Controlled Solid Propellants

  NAI Hui-min, LIU Fu-han, WANG Bin, XIA De-bin, ZHANG Jian, LIN Kai-feng, YANG Yu-lin

  2025,33(2):188-197, DOI: 10.11943/CJEM2024277

  Abstract(149) HTML(12) PDF(4.38 M)( 133)

  Abstract:To explore the effects of the degree of polymerization of polyvinyl alcohol （PVA） as a binder， the content of boronic acid （HB） as a crosslinking agent， polyaniline （PANI） and chitosan （CTS） as an additive on the mechanical properties and combustion start-stop performance of hydroxylamine nitrate based electrically controlled solid propellants （HAN-ECSP）， the samples were prepared by the swelling method. The mechanical properties and combustion start-stop performance were studied by tensile test， combustion start-stop performance test， X-ray diffraction （XRD）， Fourier transform infrared attenuated total reflection （FTIR-ATR）， and optical electron microscopy. The findings indicate that there is minimal variation in crystallinity among PVA samples with different degrees of polymerization. But when the polymerization degree of PVA is 2400±50， the ─OH content is relatively high and the mechanical properties are better. The HB content influences the crosslinking density of PVA. PVA₂₄₉₉-20% HB exhibits a relatively low density and uneven surface depressions. PVA₂₄₉₉-18% HB demonstrates superior mechanical properties， boasting a tensile strength of 2.36 MPa and a tensile elongation of 482.92%. However， it is highly corrosive to the copper electrode. Therefore， an initial network is formed by cross-linking 8% HB with PVA. The PANI or CTS is added to form a semi-interpenetrating polymer network with PVA. The tensile strength of the PVA₂₄₉₉-8% HB-2% PANI propellant is 1.03 MPa， and its elongation at break is 432.85%. The tensile strength of the PVA₂₄₉₉-8% HB-1% CTS propellant is 0.87 MPa， and its elongation at break is 302.28%. Both PANI and CTS addition to the ECSP can achieve combustion start‐stop， and the maximum burning temperature can reach 800 ℃.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64
  • 65
  • 66
  • 67
  • 68
  • 69
  • 70
  • 71
  • 72
  • 73
  • 74
  • 75
  • 76
  • 77
  • 78
  • 79
  • 80
  • 81
  • 82
  • 83
  • 84
  • 85
  • 86
  • 87
  • 88
  • 89
  • 90

>Reviews

• Research Progress on the Design of Flexible Chain Structures for Energetic Polyether Binders

  PENG Xiao-yue, WANG Yu, SUN Qian, WANG Ya, LI Ya-nan, ZHENG Wen-fang

  2025,33(2):198-211, DOI: 10.11943/CJEM2024243

  Abstract(174) HTML(15) PDF(1.71 M)( 113)

  Abstract:Energetic polyether binders， as the backbone and matrix of propellants， are the basis to improve the energy levels， mechanical properties， and processing properties of propellants. However， the polar energetic groups in polyethers hinder the movement of molecular chains and reduce the flexibility of polymer chains， leading to a decrease in mechanical properties and restricting the development of energetic solid propellants. Copolymerization modification of energetic polyethers and introduction of flexible structural units to improve flexibility are effective methods to obtain polyethers with diverse structures and adjustable properties. This paper summarizes different types of flexible chain segments introduced into energetic polyethers in recent years， and discusses from the perspectives of main-chain flexible chain segments and side-chain flexible chain segments to explain their effects on the mechanical properties and processing performance of energetic polyether binders. Future design of the flexible chain structure of energetic polyether binders is also discussed， which will provide a direction for the design and development of new types of energetic binders.

  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 32
  • 33
  • 34
  • 35
  • 36
  • 37
  • 38
  • 39
  • 40
  • 41
  • 42
  • 43
  • 44
  • 45
  • 46
  • 47
  • 48
  • 49
  • 50
  • 51
  • 52
  • 53
  • 54
  • 55
  • 56
  • 57
  • 58
  • 59
  • 60
  • 61
  • 62
  • 63
  • 64