发明名称 一种光学元件吸收损耗测量的数据处理改进方法
摘要 一种光学元件吸收损耗测量的数据处理改进方法,连续激励光束经透镜聚焦后照射光学元件表面,光学元件因吸收入射激光束能量引起内部温度上升。通过建立一种更能真实反映光学元件所处物理实际的精确温度模型,基于此精确温度模型,通过对一定尺寸不同热物理性质光学元件不同位置温升的数值模拟,得到不同于国际标准ISO11551所采用的最佳温度探测位置。再将温度传感器调至此最佳温度探测位置测量光学元件表面的温度数据,将测量得到的温度数据拟合到国际标准ISO11551中所采用的均匀温度模型得到样品的吸收损耗值。与现有国际标准ISO11551采用固定温度探测位置相比较,本发明提出针对不同待测光学元件尺寸和热物理参数,采用位置可调的温度传感器以实现更精确测量光学元件吸收损耗的目的。
申请公布号 CN102053006A 申请公布日期 2011.05.11
申请号 CN201010535210.2 申请日期 2010.11.03
申请人 中国科学院光电技术研究所 发明人 李斌成;王艳茹;高卫东
分类号 G01M11/00(2006.01)I;G01M11/02(2006.01)I;G01N25/20(2006.01)I 主分类号 G01M11/00(2006.01)I
代理机构 北京科迪生专利代理有限责任公司 11251 代理人 卢纪
主权项 1.一种光学元件吸收损耗测量的数据处理方法,其特征在于实现步骤如下:(1)连续加热激光束经透镜聚焦后入射在待测光学元件表面,所述光学元件因吸收入射激光束能量,引起光学元件内温度梯度分布;(2)建立精确理论温度模型准确描述所述光学元件因受激光照射所引起的温升分布,并且基于此精确温度模型的数值模拟确定待测光学元件的最佳温度探测位置;所述建立的精确温度模型为:一光斑半径为a的加热激光束照射在半径为b,厚度为d的圆柱型光学元件表面,引起温升分布模型为:ΔT(r,z,t)=0                    (t<t<sub>d</sub>)<maths num="0001"><![CDATA[<math><mrow><mi>&Delta;T</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>D</mi><msub><mi>K</mi><mi>th</mi></msub></mfrac><mfrac><mrow><msub><mi>A</mi><mn>0</mn></msub><mi>P</mi></mrow><mrow><mi>&pi;</mi><msup><mi>a</mi><mn>2</mn></msup></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mo>&infin;</mo></munderover><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mo>&infin;</mo></munderover><mfrac><mn>1</mn><mrow><mi>N</mi><mrow><mo>(</mo><msub><mi>&beta;</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>N</mi><mrow><mo>(</mo><msub><mi>&eta;</mi><mi>p</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>J</mi><mn>0</mn></msub><mrow><mo>(</mo><msub><mi>&beta;</mi><mi>m</mi></msub><mi>r</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>&eta;</mi><mi>p</mi></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>&eta;</mi><mi>p</mi></msub><mi>z</mi><mo>)</mo></mrow><mo>+</mo><mi>H</mi><mi>sin</mi><mrow><mo>(</mo><msub><mi>&eta;</mi><mi>p</mi></msub><mi>z</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths>(1)                           (t<sub>d</sub>≤t≤t<sub>0</sub>)<maths num="0002"><![CDATA[<math><mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>D</mi><mrow><mo>(</mo><msubsup><mi>&beta;</mi><mi>m</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&eta;</mi><mi>p</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><msub><mi>&eta;</mi><mi>p</mi></msub><mrow><mi>D</mi><mrow><mo>(</mo><msubsup><mi>&beta;</mi><mi>m</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&eta;</mi><mi>p</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>&CenterDot;</mo><mfrac><msup><mi>a</mi><mn>2</mn></msup><mn>2</mn></mfrac><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msubsup><mi>&beta;</mi><mi>m</mi><mn>2</mn></msubsup><msup><mi>a</mi><mn>2</mn></msup></mrow><mn>8</mn></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>(2)<maths num="0003"><![CDATA[<math><mrow><mi>&Delta;T</mi><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>z</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>D</mi><msub><mi>K</mi><mi>th</mi></msub></mfrac><mfrac><mrow><msub><mi>A</mi><mn>0</mn></msub><mi>P</mi></mrow><mrow><mi>&pi;</mi><msup><mi>a</mi><mn>2</mn></msup></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mo>&infin;</mo></munderover><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mo>&infin;</mo></munderover><mfrac><mn>1</mn><mrow><mi>N</mi><mrow><mo>(</mo><msub><mi>&beta;</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mi>N</mi><mrow><mo>(</mo><msub><mi>&eta;</mi><mi>p</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>J</mi><mn>0</mn></msub><mrow><mo>(</mo><msub><mi>&beta;</mi><mi>m</mi></msub><mi>r</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><mi>&eta;</mi><mi>p</mi></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>&eta;</mi><mi>p</mi></msub><mi>z</mi><mo>)</mo></mrow><mo>+</mo><mi>H</mi><mi>sin</mi><mrow><mo>(</mo><msub><mi>&eta;</mi><mi>p</mi></msub><mi>z</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><msub><mi>&eta;</mi><mi>p</mi></msub><mrow><mi>D</mi><mrow><mo>(</mo><msubsup><mi>&beta;</mi><mi>m</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&eta;</mi><mi>p</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>(t>t<sub>0</sub>)<maths num="0004"><![CDATA[<math><mrow><mo>&CenterDot;</mo><mfrac><msup><mi>a</mi><mn>2</mn></msup><mn>2</mn></mfrac><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msubsup><mi>&beta;</mi><mi>m</mi><mn>2</mn></msubsup><msup><mi>a</mi><mn>2</mn></msup></mrow><mn>8</mn></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>D</mi><mrow><mo>(</mo><msubsup><mi>&beta;</mi><mi>m</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&eta;</mi><mi>p</mi><mn>2</mn></msubsup><mo>)</mo></mrow><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>D</mi><mrow><mo>(</mo><msubsup><mi>&beta;</mi><mi>m</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&eta;</mi><mi>p</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow></math>]]></maths>(3)其中:β<sub>m</sub>、η<sub>p</sub>分别为满足方程HJ<sub>o</sub>(β<sub>m</sub>b)=β<sub>m</sub>J<sub>1</sub>(β<sub>m</sub>b)、<img file="FSA00000337557200015.GIF" wi="406" he="133" />的正根,m和p分别表示满足上述方程的第m和第p个正根,m、p=1、2、3...;N(β<sub>m</sub>)、N(η<sub>p</sub>)分别满足:<maths num="0005"><![CDATA[<math><mrow><mfrac><mn>1</mn><mrow><mi>N</mi><mrow><mo>(</mo><msub><mi>&beta;</mi><mi>m</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mn>2</mn><mrow><msubsup><mi>J</mi><mn>0</mn><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>&beta;</mi><mi>m</mi></msub><mi>b</mi><mo>)</mo></mrow></mrow></mfrac><mo>&CenterDot;</mo><mfrac><msubsup><mi>&beta;</mi><mi>m</mi><mn>2</mn></msubsup><mrow><msup><mi>b</mi><mn>2</mn></msup><mrow><mo>(</mo><msup><mi>H</mi><mn>2</mn></msup><mo>+</mo><msubsup><mi>&beta;</mi><mi>m</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mrow><mo>(</mo><msub><mi>&eta;</mi><mi>p</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mn>2</mn><mrow><mrow><mo>(</mo><mi>d</mi><mo>+</mo><mfrac><mi>H</mi><mrow><msubsup><mi>&eta;</mi><mi>p</mi><mn>2</mn></msubsup><mo>+</mo><msup><mi>H</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>&eta;</mi><mi>p</mi><mn>2</mn></msubsup><mo>+</mo><msup><mi>H</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>+</mo><mi>H</mi></mrow></mfrac></mrow></math>]]></maths>式中t<sub>d</sub>和t<sub>0</sub>分别为加热激光束照射光学元件开始时间点和结束时间点,H=h/K<sub>th</sub>,h为待测光学元件与周围环境的热交换系数。D、K<sub>th</sub>分别为待测光学元件的热扩散率和热导率;A<sub>0</sub>为待测光学元件在入射波长处的吸收损耗值,P为入射激光功率r,z分别表示光学元件上的径向和深度方向的位置坐标;t表示时间;所述的最佳温度探测位置的确定方法为:A.对于已知热导率K<sub>th</sub>的待测光学元件,通过不同的待测光学元件半径b和厚度d,在不同径向位置r处温升的精确温度模型数值模拟,用均匀温度模型进行拟合时吸收损耗值最接近假设值的位置即最佳温度探测位置;B.对于未知热导率K<sub>th</sub>的待测光学元件,通过对热导率K<sub>th</sub>在从0.2到50W/m·K范围内变化时、不同待测光学元件的半径b和厚度d,在不同r位置处温升的精确温度模型数值模拟,用均匀温度模型进行拟合时,吸收损耗拟合值均能接近假设值的温度探测位置即该尺寸光学元件的最佳温度探测位置;(3)将高灵敏度温度探测元件调至最佳温度探测位置处,根据国际标准ISO11551所规定的测量过程测量待测光学元件在激光照射前、照射过程中以及照射后的温度变化数据;(4)再利用国际标准ISO11551中所描述的均匀温度模型拟合所测温度数据得到待测光学元件的吸收损耗绝对值。
地址 610209 四川省成都市双流350信箱