Phase Retrieval and Reconstruction of Coherent Diffraction Imaging
-
摘要: 相干衍射成像是一种对材料体密度敏感的超高分辨成像技术。相较于传统表面敏感的超高分辨成像技术，相干衍射成像利用了硬X射线的强穿透能力，可以深入材料体内部进行成像，且成像分辨能力可以根据成像布局进行调整，最高达到原子级空间分辨能力。这种灵活的空间分辨调整依赖于相干衍射成像独特的相位复原技术，即通过对图像成像强度的过采样，利用含约束的迭代算法同时获得光场的强度及相位，进而对样品进行重建；同时结合图像定向及组合技术，相干衍射成像可以实现对样品的三维重建。本文主要从成像原理、复原算法和重建方法介绍相干衍射成像技术，并结合实验进展及模拟研究展示该技术在多种重建情形下具备的诊断能力，以期较为全面地给出相干衍射成像技术的发展趋势。Abstract: The coherent diffraction imaging (CDI) is an ultra-high resolution imaging technique that is sensitive to the density of the material. Compared to the surface-sensitive imaging methods with ultra-high resolution, the CDI is able to probe the interior of the sample by taking advantages of hard X-rays. According to the imaging layout, the space resolution of CDI is variable and can reach up to an atomic scale. This feature depends on the iterative phase retrieval method that almost becomes the signature of CDI. Based on oversampling a sample in a detected image, the phase and intensity of X-ray beam can be retrieved simultaneously by iterative calculations with constraints, and then are used to reconstruct the sample. Meanwhile, the three-dimensional reconstruction could be realized by combining image orientating and merging techniques. Here we present the imaging theory, phase retrieval and reconstruction methods of the CDI technique, and its diagnostic ability in a variety of reconstruction situations by experimental and simulation examples, to hopefully provide a systematic introduction of its development.
-
管道运输是近现代主要运输方式之一。在石油工业促进下管道运输逐步发展成大范围、大跨度、大容量、多介质运输网络系统。然而，复杂管网系统可能面临诸如维护成本高、检测周期长、工艺更新快、极端环境影响加大的问题，致使管路系统发生失效破坏，造成严重的经济损失和环境污染。其中由疲劳引起的失效破坏因预测难度大、事故影响范围广等特点成为研究热点。
在工程结构中疲劳是主要的失效机理之一[1]。在材料科学中疲劳是循环力作用下材料产生的失效。其机理可以描述为：当载荷超过某极限时会在应力集中处产生微裂纹，造成材料的几何不连续[2]，并逐步发展至不稳定状态而导致突然失效[3]。在ASTM标准中使用疲劳寿命(Nf)表征材料在自然失效前所能承担的应力循环次数[4]。对于一些材料，此参数即为低应力幅状态下判断材料是否失效的理论值[5]。这里，裂纹成为研究重点，前人对裂纹缺陷与疲劳失效之间的关系做了大量研究[6-9]。其中Hussain[10]在1997年对已提出的一些短裂纹疲劳分析模型和观点进行了对比和评判。在石油和天然气管道运输中，运输介质的特殊性和管道铺设的地理特殊性使工况变得异常复杂，相关研究包含：水土侵蚀后悬空管的疲劳失效及影响因素[11]，悬跨管路中波浪、海流和不同裂纹凹痕对寿命的影响[12]，从模态、位移、转角响应和动应力幅值方面对冰振下天然气管线进行分析[13]等。同时还研究了管系中非常温运输介质的动态变化对裂纹萌生的影响，并对裂纹的萌生方向进行了理论预测，依照Fatemi-Socie准则和Miner准则进行了解析计算[14]。
到目前为止，管道疲劳计算主要运用针对不同裂纹形式所建立的模型进行数值计算，同时就不同法兰管件在不同载荷和裂纹形式下进行力学模拟和分析，所涉及的疲劳计算软件有MSC.Fatigue、nCode DesignLife、CAE SAR Ⅱ、ANSYS。为模拟实际工况，软件间的耦合同样是研究热点，例如：利用Fluent和Abaqus耦合分析超音速进气道流振动[15]，通过ANSYS和Franc 3D联合分析含裂纹的核级管道[16]，使用ANSYS Workbench和CFX估算飞机液压管道的动态特性[17]。
当前研究主要根据实际管道中发生破坏的管道部件进行失效验证和失效机理模拟，对于疲劳的预测主要集中在理论推导和拟合，而对整体管道的疲劳模拟和预测比较缺乏。主要原因是模型生成软件、力学计算软件和疲劳分析软件三者之间的信息传递存在不兼容现象。ANSYS Workbench平台很好地解决了此问题，它包含自底向上的几何模型创建模块、静态和动态的力学分析模块、通用的流体计算软件和nCode DesignLife疲劳分析软件，能在一个工程中通过简单操作将不同的计算模块联系起来进行耦合计算，解决了因复杂模型造成的大量数据在不同软件之间传递的困难。
本研究针对管道系统的疲劳分析，利用ANSYS Workbench的各种模块，并结合单向流固耦合方法，研究复杂管道内流体及其激振力对整个管道系统疲劳的影响。
1. 数值模拟
1.1 模块选择与模型建立
本研究的计算模块包含：有限元基本模块、Fluent计算模块和nCode DesignLife疲劳模块。其中流体和固体部分采用单向流固耦合方法，并将计算数据与疲劳软件共享，得到理想的疲劳分析结果。
直管尺寸参考ASME B36.19M或ASME B36.10M标准，弯头(Elbows A & B)、支座三通(Supported Tee)、四通(Four-Way)、异径三通(Reducer Tee)、法兰(Flange)、法兰垫片(Spacer)参考ASME B16.1标准，选用250磅级，其中法兰片中螺栓参考ASME B18.2.1重型六角头螺栓，螺母参考ASME B18.2.2重型六角螺母。具体尺寸为：法兰直径444.5 mm，法兰厚度47.752 mm，凸台面直径357.124 mm；毂直径320.548 mm，毂面贯穿长度60.452 mm；螺栓件数量为16组，螺栓直径25.4 mm；弯头及三通中基础面到中心的距离266.7 mm，基础面尺寸254.0 mm，基础面厚度31.75 mm，肋板厚度22.352 mm，螺栓孔排列半径和孔直径200.152 mm；直管中弯头中心到基本面的距离292.1 mm，长弯头中心到基本面的距离419.1 mm。
1.2 装配与网格划分
ANSYS Workbench中通过Geometry模块完成模型装配，得到如图 1所示的管道系统。定义四通(Four-way)、三通出口(Tee-Out)及B类弯头(Elbow B)所在平面为x-z平面。
整体网格节点数为4 814 382，单元数为16 531 216，单元类型为四面体(Tet 4)和三棱柱(Wed 6)。为判定网格质量，统计了网格正交品质和倾斜度，结果如图 2和图 3所示。
1.3 参数确定
参数主要为求解方程中的常量，与模型、工况和求解方法密切相关。对于结构计算部分，根据GB 50316-2008工业金属管道设计规范和ASME B16.5标准(工业管道压力分级：10~42 MPa为高压，大于42 MPa为超高压)，结合本研究中流体为液态水，选择较低的压力12 MPa作为工况压力。剩余参数单独讨论。
(1) 流体计算参数
将管道主体模型抑制后，Fluent软件内物理模型选择Realizable k-ε湍流模型，材料属性中选择water-liquid(h20〈1〉)作为流体材料，计算域中的实体均选择fluid。求解方法选用SIMPLEC方法。速度入口边界条件需要确定雷诺数Re、湍流强度I和湍流尺度λ，由下式确定
Re=ρv0Dη,I≈0.16(ReDH)1−/8,λ=0.07L (1) 式中：ρ为密度，v0为入口速度，D为几何限度(直径)，η为流体黏度，L为特征长度(本研究中为管道直径)，ReDH表示基于水力直径的雷诺数。
(2) 疲劳计算参数
分析疲劳时，首先需要导入载荷曲线，即力-时间或振幅-时间曲线，通常由实验获得。本研究着眼于不同子系统与整体系统之间疲劳的相互影响，而通过实验获得与本系统匹配良好的载荷曲线较为困难，为此人工导入常见的白噪声载荷曲线(见图 4，其中纵轴振幅表示加载应力与最大应力的比值)，在所设计的工况下获得接近实际的疲劳分析结果。
模拟使用的算法必须与材料模型中提供的S-N曲线(一定循环特征下标准试件的疲劳强度与疲劳寿命的关系曲线，也称应力-寿命曲线)类型一致，且因系统的复杂性，应力组合方法选择带符号的von-Mises等效应力(Signed von Mises)，平均应力修正使用FKM，即使用4个系数定义平均应力的敏感程度。小循环修正选择BS7608，即最大应力循环高于S-N曲线的NC1(106)时，S-N斜率关系式变为b2=-b1/(2b1-1)，否则b2=0，且损伤为零。图 5给出了本研究所使用材料灰铸铁(Grey Cast Iron BS1452 Grade250)和结构钢(Structural Steel BS4360 Grade 40B)的标准S-N曲线，图中UTS表示抗拉极限强度，S0为疲劳极限强度，b1、b2为斜率, σ为应力，n为循环次数。
2. 结果与讨论
计算结果包含流体、结构和疲劳3部分。流体结果主要包含管路系统各部分流速云图。结构计算结果主要包含流体映射力、应力和应变最值及云图。疲劳结果包含管系寿命云图、关键部件和子系统的应力-循环次数曲线。
2.1 流场结果分析
已知工业管道系统的正常流速为1.5~3.0 m/s。首先通过流体计算结果选择分析工况。不同入口速度v0下管道子系统流速云图如图 6所示，其中v为流速。
图 6 不同入口速度下管路子系统流速Figure 6. Flow rate of pipeline subsystem for different inlet velocities可以看出：入口速度为1 m/s时子系统内流速仅0.36 m/s，远低于要求的正常流速；当入口速度为10 m/s时，子系统流速达到上限要求。因此，以下主要讨论入口速度为3、5、7和10 m/s时的疲劳。
2.2 结构计算结果分析
静力学计算包含映射力导入、管道内部压力设置、边界条件设置。本研究主要针对管道系统在随机激振力作用下的疲劳响应，为避免连接工艺造成的疲劳失效，管路连接部分设置为固支连接。映射力、结构应力最大值σmax和应变最大值εmax的计算结果见表 1和表 2，部分部件的应力云图见图 7和图 8。
表 1 管道系统内流体对结构的映射力Table 1. Mapping force of fluid in pipeline system to structureComponent v0/(m·s-1) Mapping force/N x y z Elbow A 3 1.86 124.00 -147.94 5 282.56 39.96 -332.05 7 561.72 142.27 -707.51 10 5 677.60 -2 269.10 -797.02 Elbow B 3 -0.32 -44.00 2.26 5 865.42 -1 068.40 0.24 7 1 745.20 -2 095.90 -9.68 10 18 628.00 -18 369.00 580.01 Tee-support 3 136.32 -0.23 -3.15 5 682.44 0.04 -3.56 7 1 093.30 -1.14 -0.15 10 2 715.60 -6.83 -32.04 Tee-out 3 472.57 -0.01 343.97 5 1 707.60 -0.65 1 092.30 7 3 031.30 -0.24 2 002.40 10 13 840.00 -20.37 9 474.50 Four-way 3 -8.02 0.25 -186.31 5 -7.74 1.70 -615.87 7 -19.67 3.66 -1 173.30 10 126.58 13.93 -2 455.70 表 2 计算结果Table 2. Calculated resultsv0/(m·s-1) σmax/MPa εmax/10-3 Position with maximum stress 3 285.34 2.28 Tee-out 5 285.37 2.28 Tee-out 7 285.30 2.28 Tee-out 10 287.11 2.29 Tee-out 由表 1数据可知：不同入口速度条件下异径三通出口x-z平面上的映射力远大于其他部件，说明在系统中此部分的疲劳寿命可能较其他部件低；通过法兰部件改变流体流动方向时，法兰部件所在平面两个方向上的映射力远大于另一方向，说明法兰部件较其他管道部件更容易发生疲劳失效或寿命降低。
从图 7和图 8显示的入口流速为5 m/s时四通和三通出口的静态应力计算结果可以看出，应力最大值发生在法兰管件内部的界面变化处，其余部分均未达到屈服应力。这说明所选择的应力疲劳分析方法是可行的，并且管系连接件可能是发生疲劳的重点部位。
从表 1和表 2数据不难看出，应力最大值(σmax)和应变最大值(εmax)受管系初始流体速度的影响较小，因此管路系统整体疲劳寿命结果的影响因素可能是疲劳分析时向管道内施加的随机激振力。
2.3 疲劳计算结果分析
疲劳计算结果包含整体系统和重点部件的疲劳寿命云图，以及各子系统和各部件的应力-疲劳寿命(σ-lg n)曲线，用以和标准S-N曲线(见图 5)对比。
2.3.1 疲劳云图分析
图 9为管道整体系统的疲劳云图，图 10、图 11和图 12分别为四通、出流三通和支座三通的疲劳云图。
可以看出：在整体管道系统中未发生疲劳失效破坏，说明管路在本研究所设条件下可以正常运行；但是在出流三通、四通、支座三通位置出现了疲劳寿命极低值，与结构计算中预测的可能出现寿命低值的位置和最终模拟位置相同；最小寿命为1 948，大于抗拉极限强度(UTS)寿命1 000，再次验证了之前选用应力疲劳模型的正确性。
2.3.2 应力-寿命曲线分析
通过在nCode DesignLife软件中提取一定数量模型结点，将这些结点包含的应力和寿命数据记录下来，获得对应部分的S-N曲线，其中包含A类弯头、B类弯头、支座三通、出流三通和四通，如图 13~图 17所示，并与标准S-N曲线(见图 5)进行对比。图 13~图 17所示曲线的变化规律可大致归纳为：不同部件和子系统在相同循环区间内呈现的变化趋势相同，不同入口流速条件下曲线间未发生明显偏移，当循环次数的对数值在3~7范围内时曲线整体呈现非线性变化趋势。
标准S-N曲线(见图 5)中，灰铸铁和结构钢材料曲线在103~107的循环区间内呈线性变化。而模拟结果中S-N曲线呈现下凹的非线性变化，说明疲劳寿命对应力的敏感度较标准值从常量变为变量，即在寿命区间前段应力变化引起的疲劳变化较小，而在寿命区间后半段应力变化引起的疲劳变化较大。同时也发现抗拉极限强度和疲劳极限强度所对应的寿命得到较大的增强，但在极值之间的寿命却并非等值增加。
从图 18中不难发现，在入口流速为5 m/s的条件下，5种法兰管件在不同取值区间上都满足整体S-N曲线变化趋势，且与单个部件的S-N曲线变化趋势基本相同。说明在管道系统中各部件与整体拥有基本相同的特征曲线，验证了系统对部件的影响同样能反作用于整体。
3. 结论
(1) 管道系统中只要保证各个子系统内流速处于合理区间内，入口流速对整体系统的影响很小，因此工业实际中需要对各部分进行流速监控，保证其处在合理可控范围之内。
(2) 管道系统内的连接部是发生疲劳失效的重点区域，其中连接不同子系统的三通和四通类法兰管件是影响整体疲劳寿命的关键，并且其疲劳失效多发生于相贯线处的管道内部，需要重点关注。
(3) 管路系统有助于增加内部各个法兰管件在极限状态下的疲劳寿命，降低疲劳敏感度。因此各个法兰管件在短期内可以承受更大的应力变化，在寿命的中后期能更快地趋于稳定，即疲劳极限强度得以提高，但在正常工作条件下对疲劳寿命的增强没有极限状态明显。
-
图 8 模拟的超高分辨二维CDI实验
Figure 8. Simulation of 2-dimensional CDI experiment with ultra-high resolution
-
[1] CHAO W, HARTENECK B D, LIDDLE J A, et al. Soft X-ray microscopy at a spatial resolution better than 15 nm [J]. Nature, 2005, 435(7046): 1210. doi: 10.1038/nature03719 [2] BARBER J L, BARNES C W, SANDBERG R L, et al. Diffractive imaging at large fresnel number: challenge of dynamic mesoscale imaging with hard X-rays [J]. Physical Review B, 2014, 89(18): 184105. doi: 10.1103/PhysRevB.89.184105 [3] XIAO X H, SHEN Q. Wave propagation and phase retrieval in fresnel diffraction by a distorted-object approach [J]. Physical Review B, 2005, 72(3): 033103. doi: 10.1103/PhysRevB.72.033103 [4] MIAO J W, AMONETTE J E, NISHINO Y, et al. Direct determination of the absolute electron density of nanostructured and disordered materials at sub-10-nm resolution [J]. Physical Review B, 200, 68(1): 012201. [5] SAYER D. Some implications of a theorem due to shannon [J]. Acta Crystallographica, 1952, 5: 843. [6] GERCHBERG R W, SAXTON W O. A practical algorithm for the determination of phase from image and diffraction plane pictures [J]. Optik, 1972, 35: 237. [7] FIENUP J R. Phase retrieval algorithm: a comparison [J]. Applied Optics, 1982, 21: 2758. doi: 10.1364/AO.21.002758 [8] ELSER V. Phase retrieval by iterated projections [J]. Journal of the Optical Society of America, 2003, 20(1): 40. doi: 10.1364/JOSAA.20.000040 [9] CHEN C C, MIAO J, WANG C W, et al. Application of optimization technique to noncrystalline X-ray diffraction microscopy: guided hybrid input-output method [J]. Physical Review B, 2007, 76(6): 064113. doi: 10.1103/PhysRevB.76.064113 [10] LUKE D R. Relaxed averaged alternating reflections for diffraction imaging [J]. Inverse Problems, 2005, 21: 37. doi: 10.1088/0266-5611/21/1/004 [11] MARCHESINI S, HE H, CHAPMAN H N, et al. X-ray image reconstruction from a diffraction pattern alone [J]. Physical Review B, 2003, 68: 140101. doi: 10.1103/PhysRevB.68.140101 [12] MIAO J, SYAER D, CHAPMAN H N. Phase retrieval from the magnitude of the fourier transforms of nonperiodic objects [J]. Josa A, 1998, 15: 1662. doi: 10.1364/JOSAA.15.001662 [13] 周光照, 佟亚军, 陈灿, 等. 相干X射线衍射成像的数字模拟研究 [J]. 物理学报, 2011, 60(2): 028701.ZHOU G Z, TONG Y J, CHEN C, et al. Digital simulation for coherent X-ray diffractive imaging [J]. Acta Physica Sinica, 2011, 60(2): 028701. [14] VARTANYANTS I A, ROBINSON I K. Partial coherence effects on the imaging of small crystals using coherent X-ray diffraction [J]. Journal of Physics: Condensed Matter, 2001, 13(47): 10593. doi: 10.1088/0953-8984/13/47/305 [15] MIAO J W, CHARALAMBOUS P, KIRZ J, et al. Extending the methodology of X-ray crystallography to allow imaging of micrometer-sized non-crystalline specimens [J]. Nature, 1999, 400: 342. doi: 10.1038/22498 [16] MIAO J W, NISHINO Y, KOHNURA Y, et al. Quantitative image reconstruction of GaN quantum dots from oversampled diffraction intensities alone [J]. Physical Review Letters, 2005, 95(8): 085503. doi: 10.1103/PhysRevLett.95.085503 [17] NISHINO Y, TAKAHASHI Y, IMAMOTO N, et al. Three-dimensional visualization of a human chromosome using coherent X-ray diffraction [J]. Physical Review Letters, 2009, 102(1): 018101. doi: 10.1103/PhysRevLett.102.018101 [18] EKEBERG T, SVENDA M, ABERGEL C, et al. Three-dimensional reconstruction of the giant mimivirus particle with an X-ray free-electron laser [J]. Physical Review Letters, 2015, 114(9): 098102. doi: 10.1103/PhysRevLett.114.098102 [19] DUANE N T, ELSER V. Reconstruction algorithm for single-particle diffraction imaging experiments [J]. Physical Review E, 2009, 80(2): 026705. doi: 10.1103/PhysRevE.80.026705 [20] MIAO J W, CHEN C C, SONG C, et al. Three-dimensional GaN-Ga2O3 core shell structure revealed by X-ray diffraction microscopy [J]. Physical Review Letters, 2006, 97(21): 215503. doi: 10.1103/PhysRevLett.97.215503 [21] TAKAHASHI Y, NISHINO Y, TSUTSUMI R, et al. High-resolution projection image reconstruction of thick objects by hard X-ray diffraction microscopy [J]. Physical Review B, 2010, 82(21): 214102. doi: 10.1103/PhysRevB.82.214102 [22] THIBAULT P, DIEROLF M, MENZEL A, et al. High-resolution scanning X-ray diffraction microscopy [J]. Science, 2008, 321(5887): 379. doi: 10.1126/science.1158573 [23] RODENBURG J M, HURST A C, CULLIS A G, et al. Hard-X-ray lensless imaging of extended objects [J]. Physical Review Letters, 2007, 98(3): 034801. doi: 10.1103/PhysRevLett.98.034801 [24] KLAUS G, PIERRE T, SEBASTIAN K, et al. Quantitative biological imaging by ptychographic X-ray diffraction microscopy [J]. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(2): 529. doi: 10.1073/pnas.0905846107 [25] DIEROLF M, MENZEL A, THIBAULT P, et al. Ptychographic X-ray computed tomography at the nanoscale [J]. Nature, 2010, 467(7314): 436. doi: 10.1038/nature09419 [26] ROBINSON I K, VARTANYANTS I A, WILLIAMS G J, et al. Reconstruction of the shapes of gold nanocrystals using coherent X-ray diffraction [J]. Physical Review Letters, 2001, 87(19): 195505. doi: 10.1103/PhysRevLett.87.195505 [27] WILLIAMS G J, PFEIFER M A, VARTANYANTS I A, et al. Three-dimensional imaging of microstructure in Au nanocrystals [J]. Physical Review Letters, 2003, 90(17): 175501. doi: 10.1103/PhysRevLett.90.175501 [28] PFEIFER M A, WILLIAMS G J, VARTANYANTS I A, et al. Three-dimensional mapping of a deformation field inside a nanocrystal [J]. Nature, 2006, 442(7098): 63. doi: 10.1038/nature04867 [29] NEWTON M C, LEAKE S J, HARDER R, et al. Three-dimensional imaging of strain in a single ZnO nanorod [J]. Nature Materials, 2010, 9(2): 279. [30] HARDER R, ROBINSON I. Coherent X-ray diffraction imaging of strain at the nanoscale [J]. Nature Materials, 2009, 8(4): 291. doi: 10.1038/nmat2400 [31] GANG X, OUSSAMA M, MANFRED R, et al. Coherent X-ray diffraction imaging and characterization of strain in silicon-on-insulator nanostructures [J]. Advanced Materials, 2014, 26(46): 7747. doi: 10.1002/adma.v26.46 -
下载:
下载:
点击查看大图
计量
- 文章访问数: 9566
- HTML全文浏览量: 3096
- PDF下载量: 43
