《应用feko计算入门,第一个例子》由会员分享，可在线阅读，更多相关《应用feko计算入门,第一个例子（21页珍藏版）》请在金锄头文库上搜索。
1、由一个详细的例子来了解 FEKO 是如何组织以上模块工作的。例子来源于 FEKO 安装目录下的FEKOexamplesprojec tss tar tup,是计算金属球放置在半波阵 子天线前的电磁特性。如图示。启动CADFEKO。出现如图所示的默认窗口。选择菜单中的FileSave as，存储 工程文件。说明：工程将以相同文件名称，不同文件后缀的形式保存。詁标准工具CADFEKO - Untitled - 3D vie菜单 gilg gmit 卑詐 EieMietiry 也豁h星如 世ip取档 /dpJ & K圜釘f f 氽再毛熬F C 口倉専0痔 测园分i|. VariablescO = 1
2、/iqt (tpiO*muD* tpiO = 8.854187817611 =muO = pi*4e-7pi e0 = sqrt tmuD/epsO) pointsGeometryMtshci宀引回削净诅入巴萍*巴 士同模型树结构3D模型&。卄w二改模工具儿何树结构StartingprQjctB 消息状态建模丁具定义变量。选择菜单中的Geometry Add visible，或者双击模型树区的，定义模型中所需要的8个变量：lam=3、dis=2.5、len=lam/2、rad=1、 seg_len=(lam/15)/1.2、 segrad=0. 002, 、 tri_len=(lam/6.3)
3、/1.2、f_len=0.99*seg_len 。说明：定义了包括波长、剖分精度、模型尺寸等量。定义后的变量可以在模型树区被查看到。变量只能是以字母、数字和下划线组成。创建细线模型。选择菜单中的： Geometrycurvepolyline,或者单击建模工 具中的 ， 创建一段细线模型（各点坐标如图示）。点击 豎 可以增加细线的结点，把label值修改为dipole，标识清楚。分别单击Create 和Close完成细线建模。说明：这里只创建了半波阵子的一半，另一半将在DITFEK0中以电对称方式创 建。结点1到结点2的距离只有f_len=0.99*seg_len小于线段的单元剖分精 度seg_
4、len,这是为保证馈电加载到唯一剖分单元准备的。选择细线的一段。在3D模型区，可以观察到已经完成的细线。单击3D模型区 固（最佳方式显示），让细线更容易被观察。单击选择工具区的（选择线段），在 3D 模型区中点击细线，选择细线的底部一段。说明：在3D模型区，为了方便观察，按住Shift和鼠标左键或滚动鼠标滑轮可 以缩放模型，按住 Ctrl 和鼠标左键可以移动模型，直接按住鼠标左键可以旋转 模型。重命名细线各部分的名称。在几何树区可以观察到Edges下有两个wire （这 里略去Wire后数字编号，因为读者在模仿此工程中得到的数字编号会和图示编 号不一致，但略去数字编号不影响叙述所指的唯一性。后
5、面亦如此，）其中一条 已经被选中，按下F2键或者右键选择rename，重命名这条线段为feed；另一条 线段命名为 wire1。说明：几何树区线和面的名称是EDITFEKO可以调用的。重新命名是为了在 EDITFEKO中可以调用它们。创建金属球。选择菜单中GeometrySolidSphere，或者单击建模工具区中的 Q，设置Radius值为rad，单击 Create 创建金属球面。按下3D模型区冈， 调整到最合适的观察角度。说明：这样创建是个实心金属球，而EKO中的矩量法（MOM）默认是以表面积分方程求解，以面剖分为基础的，需要把实心球体变成空心球壳。设置金属球为空心。确认球被选中（单击模型
6、树区中GeometrySphere）,右 键选择几何树中RegionsRegion的Properties,把内部介质值改为Free space。 说明：实心金属体结构球变成了面结构球面。以尸0分裂球面。确认模型树区的GeometrySphere被选中，单击改模工具区的 ，键入如图示数值，把模型分裂。说明：联合在一起的球面被分开成两部分Ori ginU vectorX Io. 0O0.O0.Y |0?0 z fi?oSgt to pl：=lTlhGlub：J. XYGlobal ZXGlub：al YZW orkfil ：=ltlhOK删除 ySplit_back 部分，右键 单击选删除。说明：
7、模型树区GeometrySplit_fron部分为剩下的y0部分球面，原来的球 面已经不存在。如果想恢复可以单击标准工具中的，也可以右键选择模型树 区GeometrySplit_frontSphere的Copy，把模型复制 出来，再把 GeometrySplit_front部分删除。这也是在以后修改模型时经常用的。以x=0分裂球面。确认模型树区的GeometrySplit_front部分被选中，再次 单击改模工具区的键入如图示数值，把模型分裂。删除xSplit_back部分, 右键单击选删除。di p ul eDele1Delfa.FLerL：dineCopyCtrl+KCojiy sp e
8、i： i ：1AsEemblyProperti es Spl i t_frurLt3 (iSpherel i： (JSpli t_fr ontU 丄；#pl i t_fruntS SpherelMeshes以z=0分裂球面。确认模型树区的GeometrySplit_front部分被选中，再次 单击改模工具区的C),键入如图示数值，把模型分裂。sldldldldldldldldl删除 zSplit_back 部分，右键单击选删除。右键单击 GeometrySplit_front, 重命名为Partial_sphere，标识清楚。固和犬确认模型。按下3D模型区默认角度观察），调整到最合适的观察角度
9、。如下图示。如果没有显示此图，请确认3D模型区的显示/隐藏模型）Create和门晔完成网格剖分。说明：网格剖分完成，在模型树区生成MeshsPartial_sphere和M eshsdipole 两个模型结构的剖分，选中它们，几何树区就有相应的bels。这些labels是 EDITFEKO 可以直接调用的，后面将用到它们。细线被剖分成线段单元，球面被 剖分成三角贴片单元。导入模型到EDITFEKO界面。选择菜单中的RunEDITFEKO，存储已经完成的模 型，进入EDITFEKO界面。说明：这个界面全部是用卡片进一步建模或控制模型参数把，光标移到文本编辑 的某一行，按F1可以编辑此卡片，Esc
10、退出卡片编辑。把光标移到空行处，单 击快速卡片区的相应卡片，可以创建卡片。*行是注释语句。先导入三角形剖分结构。把光标移到“IN 8 31 startup.cfm”行，按下F1 ， 在 快速 卡片 区出现 IN 卡片 编辑界面 ， 调整 参数 如图示 ， 按下确定。说明：默认IN卡导入了全部模型，本例为了运用对称提高计算效率，先把球面剖分（三角贴片单元）导入。模型关于x=0几何对称和y=0磁对称，把1/8球面变成1/2球面。选择菜单 中 Geometry CardModify GeometrySY ，在快速卡片区出现了 SY 卡编辑界面， 调整参数如图示。说明：几何对称会减少计算时间，不减少计
11、算内存；电和磁对称会减少待定系数 的求解个数，既减少计算时间，也减少计算内存。就半波阵子而言，两个臂是电 对称结构，围绕两臂周围是磁对称结构。再导入阵子的上臂。选择菜单中 Geome try CardModify Geome try IN,在快 速卡片区出现了 IN卡编辑界面，调整参数如图示。说明：再把半波阵子的上臂导入模型，它和已经存在的上半球共同组成电对称结 构。细线的剖分单元是segments。阵子上臂和上半球关于 z=0 电对称。 选择菜单中 Geometry CardModifyGeometrySY，在快速卡片区出现了 SY卡编辑界面，调整参数如图示。SY - Specify syi
12、nmetiy of the geoinetiySelect symmetry for the planeElectricalMagneticNoneGeometricalSelect symmetry for the plane y=Dc* NoneC ElectricalC GeometricalC MagneticSelect symmetry for the plane z=0NoneGeometrical( ElectricalMagnetic最后导入馈电部分。选择菜单中 Geome try CardModify Geome try IN,在快 速卡片区出现了 IN卡编辑界面，调整参数
13、如图示。说明：在这个过程中，读者随时可以运行菜单中的RunPREFEKO,再运行RunP0STFEK0查看模型是否建立正确。定义频率扫描的范围。在文本编辑区移动光标到EG卡后的空行中，选择菜单 中 Geometry CardPoints 、 labels 、 Segmentations# ， 编 辑 变 量#f_start=0.8*#c0/#lam ；单击接着编辑变量#f_end=1.2*#c0/#lam。说明：EG卡表明建模结束，其后是控制卡片。注意这里的变量前面都有。设置计算频率。在文本编辑区移动光标到FR卡行，按下F1到FR卡编辑界面。 调整参数如下图示。说明：从#f_start到#f_
14、end共设置了 9个线性变化的频点。设置电压源激励。在文本编辑区移动光标到“* Set source”下的空行，选 择菜单中的Control CardExcitationsAl，编辑Al卡如图示。其中“Apply source to las t segmen t with label 的值为“ dipole.feed说明：把电压源加到dipole.feed上，这里的dipole指的是前面CADFEK0中模 型树区的meshsdipole而feed指的是其下面几何模型区的abelsfeed要求输出表面电流。在文本编辑区移动光标到“ * Sol ution control下的 空行，选择菜单中的Control CardOutput controlOS，编辑OS卡如图示。OS - Output cunemtsExtract currents to output fileNo currents* Ml currentsOnly the currents on trianglesOnly the segment currentsCurrents on structures with single labelCurrents on structures
