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

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

引用第26楼robben1234于2008-11-28 10:24发表的  :
对于这个模型在Feko里面仿真，使用MLFMM方法。 }5u$/c@f1  
总共在两台机子上作了计算： XP(32位），Feko(32位） +Y6=;*j$  
.......

jq%Qc9y  
0)+F}SyyD  
ljQru ^(u  
增加内存啊，现在内存很便宜，只要你的主板能插几根内存，都给他插上，另外打听一下有没有外置内存插槽 3N[t2Y1r  
yd'>Mw  
用MLFMM算，很耗内存，算个圆球还可以，算飞机就不要奢望了

引用第30楼jimliu12于2008-11-28 17:04发表的  : Bs`mzA54  
羡慕，我的机子比拟的内存少呀

9^u}~e #(  
E~@&&dU8  
0q4PhxR`e  
呵呵，我估计你用他那个机器算300MHz的F117，也不一定能算，速度不会快到那去，你可以请他算算，看看速度和内存消耗有多少 \G#_z|'dN  
用工作站应该是王道！ &S=Qu?H  
    

引用第32楼jetwolf于2008-11-28 22:04发表的 :
呵呵，我估计你用他那个机器算300MHz的F117，也不一定能算，速度不会快到那去，你可以请他算算，看看速度和内存消耗有多少 T?
f{.a)  
用工作站应该是王道！ &+@`Si=  
.......

4U}J?EB?K  
Yb`b/BMR  

图片:RCS_300MHz.jpg

Q:Pp'[ RK  
sRZ<c  
300MHz  ... MOM+MLFMM(CFIE) Y$"
m*0  
Plane Wave 正向打入， 在 XY 平面上每15度取1點的時機器速度跟記憶體需求。 &n_f.oUc  
D&"l
u*"tg  
          Scaling of the geometry with factor SKAL =  1.00000E-03 axt
;}8  
          Surface of all triangles in m*m:     242.025892 wT- <#+L\  
- EGZ  
                           FAST MULTIPOLE METHOD (FMM) Ssz;d&93  
 Multilevel FMM is used "OYD9Q''  
 Storage of some elements with single precision (saves memory) 2]9 2J  
]NBx5m+y@i  
                                DATA FOR MEMORY USAGE ,D*bLXWh  
 Number of metallic triangles:        52614                   max. triangles:MAXNDR   =    52614 #_S]\=N(  
 Number of dielectric triangles:          0 _V1:'T8  
 Number of diel. GO triangles:            0 8@LWg d  
 Number of FEM surface triangles:         0 'dh{q`#0  
 Number of metallic segments:             0                   max. segments: MAXNSEG  =        0 [pp|*@1T  
 Number dielectr./magnet. cuboids:        0                   max. cuboids:  MAXNQUA  =        0 `?R{sNr.  
 Number of tetrahedral elements:          0                   max. tetrah.:  MAXNTETRA=        0 4v_Hh<%  
 Number of edges in PO region:            0                   max. edges:    MAXPOKA  =        0 EATVce]T  
 Number of wedges in PO region:           0                   max. wedges:   MAXPOKL  =        0 ;SjNZi)4d  
 Number of Fock regions:                  0                   max. Fock reg.:MAXFOGE  =        0 f42F@M(:  
 Number of polygonal surfaces:            0                   max. surfaces: MAXPOLYF =        0 P2f^]z  
                                                              max. corner p.:MAXPOLYP =        0
1Dc6v57  
 Number of UTD cylinders:                 0  )sdHJ  
 Number of metallic edges (MoM):        78921  unknown:    78921   max. edges:   MAXNKA    =    78921 ebJTrh<{  
                                            0  unknown:        0   (magnet.) .TdFI"Yn  
 Number of metallic edges (PO):             0  unknown:        0   (electr.)         0 (magnet.) gsEcvkj*  
 Number of dielectric edges (MoM):          0  unknown:        0   (electr.) CUG"2K9  
                                            0                  0   (magnet.) '"\n,3h  
 Number of dielectric edges (PO):           0  unknown:        0   (electr.)         0 (magnet.) @$^4Av-  
 Number of edges FEM/MoM surface:           0  unknown:        0   (electr.) F1_,V?  
                                            0                  0   (magnet.) F.a
G7  
 Number of nodes between segments:          0  unknown:        0   max. nodes:   MAXNKNO   =        0 0AffD:  
 Number of connection points:               0  unknown:        0   max. conn.:   MAXNV     =        0 hZ>m:es  
 Number of dielectric cuboids:              0  unknown:        0   max. cuboids: MAXNQUA   =        0 4FZR }e\  
 Number of magnetic cuboids:                0  unknown:        0 }ri
M-  
 Number of basis funct. for MoM:        78921  unknown:    78921   max. basisf.  MAXNZEILE =    78921 &}Wi@;G]2
 
 Number of basis funct. for PO:             0  unknown:        0   max. basisf.  MAXNKAPO  =        0 ,D 
}Ka?  
]g;+7  
Cj4Y, N  
                    SUMMARY OF REQUIRED TIMES IN SECONDS S .x>w/  
                                            CPU-time      runtime s+fxv(,"c  
 Reading and constructing the geometry         1.891        1.890 |;~2y>E  
 Checking the geometry                         1.578        1.578 s#aj5_G  
 Initialisation of the Greens function         0.000        0.000 6Y\TVRR  
 Calcul. of coupling for PO/Fock               0.000        0.000 )V>OND  
 Calcul. of the FMM transfer function          0.657        0.656 QgYt(/S  
 Fourier transform of FMM basis funct.         1.031        1.031 '{.4~:  
 Calcul. of matrix elements                   60.156       60.157 J|^XD<Y  
 Calcul. of right-hand side vector             0.047        0.047 j2IK\~W?-  
 Preconditioning system of linear eqns.       92.031       92.031 T-STM"~%  
 Solution of the system of linear eqns.       10.625       10.625 i}C%8}%  
 Determination of surface currents             0.000        0.000 k*6"!J%A  
 Calcul. of impedances/powers/losses           0.000        0.000 :w5
g!G?z  
 Calcul. of averaged SAR values                0.000        0.000 z&9MkbH1  
 Calcul. of power ideal receiving ant.         0.000        0.000 ?.F^Oi6 u  
 Calcul. of cable coupling                     0.000        0.000 e~(e&4pb  
 Calcul. of electric near field                0.000        0.000 ;qUB[Kw  
 Calcul. of magnetic near field                0.000        0.000 f1'X<VA  
 Calcul. of far field                          0.234        0.235 !LpjTMYs  
                            other              0.281        0.281 8'zl\:@N  
                                         -----------  ----------- @J5TDq @  
                         total times:        168.531      168.531 ?qSwV.l]d  
               (total times in hours:          0.047        0.047) Yl'8" \HF  
 Specified CPU-times are referring to the master process only r^msJ|k8[  
 Sum of the CPU-times of all processes:     1348.249 seconds (    0.375 hours) .-T P1C  
 On average per process:                     168.531 seconds (    0.047 hours) s:xJ }Ll  
 Peak memory usage during the whole solution:  145.734 MByte r0[<[jEh  
 (refers to the master process only) z a_0-G%C2  
 Sum of the peak memory of all processes:        1.182 GByte h x5M)8#+  
 On average per process:                       151.294 MByte UcB&pt&  
Finished: 2008-11-29 12:20:34                    

300MHz  ... MOM+MLFMM(CFIE) ]:2Ro:4Yv  
Plane Wave 正向打入， 在 XY 平面上每5度取1點的時機器速度跟記憶體需求。 baTd;`Pn  
}Ptv[{q]GE  

图片:RCS_300MHz_POSTFEKO.jpg

[hH>BEtm  
bk=;=K  
                    SUMMARY OF REQUIRED TIMES IN SECONDS Y?K?*`Pkc1  
$'d,X@}8  
                                            CPU-time      runtime jz;{,F  
Reading and constructing the geometry         1.907        1.906 ^jS1g*nrN  
Checking the geometry                         1.625        1.625 q^^R
|X1  
Initialisation of the Greens function         0.000        0.000 S[y_Ewzq  
Calcul. of coupling for PO/Fock               0.000        0.000 .L1[Rv3  
Calcul. of the FMM transfer function          0.641        0.641 FcZ)^RQ4G  
Fourier transform of FMM basis funct.         1.031        1.031 0(u}z  
Calcul. of matrix elements                   60.156       60.156 T1Z*>(M  
Calcul. of right-hand side vector             0.032        0.031 <Pe'&u  
Preconditioning system of linear eqns.      104.515      104.516 WnOYU9;%  
Solution of the system of linear eqns.       10.313       10.312 jZmL7 V  
Determination of surface currents             0.000        0.000 +.^pAz U}R  
Calcul. of impedances/powers/losses           0.000        0.000 #pW!(tfN^a  
Calcul. of averaged SAR values                0.000        0.000 Z]2z*XD  
Calcul. of power ideal receiving ant.         0.000        0.000 $K\e Pfk  
Calcul. of cable coupling                     0.000        0.000 2{l|<'  
Calcul. of electric near field                0.000        0.000 ZR)M<*$  
Calcul. of magnetic near field                0.000        0.000 L">jSZW[[  
Calcul. of far field                          0.140        0.141 1lA? 5:  
                            other              0.312        0.313 =
2rdbq6R  
                                         -----------  ----------- @sZ' --Y  
                         total times:        180.672      180.672 ]xb2W~  
               (total times in hours:          0.050        0.050) HL$7Ou  
ume70ap}m  
Specified CPU-times are referring to the master process only >
Q"3dw  
Sum of the CPU-times of all processes:     1445.374 seconds (    0.401 hours) 7tJ#0to  
On average per process:                     180.672 seconds (    0.050 hours) ?kV_!2U)'K  
O#J7GbrHO  
Peak memory usage during the whole solution:  145.734 MByte XJl 3\*  
(refers to the master process only) \0;w7tdo  
Sum of the peak memory of all processes:        1.182 GByte #7:ah  
On average per process:                       151.294 MByte J%q)6&  
`e7vSp  
Finished: 2008-11-29 12:59:46                    

谢谢webhugo网友的计算。 Ts=TaRwWf  
对于你的计算结果有两个疑问： !@T~m1L eY  
1）"Sum of the peak memory of all processes:        1.182 GByte ” 28}L.>5k  
    这个是不是说明对内存需求没有想像中那么大？ ose)\rM'  
2）从计算结果看，F117的RCS并没有达到公布数据中0.1m2，那么小的量级 l[WX77bp=  
    而且RCS最大点位置也没有在正测向，而是在尾喷口（后向）。 2`FDY3n  
    这一趋势跟我之前用PO算J10模型的结果吻合。 +^rt48${ y  
大家是这么看这个问题呢？ XA69t2J~F  
  

引用第35楼robben1234于2008-11-29 13:35发表的 : $&. rS.*  
谢谢webhugo网友的计算。 c^}DBvG,  
对于你的计算结果有两个疑问： '2ACZcjDSv  
1）"Sum of the peak memory of all processes:        1.182 GByte ” j.SE'a_  
    这个是不是说明对内存需求没有想像中那么大？ :2AlvjvjZ  
2）从计算结果看，F117的RCS并没有达到公布数据中0.1m2，那么小的量级 Md>9Daa~  
.......

\&^U9=uq  
!~ZAm3GwL  
1) YES, 在 300MHz 真的是不用很大的內存。IF, 你使用 MOM+MLFMM 的方法來做。 G5!!^p~  
是很節省記憶體的..不然就不會有人說他是電大尺寸的終結者。 u~1 ,88&U  
~Da-|FKa>  
我記得 CST 也有 MLFMM 的方法..但是跟 FEKO 比起來真的是遜很多。一樣的問題 FEKO 很快就算完了.. CST 的計算時間是 FEKO 的20倍以上。 U/>l>J
5  
個人還是認為要評估軟件還是要深入的去用，不要相信廣告，很多人都是看了廣告後就覺得可行.但是沒仔細廣告後面隱藏的東西。 UkM#uKr:  
-}"nb-RR\  
2) CFIE 的部分，僅適合封閉區面，是可以再加速計算速度。 }n?D#Pk,  
PS.. 由於我自己還不是很擅長修模型.... 不知道有誰能幫我將較詳細的圖做好，進氣口就留進氣口， RtF!(gd  
只要進氣口跟結構能摳出封閉的區域即可。我想在解這類的問題對 FEKO 不事件難事。 Ap!i-E,"J  
接下來只剩下 CPU Time 而已。 nF=[m; ~  
Example PS`F  
----------- )S|&3
\  
   AIR    ] D]+@pKb  
             ] 5_4=(?<  
----------- kQO5sX$;  
zKNac[:  
3) 個人認為只要能先確定可以算..接下來只要合理的假設，應該都可以做出您要的答案。 ~5n?=  
Q`z2SYz>  
4) 能否很明確的告知，您想從那個方向看..想看什麼..問題就是這麼簡單。

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

qnf\K
}  
^tl&FWF  
回答webhugo的问题： Sigu p#.p  
为什么我一直认为这个结果不可信呢？主要是和《雷达目标特性》（黄培康，殷红成，许小剑编著，电子工业出版社，2005） PP68，给的数据不符。 p&3~n: Fo  
文中叙述如下： ph@2[rUp  
图中横坐标方位角0度表示飞机鼻锥方向对准雷达来波，纵坐标RCS的单位为dB.m2，即相对于1m2的分贝数。 8iN@n8O  
给出2个结论： :c0 |w  
1）频率越高，RCS对姿态角越为敏感，经计算，f>10GHz时自相关姿态角小于1度；f>1GHz时，计算雷达检测频率时可以认为飞机属快起伏目标；在200MHz－800MHz范围内，典型飞机处于目标谐振区高端，具有谐振现象，且起伏缓慢。 QjyJmW("Z  
2）RCS随方位角有规律地变化，通常最小值在5度－20度范围内（非0度时最小）；方位角处于飞机后掠翼前缘法线方向时RCS曲线呈现突增，称成为“突增线”（例如图3.1在方位58度），可用于对飞机翼后掠角的识别；通常在正侧向（方位角90度）RCS最大，可达几百平方米，因此一个飞机目标的起伏动态范围可达40dB。 #oEq)Vq>g|  
#L xfE<^  
h
Pr
E  
FEKO中的仿真设置我认为是正确的。数据不符的原因可能会是没有对进气口和尾喷口等（影响RCS）关键部位进行更好的模拟。 anFl:=  
但是同一个网格模型情况下怎么对不同部位进行不同仿真方法设置呢？这是一个问题

上面我是假設 1個平面波射入在每一個方向看到的值 Fv6<Cz6L  
ndIU0kq3  
進氣口和尾噴口等... 真的是個大問題...But..我對 3D CAD熟悉度不夠..可能還是要麻煩大家的幫忙.. &% \`Lwh  
W+0VrH 0F  
但是同一個網格模型情況下怎麽對不同部位進行不同仿真方法設置呢？ yy6?16@  
<--- 我不是很清楚您的意思.. But, 假使模型都一樣.. plane wave 也一樣，可以藉由儲存 STR檔(PS Card)，只對 RCS 的部分做後處理計算即可。 U3yIONlt  
icPg<>
TQ  
U:.  
IF 平面波射入.. 自己看到自己的值喔 .. f117_monostatic_rcs_300MHz e.|_=Gd2/  
總共計算了361個點 ... 浪費很多時間..(Calcul. of far field ) *!*J5/b  
r%/*,lLO  

描述:f117_monostatic_rcs_CADFEKO

图片:f117_monostatic_rcs_CADFEKO.jpg

s~(iB{-  
|VL(#U  

描述:F117Jetwolf_Fix_f117-100f117_monostatic_rcs_300MHz

图片:f117_monostatic_rcs_300MHz.jpg

V5{^R+_)Ya  
G:2m)0bW  
                    SUMMARY OF REQUIRED TIMES IN SECONDS ^#4?v^QNh  
2nRL;[L*.  
                                            CPU-time      runtime F
hn=}7|4q  
Reading and constructing the geometry         1.844        1.843 \VX~'pkrd/  
Checking the geometry                         1.578        1.579 {~"6/L  
Initialisation of the Greens function         0.000        0.000 Vuqm{bo^  
Calcul. of coupling for PO/Fock               0.000        0.000 e/;chMCq  
Calcul. of the FMM transfer function          0.640        0.641 S:En9E  
Fourier transform of FMM basis funct.         1.031        1.031 fscAG\>8  
Calcul. of matrix elements                   60.140       60.140 O0gLu1*1v  
Calcul. of right-hand side vector            12.873       12.865 ~e<v<92Xu  
Preconditioning system of linear eqns.       97.407       97.406 6*<=(SQI  
Solution of the system of linear eqns.     3545.168     3545.183 0#rv.rJ{  
Determination of surface currents             0.000        0.000 # N.(ZP  
Calcul. of impedances/powers/losses           0.263        0.265 a!>yX ex  
Calcul. of averaged SAR values                0.000        0.000 I!ykm\<  
Calcul. of power ideal receiving ant.         0.000        0.000 6h|@Bz/A  
Calcul. of cable coupling                     0.000        0.000 r%g?.4o*b  
Calcul. of electric near field                0.000        0.000 +0Rr5^8u  
Calcul. of magnetic near field                0.000        0.000 0/."R;  
Calcul. of far field                         56.835       56.832 N@>,gm@UU  
                            other             29.877       29.871 qi\!<clv  
                                         -----------  ----------- /i]!=~\qFs  
                         total times:       3807.656     3807.656 siHS@S  
               (total times in hours:          1.058        1.058) |6pNe T[  
_ s]=g  
Specified CPU-times are referring to the master process only :_
q  
Sum of the CPU-times of all processes:    30461.250 seconds (    8.461 hours) 'V5^D<1P  
On average per process:                    3807.656 seconds (    1.058 hours) P/'~&*m-  
H3"[zg9L:a  
Peak memory usage during the whole solution:  145.736 MByte E Cx_ [|3{  
(refers to the master process only) 2>3g
C_^go  
Sum of the peak memory of all processes:        1.182 GByte aNpeePF)z  
On average per process:                       151.294 MByte ?5#Ng,8iT  
64^dy V,;  
Finished: 2008-11-29 18:22:02                    

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