zcy
/
openEMS
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

bug fix: get max freq from xml to setup conducting sheet model 

Signed-off-by: Thorsten Liebig <Thorsten.Liebig@gmx.de>
pull/12/head
Thorsten Liebig 2014-04-08 21:32:18 +02:00
parent d09931e28b
commit 865e817de7
2 changed files with 12 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
matlab/examples/transmission_lines/MSL_Losses.m
View File

@ -31,7 +31,7 @@ f_stop = 25e9;
lambda = c0/f_stop; lambda = c0/f_stop;
%% setup FDTD parameters & excitation function %%%%%%%%%%%%%%%%%%%%%%%%%%%% %% setup FDTD parameters & excitation function %%%%%%%%%%%%%%%%%%%%%%%%%%%%
FDTD = InitFDTD(1e6,1e-4); FDTD = InitFDTD('endCriteria',1e-4);
FDTD = SetGaussExcite(FDTD,0.5*(f_start+f_stop),0.5*(f_stop-f_start)); FDTD = SetGaussExcite(FDTD,0.5*(f_start+f_stop),0.5*(f_stop-f_start));
BC = {'PML_8' 'PML_8' 'PML_8' 'PML_8' 'PEC' 'PML_8'}; BC = {'PML_8' 'PML_8' 'PML_8' 'PML_8' 'PEC' 'PML_8'};
FDTD = SetBoundaryCond( FDTD, BC ); FDTD = SetBoundaryCond( FDTD, BC );
@ -56,11 +56,11 @@ CSX = AddBox( CSX, 'RO4350B', 0, start, stop );
CSX = AddConductingSheet( CSX, 'gold', MSL.conductivity, MSL.thickness ); CSX = AddConductingSheet( CSX, 'gold', MSL.conductivity, MSL.thickness );
portstart = [ mesh.x(1), -MSL.width/2, substrate.thickness]; portstart = [ mesh.x(1), -MSL.width/2, substrate.thickness];
portstop = [ mesh.x(1)+MSL.port_dist, MSL.width/2, 0]; portstop = [ mesh.x(1)+MSL.port_dist, MSL.width/2, 0];
[CSX,portstruct{1}] = AddMSLPort( CSX, 999, 1, 'gold', portstart, portstop, 0, [0 0 -1], 'ExcitePort', 'excite', 'FeedShift', 10*resolution, 'MeasPlaneShift', MSL.port_dist); [CSX, port{1}] = AddMSLPort( CSX, 999, 1, 'gold', portstart, portstop, 0, [0 0 -1], 'ExcitePort', 'excite', 'FeedShift', 10*resolution, 'MeasPlaneShift', MSL.port_dist);
portstart = [mesh.x(end), -MSL.width/2, substrate.thickness]; portstart = [mesh.x(end), -MSL.width/2, substrate.thickness];
portstop = [mesh.x(end)-MSL.port_dist, MSL.width/2, 0]; portstop = [mesh.x(end)-MSL.port_dist, MSL.width/2, 0];
[CSX,portstruct{2}] = AddMSLPort( CSX, 999, 2, 'gold', portstart, portstop, 0, [0 0 -1], 'MeasPlaneShift', MSL.port_dist ); [CSX, port{2}] = AddMSLPort( CSX, 999, 2, 'gold', portstart, portstop, 0, [0 0 -1], 'MeasPlaneShift', MSL.port_dist );
start = [mesh.x(1)+MSL.port_dist, -MSL.width/2, substrate.thickness]; start = [mesh.x(1)+MSL.port_dist, -MSL.width/2, substrate.thickness];
stop = [mesh.x(end)-MSL.port_dist, MSL.width/2, substrate.thickness]; stop = [mesh.x(end)-MSL.port_dist, MSL.width/2, substrate.thickness];
@ -80,8 +80,7 @@ RunOpenEMS( Sim_Path, Sim_CSX ,'');
%% post-processing %% post-processing
close all close all
f = linspace( f_start, f_stop, 1601 ); f = linspace( f_start, f_stop, 1601 );
port{1} = calcPort( portstruct{1}, Sim_Path, f, 'RefImpedance', 50); port = calcPort(port, Sim_Path, f, 'RefImpedance', 50);
port{2} = calcPort( portstruct{2}, Sim_Path, f, 'RefImpedance', 50);
s11 = port{1}.uf.ref./ port{1}.uf.inc; s11 = port{1}.uf.ref./ port{1}.uf.inc;
s21 = port{2}.uf.ref./ port{1}.uf.inc; s21 = port{2}.uf.ref./ port{1}.uf.inc;

9
openems.cpp
View File

@ -712,7 +712,14 @@ int openEMS::SetupFDTD(const char* file)
if ((m_CSX->GetQtyPropertyType(CSProperties::LORENTZMATERIAL)>0) || (m_CSX->GetQtyPropertyType(CSProperties::DEBYEMATERIAL)>0)) if ((m_CSX->GetQtyPropertyType(CSProperties::LORENTZMATERIAL)>0) || (m_CSX->GetQtyPropertyType(CSProperties::DEBYEMATERIAL)>0))
FDTD_Op->AddExtension(new Operator_Ext_LorentzMaterial(FDTD_Op)); FDTD_Op->AddExtension(new Operator_Ext_LorentzMaterial(FDTD_Op));
if (m_CSX->GetQtyPropertyType(CSProperties::CONDUCTINGSHEET)>0) if (m_CSX->GetQtyPropertyType(CSProperties::CONDUCTINGSHEET)>0)
FDTD_Op->AddExtension(new Operator_Ext_ConductingSheet(FDTD_Op,m_Exc->GetMaxFrequency())); {
double f_max=0;
FDTD_Opts->QueryDoubleAttribute("f_max",&f_max);
if (f_max>0)
FDTD_Op->AddExtension(new Operator_Ext_ConductingSheet(FDTD_Op,f_max));
else
cerr << __func__ << ": Error, max. frequency is <=0, disabling conducting sheet material..." << endl;
}
//check all properties to request material storage during operator creation... //check all properties to request material storage during operator creation...
SetupMaterialStorages(); SetupMaterialStorages();