CN112102892A - Method for determining temperature correction coefficient of energetic material chemical combination process - Google Patents

Method for determining temperature correction coefficient of energetic material chemical combination process Download PDF

Info

Publication number
CN112102892A
CN112102892A CN202010827564.8A CN202010827564A CN112102892A CN 112102892 A CN112102892 A CN 112102892A CN 202010827564 A CN202010827564 A CN 202010827564A CN 112102892 A CN112102892 A CN 112102892A
Authority
CN
China
Prior art keywords
temperature
correction coefficient
temperature correction
tmr
formula
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010827564.8A
Other languages
Chinese (zh)
Other versions
CN112102892B (en
Inventor
丁黎
禄旭
华栋
黄蒙
祝艳龙
安静
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Modern Chemistry Research Institute
Original Assignee
Xian Modern Chemistry Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Modern Chemistry Research Institute filed Critical Xian Modern Chemistry Research Institute
Priority to CN202010827564.8A priority Critical patent/CN112102892B/en
Publication of CN112102892A publication Critical patent/CN112102892A/en
Application granted granted Critical
Publication of CN112102892B publication Critical patent/CN112102892B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10—Analysis or design of chemical reactions, syntheses or processes
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation

Landscapes

  • Theoretical Computer Science (AREA)
  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

The invention discloses a method for determining a temperature correction coefficient of a chemical combination process of an energetic material, which obtains a thermal decomposition reaction kinetic parameter according to the thermal decomposition characteristic quantity of the material in the chemical combination process of the energetic material, further obtains the process temperature required by the maximum reaction rate reaching time of adiabatic conditions with different danger levels, and quantitatively obtains the temperature correction coefficient gamma by taking the process temperatures with different danger levels and the specific operation time as scales1. The method has universality, is suitable for evaluating the thermal stability of different energetic materials, and has the advantages of comprehensive and reliable determination of the temperature correction coefficient, strong operability, conservative result of the determined temperature correction coefficient and high safety use reliability.

Description

Method for determining temperature correction coefficient of energetic material chemical combination process
Technical Field
The invention belongs to the technical field of risk assessment of energetic materials, and mainly relates to a method for determining a temperature correction coefficient of a chemical combination process of an energetic material.
Background
The energetic material has the characteristics of easy combustion and explosion under the stimulation action of heat, static electricity, machinery and the like, and the development of the combustion and explosion risk evaluation of the combination process of the energetic material is an important content of safety production design and management and control. Materials in the energetic material combination process are poor thermal conductors, are extremely sensitive to thermal stimulation, can generate autocatalytic decomposition exothermic reaction when heated, are extremely easy to generate heat accumulation to cause thermal explosion, and are the inducing factors of a large number of safety accidents, so that the combustion explosion risk evaluation in the process needs to quantitatively determine the temperature correction coefficient in the process. The invention quantitatively determines the temperature correction coefficient in the chemical combination process aiming at energetic materials.
In order to qualitatively know and quantitatively estimate the risk of combustion and explosion, a quantitative assessment method for major hazard sources of fire and explosive and product enterprises (BZA for short) is made in China according to the characteristics of weapon fire and explosive and ammunition enterprises and by referring to the assessment experience abroad-1Method). BZA-1The hazard source evaluation mathematical model and the physical significance proposed by the method are as follows:
H＝Hinner part+HOuter cover(1)
In the formula: h is the actual risk of the explosion hazard system, HInner partIs the real risk level in the system, HOuter coverThe risk degree is the real risk degree outside the system.
Wherein HInner partThe calculation method of (2) is as follows:
Hinner part＝V+KB (2)
In the formula: v is a material risk coefficient, namely the inherent static risk degree of the explosive, K is a controllable risk behavior controlled degree coefficient in the system, and B is the controllable risk degree in the system.
The controllable risk degree B in the system is calculated by the following formula:
B＝WB·D·P (3)
in the formula: wBThe calculation formula of the material energy danger coefficient of the explosive and the device thereof is shown in the specification.
WB＝V·γ (4)
Wherein gamma is a process risk coefficient, and the calculation formula is as follows:
γ＝γ1+γ2+γ3+γ4+γ5(5)
in the formula: gamma ray1Is a temperature correction coefficient, which is positively correlated with the temperature rise, gamma2For correction of coefficients by chemical agents, gamma3As a coefficient of compression correction, gamma4Correction of the coefficient for other mechanical effects (e.g. cutting, drilling, sawing, abrading, etc.), gamma5Is the static electricity correction coefficient.
Temperature correction coefficient gamma in the current technological process1The values are as follows (see Table 1)
TABLE 1 Process hazard coefficient gamma value conditions
Figure BDA0002636774750000021
As can be seen from Table 1, the temperature correction coefficient γ1The determination method is a semi-quantitative method, although corresponding evaluation parameter values exist, the actual operability is poor, the values cannot be accurately quantized one by one, the concept of the value condition is ambiguous, the quantitative criterion is lacked, and the risk degree of the process material is difficult to evaluate by effectively combining with the actual process condition.
Because the materials in the energetic material combination process are energetic materials and can generate decomposition exothermic reaction under the action of heat to cause combustion explosion accidents, whether the materials can keep stable after being stimulated by the process temperature in the process is one of important indexes of process dangerousness, accurate evaluation of thermal stability is an important way for realizing intrinsic process safety, scientific basis is provided for quantitative evaluation of dangerousness of the whole process production line and establishment of early warning and prevention and control measures, and therefore a method for determining the temperature correction coefficient in the energetic material combination process is urgently needed.
Disclosure of Invention
Correction coefficient gamma for existing temperature1The invention provides a method capable of correcting a temperature coefficient gamma1Determination of quantitative values, to achieveThe invention adopts the following technical scheme:
the method for determining the temperature correction coefficient of the energetic material chemical combination process is characterized in that a differential scanning calorimetry is utilized to test the material thermal decomposition reaction to obtain a kinetic parameter apparent activation energy EaAnd a pre-factor A is pointed, so as to obtain the process temperature T required by the reaching time of the maximum reaction rate of adiabatic conditions with different danger levelspAnd process temperature T in different hazard classespAnd the specific operation time coefficient is used as a scale for comparing and measuring the actual process temperature T of the materials in the processpsAnd actual process temperature time tsQuantitatively obtaining the temperature correction coefficient gamma1. The flow chart for determining the temperature correction coefficient is shown in fig. 1, and specifically comprises the following steps:
the method comprises the following steps: differential scanning calorimetry test for obtaining apparent activation energy E of material thermal decomposition reaction kinetic parameters in energetic material combination processaAnd pre-finger factor a.
Obtaining the thermal decomposition reaction of materials in the energetic material combination process by using a differential scanning calorimetry test, obtaining the peak temperature of the exothermic reaction causing the rapid decomposition of the materials, and obtaining the kinetic parameter apparent activation energy E of the thermal decomposition reaction of the materials by using the Kissinger equation formula (1) regressionaAnd pre-finger factor a:
Figure BDA0002636774750000031
in the formula, betaiThe heating rate is K/s; t is decomposition peak temperature, K; a is a pre-exponential factor, s-1；EaIs apparent activation energy, J/mol; r is a gas constant, J/(mol. K). By
Figure BDA0002636774750000032
To pair
Figure BDA0002636774750000033
Drawing to obtain a straight line, and calculating E from the slope of the straight lineaA is determined from the intercept.
Step two: is obtained byMaximum reaction rate arrival time TMR under adiabatic condition of same danger classadRequired process temperature Tp。
TMRadThe formula (II) is shown below:
Figure BDA0002636774750000041
in the formula, EaJ/mol as activation energy; a is a pre-exponential factor, s-1(ii) a T is the process temperature, K; r is a gas constant, J/(mol. K); c. CpThe specific heat capacity of the material is J/(g.K); qrThe exothermic amount of the thermal decomposition of the material is J/g.
TMR adopts 6 grade criterion proposed by Zurich Hazard Analysis as judgment conditionadRespectively 100h, 50h, 24h, 12h, 8h and 1h, and respectively obtaining up to 6 TMRs by formula (II)adRequired process temperature TpSequentially comprises the following steps: t isD100、TD50、TD24、TD12、TD8、TD1。
Step three: maximum reaction rate arrival time TMR at 6 adiabatic conditions of different hazard classesadTemperature T of the dividing processpAt a process temperature TpFor a scale, the actual process temperature T of the materials is compared and measuredpsQuantitatively obtaining the temperature correction coefficient gamma1。
Due to TMRadProcess temperature T required for 100hpIs TD100Thus, for TMRadNot less than 100h, when the actual process temperature T of the material isps≤TD100Time, temperature correction coefficient gamma1＝1；
Due to TMR ad50h required process temperature TpIs TD50Thus, for a TMR of 50h ≦adLess than 100h, when the actual process temperature T of the materialD100＜Tps≤TD50Time, temperature correction coefficient gamma1＝2；
Due to TMR ad24h required process temperatureTpIs TD24Thus, for a TMR of 24h ≦adLess than 50h, when the actual process temperature T of the materialD50≤Tps＜TD24Time, temperature correction coefficient gamma1＝3；
Due to TMRadProcess temperature T required for 12hpIs TD12Thus, for a TMR of 12h ≦adLess than 24h, when the actual process temperature T of the material isD24＜Tps≤TD12Time, temperature correction coefficient gamma1＝4；
Due to TMRadProcess temperature T required for 8hpIs TD8Thus, for 8h ≦ TMRadLess than 12h, when the actual process temperature T of the material isD12＜Tps≤TD8Time, temperature correction coefficient gamma1＝5；
Due to TMR ad1h required process temperature TpIs TD1Thus:
TMR for 1hadLess than 8h, when the actual process temperature T of the material isD8＜Tps≤TD1Time, temperature correction coefficient gamma1＝6～9；
For TMRadLess than or equal to 1h, when the actual process temperature T of the material isps≥TD1Time, temperature correction coefficient gamma1＝10。
Temperature correction coefficient gamma1As shown in table 2:
TABLE 2 temperature correction coefficient values
Figure BDA0002636774750000051
Step four: temperature correction coefficient gamma1When the value range is 6-9, namely the actual process temperature TD8≤Tps＜TD1Then, the specific operation time coefficients of different danger levels are taken as a scale, and the temperature correction coefficient gamma is further quantitatively obtained1。
γ1When the value range is 6-9, TMR is utilizedadThe operation time is further consistent with the operation time in the actual process of the materialsStep (5), calculating a specific operation time coefficient by the formula (III):
Figure BDA0002636774750000052
in the formula: TMRadThe time for the maximum reaction rate to reach adiabatic conditions at the process temperature, h; t is the operation time of the material in the process, h. The physical meaning of the specific operating time coefficient is the safety margin of the process operating time when gamma is1In the case of 6 to 9, the temperature correction coefficient γ is determined based on 1.5 times, 2 times, and 4 times of safety margin, i.e., the specific operating time coefficient is 1, 2.5, and 41The determination schemes at 6-9 are shown in Table 3.
TABLE 3 temperature correction coefficient γ1Value scheme at 6-9
Figure BDA0002636774750000061
Compared with the prior art, the method for determining the temperature correction coefficient of the energetic material chemical combination process has the following beneficial effects:
(1) temperature correction coefficient gamma of the invention1The process temperatures of different danger levels are used as a measuring scale for determining, the process temperatures are obtained from thermal decomposition kinetic parameters based on the intrinsic thermal decomposition characteristics of materials in the energetic material combination process, the material attributes of different process production lines are reflected by combining process conditions, and the method has universality and is suitable for evaluating the thermal stability of different energetic material materials;
(2) temperature correction coefficient gamma of the invention1Determining the process and specific operation time coefficient required by the maximum reaction rate reaching time under the adiabatic conditions with different danger levels within the range of 6-9 as a scale, wherein the specific operation time coefficient scale is determined based on the safety margin of the process operation time, and the actual process temperature and the actual operation time are evaluated, so that the determination of the temperature correction coefficient is comprehensive and reliable, and the operability is high;
(3) quantitative determination of temperatureCorrection coefficient gamma1Providing scientific basis for the quantitative risk assessment and the establishment of early warning and prevention and control measures of the whole process production line;
(4) the invention adopts TMR (maximum reaction rate arrival time) under adiabatic conditionadDetermining a temperature correction coefficient gamma according to the corresponding temperature1The value range of (1) is that the heat generated by the decomposition reaction of the materials under the adiabatic condition is not dissipated, and the condition is the most severe, so that the temperature obtained under the adiabatic condition is the most conservative temperature, the result of the determined temperature correction coefficient is conservative, and the safety use reliability is high.
Drawings
FIG. 1 shows the temperature correction coefficient γ of the present invention1The value scheme process of (2);
fig. 2 is a DSC curve with a temperature rise rate β of 2.5, 5, 10, 20 ℃/min for an example of the present invention;
FIG. 3 kinetic parameters of the thermal decomposition reaction of the BunENNA feedstock of the example of the invention.
Detailed Description
The invention relates to a method for determining a temperature correction coefficient of a chemical combination process of an energetic material, which takes a material BuNENA of a butyl nitrooxyethyl nitramine chemical combination process as a research object, adopts a differential scanning calorimetry to test the thermal decomposition reaction of the material to obtain thermal decomposition characteristic quantities at different heating rates, and obtains the apparent activation energy E of a kinetic parameter of the thermal decomposition reaction of the material according to Kissinger equation regressionaAnd a pre-exponential factor A, thereby obtaining TMR between the arrival of maximum reaction rates of adiabatic conditions of different danger levelsadRequired process temperature TpAnd process temperature T in different hazard classespAnd the specific operation time coefficient is used as a scale for comparing and measuring the actual process temperature T of the materials in the processpsAnd actual process temperature time tsQuantitatively obtaining the temperature correction coefficient gamma1. The process is shown in figure 1 and specifically comprises the following steps:
the method comprises the following steps: differential scanning calorimetry test for obtaining apparent activation energy E of material thermal decomposition reaction kinetic parameters in Bunena chemical combination processaAnd pre-finger factor a.
Performing Differential Scanning Calorimetry (DSC) tests by using a Differential Scanning Calorimeter (DSC), and performing temperature programmed test on the material thermal decomposition reaction at a temperature rise rate β of 2.5, 5, 10, 20 ℃/min, respectively, to obtain a thermal decomposition curve as shown in fig. 2, wherein the material has two exothermic peaks in the thermal decomposition process, the first exothermic heat of decomposition of the material causes a subsequent material to rapidly decompose and cause a subsequent rapid exothermic reaction, so that kinetic parameters are calculated by using the peak temperature of the first exothermic peak, the thermal decomposition curve is as shown in fig. 2, and the obtained thermal decomposition peak temperatures are as shown in table 4:
TABLE 4 thermal decomposition parameters of Bunena materials at different heating rates
Figure BDA0002636774750000071
Obtaining the kinetic parameter apparent activation energy E of the material thermal decomposition reaction by adopting Kissinger equation formula (I) regressionaAnd pre-finger factor a:
Figure BDA0002636774750000081
in the formula, betaiThe temperature rise rate is DEG C/s; t is the decomposition peak-to-peak temperature (absolute temperature), K; a is a pre-exponential factor, s-1；EaIs apparent activation energy, J/mol; r is a gas constant, J/(mol. K). By
Figure BDA0002636774750000082
To pair
Figure BDA0002636774750000083
Plotting, a straight line can be obtained, and as shown in FIG. 3, the apparent activation energy E is obtained from the slope of the straight linea118.60kJ/mol, the pre-factor A is determined from the intercept24.50 s-1。
Step two: TMR for obtaining maximum reaction rate reaching time under adiabatic condition of different danger gradesadRequired process temperature Tp。
TMRadThe formula (II) is shown below:
Figure BDA0002636774750000084
in the formula, EaJ/mol as activation energy; a is a pre-exponential factor, s-1(ii) a T is the process temperature, K; r is a gas constant, J/(mol. K); c. CpThe specific heat capacity of the material is J/(g.K); qrThe exothermic amount of the thermal decomposition of the material is J/g.
Will have apparent activation energy Ea118.60kJ/mol, meaning that the factor A is e24.50 s-1Gas constant R is 8.314J/mol g, specific heat capacity of material cp2J/(g.K), exothermic heat Q of material thermal decompositionrThe formula (II) is substituted by 2182J/g, and TMR under 6 possible grades is obtained by calculationadThe process temperature T is required to be 100h, 50h, 24h, 12h, 8h and 1h respectivelyD100、TD50、TD24、TD12、TD8、TD1。
TMRadProcess temperature T required for 100hpIs TD100＝64.6℃，
TMRadThe required process temperature T is 50hpIs TD50＝70.6℃，
TMRadProcess temperature T required at 24hpIs TD24＝77.1℃，
TMRadProcess temperature T required for 12hpIs TD12＝83.5℃，
TMRadProcess temperature T required for 8hpIs TD8＝87.3℃，
TMRadProcess temperature T required for 1hpIs TD1＝108.4℃。
Step three: maximum reaction rate arrival time TMR at 6 adiabatic conditions of different hazard classesadTemperature T of the dividing processpAt a process temperature TpFor a scale, the actual process temperature T of the materials is compared and measuredpsQuantitatively obtaining the temperature correction coefficient gamma1。
Due to TMRadProcess temperature T required for 100hpIs TD100Thus, for TMRadNot less than 100h, when the actual process temperature T of the material isps≤TD100When is, i.e. Tps<Temperature correction coefficient gamma of 64.6 DEG C1＝1；
Due to TMR ad50h required process temperature TpIs TD50Thus, for a TMR of 50h ≦adLess than 100h, when the actual process temperature T of the materialD100＜Tps≤TD50When it is 64.6 ≦ TpsTemperature correction coefficient gamma of < 70.61＝2；
Due to TMR ad24h required process temperature TpIs TD24Thus, for a TMR of 24h ≦adLess than 50h, when the actual process temperature T of the materialD50≤Tps＜TD24When it is 64.6 ≦ TpsTemperature correction coefficient gamma of < 77.11＝3；
Due to TMRadProcess temperature T required for 12hpIs TD12Thus, for a TMR of 12h ≦adLess than 24h, when the actual process temperature T of the material isD24＜Tps≤TD12When it is 77.1 ≦ TpsTemperature correction coefficient gamma of < 83.51＝4；
Due to TMRadProcess temperature T required for 8hpIs TD8Thus, for 8h ≦ TMRadLess than 12h, when the actual process temperature T of the material isD12＜Tps≤TD8When T is 83.5 ≦ TpsTemperature correction coefficient gamma of < 87.31＝5；
Due to TMR ad1h required process temperature TpIs TD1Thus:
TMR for 1hadLess than 8h, when the actual process temperature T of the material isD8＜Tps≤TD1When it is 87.3 ≦ TpsTemperature correction coefficient gamma of < 108.41＝6～9；
For TMRadLess than or equal to 1h, when the actual process temperature T of the material isps≥TD1I.e. Tps>Temperature correction coefficient gamma of 108.4 DEG C1＝10；
Step four: temperature correction coefficient gamma1When the value range is 6-9, namely the actual process temperature TD8≤Tps＜TD1Then, the specific operation time coefficients of different danger levels are taken as a scale, and the temperature correction coefficient gamma is further quantitatively obtained1。
γ1When the value is 6-9, TMR is utilizedadAnd the operation time t of the actual material in the process is further determined, and the calculation formula of the specific operation time coefficient is shown in the formula (III):
Figure BDA0002636774750000101
in the formula: TMRadThe time for the maximum reaction rate to reach under adiabatic conditions at the process temperature, h; t is the operation time of the material in the process, h. The physical meaning of the specific operating time coefficient is the safety margin of the process operating time when gamma is1In the case of 6 to 9, the temperature correction coefficient γ is divided based on 1.5 times, 2 times, and 4 times of safety margin, i.e., the specific operating time coefficient is 1, 2.5, and 41The determination at 6-9 is shown in Table 3:
TABLE 3 temperature correction coefficient γ1Values in the range of 6 to 9
Figure BDA0002636774750000102
Based on the above results, the temperature correction coefficient of Bunena material in the procedure of butylnitrooxyethylnitramine combination is determined as shown in Table 4
TABLE 4 temperature correction coefficient values for the Bunena synthesis process
Figure BDA0002636774750000103
Figure BDA0002636774750000111
The actual process temperature of the material in the process of synthesizing the butyl nitrooxyethyl nitramine is 35 ℃, so the temperature correction coefficient gamma of the material under the process condition is determined according to the division of temperature hazard grades1＝1。

Claims (1)

1. A method for determining a temperature correction coefficient of an energetic material compounding process is characterized by comprising the following steps:
the method comprises the following steps: differential scanning calorimetry test for obtaining apparent activation energy E of material thermal decomposition reaction kinetic parameters in energetic material combination processaAnd pre-finger factor a:
utilizing a differential scanning calorimetry test to obtain the thermal decomposition reaction of materials in the chemical combination process of energetic materials, obtaining the peak temperature of the exothermic reaction causing the rapid decomposition of the materials, and utilizing Kissinger equation formula (I) regression to obtain the apparent activation energy E of kinetic parameters of the thermal decomposition reaction of the materialsaAnd pre-finger factor a:
Figure FDA0002636774740000011
in the formula, betaiThe heating rate is K/s; t is decomposition peak temperature, K; a is a pre-exponential factor, s-1；EaIs apparent activation energy, J/mol; r is a gas constant, J/(mol. K);
step two: TMR for obtaining maximum reaction rate reaching time under adiabatic condition of different danger gradesadRequired process temperature Tp；
TMRadThe formula (II) is shown below:
Figure FDA0002636774740000012
in the formula, EaIs apparent activation energy, J/mol; a is a pre-exponential factor, s-1(ii) a T is the process temperatureDegree, K; r is a gas constant, J/(mol. K); c. CpThe specific heat capacity of the material is J/(g.K); qrThe heat release of the material thermal decomposition is J/g;
TMR using the 6-class criterion of the Zurich Risk analysis method as a determination conditionadRespectively 100h, 50h, 24h, 12h, 8h and 1h, and respectively obtaining up to 6 TMRs by formula (II)adRequired process temperature TpSequentially comprises the following steps: t isD100、TD50、TD24、TD12、TD8、TD1；
Step three: maximum reaction rate arrival time TMR for 6 different hazard class adiabatic conditions at different hazard classesadTemperature T of the dividing processpAt a process temperature TpFor a scale, the actual process temperature T of the materials is compared and measuredpsQuantitatively obtaining the temperature correction coefficient gamma1Comprises the following steps:
Figure FDA0002636774740000021
step four: temperature correction coefficient gamma1The value range is 6-9, namely the actual process temperature TD8≤Tps＜TD1Then, the specific operation time coefficients of different danger levels are taken as a scale, and the temperature correction coefficient gamma is further quantitatively obtained1：
γ1When the value range is 6-9, TMR is utilizedadAnd further determining the operation time in the actual process of the material, and calculating a specific operation time coefficient by a formula (III):
Figure FDA0002636774740000022
in the formula: TMRadThe maximum reaction rate reaching time under the adiabatic condition, and t is the operation time of the material in the process; the specific operation time coefficient is the safety margin of the process operation time when the gamma is1When the content is 6-9 times, the safety is 1.5 times, 2 times and 4 timesDetermining a temperature correction coefficient of the material, namely gamma, according to the margin, namely a specific operation time coefficient of 1, 2.5 and 41Temperature correction coefficient gamma at 6-91Comprises the following steps:
Figure FDA0002636774740000023
。
CN202010827564.8A 2020-08-17 2020-08-17 Method for determining temperature correction coefficient of energetic material combination process Active CN112102892B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010827564.8A CN112102892B (en) 2020-08-17 2020-08-17 Method for determining temperature correction coefficient of energetic material combination process

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010827564.8A CN112102892B (en) 2020-08-17 2020-08-17 Method for determining temperature correction coefficient of energetic material combination process

Publications (2)

Publication Number Publication Date
CN112102892A true CN112102892A (en) 2020-12-18
CN112102892B CN112102892B (en) 2023-05-09

Family

ID=73753020

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010827564.8A Active CN112102892B (en) 2020-08-17 2020-08-17 Method for determining temperature correction coefficient of energetic material combination process

Country Status (1)

Country Link
CN (1) CN112102892B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103776776A (en) * 2013-08-12 2014-05-07 河南教育学院 Method and device for researching quick reaction mechanism of energetic material
CN107391870A (en) * 2017-08-02 2017-11-24 中国科学技术大学 The computational methods of thermal runaway critical condition
CN109949874A (en) * 2019-04-16 2019-06-28 南京工业大学 Risk grading method for safety assessment in fine chemical production process
CN110010205A (en) * 2019-03-21 2019-07-12 南京理工大学 A method of assessment self-catalysis substance thermal stability
CN110414868A (en) * 2019-08-13 2019-11-05 南京工业大学 Method for evaluating thermal runaway risk degree in chemical process
CN110967365A (en) * 2019-12-17 2020-04-07 河北科技大学 Method for analyzing pyrolysis kinetics of metal organic framework material based on combination of thermal analysis and mass spectrometry

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103776776A (en) * 2013-08-12 2014-05-07 河南教育学院 Method and device for researching quick reaction mechanism of energetic material
CN107391870A (en) * 2017-08-02 2017-11-24 中国科学技术大学 The computational methods of thermal runaway critical condition
CN110010205A (en) * 2019-03-21 2019-07-12 南京理工大学 A method of assessment self-catalysis substance thermal stability
CN109949874A (en) * 2019-04-16 2019-06-28 南京工业大学 Risk grading method for safety assessment in fine chemical production process
CN110414868A (en) * 2019-08-13 2019-11-05 南京工业大学 Method for evaluating thermal runaway risk degree in chemical process
CN110967365A (en) * 2019-12-17 2020-04-07 河北科技大学 Method for analyzing pyrolysis kinetics of metal organic framework material based on combination of thermal analysis and mass spectrometry

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YI-MING LU 等: "Evaluation of thermal hazards based on thermokinetic parameters of 2-(1-cyano-1-methylethyl)azocarboxamide by ARC and DSC", 《JOURNAL OF THERMAL ANALYSIS AND CALORIMETRY》 *
孙笑 等: "HNF的热分解动力学和热安全性", 《含能材料》 *
张引莉 等: "含能银配合物的爆炸性能和热动力学参数及催化燃烧性能研究", 《化学研究与应用》 *

Also Published As

Publication number Publication date
CN112102892B (en) 2023-05-09

Similar Documents

Publication Publication Date Title
Zhou et al. A novel method for kinetics analysis of pyrolysis of hemicellulose, cellulose, and lignin in TGA and macro-TGA
CN1849501B (en) Low thermal inerta scanning adiabatic calorimeter
Razus et al. The rate of pressure rise of gaseous propylene–air explosions in spherical and cylindrical enclosures
CN102608284B (en) Method for determining explosion limit of multi-component mixed gas
WO2021026810A1 (en) Method for assessing risk of thermal runaway in chemical process
CN107478674B (en) Microcalorimetry for evaluation of stability of single-base propellant during 10-year storage
CN108872814B (en) Method for evaluating insulation life of oil paper in high-voltage current transformer
CN109856039A (en) Inner screw channel type ethane cracking furnace pipe residue lifetime estimation method based on L-M parametric method
Fitzgerald The kinetics of coal carbonization in the plastic state
CN108548748A (en) A kind of gravitational thermal analysis method and device
CN112033998B (en) Thermal insulation acceleration calorimetry-based explosive material thermal stability grading method
CN110411864A (en) High-temperature creep life prediction analysis calculation method based on creep activation energy
CN112102892A (en) Method for determining temperature correction coefficient of energetic material chemical combination process
CN117233210B (en) Polyethylene dust explosion hazard prediction method
CN112033997B (en) Explosive thermal stability grading method based on differential scanning calorimetry
Liu et al. New findings on thermal degradation properties of fluoropolymers: Combined kinetic analysis for a completely overlapped reaction
CN112067657B (en) Method for determining temperature correction coefficient of energetic material chemical combination process
CN208043587U (en) A kind of Thermgravimetric Analysis Apparatus
CN115223666A (en) Comprehensive evaluation method for fire hazard risk of material
CN108037263A (en) A kind of thermal power plant material ageing detection method
Smith et al. THE DENSITIES AND DEGREES OF DISSOCIATION OF THE SATURATED VAPORS OF THE AMMONIUM HALIDES, AND THE RELATED THERMAL DATA.
Liu et al. Thermal Kinetic Parameters of Lead Azide and Lead Styphnate with Antistatic Additives
CN111028898A (en) Method for evaluating damage failure life of aluminum electrolysis cathode material
JP6657867B2 (en) Estimation method of coke shrinkage
Domalski et al. The heat of formation of boron trifluoride by direct combination of the elements

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant