Ｖｏ１．３４ 高等学校化学学报 ＣＨＥＭＩＣＡＬ ＪＯＵＲＮＡＬ ＯＦ ＣＨＩＮＥＳＥ ＵＮＩＶＥＲＳＩＴＩＥＳ Ｎｏ．１１ ２５５８～２５６５ ２０１３年１１月 』ｄｏｉ：１０．７５０３／ｃｊｃｕ２０１３０５６８ 不同温度下ＨＭＸ和ＲＤＸ晶体的 感度判别和力学性能预估 ——分子动力学比较研究 刘冬梅 ，赵 丽 ，肖继军 ，陈 军 ，姬广富 ， 朱伟。，赵锋 ，吴 强 ，肖鹤呜 （１．南京理工大学化工学院，分子与材料计算研究所，南京２１００９４； ２．中国工程物理研究院冲击波物理与爆轰物理国防技术重点实验室，绵阳６２１９００； ３．嘉兴学院生物与化学工程学院，嘉兴３１４００１） 摘要为探讨和比较炸药黑索金（ＲＤＸ）和奥克托今（ＨＭＸ）晶体的结构、能量和力学性质随温度的递变规 律，在ＣＯＭＰＡＳＳ力场和ＮＰＴ系综下，对其合适的等原子数超晶胞模型分别进行５个温度（１９５，２４５，２９５， ３４５和３９５ Ｋ）下的分子动力学（ＭＤ）周期性模拟研究．结果表明，随着温度的升高，ＲＤＸ和ＨＭＸ晶体引发 键Ｎ—ＮＯ 的最大键长（Ｌ…）的逐渐增大以及引发键连双原子作用能（Ｅ ）和内聚能密度（ＣＥＤ）的逐渐减 小均与感度随温度升高而增大的实验结果一致；且在各温度下，ＲＤＸ晶体的，Ｊ…均大于ＨＭＸ晶体的Ｌ…， 与ＨＭＸ相比，ＲＤＸ的Ｅ 一 和ＣＥＤ均较小，上述结果与ＲＤＸ比ＨＭＸ感度大的实验结果相符．由此表明，在 一定条件下，Ｌ…，Ｅ 和ＣＥＤ可用于高能物质的热和撞击感度的相对大小的判别．基于ＭＤ模拟原子运动 轨迹，用静态法求得２种晶体的弹性力学性能，发现拉伸、体积和剪切模量均随温度的升高而递减，与实验 结果一致． 关键词环三亚甲基三硝胺；环四亚甲基四硝胺；分子动力学模拟；感度；力学性能 ０６４１；ＴＪ５５ 文献标志码Ａ 中图分类号ＲＤＸ（环三亚甲基三硝胺，又名黑索金）和ＨＭＸ（环四亚甲基四硝胺，又名奥克托今）均为常见的 高能单体猛炸药，可作为混合炸药、传爆药、固体推进剂和发射药配方中的高能组分或氧化剂使用， 而应用在航天、国防和国民经济等诸多领域．因为二者均为环杂硝胺类化合物，分子结构相似（图１） 并具有类似的热解和起爆机理 ｌ卫Ｊ，放在一起作比较研究，可抵消非本质因素造成的系统误差，同时可 利于突出共性和差异，揭示本质和规律性．虽然关于ＲＤＸ和ＨＭＸ的分子动力学（ＭＤ）理论研究已较 多 “ ，但对二者的结构、能量、性质及其温度效应作比较的ＭＤ模拟至今仍未见报道． 感度是炸药在外界刺激下发生爆炸的难易程度，其一直是炸药化学、工艺学及爆炸理论关注的重 点．各类型感度主要由实验测定¨ ，但人们希望可以通过理论方法对其加以判别．但迄今为止，在微 观层次上，利用量子化学等理论方法判别单体炸药感度的报道较多＿】 系的报道＿ｌ ；然而， ，近期才有借助ＭＤ模拟结 果，以研究引发键最大键长（Ｌ ）、引发键连双原子作用能（Ｅ 一。）和内聚能密度（ＣＥＤ）与感度之间关 ，Ｅ 一 和ＣＥＤ等理论参量是否适用于不同炸药感度相对大小的判别和比较 等问题有待解决． 收稿日期：２０１３－０６－１８． 基金项目：国家自然科学基金委员会与中国工程物理研究院联合基金（批准号：Ｕ１２３０１２０）和国防基础科研计划（批准号 Ｂ１５２０１ １０００２）资助． 联系人简介：肖继军，男，博士，教授，博士生导师，主要从事理论和计算化学在高分子材料中的应用研究． Ｅ—ｍａｉｌ：ｘｉａｏ ｊｉｊｕｎ＠ｍａｉｌ．ｎｊｕｓｔ．ｅｄｕ．ｃｎ 肖鹤鸣，男，教授，博士生导师，主要从事量子化学、分子动力学及其对含能材料结构、性能和设计的应用基础研究 Ｅ—ｍａｉｌ：ｘｉａｏ＠ｍａｉｌ．ｎｊｕｓｔ．ｅｄｕ．ｃｎ Ｎ０．１１ 刘冬梅等：不同温度下ＨＭＸ和ＲＤＸ晶体的感度判别和……预估——分子动力学比较研究 （Ａ） ｏ ｏ （Ｂ） 。’ ２５５９ ｌ ｌ１ Ｈ２Ｃ／ ＼ｃＨ＾ ／』＼ｃ／ ＼Ｎ 。 Ｈ＼Ｃ＼Ｎ／ Ｃ ’ ０ ０ Ｆｉｇ．１ Ｍｏｌｅｃｕｌａｒ ｓｔｒｕｃｔｕｒｅｓ ｏｆ ＲＤＸ（Ａ）ａｎｄ ＨＭＸ（Ｂ） 力学性能是高能材料最重要的实用性能之一，关系到其制备、加工、生产和使用．力学性能和感 度一样，主要靠实验测定．但实验只能提供常温下的测定结果，且至今缺乏统一规范的比较标准．提 供不同温度下的感度判别方法和力学性能预估结果，不仅具有重要的理论意义，而且具有重大的实用 价值．提供物质结构及性能随温度而变化的统计规律，正是ＭＤ方法的优势和特长，将其用于有安全 隐患的高能物质则更为合适 引． 因此，本文选用Ａｃｃｅｌｒｙｓ公司的Ｍｓ程序和ＣＯＭＰＡＳＳ力场 ，构建了适合用于相互比较的 ＨＭＸ和ＲＤＸ的等原子超晶胞模型，在５个不同温度（１９５，２４５，２９５，３４５和３９５ Ｋ）和ＮＰＴ系综下，完 成ＭＤ研究，重点探讨了结构、能量和力学性能及其与感度的关系，并将结果用于不同体系（ＨＭＸ和 ＲＤＸ）感度判别的比较研究． １ ＭＤ模拟方法和细节 １．１力场选择和模型构建 选择ＣＯＭＰＡＳＳ力场 能的研究 Ｈ～１７，２８－３０］． 进行ＭＤ模拟的原因在于其特别适用于凝聚相和硝胺类炸药的结构和性 ＲＤＸ和ＨＭＸ晶体的单胞结构取自文献［３１，３２］中的中子衍射结果．每个单胞均含８个ＲＤＸ分子 和２个ＨＭＸ分子．分别建立（２ｘ２ｘ３）ＲＤＸ超晶胞和（３ｘ３ｘ４）ＨＭＸ超晶胞，使分别包含９６个ＲＤＸ分 子和７２个ＨＭＸ分子的２个模拟周期箱中均包含２０１６个原子，即模拟体系大小和质量相等，有利于模 拟结果的合理比较．上述２种晶体的单胞结构和超晶胞初始模型如图２所示． （Ａ） （Ｂ） 一 Ｆｉｇ．２ Ｐｒｉｍｉｔｉｖｅ ｃｅｌｌｓ（Ａ，Ｃ）ａｎｄ ｉｎｉｔｉａｌ ｓｕｐｅｒｃｅｌｌ ｍｏｄｅｌｓ（Ｂ，Ｄ）ｏｆ ＲＤＸ（Ａ，Ｂ）ａｎｄ ＨＭＸ（Ｃ，Ｄ） 高等学校化学学报 １．２模拟过程 将（２ｘ２ｘ３）ＲＤＸ和（３×３×４）ＨＭＸ的超晶胞分别在ＣＯＭＰＡＳＳ力场及ＮＰＴ系综下进行ＭＤ模拟． ／ａ胛Ｊ ｄ８ ，ａ 景．Ｉ 蠹 模拟过程中温度和压力的控制分别采用Ａｎｄｅｒｓｏｎ＿３。 和Ｐａｒｒｉｎｅｌｌｏ 方法，范德华（ｖｄＷ）和静电作用 ５ Ｏ ５ Ｏ ５ ０ ５ ０ ５ 如 （Ｃｏｕｌｏｍｂ）的计算则分别采用Ａｔｏｍ—ｂａｓｅｄ和Ｅｗａｌｄ＿３ 加和方法，截断半径９．５ｘ１０圳ｍ，并进行截断尾 部校正． 依次选取温度为１９５，２４５，２９５，３４５和３９５ Ｋ；初始分子运动速度按Ｍａｘｗｅｌ１．Ｂｏｌｔｚｍａｎｎ分布确定； 基于周期性边界条件和时间平均等效于系综平均等假设，求解牛顿运动方程；积分计算采用Ｖｅｒｌｅｔ方 法．时间步长取１ ｆｓ．总模拟步长取１．５ ｎｓ，其中前０．５ ｎｓ用于平衡，后１ ＢＳ用于统计分析． １．３平衡判别和平衡结构 体系的平衡可由温度和能量同时达到平衡予以确定．通常当温度和能量在５％一１０％范围内波动 即可认为体系已达平衡． 本文以ＨＭＸ和ＲＤＸ晶体在２９５ Ｋ下的ＭＤ模拟为例，从其温度一时间和能量．时间曲线（图３）可 见，二者均已达热力学平衡． ３ｌ０ ３０５ ３００ ２９５ ２９０ 要 击 ２８５ ２８０ ２７５ ０ ５００００ １０００００ １５００００ ２０００００ ｔ／ｆｓ ，＿、 要 击 ０ ５００００ １０００００ １ ５００００ ２０００００ ｔ／ｆｓ Ｆｉｇ．３ Ｔｉｍｅ－ｅｖｏｌｕｔｉｏｎｓ ｏｆ ｔｈｅ ｔｅｍｐｅｒａｔｕｒｅｓ（Ａ，Ｃ）ａｎｄ ｅｎｅｒｇｉｅｓ（Ｂ，Ｄ）ｏｆ ＲＤＸ（Ａ，Ｂ）ａｎｄ ＨＭＸ（Ｃ，Ｄ）ｃｒｙｓｔａｌｓ ａｔ ２９５ Ｋ （Ｂ）ａｎｄ（Ｄ）口．Ｔｏｔａｌ ｎｏｎｂｏｎｄ ｅｎｅｒｙ；ｂ．ｐｏｔｇｅｎｔｉｌａ ｅｎｅｒｇｙ． 经ＭＤ模拟共获得ＲＤＸ和ＨＭＸ晶体５个温度下的平衡结构．本文仅以２９５ Ｋ的结果为例，结果 示于图４中． （Ａ） Ｆｉｇ．４ Ｅｑｕｉｈ＇ｂｒｉｕｍ ｓｔｒｕｃｔｕｒｅｓ ｏｆ ＲＤＸ（Ａ）ａｎｄ ＨＭＸ（Ｂ）ａｔ ２９５ Ｋ Ｎｏ．１１ 刘冬梅等：不同温度下ＨＭＸ和ＲＤＸ晶体的感度判别和……预估——分子动力学比较研究 ２５６１ ２结果与讨论 ２．１ 引发键最大键长与感度的关系 高能化合物的引发键是其分子中所有化学键中最弱的，在外界作用下该键将优先断裂而引发分解 ４ ２ Ｏ ８ ６ ４ ２ ０ ２ 和起爆，Ｎ—Ｎ０ 键就是硝胺类高能物质的热解或起爆的引发键 ．通常实验和理论计算只给出 分子中的平均键长，而ＭＤ模拟可提供Ｎ—Ｎ０ 键长的统计分布，对于探讨高能体系的起爆、感度和 安全性具有重要意义．表１列出了基于ＭＤ模拟轨迹所得的不同温度下ＲＤＸ和ＨＭＸ晶体中Ｎ—ＮＯ： 键的最可几键长（Ｌ ）、平均键长（Ｌ 。）和最大键长（，Ｊ ）值．以２９５ Ｋ的结果为例，图５给出了二者 的Ｎ—Ｎ０ 键的键长分布情况．图６给出了ＲＤＸ及ＨＭＸ平衡体系中Ｎ—Ｎ０ 键的最大键长与温度的 关系． Ｔａｂｌｅ １ Ｌｐ ｂ，Ｌ ａｎｄ Ｌ ｆｏｒ Ｎ—ＮＯ２ ｔｒｉｇｇｅｒ ｂｏｎｄ ｉｎ ＲＤＸ ａｎｄ ＨＭＸ ｃｒｙｓｔａｌｓ ａｔ ｄｉｆｅｒｅｎｔ ｔｅｍｐｅｒａｔｕｒｅｓ 一 一一Ｎ．． ７） Ｈ ８ ６ ４ ２ ０ ０．０８ ０．１０ ０ １２ ０ １４ ０ ｌ６ Ｏ．１８ ０．２Ｏ ０．０８ Ｏ．１０ ０．１２ ０．１４ ０．１６ ０ ｌ８ Ｏ ２Ｏ 三Ｎ—Ｎ／ｎｍ 三Ｎ—Ｎ／ｎｍ Ｆｉｇ．５ Ｎ—Ｎｏ２ ｔｒｉｇｇｅｒ ｂｏｎｄ ｌｅｎｇｔｈ ｄｉｓｔｒｉｂｕｔｉｏｎｓ ｉｎ ＲＤＸ（Ａ）ａｎｄ ＨＭＸ（Ｂ）ｃｒｙｓｔａｌｓ 近．两者平均键长和最可几键长随温度升高的变化 Ｆｉｇ．６ Ｍａｘｉｍｕｍ ｂｏｎｄ ｌｅｎｇｔｈ ｏｆ Ｎ—Ｎｏ ｔｒｉｇｇｅｒ 很小，而最大键长（Ｌ ）随温度升高显著增大，与 感度随温度升高而增大的实验结果一致．由表１和ｂｏｎｄ ｅｒｓｕｓ ｔｅｍｐｅｒａｔｕｒｅ ｉｎ ＲＤＸ（ａ）ａｎｄ ＨＭＸ（西）ｃｒｙｓｔａｌｓ 图６还可见，在每个温度下，ＲＤＸ的 …均大于ＨＭＸ，这与ＲＤＸ比ＨＭＸ感度大的实验结果相吻合． 由此可见，可利用ＲＤＸ和ＨＭＸ晶体中引发键（Ｎ—Ｎ０：）的Ｌｍ＝判别这些化合物的热和撞击感度的相 对大小，该结论与之前在混合和缺陷高能体系中已发现的结果一致 ０，Ｈ， ＇ ．虽然 布中只占１０一，但仍值得特别关注；具有 生化学反应的碰撞理论可知，具有 ２．２ 引发键连双原子作用能与感度的关系 键长在统计分 的分子能量很高，非常活泼，根据只有活化的分子才能发 的活泼分子中的引发键最容易断裂，进而引发分解和起爆． 将ＲＤＸ和ＨＭＸ晶体中引发键连双原子相互作用能（Ｅ 一 ）定义为 ＥＮ—Ｎ＝（ＥＴ—ＥＴ ）／凡 （１） 式中， 和 分别为体系在所用ＣＯＭＰＡＳＳ力场框架下的总能量和固定所有Ｎ原子而求得的总能量， ２５６２ 高等学校化学学报 凡为体系中Ｎ—ＮＯ 的数量．表２和表３分别给出了ＭＤ模拟所得到的不同温度下ＲＤＸ和ＨＭＸ晶体 －中的Ｅ 一 及其相关分量． Ｔａｂｌｅ ２ Ｉｎｔｅｒａｃｔｉｏｎ ｅｎｅｒｇｉｅｓ Ｎ＿Ｎ ａｎｄ ｔｈｅｉｒ ｃｏｍｐｏｎｅｎｔｓ ａｔ ｄｉｆｆｅｒｅｎｔ ｔｅｍｐｅｒａｔｕｒｅｓ ｆｏｒ ＲＤＸ ｃｒｙｓｔａｌｓ Ｔａｂｌｅ ３ Ｉｎｔｅｒａｃｔｉｏｎ ｅｎｅｒｇｉｅｓ Ｎ Ｎ ａｎｄ ｔｈｅｉｒ ｃｏｍｐｏｎｅｎｔｓ ａｔ ｄｉｆｅｒｅｎｔ ｔｅｍｐｅｒａｔｕｒｅｓ ｆｏｒ ＨＭＸ ｃｒｙｓｔａｌｓ 由表２、表３和图７可见，随着温度的升高， ＲＤＸ和ＨＭＸ晶体中的Ｅ 一 单调递减，与感度随温 度升高而递增的实验结果相符．从表２还可见， Ｅ 主要由式（１）中 和 中的静电力之差值所 决定，其中ｖｄＷ力和内势能的贡献量很小．由图７ 还可见，在不同温度下，ＨＭＸ晶体的该双原子作用 能Ｅ 一 均比ＲＤＸ晶体的大，与ＨＭＸ感度小于 ＲＤＸ的实验结果一致．由此可见， ～ 与 的相对大小． ２．３ 内聚能密度与感度的关系 ｉ 均与 。 Ｉｎｔｅｒａｃｔｉｏｎ ｅｎｅｒｇｙ（ＥＮ—Ｎ）ＶＳ．ｔｅｍｐｅｒａｔｕｒｅ ｏｒ ｆＲＤＸ（ａ）ａｎｄ ＨＭＸ（ｂ）ｃｒｙｓｔａｌｓ 键强度相关，均可用于判别高能物质热和撞击感度 …内聚能密度（ＣＥＤ）是单位体积内１ ｔｏｏｌ凝聚体克服分子间作用力变为气态时所需的能量，可由实 验测得．在ＭＤ模拟中，ＣＥＤ是范德华（ｖｄＷ）力与静电力之和，即分子间的非键力． 表４给出了基于ＭＤ模拟轨迹求得的不同温度下ＲＤＸ和ＨＭＸ晶体的内聚能密度及其分量，图８ 为其图示． Ｔａｂｌｅ ４ Ｃｏｈｅｓｉｖｅ ｅｎｅｒｇｙ ｄｅｎｓｉｔｉｅｓ ＣＥＤ ａｎｄ ｔｈｅｉｒ ｃｏｍｐｏｎｅｎｔｓ ｏｆ ＲＤＸ ａｎｄ ＨＭＸ ｃｒｙｓｔａｌｓ ａｔ ｄｉｆｅｒｅｎｔ ｔｅｍｐｅｒａｔｕｒｅｓ ｝ＣＥＤ＝Ｅ ｄｗ＋ＥＥｌｅｃｔ １ａｌｉｃ． 由表４和图８可见，随着温度的升高，ＲＤＸ和ＨＭＸ晶体的ＣＥＤ，ｖｄＷ力和静电力均单调递减，表 明它们由晶态变为气态时所需能量变小，这与温度升高感度增大的实验结果一致；在同一温度下， ＲＤＸ晶体的ＣＥＤ均小于ＨＭＸ晶体的相应值，与ＲＤＸ感度大于ＨＭＸ的实验结果相符． 由此表明，ＲＤＸ和ＨＭＸ晶体中的ＣＥＤ在一定条件下也可用于热感度相对大小的理论判别． Ｎｏ．１１ 刘冬梅等：不同温度下ＨＭＸ和ＲＤＸ晶体的感度判别和……预估——分子动力学比较研究 ２５６３ 目 ９ Ｆｉｇ．８ ＣＥＤ（ａ），ｖｄＷ（ｂ），ａｎｄ ｅｌｅｃｔｒｏｓｔａｔｉｃ（Ｃ）ｅｎｅｒｇｉｅｓ ． ｔｅｍｐｅｒａｔｕｒｅ ｏｆ ＲＤＸ（Ａ）ａｎｄ ＨＭＸ（Ｂ）ｃｒｙｓｔａｌｓ ２．４弹性力学性能 弹性力学性能的基本参量主要有弹性系数、工程模量和泊松比．模量是评价材料刚性的指标，是 材料抵抗弹性形变能力的度量．材料硬度与其拉伸模量（Ｅ）和剪切模量（Ｇ）密切相关，其值越大，材料 的硬度越高．材料断裂强度的大小由体模量（ ）度量，Ｋ值越大，表明断裂所需能量越大，即材料断裂 强度越大．表５和表６分别列出了基于ＭＤ模拟轨迹由静态法分析所得不同温度下ＲＤＸ和ＨＭＸ晶体 的力学性能． Ｔａｂｌｅ ５ Ｍｅｃｈａｎｉｃａｌ ｐｒｏｐｅｒｔｉｅｓ ａｔ ｄｉｆｆｅｒｅｎｔ ｔｅｍｐｅｒａｔｕｒｅｓ ｆｏｒ ＲＤＸ ｃｒｙｓｔａｌ 由表５和表６可见，随着温度升高，ＲＤＸ和ＨＭＸ晶体的拉伸模量、体模量和剪切模量均呈下降 趋势，表明其刚性、硬度和断裂强度均随温度升高而下降，与实验结果相符．但随着温度升高，体系模 量减小，按通常观点 ∞Ｊ，模量减小会使体系感度下降，这与感度随温度升高而递增的实验结果不 符．分析原因可能在于影响感度的因素很多，随着温度的升高，分子运动加快，体系的结构（如引发键 最大键长）和能量（如引发键连双原子作用能）改变较大，且结构（Ｌ ）和能量（Ｅ 一 ）对感度的影响成 为主导因素，超过了力学性能变化对感度的影响，即力学性能（模量）对感度的影响已成为次要因素． ２．５理论和实验验证 本文ＭＤ模拟研究结果，得到了之前相关理论与实验工作的充分支持．首先，随着温度升高，分子 运动加快，炸药必定更易于分解起爆，感度增大；随着温度的升高，物质的弹性力学性能必然由刚变 柔，由硬变软，均符合实验结果． 其次，经量子化学动态理论计算，在密度泛函理论（ＤＦＴ）一（Ｕ）Ｂ３ＬＹＰ／６－３１Ｇ一水平下求得的ＲＤＸ 和ＨＭＸ引发键（Ｎ—ＮＯ ）的均裂反应键离解能（ＢＤＥ），即该均裂反应的活化能，分别为１５３．８０ （１５２．０９ ）和１６６．４８ ｋＪ／ｍｏｌ ，后者大于前者，证明在外界刺激下，ＲＤＸ比ＨＭＸ更易分解起爆． ＲＤＸ和ＨＭＸ晶体前沿带隙（△Ｅ ＝ＥＬＵＭ０一Ｅ ｏＭｏ）的ＤＦＴ计算值分别为３．３２和３．６２ ｅＶ ，后者大于 前者，根据第一性原理带隙判据，说明ＨＭＸ晶体比ＲＤＸ晶体钝感． 在实验方面，真空安定性实验（１００℃，４８ ｈ放气量） １８］：ＨＭＸ为０．０６ ｍＬ／ｇ，ＲＤＸ为０．０７ ｍＬ／ｇ； 热失重实验（１００℃，４８ ｈ失重）¨ ：ＨＭＸ为０．０２５％，ＲＤＸ为０．０９％．比较发现，与ＲＤＸ相比，ＨＭＸ 的放气量小，失重也小，说明ＨＭＸ的热安定性好于ＲＤＸ．ＨＭＸ和ＲＤＸ的爆发点ｕ （５ Ｓ，延滞期）分别 ７ ９ ｍ ｎ ２５６４ 高等学校化学学报 Ｖｏ１．３４ 为３２７和２３０℃，说明ＨＭＸ的热感度小于ＲＤＸ． 此外，实验上测得的５０％爆炸特性落高 Ｈ５。 （美国１２型工具，２．５ ｋｇ落锤）：ＨＭＸ为３２ ｃｍ， ＲＤＸ为２６ ａｍ，说明ＨＭＸ的撞击感度比ＲＤＸ的小． ３ 结 论 本文在５个不同温度下，以等原子数超晶胞模型对ＨＭＸ和ＲＤＸ晶体进行了ＭＤ比较研究．二者 引发键Ｎ—ＮＯ 的平均键长（Ｌ ）均接近实验值；最大键长（Ｌ…）随温度升高而单调递增，各温度下 ＲＤＸ的 …值均大于ＨＭＸ的，与感度随温度升高而增大以及ＲＤＸ感度大于ＨＭＸ的实验结果一致， 表明引发键的 …可用于判别高能物质热和撞击感度的相对大小．模拟所得ＲＤＸ和ＨＭＸ晶体中引发 键连双原子作用能（Ｅ 一 ）和内聚能密度（ＣＥＤ）均随温度升高而单调递减，且在各温度下ＨＭＸ晶体的 Ｅ 和ＣＥＤ均大于ＲＤＸ的相应值，表明Ｅ 一 和ＣＥＤ与感度之间存在负相关关系，Ｅ 一 和ＣＥＤ亦可 用于热和撞击感度的相对大小的判别．基于ＭＤ模拟原子运动轨迹和静态力学分析，求得ＨＭＸ和 ＲＤＸ晶体的弹性力学性能，结果发现，随着温度升高，二者的模量减小，表明刚性降低，强度减小，与 实验结果一致．本文对不同温度下ＨＭＸ和ＲＤＸ晶体热和撞击感度相对大小的判别以及力学性能的预 估，为高能物质安全性研究和配方设计提供了有益的信息和参考． 感谢中国工程物理研究院化工材料研究所李金山研究员提供的帮助 参１ ２ ３ ４ 考文献 Ｂｒｉｌｌ Ｔ．Ｂ．，Ｂｒｕｓｈ Ｐ．Ｊ．，Ｇｒａｙ Ｐ．，Ｋｉｎｌｏｃｈ Ｓ．Ａ．，Ｐｈｉｌｏｓ．Ｔｒａｎｓ．Ｒ．Ｓｏｃ．Ｌｏｎｄｏｎ Ａ，１９９２，３９（１６５４），３７７＿３８５ Ｚｈｏｕ Ｔ．Ｔ．，Ｈｕａｎｇ Ｆ．Ｌ．， Ｐ ．Ｃｈｅｍ．Ｂ，２０１１，１１５，２７８—２８７ Ｂｅｄｒｏｖ Ｄ．，Ｓｉｍｔｈ Ｇ．Ｄ．，Ｓｅｗｅｌｌ Ｔ．Ｄ．， Ｃｏｍｐｕｔ．，ＡｉｄｅｄＭａｔｅｒ．Ｄｅｓ．，２００１，８，７７—８５ Ｃｈｅｍ．Ｐｈｙｓ．，２００３，１１９，７４１７－－７４２６ Ｓｅｗｅｌｌ Ｔ．Ｄ．，ＭｅｎｉｋｏｆｆＲ．，Ｂｅｄｒｏｖ Ｄ．，Ｓｍｉｔｈ Ｇ．Ｄ．， ５ ６ Ｍａｎａａ Ｍ．Ｒ．，Ｆｒｉｅｄ Ｌ．Ｅ．，Ｍｅｌｉｕｓ Ｃ．Ｆ．，Ｅｌｓｔｎｅｒ Ｍ．，Ｆｒａｕｅｎｈｅｉｍ Ｔ．， Ｃｈｅｍ．Ｐｈｙｓ．Ａ，２００２，１０６（３９），９０２４－－－９０２９ Ｓｔｒａｃｈａｎ Ａ．，ｖａｎ Ｄｕｉｎ Ａ．Ｃ．，Ｃｈａｋｒａｂｏｒｔｙ Ｄ．，Ｄａｓｇｕｐｔａ Ｓ．，Ｇｏｄｄａｒｄ ＩＩＩ Ｗ．Ａ．，Ｐｈｙｓ．Ｒｅｖ．Ｌｅｔｔ．，２００３，９１（９），０９８３０１－１—— ０９８３０１－４ Ｚｈｅｎｇ Ｌ．，Ｔｈｏｍｐｓｏｎ Ｄ．Ｌ．， Ｃｈｅｍ．Ｐｈｙｓ．，２００６，１２５（８），０８４５０５—１－－－０８４５０５－９ Ａｇｒａｗａｌ Ｐ．Ｍ．，Ｒｉｃｅ Ｂ．Ｍ．，Ｚｈｅｎｇ Ｌ．，Ｔｈｏｍｐｓｏｎ Ｄ．Ｌ．，Ｊ．Ｐｈｙｓ．Ｃｈｅｍ．Ｂ，２００６，１１０（５１），２６１８５－－２６１８８ Ｘｉａｏ Ｊ．Ｊ．，Ｚｈｕ Ｗ．，Ｍａ Ｘ．Ｆ．，Ｘｉａｏ Ｈ．Ｍ．，Ｈｕａｎｇ Ｈ．，ｌＪ Ｊ．Ｓ．，Ｍｏ１．Ｓｉｍｕ１．，２００８，３４（８），７７５一＿７７９ Ｚｈｕ Ｗ．，Ｘｉａｏ Ｊ．Ｊ．，Ｚｈｕ Ｗ．Ｈ．，Ｘｉａｏ Ｈ．Ｍ．， Ｈａｚａｒｄ．Ｍａｔｅｒ．，２００９，１６４，１Ｏ８２—１０８８ Ｃｕｉ Ｈ．Ｌ．，Ｊｉ Ｇ．Ｆ．，Ｃｈｅｎ Ｘ．Ｒ．，Ｚｈａｎｇ Ｑ．Ｍ．，Ｗｅｉ Ｄ．Ｑ．，Ｚｈａｏ Ｆ．， Ｃｈｅｍ．Ｅｎｇ．Ｄａｔａ，２０１０，５５，３１２ｌ＿３１２９ Ｇｅ Ｎ．Ｎ．，Ｗｅｉ Ｙ．Ｋ．，Ｊｉ Ｇ．Ｆ．，Ｃｈｅｎ Ｘ．Ｒ．，Ｚｈａｏ Ｆ．，Ｗｅｉ Ｄ．Ｑ．，＿，＿Ｐｈｙｓ．Ｃｈｅｍ．Ｂ，２０１２，１１６，１３６９６－－１３７０４ Ａｋｋｂａｒｚａｄｅ Ｈ．，Ｐａｒｓａｆａｒ Ｇ．Ａ．，Ｂａｙａｔ Ｙ．，Ａｐｐ１．Ｓｕ Ｓｃｉ．，２０１２，２５８，２２２６－－２２３０ Ｘｉａｏ Ｊ．Ｊ．，Ｚｈａｏ Ｌ．，Ｚｈｕ Ｗ．，Ｃｈｅｎ Ｊ．，Ｊｉ Ｇ．Ｆ．，Ｚｈａｏ Ｆ．，Ｗｕ Ｑ．，Ｘｉａｏ Ｈ．Ｍ．，Ｓｃｉ．Ｃｈｉｎａ Ｓｅｒ．Ｂ，２０１２，５５，２５８７－－２５９４ Ｘｉａｏ Ｊ．Ｊ．，Ｗａｎｇ Ｗ．Ｒ．，Ｃｈｅｎ Ｊ．，Ｊｉ Ｇ．Ｆ．，Ｚｈｕ Ｗ．，Ｘｉａｏ Ｈ．Ｍ．，Ｐｈｙｓｉｃａ Ｂ，２０１２，４０７，３５０４－－３５０９ Ｘｉａｏ Ｊ．Ｊ．，Ｗａｎｇ Ｗ．Ｒ．，Ｃｈｅｎ Ｊ．，Ｊｉ Ｇ．Ｆ．，Ｚｈｕ Ｗ．，Ｘｉａｏ Ｈ．Ｍ．，Ｃｏｍｐｕｔ．Ｔｈｅｏｒ．Ｃｈｅｍ．，２０１２，９９９，２１—２７ Ｘｉａｏ Ｊ．Ｊ．，“ｓ．Ｙ．，Ｃｈｅｎ Ｊ，，Ｊｉ Ｇ．Ｆ．，Ｚｈｕ Ｗ．，Ｚｈａｏ Ｆ．，Ｗｕ Ｑ．，Ｘｉａｏ Ｈ．Ｍ．， Ｍｏ１．Ｍｏｄｅｌ，２０１３，１９，８Ｏ３＿８０９ Ｄａｎｇ Ｈ．Ｓ．，Ｚｈｏｕ Ｆ．Ｆ．，Ｈｉｇｈ Ｅｘｐｌｏｓｉｖｅ ａｎｄ Ｒｅｌａｔｅｄ Ｍａｔｅｒｉａｌｓ ａｎｄ Ｐｅｒｆｏｒｍａｎｃｅ，Ｓｃｉｅｎｃｅ Ｐｒｅｓｓ，Ｂｅｉｊｉｎｇ，１９８９，２５２—２５３，２６６—－２６７ （董海山，周芬芬．高能炸药及其相关性能，北京：科学出版社，１９８９，２５２—２５３，２６６—２６７） ［１９］ Ｘｉａｏ Ｈ．Ｍ．，ＺｈｕＷ．Ｈ．，Ｘｉａｏ Ｊ．Ｊ．，ＰｅｉＸ．Ｑ．，Ｃｈｉｎ． Ｅｎｅｒｇ．Ｍａｔｅｒ．，２０１２，２０（５），５ｌ４＿５２７（肖鹤鸣，朱卫华，肖继军，王 桂香，裴晓琴．含能材料，２０１２，２０（５），５ｌ４—５２７） ［２０］Ｈａｒｒｉｓ Ｎ．Ｊ．，Ｌａｍｍｅｒｔｓｍａ Ｋ．，ｌ，．Ａｍ．Ｃｈｅｍ．Ｓｏｃ．，１９９７，１１９，６５８３－－－６５８９ ［２１］Ｘｉａｏ Ｈ．Ｍ．，Ｘｕ Ｘ．Ｊ．，Ｑｉｕ Ｌ．，Ｔｈｅｏｒｅｔｉｃａｌ Ｄｅｓｉｇｎ ｏｆＨｉｇｈ Ｅｎｅｒｇｙ Ｄｅｎｓｉｔｙ Ｍａｔｅｒｉａｌｓ，Ｓｃｉｅｎｃｅ Ｐｒｅｓｓ，Ｂｅｉｊｉｎｇ，２００８，１８７（肖鹤鸣，许 晓娟，邱玲．高能量密度材料的理论设计，北京：科学出版社，２００８，１８７） ａｏ Ｈ．Ｍ．，Ｑｕａｎｔｕｍ Ｃｈｅｍｉｓｔｒｙ ｏｆＥｎｅｒｇｅｔｉｃ Ｃｒｙｓｔａｌｓ，Ｓｃｉｅｎｃｅ Ｐｒｅｓｓ，Ｂｅｉｊｉｎｇ，２０１２，１９９，２１２（朱卫华，肖鹤鸣．高能晶 ［２２］ Ｚｈｕ Ｗ．Ｈ．，Ｘｉ体量子化学，北京：科学出版社，２０１２，１９９，２１２） ａｏ Ｊ．Ｊ．，Ｚｈｕ Ｗ．Ｈ．，Ｚｈｕ Ｗ．，Ｘｉａｏ Ｈ．Ｍ．，Ｍｏｌｅｃｕｌａｒ Ｄｙｎａｍｉｃｓ ｆＥｎｅｒｏｇｅｔｉｃ Ｍａｔｅｒｉａｌｓ，Ｓｃｉｅｎｃｅ Ｐｒｅｓｓ，Ｂｅｉｊｉｎｇ，２０１３，５４，６７（肖 ［２３］ Ｘｉ继军，朱卫华，朱伟，肖鹤鸣．高能材料分子动力学，北京：科学出版社，２０１３，５４，６７） ［２４］ Ｓｕｎ Ｈ．， Ｐｈｙｓ．Ｃｈｅｍ．Ｂ，１９９８，１０２，７３３８—＿７３６４ 儿 ＂ ”粥 Ｎｏ．１１ ２５ ２６ 刘冬梅等：不同温度下ＨＭＸ和ＲＤＸ晶体的感度判别和……预估——分子动力学比较研究 ３Ｏ ２５６５ Ｒｉｇｂｙ Ｄ．，Ｓｕｎ Ｈ．，Ｅｉｃｈｉｎｇｅｒ Ｂ．Ｅ．，Ｐｏｌｙｍ．Ｉｎｔ．，１９９７，４４，３１ｌ Ｓｕｎ Ｈ．．Ｒｅｎ Ｐ．，Ｆｒｉｅｄ Ｊ．Ｒ．，Ｃｏｍｐｕｔ．Ｔｈｅｏｒ．Ｐｏｌｙｍ．Ｓｃｉ．，１９９８，８，２２９—２４６ ２７ ２８ Ｂｕｎｔｅ Ｓ．Ｗ．．Ｓｕｎ Ｈ．．Ｊ．Ｐｈｙｓ．Ｃｈｅｍ．Ｂ，２０００，１０４，２４７７－－２４８９ Ｘｉａｏ Ｊ．Ｊ．，Ｆａｎｇ Ｇ．Ｙ．，Ｊｉ Ｇ．Ｆ．，Ｘｉａｏ Ｈ．Ｍ．，Ｃｈｉｎ．Ｓｃｉ．Ｂｕｌ１．，２００４，４９（２４），２５２０－－２５２４（肖继军，方国勇，姬广富，肖鹤鸣 科学通报，２００４，４９（２４），２５２０－－２５２４） Ｊ．Ｓ．， ［２９］ Ｘｕ Ｘ．Ｊ．，ＸｉａｏＨ．Ｍ．，Ｘｉａ０ Ｊ．Ｊ．，ＺｈｕＷ．，ＨｕａｎｇＨ．，Ｌｉｕ Ｌ．，Ｘｉａｏ Ｈ．Ｍ．，Ｚｈｕ Ｗ．Ｈ．，Ｘｉａｏ Ｊ．Ｊ．， ［３Ｏ］ ＯｉＰｈｙｓ．Ｃｈｅｍ．Ｂ，２００６，Ｉ１０，７２０３－－７２０７ Ｐｈｙｓ．Ｃｈｅｍ．Ｂ，２００６，１１０，１０６５１－－１０６６１ Ｃｈｏｉ Ｃ．Ｓ．．Ｐｒｉｎｃｅ Ｅ．．Ａｃｔａ Ｃｒｙｓｔａｌｌｏｇｒ．Ｓｅｃｔ．Ｂ：Ｓｔｒｕｃｔ．Ｓｃｉ．，１９７２，２８（９），２８５７－－２８６２ Ｃｈａｎｇ Ｓ．Ｃ．，Ｈｅｎｒｙ Ｐ．Ｂ．，Ａｃｔａ Ｃｒｙｓｔａｌｌｏｇｒ．Ｓｅｃｔ．Ｂ：Ｓｔｒｕｃｔ．Ｓｃｉ．，１９７０，２６，１２３５—１２４０ Ａｎｄｅｒｓｅｎ Ｈ．Ｃ．， Ｃｈｅｍ．Ｐｈｙｓ．，１９８０，７２（４），２３８４－－２３９３ ＰａｒｒｉｎｅｌｌｏＭ．，Ｒａｈｍａｎ Ａ．， Ａｐｐ１．Ｐｈｙｓ．，１９８１，５２（１２），７１８２－－７１９０ Ａｌ１ｅｎ Ｍ．Ｐ．．ＴｉｌｄｅｓｌｅＹ Ｄ．Ｊ．．Ｃｏｍｐｕｔｅｒ Ｓｉｍｕｌａｔｉｏｎ ｏｆ Ｌｉｑｕｉｄｓ，Ｏｘｆｏｒｄ Ｕｎｉｖｅｒｓｉｔｙ Ｐｒｅｓｓ，Ｏｘｆｏｒｄ，１９８７，１５６ Ｋａｍｌｅｔ Ｍ．Ｊ．．Ａｄｏｌｐｈ Ｈ．Ｇ．，Ｐｒｏｐｅｌｌａｎｔｓ Ｅｘｐｌｏｓ．Ｐｙｒｏｔｅｃｈ．，１９７９，４，３０—３４ Ｏｗｅｎ￥Ｆ．Ｊ．，Ｊａｙａｓｕｒｉｙａ Ｋ．，Ａｂｒａｈｍｓｅｎ Ｌ．，Ｃｈｅｍ．Ｐｈｙｓ．Ｌｅｔｔ．，１９８５，１１６，４３４－－－４３８ Ｍｕｌｌａｙ Ｊ．．Ｐｒｏｐｅｌｌａｎｔｓ Ｅｘｐｌｏｓ．Ｐｙｒｏｔｅｃｈ．，１９８７，１２，１２１—１２４ Ｓｅｎｓｉｔｉｖｉｔｙ Ｃｒｉｔｅｒｉｏｎ ａｎｄ Ｍｅｃｈａｎｉｃａｌ Ｐｒｏｐｅｒｔｉｅｓ Ｐｒｅｄｉｃｔｉｏｎ ｏｆ ＨＭＸ ａｎｄ ＲＤＸ Ｃｒｙｓｔａｌｓ ａｔ Ｄｉｆｆｅｒｅｎｔ Ｔｅｍｐｅｒａｔｕｒｅｓ ———Ｃｏｍｐａｒａｔｉｖｅ Ｓｔｕｄｙ ｗｉｔｈ Ｍｏｌｅｃｕｌａｒ Ｄｙｍａｍｉｃｓ Ｓｉｍｕｌａｔｉｏｎ ＬＩＵ Ｄｏｎｇ．Ｍｅｉ。ＺＨＡＯ Ｌｉ ，ＸＩＡＯ Ｊｉ．Ｊｕｎ ，ＣＨＥＮ Ｊｕｎ ，ＪＩ Ｇｕａｎｇ－Ｆｕ ，ＺＨＵ Ｗｅｉ ， 。ＺＨＡＯ Ｆｅｎｇ ，ＷＵ Ｑｉａｎｇ ，ＸＩＡＯ Ｈｅ—Ｍｉｎｇ” （１．Ｍｏｌｅｃｕｌｅ ａｎｄ Ｍａｔｅｒｉａｌ Ｃｏｍｐｕｔａｔｉｏｎ Ｉｎｓｔｉｔｕｔｅ，Ｓｃｈｏｏｌ ｏｆ Ｃｈｅｍｉｃａｌ Ｅｎｇｉｎｅｅｒｉｎｇ， №研ｎｇ Ｕｎｉｖｅｒｓｉｔｙ ｏｆＳｃｉｅｎｃｅ ａｎｄ Ｔｅｃｈｎｏｌｏｇｙ，Ｎａ ｎｇ ２１００９４，Ｃｈｉｎａ； ２．Ｎａｔｉｏｎａｌ Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｓｈｏｃｋ Ｗａｖｅ ａｎｄ Ｄｅｔｏｎａｔｉｏｎ Ｐｈｙｓｉｃｓ，Ｉｓｔｎｉｔｕｔｅ ｏｆ Ｆｌｕｉｄ Ｐｈｙｓｉｃｓ， Ｃｈｉｎａ Ａｃａｄｅｍｙ ｆ ｏＥｎｇｉｎｅｅｒｉｎｇ ｐ ｉｃｓ，Ｍｉａｎｙａｎｇ ６２１９００，Ｃｈｉｎａ； ３．Ｃｏｌｌｅｇｅ ｏｆ Ｂｉｏｌｏｇｉｃａｌ，Ｃｈｅｍｉｃａｌ Ｓｃｉｅｎｃｅｓ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ，Ｊｉａｘｉｎｇ Ｕｎｉｖｅｒｓｉｔｙ，Ｊｉａｘｉｎｇ ３ １４００１，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ Ｉｎ ｏｒｄｅｒ ｔｏ ｅｘｐｌｏｒｅ ａｎｄ ｃｏｍｐａｒｅ ｔｈｅ ｐｒｏｇｒｅｓｓｉｖｅ ｃｈａｎｇｅｓ ｏｆ ｓｔｒｕｃｔｕｒｅ，ｅｎｅｒｇｅｔｉｃ ａｎｄ ｍｅｃｈａｎｉｃａｌ ｐｒｏｐｅｒｔｉｅｓ ｏｆ ＲＤＸ（ｃｙｅｌｏｔｒｉｍｅｔｈｙｌｅｎｅ ｔｒｉｎｉｔｒａｍｉｎｅ）ａｎｄ ＨＭＸ（ｃｙｃ１ｏｔｅｔｒａｍｅｔｈｙ１ｅｎｅ ｔｅｔｒａｎｉｔｒａｍｉｎｅ）ｃｒｙｓｔａｌｓ ｗｉｔｈ ｔｅｍｐｅｒａｔｕｒｅ．ｍｏｌｅｃｕｌａｒ ｄｙｎａｍｉｃｓ ｓｉｍｕｌａｔｉｏｎ ｗａｓ ｃｏｎｄｕｃｔｅｄ ｆｏｒ ｔｈｅｉｒ ｓｕｉｔａｂｌｅ ｍｏｄｅｌｓ ｃｏｎｔａｉｎｉｎｇ ｔｈｅ ｓａｍｅ Ｂｕｍ— ｂｅｒ ｏｆ ａｔｏｍｓ．ＣＯＭＰＡＳＳ ｆｏｒｃｅ ｆｉｅｌｄ ｗａｓ ａｐｐｌｉｅｄ ｉｎ ｔｈｅ ｓｉｍｕｌａｔｉｏｎ ｉｎ Ｎ ｅｎｓｅｍｂｌｅ．Ｆｉｖｅ ｄｉｆｆｅｒｅｎｔ ｔｅｍｐｅｒａ— ｔｕｒｅｓ。ｉ．ｅ．，１９５，２４５，２９５，３４５ ａｎｄ ３９５ Ｋ，ｗｅｒｅ ｃｏｎｓｉｄｅｒｅｄ．Ｔｈｅ ｓｉｍｕｌａｔｉｏｎ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｍａｘｉｍｕｍ ｂｏｎｄ ｌｅｎｇｔｈ（Ｌ ）ｏｆ ｔｈｅｉｒ Ｎ—Ｎ ｔｒｉｇｇｅｒ ｂｏｎｄ ｉｎｃｒｅａｓｅｓ，ｔｈｅ ｉｎｔｅｒａｃｔｉｏｎ ｅｎｅｒｇｙ（ＥＮ—Ｎ）ｂｅｔｗｅｅｎ ｔｗｏ Ｎ ａｔｏｍｓ ｏｆ ｔｈｅ Ｎ—Ｎ ｔｉｒｇｇｅｒ ｂｏｎｄ ｄｅｃｒｅａｓｅｓ ａｎｄ ｔｈｅ ｃｏｈｅｓｉｖｅ ｅｎｅｒｇｙ ｄｅｎｓｉｔｙ（ＣＥＤ）ｄｅｃｒｅａｓｅｓ ａｓ ｔｈｅ ｔｅｍｐｅｒａｔｕｒｅ ｉｎ— ｃｒｅａｓｅｓ．Ｔｈｅｓｅ ｒｅｓｕｉｔｓ ａｇｒｅｅ ｗｉｔｈ ｔｈｅ ｅｘｐｅｒｉｍｅｎｔａｌ ｆａｃｔ ｔｈａｔ ｂｏｔｈ ＲＤＸ ａｎｄ ＨＭＸ ｂｅｃｏｍｅ ｍｏｒｅ ｓｅｎｓｉｔｉｖｅ ａｓ ｔｈｅ ｔｅｍｐｅｒａｔｕｒｅ ｉｎｃｒｅａｓｅｓ．Ｍｏｒｅｏｖｅｒ，ｔｈｅ Ｌ ａｘ ｖａｌｕｅ ｏｆ ＲＤＸ ｉｓ ｌａｒｇｅｒ ｔｈａｎ ｔｈａｔ ｏｆ ＨＭＸ，ａｎｄ ｔｈｅ ｖａｌｕｅｓ ｏｆ ＥＮＮ ａｎｄ ＣＥＤ ｏｆ ＲＤＸ ａｒｅ ｌｏｗｅｒ ｔｈａｎ ｔｈｏｓｅ ｏｆ ＨＭＸ ａｔ ａｌｌ ｔｅｍｐｅｒａｔｕｒｅｓ．Ｔｈｅｓｅ ａｌｓｏ ａｇｒｅｅ ｗｉｔｈ ｔｈｅ ｆａｃｔ ｔｈａｔ ＲＤＸ ｉｓ ｍｏｒｅ ｓｅｎｓｉｔｉｖｅ ｔｈａｎ ＨＭＸ．Ｉｔ ｉｓ ｔｈｕｓ ｓｕｇｇｅｓｔｅｄ ｔｈａｔ ｔｈｅ Ｌ ，ＥＮ—Ｈ ａｎｄ ＣＥＤ ｉｎｄｅｅｄ ｃａｎ ｂｅ ｕｓｅｄ ａｓ ｔｈｅｏｒｅｔｉｃａｌ ｃｒｉｔｅｒｉａ ｔｏ ｐｒｅｄｉｃｔ ｔｈｅ ｒｅｌａｔｉｖｅ ｄｅｇｒｅｅ ｏｆ ｈｅａｔ ａｎｄ ｉｍｐａｃｔ ｓｅｎｓｉｔｉｖｉｔｙ ｏｆ ｈｉｇｈ ｅｎｅｒｇｙ ｍａｔｅｒｉａｌｓ ｕｎｄｅｒ ｃｅｒｔａｉｎ ｃｏｎ— ｄｉｔｉｏｎｓ．Ｔｈｅ ｍｅｃｈａｎｉｃａｌ ｐｒｏｐｅｒｔｉｅｓ ｗｅｒｅ ｏｂｔａｉｎｅｄ ｗｉｔｈ ｓｔａｔｉｃ ａｎａｌｙｓｉｓ ｍｅｔｈｏｄ．Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｅｎｇｉ－ ｎｅｅｒｉｎｇ ｍｏｄｕｌｕｓ ｓｕｃｈ ａｓ ｔｅｎｓｉｌｅ ｍｏｄｕｌｕｓ，ｂｕｌｋ ｍｏｄｕｌｕｓ ａｎｄ ｓｈｅａｒ ｍｏｄｕｌｕｓ ｄｅｃｒｅａｓｅ ｗｉｔｈ ｔｈｅ ｉｎｃｒｅａｓｉｎｇ ｔｅｒｎ— ｐｅｒａｔｕｒｅ，ｗｈｉｃｈ ａｇｒｅｅｓ ｗｉｔｈ ｔｈｅ ｅｘｐｅｒｉｍｅｎｔｌ ｒｅｓａｕｌｔｓ． Ｋｅｙｗｏｒｄｓ Ｃｙｃｌｏｔｒｉｍｅｔｈｙｌｅｎｅ ｔｉｒｎｉｔｒａｍｉｎｅ；Ｃｙｃｌｏｔｅｔｒａｍｅｔｈｙｌｅｎｅ ｔｅｔｒａｎｉｔｒａｍｉｎｅ；Ｍｏｌｅｃｕｌａｒ ｄｙｎａｍｉｃｓ（ＭＤ） ｓｉｍｕｌａｔｉｏｎ；Ｓｅｎｓｉｔｉｖｉｔｙ；Ｍｅｃｈａｎｉｃａｌ ｐｒｏｐｅｒｔｙ （Ｅｄ．：Ｙ，Ｚ，Ａ）
