基于自适应遗传算法的爆炸冲击响应谱时域重构优化方法

孙文娟 陈海波 黄颖青

孙文娟, 陈海波, 黄颖青. 基于自适应遗传算法的爆炸冲击响应谱时域重构优化方法[J]. 高压物理学报, 2019, 33(5): 052301. doi: 10.11858/gywlxb.20180681
引用本文: 孙文娟, 陈海波, 黄颖青. 基于自适应遗传算法的爆炸冲击响应谱时域重构优化方法[J]. 高压物理学报, 2019, 33(5): 052301. doi: 10.11858/gywlxb.20180681
SUN Wenjuan, CHEN Haibo, HUANG Yingqing. Time Domain Reconstruction Optimization of Pyrotechnic Shock ResponseSpectrum via Adaptive Genetic Algorithm[J]. Chinese Journal of High Pressure Physics, 2019, 33(5): 052301. doi: 10.11858/gywlxb.20180681
Citation: SUN Wenjuan, CHEN Haibo, HUANG Yingqing. Time Domain Reconstruction Optimization of Pyrotechnic Shock ResponseSpectrum via Adaptive Genetic Algorithm[J]. Chinese Journal of High Pressure Physics, 2019, 33(5): 052301. doi: 10.11858/gywlxb.20180681

基于自适应遗传算法的爆炸冲击响应谱时域重构优化方法

doi: 10.11858/gywlxb.20180681
基金项目: 国家自然科学基金(11772322);军委装备发展部预研领域基金(6140246030216ZK01001);中国科学院战略性先导科技专项(B类)子课题(XDB22040502)
详细信息
    作者简介:

    孙文娟(1986-),女,博士研究生,讲师,主要从事爆炸冲击效应研究. E-mail: sunwenj@mail.ustc.edu.cn

    通讯作者:

    陈海波(1968-),男,博士,教授,主要从事计算力学及工程应用、振动工程研究. E-mail: hbchen@ustc.edu.cn

  • 中图分类号: V415.4

Time Domain Reconstruction Optimization of Pyrotechnic Shock ResponseSpectrum via Adaptive Genetic Algorithm

  • 摘要: 为解决现有爆炸冲击响应谱(Shock Response Spectrum,SRS)加速度重构方法依赖于大量试验数据的问题,对比了阻尼正弦与小波两种不同加速度重构方法在合成爆炸冲击响应谱时的性能。将对重构SRS质量的评估转化为与目标谱匹配度的最小值优化问题,并首次将自适应遗传算法(Adaptive Genetic Algorithm, AGA)应用于SRS重构的优化问题中。对比了交叉先行、变异先行和不定向3种不同的AGA在爆炸冲击响应谱时域重构优化中的性能,并与基本遗传算法(Genetic Algorithm, GA)进行对比。结果表明,AGA的优化结果比GA有较大幅度的改善,且不定向AGA所得结果是3种AGA方法中最好的,其SRS各频点数值均在(–3/+6)dB容差范围之内,与目标谱的匹配度更好。仿真对比算例验证了该方法在冲击响应谱的时域重构应用中具有较高的准确性和实用性,为进一步提高航天器结构在爆炸冲击载荷下响应的计算精度提供了支撑。

     

  • Metal nitrides have attracted increasing attention in a wide range application due to their fascinating properties, such as hardness, superconductivity, various types of magnetism, etc.[12] Recently, extensive efforts have focused on the synthesis and design of metal polynitrides with single and double nitrogen bonds under high pressure, as metal polynitrides are proposed to be promissing candidates for high-energy density materials (HEDMs)[35]. The decomposition of these polynitrides is expected to release huge energy due to the conversion from a single/double-bonded polymeric nitrogen (160/418 kJ/mol) into a triple-bonded nitrogen (954 kJ/mol) molecule without producing pollutants. The atomic polymeric nitrogen with cubic gauche structure (cg-N) at a high pressure and high temperature (HPHT) condition (pressure above 110 GPa and temperature greater than 2 000 K) was successfully synthetized by Eremets et al.[68]. Motivated by that synthesis, alkali metal azides AN3, constructed by spherical cations and linear molecular N3 anions with N=N double bonds, was proposed to be suitable precursor in the formation of polymeric nitrogen under pressure[9]. Experimentally, the polymeric evolutionary behaviors of N3 anions in LiN3[10], NaN3[1113], KN3[1416], RbN3[1718], and CsN3[1920] have been fully investigated under high pressure. One of the most interesting findings is the reported polymeric nitrogen nets in NaN3[11] at 160 GPa. More recently, a fully sp2-hybridized layered polymeric nitrogen structure, featuring fused 18-membered rings, in potassium supernitride (K2N16) was successfully synthesized under HPHT condition using a laser heating diamond anvil cell (DAC) technique[21]. In our previous works[2225], the N3 anions in LiN3, NaN3, and KN3 all transform firstly to “N6” molecular clusters under compression, and then to a polymerized nitrogen phase at high pressures.

    Besides the well-known alkali metal azides AN3, alkali metal diazenides A2N2 containing the N=N double bonds, were also proposed as precursors of the polymeric atomic nitrogen under pressure. In 2010, two alkali metal diazenides, Li2N2 and Na2N2, were proposed to have orthorhombic Pmmm structure[26]. However, an orthorhombic Immm structure of Li2N2[27] was then synthesized under HPHT condition. Immediately after that, a serial pressure-induced structural phase transitions of Li2N2 were reported in two independent works[2829]. Moreover, a high-pressure tetragonal I41/acd phase containing the spiral nitrogen chain was uncovered at 242 GPa[29], indicating that the polymerization of nitrogen is realized from Li2N2.

    Compared to Li2N2 and Na2N2, K2N2 have not been reported so far, neither for ambient pressure nor for high pressure. Here, we will extensively explore the high-pressure structures of K2N2 up to 150 GPa by using the first principles swarm structure searching method. The structure evolutions and chemical bonding behaviors of K2N2 within different structures were then fully studied to provide an insight into the formation of polymeric nitrogen in alkali metal diazenides.

    The variable-cell crystal structure predictions for K2N2 (1–4 formula unit (f. u.) in the simulation cell) up to 150 GPa were conducted by CALYPSO code[3031]. The effectiveness of this method has been fully confirmed in predicting high-pressure structures of various systems[3235]. During the structure search process, the 60% structures of each generation with lower enthalpies were selected to generate the structures for the next generation by particle swarm optimization (PSO) technique, and the other structures in new generation were randomly generated to increase the structural diversity. Usually, the structure searching simulation stopped at the 20th generation and 40 structures per generation were generated at each pressure point. After finishing the structure search, the most stable structures were used for further structural relaxation and property calculation, a process that is implemented in the Vienna ab Initio Simulation Package[36]. The Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional under the generalized gradient approximation[37] and the projector-augmented wave (PAW)[38] for ion-electron potentials were adopted herein. The van-der-Waals (vdW) correction was applied in the structural relaxation using Grimme’s method (DFT-D3)[39]. The cutoff energy of 550 eV and Monkhorst-Pack k meshes[40] of 2π×0.032 Å−1 were used to ensure the good convergences of total energy and maximum force on each atom. The lattice dynamic stability of each structure was checked by phonon spectrum obtained by PHONOPY code[41]. The charge transfer behavior of K2N2 was performed using Bader atoms-in-molecules (AIM) method[42].

    The ground-state structure for K2N2 at ambient pressure in our fully structural search (in Fig. 1(a)) is a monoclinic structure with the space group C2/m, which is different from those reported for Li2N2 and Na2N2[26]. In Fig. 1(a), the N2 quasimolecule dimers in C2/m structure are parallel to each other, and are tilted by a certain angle within a-b plane. The N-N bond length is 1.192 Å, which is a typical value for double bond. Also, the calculated Mayer bond order[43] for the N2 anion in the C2/m is 2.18, being consistent with double bond. Due to the weak interaction between the N2 quasimolecules in the C2/m phase, the application of gigapascal pressure leads to the orientation change of N2 quasimolecules and the structure transformation. Indeed, two stable orthorhombic structures, i.e. Pmmm and Cmmm, for K2N2 were uncovered at 2 and 5 GPa, respectively, as shown in Fig. 1(b) and Fig. 1(c). Remarkably, the unveiled Pmmm phase and the previously reported ground-state structure of Na2N2 are exactly isostructrual (the Pmmm phase denoted as the Na2N2-type phase hereafter), and the Cmmm phase of K2N2 and the high-pressure metastable structures of Na2N2 are isostructrual as well[26]. Furthermore, both Na2N2-type and Cmmm orthorhombic phases have the similar atomic arrangement that the parallel N2 dimers are perpendicular to the K atomic layer. In the Na2N2-type and Cmmm phases, the calculated N-N bond lengths are 1.220 and 1.254 Å, while the characteristics of the N=N double bond are as that of the C2/m structure in the ground-state. As the pressures rises to 150 GPa, a stable monoclinic C2/c structure, containing a folded one-dimensional nitrogen chains (see Fig. 1(d)), is identified for K2N2 for the first time. The lengths of two nonequivalent N-N bond in C2/c structure are 1.474 and 1.373 Å at 135 GPa, indicating genuine N―N single bond, and thus the nature of the structure is similar to that of the high-pressure phase of NaN3[24]. Table 1 lists the structural details of the lattice parameters and the atomic positions of various K2N2 structures at selected pressure points.

    Figure  1.  Crystal structure of C2/m, Na2N2-type, Cmmm, and C2/c phases (The large and small spheres represent K and N atoms, respectively.)
    Table  1.  Optimized structural parameters of the C2/m, Na2N2-type, Cmmm, and C2/c phases of K2N2
    Phase Pressure/GPa Lattice parameters dN-N Atomic fractional coordinates
    C2/m0a = 7.562 Å, b = 3.912 Å, c = 10.923 Å,
    α = γ = 90°, β= 134.955°
    1.192K 4i (0.613, 0, 0.260)
    N 4i (0.519, 0, 0.456)
    Na2N2-type2.5a = 3.326 Å, b = 4.375 Å, c = 5.694 Å,
    α = β = γ = 90°
    1.220K1 1e (0, 0.500, 0)
    K2 1c (0, 0, 0.500)
    N 2t (0.500, 0.500, 0.393)
    Cmmm20.0a = 6.896 Å, b = 5.104 Å, c = 2.855 Å,
    α = β = γ = 90°
    1.254K 4g (0.694, 0, 0)
    N 4j (0, 0.877, 0.500)
    C2/c135.0a = 6.971 Å, b = 4.099 Å, c = 4.510 Å,
    α = γ = 90°, β= 81.342°
    1.474, 1.373K 8f (0.669, 0.909, 0.198)
    N 8f (0.041, 0.887, 0.102)
     | Show Table
    DownLoad: CSV

    The phonon dispersion curves of the four predicted structures are illustrated in Fig. 2 in order to check their dynamic stabilities. Fig. 2 shows that all the K2N2 phases are dynamically stable as the eigen frequencies of their lattice vibrations are positive for all wavevectors in the Brillouin zone. It should be noted that the predicted C2/m, Na2N2-type, and Cmmm phases exhibit higher eigen frequencies of around 50 THz, which correspond to the vibrations of N2 quasimolecules with stronger N=N double bonds. By contrast, the C2/c phase shows relative lower vibration frequencies of around 30−35 THz, corresponding to the vibration frequencies of weaker N―N single bonds.

    Figure  2.  Calculated phonon curves of the C2/m (a), Na2N2-type (b), Cmmm (c), and C2/c (d) phases at selected pressure points

    Herein, we evaluated the energy involved in the decomposition of the C2/c phase with polymerized N―N single bond chains into metal K and gaseous N2. The energy involved in the C2/c-K2N2 decompostion was determined to be 0.344 kJ/g, which is about three times of that of conventional energetic materials (0.104−0.105 kJ/g)[44], but much smaller than those reported for high-energy density materials (3.61 kJ/g for BeN4, 5.30 kJ/g for GaN5, 4.81 kJ/g for GaN10)[4546]. This may be due to the low nitrogen content in K2N2, so the nitrogen-rich KNx (x> 1) is expected to be potential HEDM under high pressure.

    To determine the phase transition pressure of K2N2, we plotted in Figs. 3(a) and 3(b) the enthalpy difference of the predicted C2/m, Cmmm, and C2/c structures relative to the Na2N2-type phase in the pressure range of 0−150 GPa. The inset of Fig. 3(a) shows that above 1.7 GPa the Na2N2-type phase becomes energetically preferable with respect to the ground-state C2/m structure and is enthalpically stable up to about 3.6 GPa, above which it transforms into the Cmmm structure. From Fig. 3(b), it can be seen that the Cmmm structure persists up to 122 GPa and then transforms to the C2/c structure according to our structural predictions and enthalpy difference calculations. The wide pressure range of stable Cmmm phase may be attributed to the relative smaller ionic radius and atomic mass of K atom,in contrast with the heavier Rb and Cs atoms. It is thus expected that the chemical precompression exerted by heavier atoms, such as Rb and Cs, in A2N2 (as for Cs2N2[43]) can yield lower pressures for the polymerization of one-dimension N chains. From the calculated pressure dependence of the volume per f. u. in Figs. 3(c) and 3(d), the volume collapses of the first-order transformations C2/m→Na2N2-type, Na2N2-type→Cmmm, and CmmmC2/c are 14.4%, 22.5%, and 4.0% respectively. Thus, the Cmmm phase is much stiffer than low-pressure C2/m and Na2N2-type phases under compression. This is consistent with the increasing coordination number of the K atom in the different K2N2 phases under compression in Figs. 3(c) and 3(d). The densification of K2N2 under pressure is realized by increasing the cation coordination number from 5 in ground-state C2/m phase to 8 in Na2N2-type (K2 atoms are coordinated by eight N atoms and K1 atoms are fourfold coordinated by four N atoms) and Cmmm phases, and finally to 10 in C2/c phase.

    Figure  3.  (a)–(b) Enthalpy differences of the different predicted structures relative to the Na2N2-type structure under pressure;(c)–(d) pressure dependence of the volume per f. u. of each structure for K2N2

    To get deeper insight into the electronic structure evolutions of these K2N2 phases under pressure, we calculated the total density of states (DOS) and the site projected DOS, as shown in Fig. 4, where the Fermi level is the vertical dotted line. In contrast with the semiconducting nature of C2/c phase (Fig. 4(d)), the C2/m, Na2N2-type, and Cmmm phases are all metallic with the N-p electron crossing the Fermi level. Thus, the K2N2 possesses a metal-semiconductor transition, that is similar to that of Li-N system in recent work[28]. From Fig. 4(a)–Fig. 4(c), the C2/m, Na2N2-type, and Cmmm phases exhibit similar projected DOS profiles in the valence band region, the K-s and K-p atomic orbitals are located at the deep energy levels (around −32 and −16 eV, respectively), and the N-s atomic orbitals mainly lie around −22 − −10 eV. Moreover, the interaction between K and N atoms is extremely weak due to the absence of overlapping regions in their atomic projected DOSs. From the inspection of the DOS profiles, a clear broadening of different atomic decomposed states is seen in the high-pressure C2/c phase (Fig. 4(d)), leading to clear hybridizations between K-p and N-p atomic orbitals. In contrast to the first three low-pressure K2N2 phases, this distinction is partially due to the increasing coordination number of K atoms in C2/c phase and to the resulting chemical bonding change that will be illustrated in the following Bader charge transfer analysis.

    Figure  4.  Total and site projected DOSs of C2/m (a), Na2N2-type (b), Cmmm (c), and C2/c (d) phases

    Fig. 5 presents the PBE-resulted energy band structures of the four K2N2 phases, the projected weights of K-s, K-p, and N-s orbital electrons are not taken into account herein because of their exceptionally tiny contributions around the Fermi level in Fig. 4. Compared to the strong metallic nature of the Cmmm phase, which is mainly dominated by N-px/N-pz at the considered energy range of −2 − 2 eV (Fig. 5(c)), the C2/m ( Na2N2-type) phase exhibits a relative weaker metallic behavior that two bands of N-py/N-pz (N-px/N-py) crossing the Fermi level in Fig. 5(a) (Fig. 5(b)). The indirect semiconducting nature of C2/c phase (band gap of 2.0 eV) is demonstrated in Fig. 5(d), with the valence band maximum (VBM) and conduction band minimum (CBM) is located at A and E points, respectively. It is further found that the VBM and CBM of the C2/c phase are principally composed of N-px and N-pz orbitals, respectively.

    Figure  5.  Projected weights of N-p orbitals in the band structures of C2/m (a), Na2N2-type (b), Cmmm (c), and C2/c phases (d) (The Fermi level is indicated by horizontal lines and the black solid lines denote the energy structures of each phase of K2N2)

    Fig. 6 shows the projected two-dimensional (2D) electron localization functions (ELF) of the four K2N2 phases to further reveal the atomic chemical bonding. The ELF is used to characterize the localized distribution of electrons of different atoms. The ELF values of 1, 0.5, and close to 0 indicate strong electron localization, electron gas, and nonelectron localization, respectively. From Fig. 6(a)–Fig. 6(c), the high-ELF regions are distributed between two N atoms, indicating the strong covalent N-N bonding in the C2/m, Na2N2-type, and Cmmm phases. Lone-pair electrons of K remain from the C2/m to Cmmm phases under pressure, as illustrated by the ELF value of 0 around K atoms, which is in conformity with the projected atomic orbital DOSs described in Fig. 4(a)–Fig. 4(c). For the compressed C2/c phase, some lone-pair electrons of K still remain, yet some of them are transferred to the regions of the K-N bond. Hence, the pressure compels the K-p lone-pair electrons to participate in the chemical bonding with N-p electrons under compression, realizing the high pressure C2/c phase with increased coordination number of K atoms. The electron localization on N-N bond in the C2/c phase (Fig. 6(d)) is weaker than that in the other three low-pressure K2N2 phases, confirming that the N―N single bond of the polymerized one-dimensional nitrogen chains is weaker than the N=N double bond of the N2 quasimolecules.

    Figure  6.  Contours of ELF for the C2/m on the (010) plane (a), Na2N2-type on the (010) plane (b), Cmmm on the (100) plane (c), and C2/c on the (010) plane (d)

    Finally, the Bader charges in the four phases are listed in Table 2 by using the Bader AIM method. For the ground-state C2/m phase, each N atom accepts 0.491e which corresponds to the electron loss per K atom. When phase transits into the Na2N2-type and Cmmm phases, a relatively stronger ionic bonding feature of K-N bond is disclosed because that the charge transfer is raised to approximately 0.645e − 0.690e as a result of the electron gain of N. Contrarily, it is noted that the atomic charge transfer from the K atom to N atom decrease to 0.623e in the compressed C2/c phase. This variation may be originated from that partial K atoms participate in covalent bonding with the N atom under compression, as indicated by the orbital hybridizations of K-p and N-p in Fig. 4(d). It then can be concluded that the chemical bonding in the C2/c phase containing one-dimensional polymerized nitrogen chains is a complex mixture of covalent and ionic characteristics.

    Table  2.  Calculated Bader charges of K and N atoms in C2/m, Na2N2-type, Cmmm, and C2/c phases
    Phase Pressure/GPa Atom Charge value/e Charge transfer/e
    C2/m 0 K
    N
    8.509 (×2)
    5.491 (×2)
    +0.491
    −0.491
    Na2N2-type 2.5 K1
    K2
    N
    8.407
    8.303
    5.645 (×2)
    +0.593
    +0.697
    −0.645
    Cmmm 20.0 K
    N
    8.310 (×2)
    5.690 (×2)
    +0.690
    −0.690
    C2/c 135.0 K
    N
    8.377 (×2)
    5.623 (×2)
    +0.623
    −0.623
     | Show Table
    DownLoad: CSV

    In summary, we have performed a systematic crystal structure search for alkali metal diazenide K2N2 under pressure using the advanced swarm intelligence structure simulations. The ground-state phase of K2N2 is determined to be a new monoclinic C2/m structure with N2 quasimolecules, and three high-pressure structures of K2N2, including Na2N2-type (1.7−3.6 GPa), Cmmm (3.6−122.0 GPa), and C2/c (above 122.0 GPa), are identified. Remarkably, the pressure-induced polymerization of one-dimensional nitrogen chains in the high-pressure C2/c semiconducting phase is firstly reported, attributed to the participation of K-p lone-pair electrons in the chemical bonding under pressure. The occurrences of the high-pressure phases of Na2N2-type, Cmmm, and C2/c follow the increased coordination numbers (from 5 to 10) of K atoms with the incresed pressures, accompanied by significant N-N bonding modification from N=N dimers to polymerized N―N single bond one-dimensional chains. Furthermore, a metal-semiconductor transition of K2N2 under pressure is demonstrated by the electronic structures calculations. The present findings provide fundamental insights into the high-pressure structures and the polymerization of N atoms of alkali metal diazenides, supplying useful guidance for the syntheses of K2N2 in further experimental efforts.

  • 图  冲击响应谱概念示意图

    Figure  1.  Graphical representation of the shock response spectrum

    图  阻尼正弦与小波合成SRS计算结果对比

    Figure  2.  Comparison results of damped sine and wavelet

    图  交叉概率和变异概率取值分析结果

    Figure  3.  Average value of Pc and Pm

    图  参数灵敏度分析前后远场SRS结果对比

    Figure  4.  Comparison results of empirical parameters and optimized parameters for far-field SRS

    图  参数灵敏度分析前后中场和近场SRS结果对比

    Figure  5.  Comparison results of empirical parameters and optimized parameters for mid-field and near-field SRS

    图  GA与AGA优化结果对比

    Figure  6.  Comparison results of GA and AGA

    图  加速度时间历程曲线

    Figure  7.  Acceleration time-history curves

    图  不定向AGA优化结果与文献[19]结果对比

    Figure  8.  Comparison results of uncertain-order AGA and Ref. [19]

    图  不同种群数目下中场目标SRS结果对比

    Figure  9.  Comparison results of mid-field SRS under different population numbers

    表  1  决策变量的取值范围

    Table  1.   Variation ranges of the decision variables

    Optimization variableVariation range
    Am(1/4 to 1/3)A0 (g)
    tdm[0.0001, 0.015] (s)
    ξm[0.001, 0.1]
    Nm[5, 27] (odd number)
    下载: 导出CSV

    表  2  典型爆炸冲击响应谱规范

    Table  2.   Specification of SRS

    Far fieldMedium fieldNear field
    Frequency/HzAmplitude/gFrequency/HzAmplitude/gFrequency/HzAmplitude/g
    100 80 100 150 200 250
    450 600 300 200 1 0004 000
    9001 000 1 5003 000 1 2005 000
    10 0001 00010 0003 000 10 0005 000
    下载: 导出CSV

    表  3  阻尼正弦与小波合成SRS计算结果对比

    Table  3.   Comparison results of damped sine and wavelet

    ParameterFar fieldMedium fieldNear field
    Damped sineWaveletDamped sineWaveletDamped sineWavelet
    Objective function value/g 44.2294.9103.7890.1155.51491
    Time/s142.2142.1139.7139.2 81.9 80.2
    下载: 导出CSV

    表  4  AGA选用参数

    Table  4.   Parameters of AGA

    ParameterValueParameterValueParameterValue
    Population40Pm10.1 Pc00.75
    Maximum evolutionary generation200Pm20.05 PLBm0.01
    Pc10.9Pm30.005PUBm0.1
    Pc20.5PLBc0.5 Pm00.05
    Pc30.1PUBc0.9
    下载: 导出CSV

    表  5  GA与AGA优化结果对比

    Table  5.   Comparison results of GA and AGA (OFV: objective function value)

    AlgorithmFar fieldMedium fieldNear field
    OFV/gCurrent generationTotal time/sOFV/gCurrent generationTotal time/sOFV/gCurrent generationTotal Time/s
    GA44.20200142.2103.7200139.7155.5200 81.9
    Crossover first AGA44.05115147.1102.9102135.1156.8127 82.9
    Mutation first AGA44.29128139.8103.9137136.1155.4158 83.6
    Uncertain-order AGA44.2589172.1103.187143.0155.398102.2
    下载: 导出CSV

    表  6  不同种群数目下中场目标SRS优化值与计算时间对比

    Table  6.   Comparison of optimization values and calculation time for mid-field SRS under different population numbers

    PopulationFinal optimization value/gCalculation time of
    200 generations/s
    2077.94 53.34
    4048.48139.70
    6044.53148.80
    8037.93159.70
    10027.88221.80
    下载: 导出CSV
  • [1] LEE J R, CHIA C C, KONG C W. Review of pyroshock wave measurement and simulation for space systems [J]. Measurement, 2012, 45(4): 631–642. doi: 10.1016/j.measurement.2011.12.011
    [2] HIMELBLAU H, KERN D L, MANNING J E, et al. Dynamic environment critical: NASA-HDBK-7005 [R]. Washington, DC: National Aeronautics and Space Administration, 2001.
    [3] ECSS. Space engineering-mechanical shock design and verification handbook: ECSS-E-HB-32-25A [R]. The Netherlands: European Space Agency, 2015.
    [4] ALEXANDER J E. A new method to synthesize a shock response spectrum compatible base acceleration to improve multi-degree of freedom system response [D]. Minneapolis, MN: University of Minnesota, 2015.
    [5] HWANG J H, DURAN A. Stochastic shock response spectrum decomposition method based on probabilistic definitions of temporal peak acceleration, spectral energy, and phase lag distributions of mechanical impact pyrotechnic shock test data [J]. Mechanical Systems and Signal Processing, 2016, 76: 424–440.
    [6] CHONG S Y, LEE J R, KONG C W. Shock response spectra reconstruction of pointwise explosive-induced pyroshock based on signal processing of laser shocks [J]. Shock and Vibration, 2014, 2014: 1–14.
    [7] 杜志鹏, 汪玉, 杨洋, 等. 舰艇水下爆炸冲击信号拟合及应用 [J]. 振动与冲击, 2010, 29(3): 182–184. doi: 10.3969/j.issn.1000-3835.2010.03.044

    DU Z P, WANG Y, YANG Y, et al. Curve fit method for naval underwater explosion shock and its application [J]. Journal of Vibration and Shock, 2010, 29(3): 182–184. doi: 10.3969/j.issn.1000-3835.2010.03.044
    [8] 马道远, 庄方方, 徐振亮. 基于遗传算法的冲击响应谱时域合成方法 [J]. 强度与环境, 2015, 42(5): 49–53.

    MA D Y, ZHUANG F F, XU Z L. Time-domain synthesis method for shock response spectrum based on genetic algorithm [J]. Structure and Enviroment Engineering, 2015, 42(5): 49–53.
    [9] YE Z, LI Z, XIE M. Some improvements on adaptive genetic algorithms for reliability-related applications [J]. Reliability Engineering & System Safety, 2010, 95(2): 120–126.
    [10] SRINIVAS M, PATNAIK L M. Adaptive probabilities of crossover and mutation in genetic algorithms [J]. IEEE Transactions on Systems, Man, and Cybernetics, 1994, 24(4): 656–667. doi: 10.1109/21.286385
    [11] YAN M, HU H, OTAKE Y, et al. Improved adaptive genetic algorithm with sparsity constraint applied to thermal neutron CT reconstruction of two-phase flow [J]. Measurement Science and Technology, 2018, 29(5): 1–14.
    [12] BIOT M A. Transient oscillations in elastic systems [D]. Pasadena, CA: California Institute of Technology, 1932.
    [13] LI B W, LI Q M. Damage boundary of structural components under shock environment [J]. International Journal of Impact Engineering, 2018, 118: 67–77. doi: 10.1016/j.ijimpeng.2018.04.002
    [14] SMALLWOOD D O. An improved recursive formula for calculating shock response spectra [J]. Shock and Vibration Bulletin, 1981, 51(2): 211–217.
    [15] ISO. Mechanical vibration and shock—signal processing—part 4: shock response spectrum analysis: ISO/WD18431-4-2007 [S]. Geneva, Switzerland: ISO, 2007.
    [16] SMALLWOOD D O. A family of transients suitable for reproduction on a shaker based on the cos m(x) window [J]. Journal of the Institute of Environmental Sciences and Technology, 2002, 45(1): 178–184.
    [17] SIAM N. Development of an efficient analysis method for prediction and structural dimensioning of space structures subjected to shock loading [D]. Luleå, Sweden: Luleå University of Technology, 2010.
    [18] DILHAN D, CIPOLLA V, GRZESKOWIAK H, et al. Pyroshock generation [C]//European Conference on Spacecraft Structures, Materials and Mechanical Testing. The Netherlands, 2005: 1–10.
    [19] MONTI R, PAOLO G. Dynamic load synthesis for shock numerical simulation in space structure design [J]. Acta Astronautica, 2017, 137: 222–231. doi: 10.1016/j.actaastro.2017.04.023
    [20] BAI X S, YAN W S, GE S S, et al. An integrated multi-population genetic algorithm for multi-vehicle task assignment in a drift field [J]. Information Sciences, 2018, 453: 227–238. doi: 10.1016/j.ins.2018.04.044
  • 加载中
图(9) / 表(6)
计量
  • 文章访问数:  12822
  • HTML全文浏览量:  3500
  • PDF下载量:  28
出版历程
  • 收稿日期:  2018-11-09
  • 修回日期:  2018-11-27

目录

/

返回文章
返回