[互助]仿真F117的RCS时网格出错，寻求大侠解决

羡慕，我的机子比拟的内存少呀

引用第26楼robben1234于2008-11-28 10:24发表的  :
对于这个模型在Feko里面仿真，使用MLFMM方法。 4gyC?#Ed
e  
总共在两台机子上作了计算： XP(32位），Feko(32位） _SVIY@K|/  
.......

O$ p  
\L6kCY  
k 5~#_D>  
增加内存啊，现在内存很便宜，只要你的主板能插几根内存，都给他插上，另外打听一下有没有外置内存插槽 "#<P--E9  
\&tv *  
用MLFMM算，很耗内存，算个圆球还可以，算飞机就不要奢望了

引用第30楼jimliu12于2008-11-28 17:04发表的  : ;Q8LA",5d  
羡慕，我的机子比拟的内存少呀

#*)X+*  
%D $+Z(  
am(jmf::  
呵呵，我估计你用他那个机器算300MHz的F117，也不一定能算，速度不会快到那去，你可以请他算算，看看速度和内存消耗有多少 y{O817 \  
用工作站应该是王道！ nAY'1!Oi  
    

引用第32楼jetwolf于2008-11-28 22:04发表的 :
呵呵，我估计你用他那个机器算300MHz的F117，也不一定能算，速度不会快到那去，你可以请他算算，看看速度和内存消耗有多少 j@u]( nf  
用工作站应该是王道！ E*AI}:or;  
.......

hZ`<ID  

U-U"RC>  

图片:RCS_300MHz.jpg

38E %]*5F  
7z;X@+O}s  
300MHz  ... MOM+MLFMM(CFIE) !*e1F9k  
Plane Wave 正向打入， 在 XY 平面上每15度取1點的時機器速度跟記憶體需求。 Rl{e<>O\^  
ISHzlEY  
          Scaling of the geometry with factor SKAL =  1.00000E-03 yQ!I`T>a  
          Surface of all triangles in m*m:     242.025892 cXo
d43  
5RsO^2V:  
                           FAST MULTIPOLE METHOD (FMM) 9T#${NK  
 Multilevel FMM is used 7SJR_G6,{  
 Storage of some elements with single precision (saves memory) U[EZ,7n8  
CC B'  
                                DATA FOR MEMORY USAGE SxI-pH'  
 Number of metallic triangles:        52614                   max. triangles:MAXNDR   =    52614 38#Zlcf  
 Number of dielectric triangles:          0 Y?v{V>
;*A  
 Number of diel. GO triangles:            0 |R&cQKaQ`  
 Number of FEM surface triangles:         0 572{DC&T  
 Number of metallic segments:             0                   max. segments: MAXNSEG  =        0 /Os6i&;  
 Number dielectr./magnet. cuboids:        0                   max. cuboids:  MAXNQUA  =        0 _)kTlX:,  
 Number of tetrahedral elements:          0                   max. tetrah.:  MAXNTETRA=        0 H):(8/>(  
 Number of edges in PO region:            0                   max. edges:    MAXPOKA  =        0 l0w<NZF  
 Number of wedges in PO region:           0                   max. wedges:   MAXPOKL  =        0 WCD)yTg:ES  
 Number of Fock regions:                  0                   max. Fock reg.:MAXFOGE  =        0 G _-JR  
 Number of polygonal surfaces:            0                   max. surfaces: MAXPOLYF =        0 XY^]nm-{I  
                                                              max. corner p.:MAXPOLYP =        0 aY
-7K._</  
 Number of UTD cylinders:                 0 6o d^+>U  
 Number of metallic edges (MoM):        78921  unknown:    78921   max. edges:   MAXNKA    =    78921 ~JQ6V?fucD  
                                            0  unknown:        0   (magnet.) Mz6PH)e;  
 Number of metallic edges (PO):             0  unknown:        0   (electr.)         0 (magnet.) `Kbf]"4q  
 Number of dielectric edges (MoM):          0  unknown:        0   (electr.) 8+@j %l j  
                                            0                  0   (magnet.) Ji4JP0  
 Number of dielectric edges (PO):           0  unknown:        0   (electr.)         0 (magnet.) 8I[=iU7]l  
 Number of edges FEM/MoM surface:           0  unknown:        0   (electr.) Ef$a&*)PH  
                                            0                  0   (magnet.) Qo~|[]GE  
 Number of nodes between segments:          0  unknown:        0   max. nodes:   MAXNKNO   =        0 +Ly@5y"  
 Number of connection points:               0  unknown:        0   max. conn.:   MAXNV     =        0 19b@QgfWpb  
 Number of dielectric cuboids:              0  unknown:        0   max. cuboids: MAXNQUA   =        0 Nsn~mY%  
 Number of magnetic cuboids:                0  unknown:        0 *3\*GatJ  
 Number of basis funct. for MoM:        78921  unknown:    78921   max. basisf.  MAXNZEILE =    78921 jS ?#c+9  
 Number of basis funct. for PO:             0  unknown:        0   max. basisf.  MAXNKAPO  =        0 ShesJj  
wZfY~  
Z`<5SHQd  
                    SUMMARY OF REQUIRED TIMES IN SECONDS %uw7sG
z\  
                                            CPU-time      runtime _9C,N2a{C  
 Reading and constructing the geometry         1.891        1.890 \q@Co42n\  
 Checking the geometry                         1.578        1.578 $p|Im,  
 Initialisation of the Greens function         0.000        0.000 l;d4Le  
 Calcul. of coupling for PO/Fock               0.000        0.000 7!dj&?  
 Calcul. of the FMM transfer function          0.657        0.656 }jg1..)"<  
 Fourier transform of FMM basis funct.         1.031        1.031 :cvT/xhO  
 Calcul. of matrix elements                   60.156       60.157 <?>tjCg'  
 Calcul. of right-hand side vector             0.047        0.047 \`;1[m  
 Preconditioning system of linear eqns.       92.031       92.031 Ggry,3X3  
 Solution of the system of linear eqns.       10.625       10.625 I(H9-!&  
 Determination of surface currents             0.000        0.000 goe%'k,  
 Calcul. of impedances/powers/losses           0.000        0.000 {l"(EeW6)  
 Calcul. of averaged SAR values                0.000        0.000 -}u1ZEND  
 Calcul. of power ideal receiving ant.         0.000        0.000 #>M^BOR8  
 Calcul. of cable coupling                     0.000        0.000 E#R1
 
 Calcul. of electric near field                0.000        0.000 2m^qXE$  
 Calcul. of magnetic near field                0.000        0.000 eLIZ<zzW0}  
 Calcul. of far field                          0.234        0.235 ik#ti=.  
                            other              0.281        0.281 ,N1pww?  
                                         -----------  ----------- :nOI|\rC  
                         total times:        168.531      168.531 lVCnu>8  
               (total times in hours:          0.047        0.047) y* :C~  
 Specified CPU-times are referring to the master process only l >~Rzw  
 Sum of the CPU-times of all processes:     1348.249 seconds (    0.375 hours) IIN,Da;hD  
 On average per process:                     168.531 seconds (    0.047 hours) 21O@yNpS$  
 Peak memory usage during the whole solution:  145.734 MByte WjguM  
 (refers to the master process only) D"gv:RojD  
 Sum of the peak memory of all processes:        1.182 GByte [y~kF?a  
 On average per process:                       151.294 MByte I? ="Er[g}  
Finished: 2008-11-29 12:20:34                    

300MHz  ... MOM+MLFMM(CFIE) K`j:F>b  
Plane Wave 正向打入， 在 XY 平面上每5度取1點的時機器速度跟記憶體需求。 #3{{[i(;i  
vT @25  

图片:RCS_300MHz_POSTFEKO.jpg

)/Gi-::  
#_|O93HN'  
                    SUMMARY OF REQUIRED TIMES IN SECONDS CJDNS21m  
]g }5p4*&  
                                            CPU-time      runtime sl(go^  
Reading and constructing the geometry         1.907        1.906 nywC]T  
Checking the geometry                         1.625        1.625 j,SZJ{ebXg  
Initialisation of the Greens function         0.000        0.000 us8HXvvp{  
Calcul. of coupling for PO/Fock               0.000        0.000 V6h8+|hK  
Calcul. of the FMM transfer function          0.641        0.641 ];5Auh0o  
Fourier transform of FMM basis funct.         1.031        1.031 AX'-}5T=  
Calcul. of matrix elements                   60.156       60.156 ^t*Ba>A  
Calcul. of right-hand side vector             0.032        0.031 y&eU\>M  
Preconditioning system of linear eqns.      104.515      104.516 Oqt{ uTI~  
Solution of the system of linear eqns.       10.313       10.312 G'';VoW=  
Determination of surface currents             0.000        0.000 .,U4 ATO  
Calcul. of impedances/powers/losses           0.000        0.000 BjfTt:kY  
Calcul. of averaged SAR values                0.000        0.000 Rlm28  
Calcul. of power ideal receiving ant.         0.000        0.000 _ysakn  
Calcul. of cable coupling                     0.000        0.000 LU8:]zOY  
Calcul. of electric near field                0.000        0.000 ^QG<_Dm]  
Calcul. of magnetic near field                0.000        0.000 fU\;\  
Calcul. of far field                          0.140        0.141 3xKgj5M  
                            other              0.312        0.313 A0,e3gb  
                                         -----------  ----------- P2 qC[1hYH  
                         total times:        180.672      180.672 cct/mX2&~  
               (total times in hours:          0.050        0.050) 86!$<!I  
YXh!+}  
Specified CPU-times are referring to the master process only 8J
- ?bo  
Sum of the CPU-times of all processes:     1445.374 seconds (    0.401 hours) ]0SqLe  
On average per process:                     180.672 seconds (    0.050 hours) B2KBJ4rI[1  
O(9*
VoD  
Peak memory usage during the whole solution:  145.734 MByte yq<W+b/  
(refers to the master process only) (_+ux1h6^  
Sum of the peak memory of all processes:        1.182 GByte JoZzX{eu"  
On average per process:                       151.294 MByte l8 $.k5X  
:XoR~syT  
Finished: 2008-11-29 12:59:46                    

谢谢webhugo网友的计算。 <Lt"e8Z>x  
对于你的计算结果有两个疑问： _\= /~>Xl  
1）"Sum of the peak memory of all processes:        1.182 GByte ” Ol>/^3a=  
    这个是不是说明对内存需求没有想像中那么大？ Ir*{IV
vej  
2）从计算结果看，F117的RCS并没有达到公布数据中0.1m2，那么小的量级 -%K}~4J  
    而且RCS最大点位置也没有在正测向，而是在尾喷口（后向）。 'WBhW5@  
    这一趋势跟我之前用PO算J10模型的结果吻合。 :S!!J*0  
大家是这么看这个问题呢？ hstGe>f[6  
  

引用第35楼robben1234于2008-11-29 13:35发表的 : m@",Zr`f=  
谢谢webhugo网友的计算。 {9cjitl  
对于你的计算结果有两个疑问： X]=8Oa  
1）"Sum of the peak memory of all processes:        1.182 GByte ” 7s Gf_`Z  
    这个是不是说明对内存需求没有想像中那么大？ WM ]eb, 8q  
2）从计算结果看，F117的RCS并没有达到公布数据中0.1m2，那么小的量级 tCv}+7)  
.......

5!Ovd O}g  
N%)q.'M  
1) YES, 在 300MHz 真的是不用很大的內存。IF, 你使用 MOM+MLFMM 的方法來做。 <driD'=F  
是很節省記憶體的..不然就不會有人說他是電大尺寸的終結者。 Q9&H/]"v  
I2,AT+O<  
我記得 CST 也有 MLFMM 的方法..但是跟 FEKO 比起來真的是遜很多。一樣的問題 FEKO 很快就算完了.. CST 的計算時間是 FEKO 的20倍以上。 '*<I<? z;  
個人還是認為要評估軟件還是要深入的去用，不要相信廣告，很多人都是看了廣告後就覺得可行.但是沒仔細廣告後面隱藏的東西。 a&7uRR26  
4v9d& m!<  
2) CFIE 的部分，僅適合封閉區面，是可以再加速計算速度。 gOm8 O,  
PS.. 由於我自己還不是很擅長修模型.... 不知道有誰能幫我將較詳細的圖做好，進氣口就留進氣口， O-3aU!L  
只要進氣口跟結構能摳出封閉的區域即可。我想在解這類的問題對 FEKO 不事件難事。 zu 7Fq]zD  
接下來只剩下 CPU Time 而已。 OMBH[_  
Example z:&/O&?  
----------- /E:BEm!  
   AIR    ] cNye@}$lu  
             ] :BB=E'293  
----------- Q!r` G  
((=T E  
3) 個人認為只要能先確定可以算..接下來只要合理的假設，應該都可以做出您要的答案。 yA';~V\V{>  
ryb81.|  
4) 能否很明確的告知，您想從那個方向看..想看什麼..問題就是這麼簡單。

引用第36楼webhugo于2008-11-29 15:29发表的 :
1) YES, 在 300MHz 真的是不用很大的內存。IF, 你使用 MOM+MLFMM 的方法來做。 ",J&UTUh  
是很節省記憶體的..不然就不會有人說他是電大尺寸的終結者。
.......

N)AlQ'Lwx  

!H[01  
回答webhugo的问题： 6H5o/)Q~  
为什么我一直认为这个结果不可信呢？主要是和《雷达目标特性》（黄培康，殷红成，许小剑编著，电子工业出版社，2005） PP68，给的数据不符。 `JpFqZ'58  
文中叙述如下： dr+(C[=  
图中横坐标方位角0度表示飞机鼻锥方向对准雷达来波，纵坐标RCS的单位为dB.m2，即相对于1m2的分贝数。 
C"B'Dj  
给出2个结论： XWQ
`]m)  
1）频率越高，RCS对姿态角越为敏感，经计算，f>10GHz时自相关姿态角小于1度；f>1GHz时，计算雷达检测频率时可以认为飞机属快起伏目标；在200MHz－800MHz范围内，典型飞机处于目标谐振区高端，具有谐振现象，且起伏缓慢。 {"%a-*@%  
2）RCS随方位角有规律地变化，通常最小值在5度－20度范围内（非0度时最小）；方位角处于飞机后掠翼前缘法线方向时RCS曲线呈现突增，称成为“突增线”（例如图3.1在方位58度），可用于对飞机翼后掠角的识别；通常在正侧向（方位角90度）RCS最大，可达几百平方米，因此一个飞机目标的起伏动态范围可达40dB。 LExm#T`  
E2YVl%.
 
0I<L<^s3^U  
FEKO中的仿真设置我认为是正确的。数据不符的原因可能会是没有对进气口和尾喷口等（影响RCS）关键部位进行更好的模拟。 ('UTjV  
但是同一个网格模型情况下怎么对不同部位进行不同仿真方法设置呢？这是一个问题

上面我是假設 1個平面波射入在每一個方向看到的值 /tj]^QspS  
T[1iZ  
進氣口和尾噴口等... 真的是個大問題...But..我對 3D CAD熟悉度不夠..可能還是要麻煩大家的幫忙.. V.*M;T\i  
D|)_c1g  
但是同一個網格模型情況下怎麽對不同部位進行不同仿真方法設置呢？ p:y\{k"  
<--- 我不是很清楚您的意思.. But, 假使模型都一樣.. plane wave 也一樣，可以藉由儲存 STR檔(PS Card)，只對 RCS 的部分做後處理計算即可。 [.xY>\e
 
T\.(e*hC  
jGz~}&B  
IF 平面波射入.. 自己看到自己的值喔 .. f117_monostatic_rcs_300MHz 99ZWB  
總共計算了361個點 ... 浪費很多時間..(Calcul. of far field ) ;xz_H$g  
TY %zw6 #p  

描述:f117_monostatic_rcs_CADFEKO

图片:f117_monostatic_rcs_CADFEKO.jpg

lkTA"8d  
t*H2
;|zn_  

描述:F117Jetwolf_Fix_f117-100f117_monostatic_rcs_300MHz

图片:f117_monostatic_rcs_300MHz.jpg

;6pB7N  
^`id/  
                    SUMMARY OF REQUIRED TIMES IN SECONDS I]W7FZ=o  
;k`51=Wi  
                                            CPU-time      runtime o^X3YaS)  
Reading and constructing the geometry         1.844        1.843 6}I X{nQI  
Checking the geometry                         1.578        1.579 TBPu&+3  
Initialisation of the Greens function         0.000        0.000 vk
hPE(f  
Calcul. of coupling for PO/Fock               0.000        0.000 v5a\}S<(  
Calcul. of the FMM transfer function          0.640        0.641 f<wYJGI  
Fourier transform of FMM basis funct.         1.031        1.031 Cfmd*,  
Calcul. of matrix elements                   60.140       60.140 J(d+EjC  
Calcul. of right-hand side vector            12.873       12.865 6}N`YOJ.  
Preconditioning system of linear eqns.       97.407       97.406 I@Hx LEGj  
Solution of the system of linear eqns.     3545.168     3545.183 ( `bb1g
z  
Determination of surface currents             0.000        0.000 >JE+g[$@  
Calcul. of impedances/powers/losses           0.263        0.265 Sxc)~y  
Calcul. of averaged SAR values                0.000        0.000 2?q>yL!Gz  
Calcul. of power ideal receiving ant.         0.000        0.000 T})q/oUqK  
Calcul. of cable coupling                     0.000        0.000 -F`GZ  
Calcul. of electric near field                0.000        0.000 uCB9;+ Hjw  
Calcul. of magnetic near field                0.000        0.000 WJONk_WAc  
Calcul. of far field                         56.835       56.832 {v]L|e%{  
                            other             29.877       29.871 B3&C&o.h  
                                         -----------  ----------- 81y<Uz 6  
                         total times:       3807.656     3807.656 (["kbPma  
               (total times in hours:          1.058        1.058) 4PxP*j  
/mz.HCs  
Specified CPU-times are referring to the master process only +&VY6(Zj+*  
Sum of the CPU-times of all processes:    30461.250 seconds (    8.461 hours) ,NZllnW  
On average per process:                    3807.656 seconds (    1.058 hours) 'a+^= c  
{Dl@/fz  
Peak memory usage during the whole solution:  145.736 MByte fiWN^sTM  
(refers to the master process only) 5)XUT`;'){  
Sum of the peak memory of all processes:        1.182 GByte ynM~&]fk#k  
On average per process:                       151.294 MByte 4y3c=L No  
& tT6.@kH  
Finished: 2008-11-29 18:22:02                    

进气道的问题就是腔体的问题，存在着多次反射，目前还不能很好的而通用的计算，东南崔铁军曾经说过