[转帖]微波天线辐射的命令流例子
本帖最后由 wdhd 于 2016-3-7 14:19 编辑给大家一篇微波天线辐射的命令流参考，呵呵，谁能告诉我pml层的定义方法和载荷的加载方式，gui格式的，谢拉
/batch,list/title, Half Wavelength Dipole Antenna /nopr/prep7! structure dimensionsfreq=300e6wavel=3.e8/freqbeta=2.*3.1415926535/wavelcurr=1.e-2 ! 1/8 of total current density_l=wavel/2a=0.5*wavel! define elements and materialet,11,200,7et,1,120,1et,2,120,1,,,1mp,murx,1,1.mp,perx,1,1.! define mesh sizeh=wavel/20_n1=4_nz1=nint(_l/h)+1_nz2=4_nz3=4_npml=6! set up computational domain*dim,_a,array,4_a(1)=0_a(2)=_a(1)+a_a(3)=_a(2)+_n1*h_a(4)=_a(3)+_npml*h_c1=0_c2=_c1+_l/2_c3=_c2+_nz2*h_c4=_c3+_nz3*h_c5=_c4+_npml*h*do,_i,1,3 *do,_j,1,3 rect,_a(_i),_a(_i+1),_a(_j),_a(_j+1) *enddo*enddoaglue,allesize,hamesh,allasel,s,loc,x,_a(1),_a(3)asel,r,loc,y,_a(1),_a(3)cm,_area1,areatype,1mat,1esize,,_nz1asel,s,area,,_area1vext,all,,,0,0,(_c2-_c1)esize,,_nz2asel,s,loc,z,_c2vext,all,,,0,0,(_c3-_c2)asel,s,loc,z,_c3vext,all,,,0,0,(_c4-_c3)type,2esize,,_nz1asel,s,loc,z,0asel,u,,,_area1vext,all,,,0,0,(_c2-_c1)esize,,_nz2asel,s,loc,z,_c2asel,r,loc,x,_a(1),_a(3)asel,r,loc,y,_a(1),_a(3)cm,_area1,areaasel,s,loc,z,_c2asel,u,,,_area1vext,all,,,0,0,(_c3-_c2)esize,,_nz3asel,s,loc,z,_c3asel,r,loc,x,_a(1),_a(3)asel,r,loc,y,_a(1),_a(3)cm,_area1,areaasel,s,loc,z,_c3asel,u,,,_area1vext,all,,,0,0,(_c4-_c3)esize,,_npmlasel,s,loc,z,_c4vext,all,,,0,0,(_c5-_c4)asel,a,loc,z,_c1aclear,allesel,s,type,,11edel,allallsnummrg,all! define equivalent source surfacensel,s,loc,x,_a(1),_a(2)nsel,r,loc,y,_a(1),_a(2)nsel,r,loc,z,_c1,_c3esln,s,1,allnsel,s,loc,x,_a(2)nsel,a,loc,y,_a(2)nsel,a,loc,z,_c3sf,all,mxwfnsel,allesel,all! define boundary conditionnsel,s,loc,x,_a(4)nsel,a,loc,y,_a(4)nsel,a,loc,z,_c1nsel,a,loc,z,_c5d,all,ax,0.allsel,all_ll=_c2-_c1_hz=_ll/_nz1*do,_i,0,_nz1 !define sinusoid line current source nsel,s,loc,x,0 nsel,r,loc,y,0 nsel,r,loc,z,_i*_hz cc= sin(beta*(_ll-_i*_hz))*curr !current density distribution *if,cc,eq,0,then cc=1.e-9 *endif bf,all,js,0,0,cc*enddoallsfini! perform solution/solueqslv,sparseantype,harmicharfrq,freqsolvefini/post1set,1,1hfsys,11,pmc,pmc,pec !set up symmetryhfang,,0,360,0,180 !set up radiation solid angleplhffar,field,EF,0,0,,0,180,18,10.,2,Z !plot theta component of E on E-planeplhffar,patt,polar,0,0,,0,180,18 !radiation pattern on E-planeplhffar,dgain,polar,0,0,,0,180,18prhffar,dgain,max !print out directivityprhffar,prad !print out radiation powerfini
又是一大坨
页:
[1]