输入超高斯指数,光分布沿x轴是通常的高斯分布,而y轴呈平顶分布。从左至右图形拖动窗口下方的滑动条,能看到由于吸收而导致的吸收泵浦能量强度的减少。单击图形可以得到最大高度的图形。
2.3. 定义棒的冷却
选定“Boundaries”标签打开如图5所示的窗口。
在对话框中可以对棒的每个表面分别定义冷却条件,可以选择用固体冷却也可以选择用液体。对于液体冷却需要选定 “Fluid Cooling”选项。
对于固体冷却情况,其保持常温的表面温度可以在“Temperature, K”对话框中定义。而对于液体冷却情况,定义液体的体积温度,另外,在底部的膜系数(它描述与界面处的固体或液体的热传导)对话框必须定义。更详细的描述请参阅LASCAD手册6.10.3节。
对话框“Reference temperature”可以用来计算畸变和相应的晶体在受热前的温度。当边界温度是用开氏温度定义时,注意输入正确的温度值。
图6
2.4. 定义工作物质参数
选定“Material Param”标签打开如图6所示的窗口。
“absorption coefficient”描述由于泵浦光的吸收,泵浦光按 呈指数衰减,它也与掺杂浓度有关系。各向异性的晶体对泵浦光的吸收还与泵浦光的偏振态有关,这种情况下的吸收因子可以定义为两个独立的偏振方向。更详细的描述请参阅LASCAD 手册的附录部分。对Nd:YAG晶体,只需要输入相等的值到相应的对话框。
另外,更复杂的内容(比如,温度对物质参数的依赖,由两种不同物质组成的固体物质的情况)请参阅使用手册。
还可以将物质参数保存为一个文件便于以后在新项目中使用这些参数。
2.5. 定义合成物质
选择“Doping &Mats”标签打开如图7所示的窗口。
在现在这个例子中,我们用输入框定义端盖未掺杂的棒的泵浦区域范围[mm]from z=1.0 to z=5.0意思是说,泵浦光只有在1 ≤ z ≤ 5 mm这个掺杂区域才被吸收。由于我们例子里使用的晶体是6mm长,我们定义两非掺杂端盖为1mm长。
图7
这个窗口显示的参数与在“Material Param”窗口的选择有关。
2.6. 定义控制FEA计算程序的选项
选定“FEA Options”标签打开如图8所示的窗口。
有限元代码(FEA code)使用晶体中一个规则长方形网格来计算,这个长方形与它所相邻的那些小的不规则的网格有关。
“Mesh size in x- and y-direction”是指与晶轴垂直的体积元的边界长。
“Mesh size in z-direction”指与晶轴平行的体积元的边长。
单击“Estimated number of elements”按钮可以得到进行计算所需的随机存储内存(RAM)的大小。可以按照下面的规则计算所需RAM:如果有n个元素,则至少需要2n/1000MB RAM,否则计算机会因为内存与硬盘的数据交换速度太慢而使得计算速度很慢。(你可以这样在任务管理窗口来确定当前空RAM的大小:同时按下WinNt/2K/XP操作系统的Ctrl-Alt-Del 三键,从弹出的窗口中选定任务管理按钮,打开“System Performance”窗口,在物理内存窗口(Physical Memory)显示了以K计算的当前RAM。
对话框“Convergence limits”用来控制迭代计算程序的收敛性。对热分析情况,默认值为1.0E-7时如果最高温度的前7位没有变化,代码将停止运行。对结构分析的限制是最大节点位移的绝对温度值。如果最大迭代次数超过“Maximum number of iterations”对话窗口中设置的次数,不论其收敛性如何,迭代程序将停止。在对话框“Directory for output of FEA results”中可以定义生成的FEA代码文件的保存路径,默认的路径是工作路径的子路径LascadFEA。
在对话框“Position of cutting plane perpendicular to z-axis”中可以定义放置一个垂直于z轴的切平面并以3D可视化方式显示晶体中物理量的分布,下一节会具体描述3D可视化的问题。如果输入的值为0,切平面将沿着棒的轴线和垂直于y轴方向。如果输入的值大于0但小于棒长,切平面垂直于棒的轴线。如果输入值大于棒的长度,整个晶体都会显示出来。
按下“Apply and run FEA”按钮可以将输入的数据传送到内部变量后运行FEA代码。弹出的“Finite Element Analysis”窗口显示当前迭代计算的次数。此外,它还在热分析时显示最高温度,在结构分析时显示最大节点位移的值。计算完成后将出现“FEA finished successfully”消息框,单击“ok”按钮关闭对话框。
需要注意的是,对于元素个数较多的情况,FEA的初始化和网格的生成需要一段时间。
图8
图9
2.7. FEA结果的可视化
LASCAD提供两种方式来可视化FEA结果。
2.7.1. 3D观察器
点击LASCAD主菜单选择“FEA/3D Visualizer”命令来显示FEA结果的3D图形,这是一个基于OpenGL的易于操控的三维观察器。图9显示了一个温度分布的3D图形。由于在2.5节中定义了端面非掺杂,因此温度上升限制在棒的内区域。
2.7.2. 2D数据分布和抛物线拟合
选定主菜单中的“FEA/2D Data Profiles and Parabolic Fit”命令,可以沿z轴显示不同位置处的数据的2D分布。
一个路径对话框弹出来,可以选择包含FEA结果数据的路径。保持默认设置,单击“Open Fit Window”按钮,如图10所示的“2D Data Profiles and Parabolic Fit”窗口将打开。这窗口显示FEA结果数据的横向分布。使用窗口右上角的下拉框,选择沿着晶体轴的位置来显示横向温度分布,与相应位置对应的网格子域被FEA生成。“FEA Data”对话框中的“Show Item”选项可以用来显示其他量的分布,比如泵浦光的分布。
2.8. 高斯模的计算
为了对FEA结果使用高斯模算法,热致折射率分布和热致晶体端面畸变分布必须沿光轴横向按抛物线型拟合。
首先,如果模的预期光斑大小比晶体直径小得多,对话框“Range of Fit”就必须定义。现在这个例子我们就使用默认设置。
单击“2D Data Profiles and Parabolic Fit”窗口中的“Compute Fit”命令按钮,显示如图11所示的抛物线拟合曲线。
红线和蓝线分别代表抛物曲线和FEA结果。拟合是对沿晶轴的所有由FEA网格程序生成的小单元得到的。也就是说,晶体被分为一系列的自聚焦(GRIN)透镜,每个透镜的折射率都成抛物线型分布。最左端小区域的左端面与晶体的左端面一致。同样,最右边的小区域的右端面与晶体的右端面一致。拟合端面的象散球面曲率半径时考虑到了这些端面畸变。
图10
再次,可以使用窗口右上角的下拉框显示沿z轴的拟合位置,拟合区域的抛物线因子在“Parab.Coeff”对话框显示。另外,这些因子被写到“FIT.dat”文件中并保存在FEA子目录下。
2.9. 在模式图中插入晶体
按下ALT键,单击两个元件的中间插入晶体,这时一个代表热透镜晶体的黄颜色的元件出现在模式图中。元件1和2所在的面分别变成了晶体的左右端面,两者的距离调整为晶体长度。为了计算模形状,所有FEA的小区域和晶体的畸变端面的ABCD矩阵通过计算抛物线因子而得到,并且与腔的端面镜的矩阵组合。所有的相邻矩阵相乘起来计算全程的ABCD矩阵,最后得出如图12所示的基模高斯形状。
图11
图13
|
|