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

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

引用第26楼robben1234于2008-11-28 10:24发表的  :
对于这个模型在Feko里面仿真，使用MLFMM方法。 5SZa,+]  
总共在两台机子上作了计算： XP(32位），Feko(32位） sT "q]  
.......

.Z#/%y3S  
'Ebjn>"  
{>
8?6m-  
增加内存啊，现在内存很便宜，只要你的主板能插几根内存，都给他插上，另外打听一下有没有外置内存插槽 JVxja<43  
aGfp"NtL  
用MLFMM算，很耗内存，算个圆球还可以，算飞机就不要奢望了

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

!ka* rd  
*(?Wzanh  
Szgo@x$^  
呵呵，我估计你用他那个机器算300MHz的F117，也不一定能算，速度不会快到那去，你可以请他算算，看看速度和内存消耗有多少 F#sm^%_2  
用工作站应该是王道！ SXm%X(JU  
    

引用第32楼jetwolf于2008-11-28 22:04发表的 :
呵呵，我估计你用他那个机器算300MHz的F117，也不一定能算，速度不会快到那去，你可以请他算算，看看速度和内存消耗有多少 G;Jqby8d  
用工作站应该是王道！ p}I,!~}  
.......

LVt{`   
(CDwl,  

图片:RCS_300MHz.jpg

IWc?E  
M`{~AIqd(  
300MHz  ... MOM+MLFMM(CFIE) |:Q`9;  
Plane Wave 正向打入， 在 XY 平面上每15度取1點的時機器速度跟記憶體需求。 XpS].P9  
`0'Bg2'  
          Scaling of the geometry with factor SKAL =  1.00000E-03 bLSXQStB  
          Surface of all triangles in m*m:     242.025892 bruM#T@}  
NFxs4:] RT  
                           FAST MULTIPOLE METHOD (FMM) C 0@tMB7  
 Multilevel FMM is used gK7bP'S8H  
 Storage of some elements with single precision (saves memory) w=Cqv~  
/1O
zX'5f  
                                DATA FOR MEMORY USAGE *%6NuZ  
 Number of metallic triangles:        52614                   max. triangles:MAXNDR   =    52614 vOMmsU F  
 Number of dielectric triangles:          0 $=) i{kGS@  
 Number of diel. GO triangles:            0 o$disJ  
 Number of FEM surface triangles:         0 bUJ5jkZ)  
 Number of metallic segments:             0                   max. segments: MAXNSEG  =        0 UM[<v9NWE  
 Number dielectr./magnet. cuboids:        0                   max. cuboids:  MAXNQUA  =        0
n^kszIu~  
 Number of tetrahedral elements:          0                   max. tetrah.:  MAXNTETRA=        0 %
uA\Le  
 Number of edges in PO region:            0                   max. edges:    MAXPOKA  =        0 (?[^##03MN  
 Number of wedges in PO region:           0                   max. wedges:   MAXPOKL  =        0 c:>&iB-Yu  
 Number of Fock regions:                  0                   max. Fock reg.:MAXFOGE  =        0 m 3k}iIU7  
 Number of polygonal surfaces:            0                   max. surfaces: MAXPOLYF =        0 r92C^h0  
                                                              max. corner p.:MAXPOLYP =        0 y<#?z 8P  
 Number of UTD cylinders:                 0 GP|G[  
 Number of metallic edges (MoM):        78921  unknown:    78921   max. edges:   MAXNKA    =    78921 gpt98:w:  
                                            0  unknown:        0   (magnet.) 35>VCjCw0  
 Number of metallic edges (PO):             0  unknown:        0   (electr.)         0 (magnet.) >M1m(u84#  
 Number of dielectric edges (MoM):          0  unknown:        0   (electr.) Qqju6}+  
                                            0                  0   (magnet.) X"'c2gaa_  
 Number of dielectric edges (PO):           0  unknown:        0   (electr.)         0 (magnet.) 8~q%H1[I\N  
 Number of edges FEM/MoM surface:           0  unknown:        0   (electr.) ZMFV iE;8  
                                            0                  0   (magnet.) ~x}/>-d  
 Number of nodes between segments:          0  unknown:        0   max. nodes:   MAXNKNO   =        0 vbXZZ  
 Number of connection points:               0  unknown:        0   max. conn.:   MAXNV     =        0 02|f@bP.  
 Number of dielectric cuboids:              0  unknown:        0   max. cuboids: MAXNQUA   =        0 ? dD<KCbP,  
 Number of magnetic cuboids:                0  unknown:        0 iobL6SUZ  
 Number of basis funct. for MoM:        78921  unknown:    78921   max. basisf.  MAXNZEILE =    78921 4|PNsHXt
 
 Number of basis funct. for PO:             0  unknown:        0   max. basisf.  MAXNKAPO  =        0 ah1DuTT/G  
{Vg8pt  
lL6bIjf  
                    SUMMARY OF REQUIRED TIMES IN SECONDS r@j$$Pk`  
                                            CPU-time      runtime n4*
hQi+d  
 Reading and constructing the geometry         1.891        1.890 +vR$%  
 Checking the geometry                         1.578        1.578 Xn<~ln  
 Initialisation of the Greens function         0.000        0.000 H?ZlJ|/c  
 Calcul. of coupling for PO/Fock               0.000        0.000 e@B
+\1  
 Calcul. of the FMM transfer function          0.657        0.656  3)5Gzn  
 Fourier transform of FMM basis funct.         1.031        1.031 X}!r4<;(  
 Calcul. of matrix elements                   60.156       60.157 C9?R*2L>  
 Calcul. of right-hand side vector             0.047        0.047 g(9\r  
 Preconditioning system of linear eqns.       92.031       92.031
j9sK P]w  
 Solution of the system of linear eqns.       10.625       10.625 ,VtrQb)Yf  
 Determination of surface currents             0.000        0.000 k{fTqKS%h  
 Calcul. of impedances/powers/losses           0.000        0.000 
#+X|,0p  
 Calcul. of averaged SAR values                0.000        0.000 G|WO  
 Calcul. of power ideal receiving ant.         0.000        0.000 }M9R5!=q  
 Calcul. of cable coupling                     0.000        0.000 Q"2t:  
 Calcul. of electric near field                0.000        0.000 Km)5;BQxg  
 Calcul. of magnetic near field                0.000        0.000 g\Gx
oR  
 Calcul. of far field                          0.234        0.235 mQV
c ZV  
                            other              0.281        0.281 Ib$*w)4:  
                                         -----------  ----------- ?B&Z x-krd  
                         total times:        168.531      168.531 ~hJ/&,vH!  
               (total times in hours:          0.047        0.047) AYB =iLa  
 Specified CPU-times are referring to the master process only Yn0l}=, n  
 Sum of the CPU-times of all processes:     1348.249 seconds (    0.375 hours) @ P[o  
 On average per process:                     168.531 seconds (    0.047 hours) ,(%?j]_P2  
 Peak memory usage during the whole solution:  145.734 MByte ]6PX4oK_t  
 (refers to the master process only) b?h
)~j5  
 Sum of the peak memory of all processes:        1.182 GByte cyR K&J  
 On average per process:                       151.294 MByte  pt`^4}  
Finished: 2008-11-29 12:20:34                    

300MHz  ... MOM+MLFMM(CFIE) #(-V^T  
Plane Wave 正向打入， 在 XY 平面上每5度取1點的時機器速度跟記憶體需求。 b Mi,z3z  
x
N0n0  

图片:RCS_300MHz_POSTFEKO.jpg

{E!ie{~  
iZGc'y  
                    SUMMARY OF REQUIRED TIMES IN SECONDS DLe?@R5  
oqOv"yLJ:  
                                            CPU-time      runtime ! FbW7"yE  
Reading and constructing the geometry         1.907        1.906 HMJx[ yD  
Checking the geometry                         1.625        1.625 ; k{w@L.@  
Initialisation of the Greens function         0.000        0.000 ,tv P"@d  
Calcul. of coupling for PO/Fock               0.000        0.000 lN*O</L,"  
Calcul. of the FMM transfer function          0.641        0.641 s2)a8<  
Fourier transform of FMM basis funct.         1.031        1.031 Y~(
#_K  
Calcul. of matrix elements                   60.156       60.156 `s(T(l  
Calcul. of right-hand side vector             0.032        0.031 !vHUe*1a{  
Preconditioning system of linear eqns.      104.515      104.516 i]gF 6:&  
Solution of the system of linear eqns.       10.313       10.312
~]BxM9  
Determination of surface currents             0.000        0.000 61>@-55k9  
Calcul. of impedances/powers/losses           0.000        0.000 .ICGGC`O  
Calcul. of averaged SAR values                0.000        0.000 Q~<$'j  
Calcul. of power ideal receiving ant.         0.000        0.000 -qz;  
Calcul. of cable coupling                     0.000        0.000 }5-^:}gL
 
Calcul. of electric near field                0.000        0.000 74ma   
Calcul. of magnetic near field                0.000        0.000 "WR)a`$UR  
Calcul. of far field                          0.140        0.141 i>i@r ;:|  
                            other              0.312        0.313 L,V\g^4$K  
                                         -----------  ----------- |H 0+.f;  
                         total times:        180.672      180.672 ST1PSuC~  
               (total times in hours:          0.050        0.050) '0D2e  
%YLdie6c  
Specified CPU-times are referring to the master process only tKgPKWP   
Sum of the CPU-times of all processes:     1445.374 seconds (    0.401 hours) v+A$CGH96  
On average per process:                     180.672 seconds (    0.050 hours) i]!CH2\  
2[: *0 DV#  
Peak memory usage during the whole solution:  145.734 MByte C4e3Itc9X  
(refers to the master process only) !Vl>?U?AN  
Sum of the peak memory of all processes:        1.182 GByte H Yt&MK  
On average per process:                       151.294 MByte x0B|CO  
=7pLU+ u  
Finished: 2008-11-29 12:59:46                    

谢谢webhugo网友的计算。 <*u[<  
对于你的计算结果有两个疑问： H@%Y!z@\  
1）"Sum of the peak memory of all processes:        1.182 GByte ” 6Ouy%]0$I3  
    这个是不是说明对内存需求没有想像中那么大？ O@St^o*A}  
2）从计算结果看，F117的RCS并没有达到公布数据中0.1m2，那么小的量级 _
 h\wH;  
    而且RCS最大点位置也没有在正测向，而是在尾喷口（后向）。 * Zb-
YA  
    这一趋势跟我之前用PO算J10模型的结果吻合。 Zn&S7a>7  
大家是这么看这个问题呢？ ) 9h5a+Z  
  

引用第35楼robben1234于2008-11-29 13:35发表的 : E41ay:duAl  
谢谢webhugo网友的计算。 l\Q--

 
对于你的计算结果有两个疑问： <Mt>v2a3Y  
1）"Sum of the peak memory of all processes:        1.182 GByte ” !=-{$& {  
    这个是不是说明对内存需求没有想像中那么大？ ;ui=7[Us  
2）从计算结果看，F117的RCS并没有达到公布数据中0.1m2，那么小的量级 /t4#-vz  
.......

&jbZL5  
t08E 2sI  
1) YES, 在 300MHz 真的是不用很大的內存。IF, 你使用 MOM+MLFMM 的方法來做。 p3Ey[kURp  
是很節省記憶體的..不然就不會有人說他是電大尺寸的終結者。 h$[tEmD%  
Eku+&f@RB  
我記得 CST 也有 MLFMM 的方法..但是跟 FEKO 比起來真的是遜很多。一樣的問題 FEKO 很快就算完了.. CST 的計算時間是 FEKO 的20倍以上。 \&KfIh8  
個人還是認為要評估軟件還是要深入的去用，不要相信廣告，很多人都是看了廣告後就覺得可行.但是沒仔細廣告後面隱藏的東西。 bhqV2y*'  
\$Jz26 -n  
2) CFIE 的部分，僅適合封閉區面，是可以再加速計算速度。 :u ruC
  
PS.. 由於我自己還不是很擅長修模型.... 不知道有誰能幫我將較詳細的圖做好，進氣口就留進氣口， EG7.FjnVu  
只要進氣口跟結構能摳出封閉的區域即可。我想在解這類的問題對 FEKO 不事件難事。 ['`Vg=O.{  
接下來只剩下 CPU Time 而已。 ,MmX(O0  
Example SU8vz/\%y  
----------- rV5QKz6'  
   AIR    ] :iKk"r,2P[  
             ] PQ}owEJ2eM  
----------- F\)?Ntj)>@  
B%L dH  
3) 個人認為只要能先確定可以算..接下來只要合理的假設，應該都可以做出您要的答案。 t XbMP  
7uI~Xo?N  
4) 能否很明確的告知，您想從那個方向看..想看什麼..問題就是這麼簡單。

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

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

上面我是假設 1個平面波射入在每一個方向看到的值 B\\M%!a>  
MuF{STE>->  
進氣口和尾噴口等... 真的是個大問題...But..我對 3D CAD熟悉度不夠..可能還是要麻煩大家的幫忙.. COH9E
\ZGF  
FJ] ?45  
但是同一個網格模型情況下怎麽對不同部位進行不同仿真方法設置呢？ tvcM< e20  
<--- 我不是很清楚您的意思.. But, 假使模型都一樣.. plane wave 也一樣，可以藉由儲存 STR檔(PS Card)，只對 RCS 的部分做後處理計算即可。 y`a]##1j$M  
W.nr&yiQ  
_g%,/y 9y  
IF 平面波射入.. 自己看到自己的值喔 .. f117_monostatic_rcs_300MHz 1eOQ;#OV  
總共計算了361個點 ... 浪費很多時間..(Calcul. of far field ) dmPAPCm%y  
j2s{rQQ  

描述:f117_monostatic_rcs_CADFEKO

图片:f117_monostatic_rcs_CADFEKO.jpg

z<OfSS_]R  
-0>s`ruor  

描述:F117Jetwolf_Fix_f117-100f117_monostatic_rcs_300MHz

图片:f117_monostatic_rcs_300MHz.jpg

pM}n)Q!{3"  
Dd'J"|jF38  
                    SUMMARY OF REQUIRED TIMES IN SECONDS cv"Bhql  
j7yUya&
 
                                            CPU-time      runtime KJi8LM  
Reading and constructing the geometry         1.844        1.843 NeR1}W  
Checking the geometry                         1.578        1.579 ?fX`z(Z  
Initialisation of the Greens function         0.000        0.000 -\~HAnh  
Calcul. of coupling for PO/Fock               0.000        0.000 q$kx/6=k  
Calcul. of the FMM transfer function          0.640        0.641 h~5gHx/a  
Fourier transform of FMM basis funct.         1.031        1.031 "{x+ \Z\  
Calcul. of matrix elements                   60.140       60.140 md;jj^8zj  
Calcul. of right-hand side vector            12.873       12.865 :aR_f`KMm  
Preconditioning system of linear eqns.       97.407       97.406 ,SH^L|I  
Solution of the system of linear eqns.     3545.168     3545.183 gB])@O%/  
Determination of surface currents             0.000        0.000 =3|5=ZU034  
Calcul. of impedances/powers/losses           0.263        0.265 n'83P%x  
Calcul. of averaged SAR values                0.000        0.000 4P?R "Lk  
Calcul. of power ideal receiving ant.         0.000        0.000 R c  
Calcul. of cable coupling                     0.000        0.000 <
Y)Aez  
Calcul. of electric near field                0.000        0.000 ^_t7{z%sA[  
Calcul. of magnetic near field                0.000        0.000 0?*":o30  
Calcul. of far field                         56.835       56.832 hVW1l&s  
                            other             29.877       29.871 sI`oz
|$  
                                         -----------  ----------- uJU;C.LX  
                         total times:       3807.656     3807.656 D-Q54"^3  
               (total times in hours:          1.058        1.058) oT i$@q  
7`
HKa@  
Specified CPU-times are referring to the master process only B8>3GZi  
Sum of the CPU-times of all processes:    30461.250 seconds (    8.461 hours) Zp+orc7  
On average per process:                    3807.656 seconds (    1.058 hours) 4VP$,|a  
V|)nUsU  
Peak memory usage during the whole solution:  145.736 MByte 4^7*R  
(refers to the master process only) >6gduD!6I  
Sum of the peak memory of all processes:        1.182 GByte #{5h6IC  
On average per process:                       151.294 MByte (97&mhs3  
~\u~>mtchu  
Finished: 2008-11-29 18:22:02                    

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