简介：本文是以十字元件为背景光源，经过一个透镜元件成像在探测器上，并显示其热成像图。
成像示意图
首先我们建立十字元件命名为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 |
