shell.gnu

这是一个 York大学开发的FDTD代码。很好用的工具

set output 
set noclip points
set clip one
set noclip two
set border
set boxwidth
set dummy x,y
set format x "%g"
set format y "%g"
set format z "%g"
set nogrid
set key
set nolabel
set noarrow
set nologscale
set logscale x 10
set offsets 0, 0, 0, 0
set nopolar
set angles radians
set noparametric
set view 60, 30, 1, 1
set samples 100, 100
set isosamples 10, 10
set surface
set nocontour
set clabel
set nohidden3d
set cntrparam order 4
set cntrparam linear
set cntrparam levels auto 5
set cntrparam points 5
set size 1,1
set data style lines
set function style lines
set xzeroaxis
set yzeroaxis
set tics in
set ticslevel 0.5
set xtics
set ytics
set ztics
set title "SE of spherical shell radius R thickness d" 0,0
set notime
set rrange [0 : 10]
set trange [-5 : 5]
set urange [-5 : 5]
set vrange [-5 : 5]
set xlabel "frequency (Hz)" 0,0
set xrange [10 : 1e+09]
set ylabel "SE (dB)" 2,0
set yrange [0 : 300]
set zlabel "" 0,0
set zrange [-10 : 10]
set noautoscale r
set noautoscale t
set noautoscale   
set autoscale z
set zero 1e-08
#
dB(a)=20*log10(a)
antidB(a)=10.0**(0.05*a)
parallel(a,b)=a*b/(a+b)
delta(x,sigma)=503.0/sqrt(x*mur*sigma)
R(x,sigma)=106+10*log10(sigma/(mur*x))
A(x,sigma)=8.69*d/delta(x,sigma)
B(x,sigma)=dB(1-exp(-2*d/delta(x,sigma)))
SEPW(x,sigma)=R(x,sigma)+A(x,sigma)+B(x,sigma)
SEthin(x)=dB(1+0.5*Z0/x)
g(x,sigma)=sqrt(j*2*pi*x*sigma*mur*mu0)
denom(x,sigma)=g(x,sigma)*sinh(g(x,sigma)*d)
num(x)=6*pi*pi*x*x*eps0*mur*mu0*R
SEE(x,sigma)=-dB(abs(num(x)/denom(x,sigma)))
SEELF(x,sigma)=-dB(3*pi*x*eps0*R/(sigma*d))
SEEHF(x,sigma)=-dB(6.0*sqrt(2.0)*pi*x*eps0*R*exp(-d/delta(x,sigma))/(sigma*delta(x,sigma)))
SEM(x,sigma)=dB(abs(cosh(g(x,sigma)*d)+0.3333*g(x,sigma)*R*sinh(g(x,sigma)*d)))
SEMLF(x,sigma)=dB(abs(1+0.3333*j*2*pi*x*mu0*R*sigma*d))
SEMHF(x,sigma)=-dB(3.0*sqrt(2.0)*delta(x,sigma)*exp(-d/delta(x,sigma))/R)
#
j = {0.0, 1.0}
c0 = 2.99792e+08
mu0 = 1.25664e-06
eps0 = 8.85419e-12
Z0 = 376.73
# Variables you can adjust for a particular materials
mur = 1.0
d = 0.003
R = 0.238
CFRP = 30000.0
Cu = 5.8e+07
Ni_paint = 800000.0
brass = 1e+07
#
pl SEM(x,CFRP),SEE(x,CFRP),SEPW(x,CFRP)