static void Main()
{
// PRACTICAL EXAMPLE: CREATE POST-TENSIONED CABLES
// This example will show you how to add post-tensioned cables to your concrete beam.
// This example was last updated 2022-05-03, using the ver. 21.1.0 FEM-Design API.
// DEFINE GEOMETRY
var p1 = new Geometry.FdPoint3d(0.0, 2.0, 0);
var p2 = new Geometry.FdPoint3d(10, 2.0, 0);
var edge = new Geometry.Edge(p1, p2, Geometry.FdVector3d.UnitY());
// CREATE BEAM
Materials.MaterialDatabase materialsDB = Materials.MaterialDatabase.DeserializeStruxml("materials.struxml");
Sections.SectionDatabase sectionsDB = Sections.SectionDatabase.DeserializeStruxml("sections.struxml");
var material = materialsDB.MaterialByName("C35/45");
var section = sectionsDB.SectionByName("Concrete sections, Rectangle, 300x900");
var bar = new Bars.Bar(
edge,
Bars.BarType.Beam,
material,
sections: new Sections.Section[] { section },
connectivities: new Bars.Connectivity[] { Bars.Connectivity.GetRigid() },
eccentricities: new Bars.Eccentricity[] { Bars.Eccentricity.GetDefault() },
identifier: "B");
// CREATE POST-TENSIONED CABLE
var shape = new Reinforcement.PtcShapeType(
start: new Reinforcement.PtcShapeStart(0.0, 0.0),
intermediates: new List <Reinforcement.PtcShapeInner>()
{
new Reinforcement.PtcShapeInner(0.4, -0.20, 0.0, 0.1),
new Reinforcement.PtcShapeInner(0.6, -0.20, 0.0),
},
end: new Reinforcement.PtcShapeEnd(0.0, 0.0, 0.9)
);
var losses = new Reinforcement.PtcLosses(
curvatureCoefficient: 0.05,
wobbleCoefficient: 0.007,
anchorageSetSlip: 0.0,
elasticShortening: 0.0,
creepStress: 0.0,
shrinkageStress: 0.0,
relaxationStress: 0.0);
var manufacturing = new Reinforcement.PtcManufacturingType(
positions: new List <double>()
{
0.3, 0.7
},
shiftX: 0.0,
shiftZ: 0.1);
var strandData = new Reinforcement.PtcStrandLibType(
name: "Custom ptc material",
f_pk: 1860.0,
a_p: 150.0,
e_p: 195000.0,
rho: 7.810,
relaxationClass: 2,
rho_1000: 0.1);
var ptc = new Reinforcement.Ptc(
bar,
shape,
losses,
manufacturing,
strandData,
jackingSide: Reinforcement.JackingSide.Start,
jackingStress: 1000.0,
numberOfStrands: 3,
identifier: "PTC");
var elements = new List <GenericClasses.IStructureElement>()
{
bar,
ptc
};
// CREATE AND OPEN MODEL
Model model = new Model(Country.S, elements);
string path = "output/post_tensioned_cable.struxml";
if (!Directory.Exists("output"))
{
Directory.CreateDirectory("output");
}
model.SerializeModel(path);
var app = new Calculate.Application();
app.OpenStruxml(Path.GetFullPath(path), false);
}
static void Main()
{
// EXAMPLE 9: READ RESULTS
// This example will show you how to model a simple supported beam,
// and read some of the results.
// This example was last updated 2022-07-01, using the ver. 21.2.0 FEM-Design API.
#region DEFINE GEOMETRY
// Define geometry
var p1 = new Geometry.FdPoint3d(2.0, 2.0, 0);
var p2 = new Geometry.FdPoint3d(10, 2.0, 0);
var mid = p1 + (p2 - p1) * 0.5;
// Create elements
var edge = new Geometry.Edge(p1, p2, Geometry.FdVector3d.UnitZ());
Materials.MaterialDatabase materialsDB = Materials.MaterialDatabase.DeserializeStruxml("materials.struxml");
Sections.SectionDatabase sectionsDB = Sections.SectionDatabase.DeserializeStruxml("sections.struxml");
var material = materialsDB.MaterialByName("C35/45");
var section = sectionsDB.SectionByName("Concrete sections, Rectangle, 300x900");
var bar = new Bars.Bar(
edge,
Bars.BarType.Beam,
material,
sections: new Sections.Section[] { section },
connectivities: new Bars.Connectivity[] { Bars.Connectivity.GetRigid() },
eccentricities: new Bars.Eccentricity[] { Bars.Eccentricity.GetDefault() },
identifier: "B");
bar.BarPart.LocalY = Geometry.FdVector3d.UnitY();
var elements = new List <GenericClasses.IStructureElement>()
{
bar
};
#endregion
#region DEFINE SUPPORTS
// Create supports
var s1 = new Supports.PointSupport(
point: p1,
motions: Releases.Motions.RigidPoint(),
rotations: Releases.Rotations.Free()
);
var s2 = new Supports.PointSupport(
point: p2,
motions: new Releases.Motions(yNeg: 1e10, yPos: 1e10, zNeg: 1e10, zPos: 1e10),
rotations: Releases.Rotations.Free()
);
var supports = new List <GenericClasses.ISupportElement>()
{
s1, s2
};
#endregion
#region DEFINE LOAD CASES/COMBINATIONS
// Create load cases
var deadload = new Loads.LoadCase("Deadload", Loads.LoadCaseType.DeadLoad, Loads.LoadCaseDuration.Permanent);
var liveload = new Loads.LoadCase("Liveload", Loads.LoadCaseType.Static, Loads.LoadCaseDuration.Permanent);
var loadcases = new List <Loads.LoadCase>()
{
deadload, liveload
};
// Create load combinations
var slsFactors = new List <double>()
{
1.0, 1.0
};
var SLS = new Loads.LoadCombination("SLS", Loads.LoadCombType.ServiceabilityCharacteristic, loadcases, slsFactors);
var ulsFactors = new List <double>()
{
1.35, 1.5
};
var ULS = new Loads.LoadCombination("ULS", Loads.LoadCombType.UltimateOrdinary, loadcases, ulsFactors);
var loadCombinations = new List <Loads.LoadCombination>()
{
SLS, ULS
};
// Create loads
var pointForce = new Loads.PointLoad(mid, new Geometry.FdVector3d(0.0, 0.0, -5.0), liveload, null, Loads.ForceLoadType.Force);
var pointMoment = new Loads.PointLoad(p2, new Geometry.FdVector3d(0.0, 5.0, 0.0), liveload, null, Loads.ForceLoadType.Moment);
var lineLoadStart = new Geometry.FdVector3d(0.0, 0.0, -2.0);
var lineLoadEnd = new Geometry.FdVector3d(0.0, 0.0, -4.0);
var lineLoad = new Loads.LineLoad(edge, lineLoadStart, lineLoadEnd, liveload, Loads.ForceLoadType.Force, "", constLoadDir: true, loadProjection: true);
var loads = new List <GenericClasses.ILoadElement>()
{
pointForce,
pointMoment,
lineLoad
};
#endregion
#region ASSEMBLE
// Add to model
Model model = new Model(Country.S);
model.AddElements(elements);
model.AddSupports(supports);
model.AddLoadCases(loadcases);
model.AddLoadCombinations(loadCombinations);
model.AddLoads(loads);
#endregion
#region SETTINGS
// define the file name
string fileName = "SimpleBeam.struxml";
fileName = Path.GetFullPath(fileName);
// Define the Units
// it is an optional operation and it can be omitted
// Default Units can be seen looking at FemDesign.Results.UnitResults.Default()
var units = new FemDesign.Results.UnitResults(Results.Length.m, Results.Angle.deg, Results.SectionalData.mm, Results.Force.daN, Results.Mass.kg, Results.Displacement.cm, Results.Stress.MPa);
// Select the results to extract
var resultTypes = new List <Results.ResultType>
{
Results.ResultType.PointSupportReaction,
Results.ResultType.NodalDisplacement
};
var bscPathsFromResultTypes = Calculate.Bsc.BscPathFromResultTypes(resultTypes, fileName, units);
#endregion
#region ANALYSIS
// Running the analysis
FemDesign.Calculate.Analysis analysisSettings = new FemDesign.Calculate.Analysis(null, null, null, null, calcCase: true, false, false, calcComb: true, false, false, false, false, false, false, false, false, false);
var fdScript = FemDesign.Calculate.FdScript.Analysis(fileName, analysisSettings, bscPathsFromResultTypes, null, true);
var app = new FemDesign.Calculate.Application();
app.RunFdScript(fdScript, false, true);
model.SerializeModel(fileName);
// Read model and results
model = Model.DeserializeFromFilePath(fdScript.StruxmlPath);
#endregion
#region EXTRACT RESULTS
IEnumerable <Results.IResult> results = Enumerable.Empty <Results.IResult>();
foreach (var cmd in fdScript.CmdListGen)
{
string path = cmd.OutFile;
var _results = Results.ResultsReader.Parse(path);
results = results.Concat(_results);
}
#endregion
#region DO SOMETHING WITH RESULTS
// Display Results on Screen
// The results are grouped by their type
var resultGroups = results.GroupBy(t => t.GetType()).ToList();
foreach (var resultGroup in resultGroups)
{
Console.WriteLine(resultGroup.Key.Name);
Console.WriteLine();
foreach (var result in resultGroup)
{
Console.WriteLine(result);
}
Console.WriteLine();
Console.WriteLine();
}
// Select a specific result
Console.WriteLine("Vertical Reaction Forces");
var zReactions = results.Where(t => t.GetType() == typeof(Results.PointSupportReaction)).Cast <Results.PointSupportReaction>();
foreach (var zReaction in zReactions)
{
Console.WriteLine($"Node {zReaction.Id}: {zReaction.Fz} {units.Force}");
}
#endregion
// ENDING THE PROGRAM
Console.WriteLine("\nPress any key to close console.");
Console.ReadKey();
}
static void Main()
{
// EXAMPLE 1: CREATING A SIMPLE BEAM
// This example will show you how to model a simple supported beam,
// and how to save it for export to FEM-Design.Before running,
// make sure you have a window with FEM-Design open.
// This example was last updated 2022-04-27, using the ver. 21.1.0 FEM-Design API.
// Define geometry
var p1 = new Geometry.FdPoint3d(2.0, 2.0, 0);
var p2 = new Geometry.FdPoint3d(10, 2.0, 0);
var mid = p1 + (p2 - p1) * 0.5;
// Create elements
var edge = new Geometry.Edge(p1, p2, Geometry.FdVector3d.UnitZ());
Materials.MaterialDatabase materialsDB = Materials.MaterialDatabase.DeserializeStruxml("materials.struxml");
Sections.SectionDatabase sectionsDB = Sections.SectionDatabase.DeserializeStruxml("sections.struxml");
var material = materialsDB.MaterialByName("C35/45");
var section = sectionsDB.SectionByName("Concrete sections, Rectangle, 300x900");
var bar = new Bars.Bar(
edge,
Bars.BarType.Beam,
material,
sections: new Sections.Section[] { section },
connectivities: new Bars.Connectivity[] { Bars.Connectivity.GetRigid() },
eccentricities: new Bars.Eccentricity[] { Bars.Eccentricity.GetDefault() },
identifier: "B");
bar.BarPart.LocalY = Geometry.FdVector3d.UnitY();
var elements = new List <GenericClasses.IStructureElement>()
{
bar
};
// Create supports
var s1 = new Supports.PointSupport(
point: p1,
motions: Releases.Motions.RigidPoint(),
rotations: Releases.Rotations.Free()
);
var s2 = new Supports.PointSupport(
point: p2,
motions: new Releases.Motions(yNeg: 1e10, yPos: 1e10, zNeg: 1e10, zPos: 1e10),
rotations: Releases.Rotations.Free()
);
var supports = new List <GenericClasses.ISupportElement>()
{
s1, s2
};
// Create load cases
var deadload = new Loads.LoadCase("Deadload", Loads.LoadCaseType.DeadLoad, Loads.LoadCaseDuration.Permanent);
var liveload = new Loads.LoadCase("Liveload", Loads.LoadCaseType.Static, Loads.LoadCaseDuration.Permanent);
var loadcases = new List <Loads.LoadCase>()
{
deadload, liveload
};
// Create load combinations
var slsFactors = new List <double>()
{
1.0, 1.0
};
var SLS = new Loads.LoadCombination("SLS", Loads.LoadCombType.ServiceabilityCharacteristic, loadcases, slsFactors);
var ulsFactors = new List <double>()
{
1.35, 1.5
};
var ULS = new Loads.LoadCombination("ULS", Loads.LoadCombType.UltimateOrdinary, loadcases, ulsFactors);
var loadCombinations = new List <Loads.LoadCombination>()
{
SLS, ULS
};
// Create loads
var pointForce = new Loads.PointLoad(mid, new Geometry.FdVector3d(0.0, 0.0, -5.0), liveload, null, Loads.ForceLoadType.Force);
var pointMoment = new Loads.PointLoad(p2, new Geometry.FdVector3d(0.0, 5.0, 0.0), liveload, null, Loads.ForceLoadType.Moment);
var lineLoadStart = new Geometry.FdVector3d(0.0, 0.0, -2.0);
var lineLoadEnd = new Geometry.FdVector3d(0.0, 0.0, -4.0);
var lineLoad = new Loads.LineLoad(edge, lineLoadStart, lineLoadEnd, liveload, Loads.ForceLoadType.Force, "", constLoadDir: true, loadProjection: true);
var loads = new List <GenericClasses.ILoadElement>()
{
pointForce,
pointMoment,
lineLoad
};
// Add to model
Model model = new Model(Country.S);
model.AddElements(elements);
model.AddSupports(supports);
model.AddLoadCases(loadcases);
model.AddLoadCombinations(loadCombinations);
model.AddLoads(loads);
// Save model then open in FEM-Design
string path = System.IO.Path.GetFullPath("simple_model.struxml");
model.SerializeModel(path);
var app = new Calculate.Application();
app.OpenStruxml(path, true);
}