public void GenerateSarSimulationInput(string outputDirectory)
{
Compound simulationSpace = new Compound("space");
Compound solidSpace = new Compound("solid-space");
simulationSpace.Add(solidSpace);
solidSpace.Add(endo);
var headPhantom = new CSXCAD.Ara.HeadPhantom();
headPhantom.Transformations.Add(new TRotateX(Math.PI / 2));
headPhantom.Transformations.Add(new TTranslate(32.0, 80.0, -headPhantom.Width / 2 - 7.0)); // TODO: Make endo width/height accessibles
solidSpace.Add(headPhantom);
double airBox = 40;
double envResolution = Math.Round(lambdaMin / 20 / unit);
double maxRatio = 1.5;
RectilinearGrid grid = new BoundingGrid_6x3();
Vector3D dutPosition = new Vector3D(
antenna.AbsoluteTransformation.X,
antenna.AbsoluteTransformation.Y,
antenna.AbsoluteTransformation.Z);
solidSpace.Transformations.Add(new TTranslate(-dutPosition));
grid.Move(-dutPosition);
grid.Add(new Vector3D(0, 0, 0));
grid.ZLines.Add(-(antenna.Parent as CSXCAD.Ara.PCB).Thickness);
grid.Sort();
grid.SmoothMesh(dutResolution, maxRatio);
grid.Add(headPhantom.BoundingBox.P1);
grid.Add(headPhantom.BoundingBox.P2);
grid.SmoothMesh(sarResolution, maxRatio);
grid.AddAirbox(airBox);
grid.SmoothMesh(envResolution, maxRatio);
simulationSpace.Add(new SARBox("SAR", frequency,
new Vector3D(headPhantom.XGridPoints.First(), headPhantom.YGridPoints.First(), headPhantom.ZGridPoints.First()),
new Vector3D(headPhantom.XGridPoints.Last(), headPhantom.YGridPoints.Last(), headPhantom.ZGridPoints.Last())));
simulationSpace.Add(new NF2FFBox("nf2ff",
new Vector3D(grid.XLines.First(), grid.YLines.First(), grid.ZLines.First()),
new Vector3D(grid.XLines.Last(), grid.YLines.Last(), grid.ZLines.Last()),
lambdaMin / 15 / unit));
grid.AddPML(8);
Simulation fdtd = new Simulation();
fdtd.Excitation = new GaussExcitation(frequency, bandwidth);
// Export
XDocument doc = new XDocument(
new XDeclaration("1.0", "utf-8", "yes"),
new XComment("CyPhy generated openEMS simulation file"),
new XElement("openEMS",
fdtd.ToXElement(),
new XElement("ContinuousStructure",
new XAttribute("CoordSystem", 0),
simulationSpace.ToXElement(),
grid.ToXElement()
)
)
);
System.IO.Directory.CreateDirectory(Path.Combine(outputDirectory));
string openEmsInput = Path.Combine(outputDirectory, "openEMS_input.xml");
doc.Save(openEmsInput);
string nf2ffInput = Path.Combine(outputDirectory, "nf2ff_input.xml");
var nf2ff = new Postprocess.NF2FF(frequency);
nf2ff.ToXDocument().Save(nf2ffInput);
File.WriteAllText(Path.Combine(outputDirectory, "run_sar.cmd"), AraRFAnalysis.Properties.Resources.run_sar);
}
public void GenerateSarSimulationInput(string outputDirectory)
{
Compound simulationSpace = new Compound("space");
Compound solidSpace = new Compound("solid-space");
simulationSpace.Add(solidSpace);
solidSpace.Add(endo);
var headPhantom = new CSXCAD.Ara.HeadPhantom();
headPhantom.Transformations.Add(new TRotateX(Math.PI / 2));
headPhantom.Transformations.Add(new TTranslate(32.0, 80.0, -headPhantom.Width / 2 - 7.0)); // TODO: Make endo width/height accessibles
solidSpace.Add(headPhantom);
double airBox = 40;
double envResolution = Math.Round(lambdaMin / 20 / unit);
double maxRatio = 1.5;
RectilinearGrid grid = new BoundingGrid_6x3();
Vector3D dutPosition = new Vector3D(
antenna.AbsoluteTransformation.X,
antenna.AbsoluteTransformation.Y,
antenna.AbsoluteTransformation.Z);
solidSpace.Transformations.Add(new TTranslate(-dutPosition));
grid.Move(-dutPosition);
grid.Add(new Vector3D(0, 0, 0));
grid.ZLines.Add(-(antenna.Parent as CSXCAD.Ara.PCB).Thickness);
grid.Sort();
grid.SmoothMesh(dutResolution, maxRatio);
grid.Add(headPhantom.BoundingBox.P1);
grid.Add(headPhantom.BoundingBox.P2);
grid.SmoothMesh(sarResolution, maxRatio);
grid.AddAirbox(airBox);
grid.SmoothMesh(envResolution, maxRatio);
simulationSpace.Add(new SARBox("SAR", frequency,
new Vector3D(headPhantom.XGridPoints.First(), headPhantom.YGridPoints.First(), headPhantom.ZGridPoints.First()),
new Vector3D(headPhantom.XGridPoints.Last(), headPhantom.YGridPoints.Last(), headPhantom.ZGridPoints.Last())));
simulationSpace.Add(new NF2FFBox("nf2ff",
new Vector3D(grid.XLines.First(), grid.YLines.First(), grid.ZLines.First()),
new Vector3D(grid.XLines.Last(), grid.YLines.Last(), grid.ZLines.Last()),
lambdaMin / 15 / unit));
grid.AddPML(8);
Simulation fdtd = new Simulation();
fdtd.Excitation = new GaussExcitation(frequency, bandwidth);
// Export
XDocument doc = new XDocument(
new XDeclaration("1.0", "utf-8", "yes"),
new XComment("CyPhy generated openEMS simulation file"),
new XElement("openEMS",
fdtd.ToXElement(),
new XElement("ContinuousStructure",
new XAttribute("CoordSystem", 0),
simulationSpace.ToXElement(),
grid.ToXElement()
)
)
);
System.IO.Directory.CreateDirectory(Path.Combine(outputDirectory));
string openEmsInput = Path.Combine(outputDirectory, "openEMS_input.xml");
doc.Save(openEmsInput);
string nf2ffInput = Path.Combine(outputDirectory, "nf2ff_input.xml");
var nf2ff = new Postprocess.NF2FF(frequency);
nf2ff.ToXDocument().Save(nf2ffInput);
File.WriteAllText(Path.Combine(outputDirectory, "run_sar.cmd"), AraRFAnalysis.Properties.Resources.run_sar);
}
public void GenerateFarFieldSimulationInput(string outputDirectory)
{
Compound simulationSpace = new Compound("space");
Compound solidSpace = new Compound("solid-space");
simulationSpace.Add(solidSpace);
Compound dut;
if (excludeEndo == true)
{
dut = endo.GetModule(slotIndex);
}
else
{
dut = endo;
}
solidSpace.Add(dut); // modifies dut parent (!)
double airBox = 40;
double maxRes = Math.Round(lambdaMin / 20 / unit);
double maxRatio = 1.5;
RectilinearGrid grid = new RectilinearGrid();
grid.Add(dut.BoundingBox.P1);
grid.Add(dut.BoundingBox.P2);
Vector3D antennaPosition = new Vector3D(
antenna.AbsoluteTransformation.X,
antenna.AbsoluteTransformation.Y,
antenna.AbsoluteTransformation.Z);
solidSpace.Transformations.Add(new TTranslate(-antennaPosition));
grid.Move(-antennaPosition);
grid.Add(new Vector3D(0, 0, 0));
grid.ZLines.Add(-(antenna.Parent as CSXCAD.Ara.PCB).Thickness);
grid.Sort();
grid.SmoothMesh(dutResolution, maxRatio);
grid.AddAirbox(airBox);
grid.SmoothMesh(maxRes, maxRatio);
simulationSpace.Add(new NF2FFBox("nf2ff",
new Vector3D(grid.XLines.First(), grid.YLines.First(), grid.ZLines.First()),
new Vector3D(grid.XLines.Last(), grid.YLines.Last(), grid.ZLines.Last()),
lambdaMin / 15 / unit));
grid.AddPML(8);
Simulation fdtd = new Simulation();
fdtd.Excitation = new GaussExcitation(frequency, bandwidth);
// Export
XDocument doc = new XDocument(
new XDeclaration("1.0", "utf-8", "yes"),
new XComment("CyPhy generated openEMS simulation file"),
new XElement("openEMS",
fdtd.ToXElement(),
new XElement("ContinuousStructure",
new XAttribute("CoordSystem", 0),
simulationSpace.ToXElement(),
grid.ToXElement()
)
)
);
if (dut is CSXCAD.Ara.Module)
{
dut.Parent = endo;
}
System.IO.Directory.CreateDirectory(outputDirectory);
string openEmsInput = Path.Combine(outputDirectory, "openEMS_input.xml");
doc.Save(openEmsInput);
string nf2ffInput = Path.Combine(outputDirectory, "nf2ff_input.xml");
var nf2ff = new Postprocess.NF2FF(frequency);
nf2ff.ToXDocument().Save(nf2ffInput);
File.WriteAllText(Path.Combine(outputDirectory, "run_farfield.cmd"), AraRFAnalysis.Properties.Resources.run_farfield);
}
public void GenerateFarFieldSimulationInput(string outputDirectory)
{
Compound simulationSpace = new Compound("space");
Compound solidSpace = new Compound("solid-space");
simulationSpace.Add(solidSpace);
Compound dut;
if (excludeEndo == true)
{
dut = endo.GetModule(slotIndex);
}
else
{
dut = endo;
}
solidSpace.Add(dut); // modifies dut parent (!)
double airBox = 40;
double maxRes = Math.Round(lambdaMin / 20 / unit);
double maxRatio = 1.5;
RectilinearGrid grid = new RectilinearGrid();
grid.Add(dut.BoundingBox.P1);
grid.Add(dut.BoundingBox.P2);
Vector3D antennaPosition = new Vector3D(
antenna.AbsoluteTransformation.X,
antenna.AbsoluteTransformation.Y,
antenna.AbsoluteTransformation.Z);
solidSpace.Transformations.Add(new TTranslate(-antennaPosition));
grid.Move(-antennaPosition);
grid.Add(new Vector3D(0, 0, 0));
grid.ZLines.Add(-(antenna.Parent as CSXCAD.Ara.PCB).Thickness);
grid.Sort();
grid.SmoothMesh(dutResolution, maxRatio);
grid.AddAirbox(airBox);
grid.SmoothMesh(maxRes, maxRatio);
simulationSpace.Add(new NF2FFBox("nf2ff",
new Vector3D(grid.XLines.First(), grid.YLines.First(), grid.ZLines.First()),
new Vector3D(grid.XLines.Last(), grid.YLines.Last(), grid.ZLines.Last()),
lambdaMin / 15 / unit));
grid.AddPML(8);
Simulation fdtd = new Simulation();
fdtd.Excitation = new GaussExcitation(frequency, bandwidth);
// Export
XDocument doc = new XDocument(
new XDeclaration("1.0", "utf-8", "yes"),
new XComment("CyPhy generated openEMS simulation file"),
new XElement("openEMS",
fdtd.ToXElement(),
new XElement("ContinuousStructure",
new XAttribute("CoordSystem", 0),
simulationSpace.ToXElement(),
grid.ToXElement()
)
)
);
if (dut is CSXCAD.Ara.Module)
{
dut.Parent = endo;
}
System.IO.Directory.CreateDirectory(outputDirectory);
string openEmsInput = Path.Combine(outputDirectory, "openEMS_input.xml");
doc.Save(openEmsInput);
string nf2ffInput = Path.Combine(outputDirectory, "nf2ff_input.xml");
var nf2ff = new Postprocess.NF2FF(frequency);
nf2ff.ToXDocument().Save(nf2ffInput);
File.WriteAllText(Path.Combine(outputDirectory, "run_farfield.cmd"), AraRFAnalysis.Properties.Resources.run_farfield);
}
