稿件修改说明稿件号201203022稿件名称Al包覆层厚度对Al-NiTi复合.doc

稿件修改说明 稿件号：201203022 稿件名称：Al包覆层厚度对Al-NiTi复合垫片密封面压紧力影响的有限元分析 感谢编辑和审稿人仔细阅读所投稿件，并提出中肯意见。 现对“退改稿”修改做出说明如下： 修改稿件时在“修订”状态进行，对“批注”的回应在修改稿中同样采用“批注”列出。 对“退改稿”中“批注”的详细回应见下（按“退改稿”中出现顺序）： 批注[z1]：尺寸应为φ111.15mm。 答：原文中“图1”此处标注错误，应为φ111.15mm。已修改。 批注[z5]：图2。 答：原文中“图2”误标为“图1”。已修改。 批注[z3]：Fig.2中图c是铝材的压缩性能还是Al-NiTi复合垫片的压缩性能？在考虑弯矩时，计算应采用Al-NiTi复合垫片的压缩-回弹性能曲线。 答：第一问：图2中c是铝材A1050的压缩性能：本文中无试验所得到的Al-NiTi复合垫片的压缩-回弹性能曲线。 第二问：在考虑弯矩时，若将Al-NiTi复合垫片视为整体，则应采用Al-NiTi复合垫片压缩性能曲线进行计算：本文中将Al-NiTi复合垫片视为A1050包覆层和NiTi合金的“组合体”，因此利用A1050压缩性能曲线（图2.c）和NiTi合金压缩性能曲线（图2.a），在Abaqus有限元软件中进行“组装”，从而进行数值分析。 批注[z9]：与金属石墨缠绕垫片相近。 答：接受文字修改。已修改。 批注[z10]：在接触的边缘部位，可能由于加载方式和网格的大小，计算结果有较大差异，能否作为结论的依据，请作者仔细考虑。可以选用不同的单元，同时对垫片进行网格加密，然后和之前的计算结果进行对比分析，排除计算方法和计算模型带来的误差。为什么Al-NiTi复合垫片压紧力呈现波纹状分布？这和单元选择有关吗？ 答：这个问题牵涉到对“边缘效应”和“波纹状”现象的解释和说明。 关于“边缘效应”的解释： “从压紧力数值来看，A1050包覆层沿径向的分布有局部区域超过了材料的压缩屈服强度（σ0.2=30MPa），这些区域的存在导致材料沿径向有“流动”的趋势，在法兰密封面和A1050包覆层表面摩擦属性的影响下，形成沿径向的应力，这些径向应力改变了垂直于包覆层表面的压紧力数值，从而形成波纹现象。”此段文字已添加至修改稿对应段落。 关于“边缘效应”现象的解释和说明： “这种边缘效应的产生源于在δ区域内，仅为A1050材料承受来自于法兰密封面的载荷，使得这一区域的材料表现出图2中（c）图A1050的受压性能。然而这种边缘效应在实际垫片使用中，受包覆工艺影响将会比较大。因为在本文计算模拟中，未考虑法兰密封面载荷作用下δ区域在厚度方向上受载会有沿径向压缩“失稳”现象产生的可能：当包覆层与NiTi合金垫片侧面（垫片环内侧面或外侧面）贴合较好，受载后不分开时，这种边缘效应会表现出较好效果；然而当包覆工艺不能保证较好的贴合效果，δ区域受载后可能出现“压溃”现象，不出现边缘效应。”此段文字已添加至修改稿相应段落。 关于“边缘效应”的进一步解释参见下图： Al-NiTi复合垫片在受载后压缩变形，可能出现两种情况，如上图所示。换句话说，所投稿件中的数值模拟是建立在“受载后-情况2”的条件下进行分析的。 另1：关于“边缘效应”可否作为结论依据：这个问题的确很重要，因此我在“修改稿”相应地方加了条件说明，即“当包覆工艺较为理想时，……”。 另2：关于网格加密：考虑到垫片是重点研究对象，因此在建立计算模型时，已经进行了加密处理。 另3：关于单元类型：这个问题非常好。在Abaqus有限元软件分析中，垫片的属性定义有很多种，最为常见的就是应用gasket单元，输入垫片的压缩-回弹曲线性能参数，就可以进行数值计算。然而gasket单元在Abaqus有限元软件中也有自己的定义范围，比如在所投稿件中应用的GK3D8单元时考虑的是垫片的厚度方向行为（Gasket Thickness Behavior），显然这是做了简化处理的。所投稿件中，A1050包覆层用的C3D8R实体单元，可以反映材料在受载时较为真实的三向力学行为。那么问题是，可不可以在所有计算模型中都采用实体单元（如C3D8R），答案是：比较困难。因为若采用实体单元对金属石墨缠绕垫片进行分析，涉及到所用金属缠绕带（可能还有内外支撑环）和石墨的性能数据，从目前国内外文献来看，仅对金属石墨缠绕带垫片细致的数值模拟分析（不考虑法兰、螺栓）都显得异常困难，最新的研究文献也仅是采用模糊化的方法进行处理，这里可以参见文献。NiTi合金是一种形状记忆合金，在应力作用下的变形行为涉及到材料本身的应力诱发相变行为（马氏体相和奥氏体相的相互转变），若需要细致研究NiTi合金垫片的压缩行为，在Abaqus有限元分析软件中，已不是实体单元能解决的，需要写入UMAT程序。鉴