这个好办:
用透明:
例子:- /prep7
- ! frequency, material and structure dimensions
- freq=15.e9
- lambda=3.e8/freq
- ra=8.e-3
- rb=ra+2.e-3
- epsr=4.
- loss=1./epsr
- wave1=3.e8/(sqrt(epsr)*freq)
- wave2=3.e8/freq
- ! define the computational domain
- h=wave2/15
- a=rb+3*h
- b=a+3*h
- c=b+4*h
- ! define the elements and materials
- et,1,HF119,1
- et,2,HF119,1,,,1
- mp,murx,1,1.
- mp,perx,1,epsr
- mp,lsst,1,loss
- mp,murx,2,1.
- mp,perx,2,1.
- ! set up numerical model
- sphere,ra,rb,0,180
- block,-a,a,0,a,-a,a
- block,-b,b,0,b,-b,b
- block,-c,c,0,c,-c,c
- vsbv,4,3,,delete,keep
- vsbv,3,2,,delete,keep
- vsbv,2,1,,delete,keep
- csys,2
- vsel,s,loc,x,0,ra/2
- vdel,all
- alls
- csys,0
- vglue,all
- ! meshing
- esize,wave1/10
- type,1
- mat,1
- vmesh,1
- esize,h
- mat,2
- vmesh,4,6,2
- type,2
- vmesh,5
- alls
- nummrg,all
- ! define equivalent source surface
- nsel,s,loc,x,-a,a
- nsel,r,loc,y,0,a
- nsel,r,loc,z,-a,a
- esln,s,1,all
- nsel,s,loc,x,a
- nsel,a,loc,x,-a
- nsel,a,loc,y,a
- nsel,a,loc,z,-a,
- nsel,a,loc,z,a
- sf,all,mxwf
- alls
- ! define boundary condition
- nsel,s,loc,x,c
- nsel,a,loc,x,-c
- nsel,a,loc,y,c
- nsel,a,loc,z,-c
- nsel,a,loc,z, c
- d,all,ax,0.
- nsel,all
- csys,2
- nsel,s,loc,x,ra
- d,all,ax,0.
- csys,0
- allsel,all
- ! define incident plane wave
- plwave,1,0,0,0,180
- fini
- ! perform solution
- /solu
- hfscat,scat
- eqslv,sparse
- antype,harmic
- harfrq,freq
- solve
- fini
- ! calculate normalized RCS
- /post1
- set,1,1
- hfsym,11,,pmc
- prhffar,rcsn,,0,0,,0,180,36,,
- prhffar,rcsn,,90,90,,0,180,36,,
- fini
复制代码
[ 本帖最后由 北极熊甲 于 2007-12-18 19:17 编辑 ] |