十字元件热成像分析

简介：本文是以十字元件为背景光源，经过一个透镜元件成像在探测器上，并显示其热成像图。

成像示意图

首先我们建立十字元件命名为Target

创建方法：

面1 ：
面型：plane
材料：Air
孔径：X=1.5, Y=6,Z=0.075,形状选择Box

辅助数据：
首先在第一行输入temperature :300K，
emissivity:0.1;

面2 ：
面型：plane
材料：Air
孔径：X=1.5, Y=6,Z=0.075,形状选择Box

位置坐标：绕Z轴旋转90度，

辅助数据：

首先在第一行输入temperature :300K，emissivity: 0.1;

Target 元件距离坐标原点-161mm;

单透镜参数设定：F=100, bend=0, 位置位于坐标原点

探测器参数设定：

在菜单栏中选择Create/Element Primitive /plane

元件半径为20mm*20,mm,距离坐标原点200mm。

光源创建：

光源类型选择为任意平面，光源半角设定为15度。

我们将光源设定在探测器位置上，具体的原理解释请见本章第二部分。

我们在位置选项又设定一行的目的是通过脚本自动控制光源在探测器平面不同划分区域内不同位置处追迹光线。

功率数值设定为：P=sin2（theta） theta为光源半角15度。我们为什么要这么设定，在第二部分会给出详细的公式推导。

创建分析面：

到这里元件参数设定完成，现在我们设定元件的光学属性，在前面我们分别对第一和第二面设定的温度和发射系数，散射属性我们设定为黑朗伯，4%的散射。并分别赋予到面一和面二。

到此，所有的光学结构和属性设定完成，通过光线追迹我们可以查看光线是否可以穿过元件。

FRED在探测器上穿过多个像素点迭代来创建热图

FRED具有一个内置的可编译的Basic脚本语言。从Visual Basic脚本语言里，几乎所有用户图形界面（GUI）命令是可用这里的。FRED同样具有自动的客户端和服务器能力，它可以被调用和并调用其他可启动程序，如Excel。因此可以在探测器像素点上定义多个离轴光源，及在FRED Basic脚本语言里的For Next loops语句沿着探测器像素点向上和向下扫描来反向追迹光线,这样可以使用三维图表查看器(Tools/Open plot files in 3D chart)调用和查看数据。
将如下的代码放置在树形文件夹 Embedded Scripts，

打开后清空里面的内容，此脚本为通用脚本适用于一切可热成像的应用。

绿色字体为说明文字，

'#Language "WWB-COM"
'script for calculating thermal image map
'edited rnp 4 november 2005

'declarations
Dim op As T_OPERATION
Dim trm As T_TRIMVOLUME
Dim irrad(32,32) As Double 'make consistent with sampling
Dim temp As Double
Dim emiss As Double
Dim fname As String, fullfilepath As String

'Option Explicit

Sub Main
    'USER INPUTS
    nx = 31
    ny = 31
    numRays = 1000
    minWave = 7    'microns
    maxWave = 11   'microns
    sigma = 5.67e-14 'watts/mm^2/deg k^4
    fname = "teapotimage.dat"

    Print ""
    Print "THERMAL IMAGE CALCULATION"

    detnode = FindFullName( "Geometry.Detector.Surface" ) '找到探测器平面节点

    Print "found detector array at node " & detnode

    srcnode = FindFullName( "Optical Sources.Source 1" ) '找到光源节点

    Print "found differential detector area at node " & srcnode

    GetTrimVolume detnode, trm
    detx = trm.xSemiApe
    dety = trm.ySemiApe
    area = 4 * detx * dety
    Print "detector array semiaperture dimensions are " & detx & " by " & dety
    Print "sampling is " & nx & " by " & ny

    'reset differential detector area dimensions to be consistent with sampling
    pixelx = 2 * detx / nx
    pixely = 2 * dety / ny
    SetSourcePosGridRandom srcnode, pixelx / 2, pixely / 2, numRays, False
    Print "resetting source dimensions to " & pixelx / 2 & " by " & pixely / 2

    'reset the source power
    SetSourcePower( srcnode, Sin(DegToRad(15))^2 )
    Print "resetting the source power to " & GetSourcePower( srcnode ) & " units"

    'zero out irradiance array
    For i = 0 To ny - 1
        For j = 0 To nx - 1
            irrad(i,j) = 0.0
        Next j
    Next i

    'main loop
    EnableTextPrinting( False )

    ypos =  dety + pixely / 2
    For i = 0 To ny - 1
        xpos = -detx - pixelx / 2
        ypos = ypos - pixely

        EnableTextPrinting( True )
        Print i
        EnableTextPrinting( False )

        For j = 0 To nx - 1

            xpos = xpos + pixelx

            'shift source
            LockOperationUpdates srcnode, True
            GetOperation srcnode, 1, op
            op.val1 = xpos
            op.val2 = ypos
            SetOperation srcnode, 1, op
            LockOperationUpdates srcnode, False

'raytrace
            DeleteRays
            CreateSource srcnode
            TraceExisting 'draw

            'radiometry
            For k = 0 To GetEntityCount()-1
                If IsSurface( k ) Then
                    temp = AuxDataGetData( k, "temperature" )
                    emiss = AuxDataGetData( k, "emissivity" )
                    If ( temp <> 0 And emiss <> 0 ) Then
                        ProjSolidAngleByPi = GetSurfIncidentPower( k )
                        frac = BlackBodyFractionalEnergy ( minWave, maxWave, temp )
                        irrad(i,j) = irrad(i,j) + frac * emiss * sigma * temp^4 * ProjSolidAngleByPi
                    End If

                End If

            Next k

        Next j

    Next i
    EnableTextPrinting( True )

    'write out file
    fullfilepath = CurDir() & "\" & fname
    Open fullfilepath For Output As #1
    Print #1, "GRID " & nx & " " & ny
    Print #1, "1e+308"
    Print #1, pixelx & " " & pixely
    Print #1, -detx+pixelx/2 & " " & -dety+pixely/2

    maxRow = nx - 1
    maxCol = ny - 1
    For rowNum = 0 To maxRow                    ' begin loop over rows (constant X)
            row = ""
        For colNum = maxCol To 0 Step -1            ' begin loop over columns (constant Y)
            row = row & irrad(colNum,rowNum) & " "     ' append column data to row string
        Next colNum                     ' end loop over columns

            Print #1, row

    Next rowNum                         ' end loop over rows
    Close #1

    Print "File written: " & fullfilepath
    Print "All done!!"
End Sub

在输出报告中，我们会看到脚本对光源的孔径和功率做了修改，并最终经过31次迭代，将所有的热成像数据以dat的格式放置于：

找到Tools工具，点击Open plot files in 3D chart并找到该文件

打开后，选择二维平面图：

谭健	2023-04-07 08:20
感谢分享 v^zu:Z*

查看本帖完整版本: [-- 十字元件热成像分析 --] [-- top --]