public static List <FEMResult> Results(List <Bar> bars, List <PointLoad> loads, Vector <double> es, Matrix <double> ed)
{
int nEL = bars.Count;
List <FEMResult> outResults = new List <FEMResult>();
for (int i = 0; i < nEL; i++)
{
FEMResult outResult1 = new FEMResult();
BarDeformation barDef1 = new BarDeformation();
BarDisplacement barDisp1 = new BarDisplacement();
BarForce barForce1 = new BarForce();
BarStrain barStrain1 = new BarStrain();
BarStress barStress1 = new BarStress();
barDisp1.UX = ed[i, 0];
barDisp1.UY = ed[i, 1];
barDisp1.UZ = ed[i, 2];
barDisp1.ObjectId = bars[i].Name;
barDisp1.ResultCase = loads[0].Loadcase.Number;
barDisp1.Position = 0;
barDisp1.Divisions = 2;
barForce1.FX = es[i];
barForce1.ObjectId = bars[i].Name;
barForce1.ResultCase = loads[0].Loadcase.Number;
barForce1.Position = 0;
barForce1.Divisions = 2;
outResult1.barDisplacement = barDisp1;
outResult1.barForce = barForce1;
FEMResult outResult2 = new FEMResult();
BarDeformation barDef2 = new BarDeformation();
BarDisplacement barDisp2 = new BarDisplacement();
BarForce barForce2 = new BarForce();
BarStrain barStrain2 = new BarStrain();
BarStress barStress2 = new BarStress();
barDisp2.UX = ed[i, 3];
barDisp2.UY = ed[i, 4];
barDisp2.UZ = ed[i, 5];
barDisp2.ObjectId = bars[i].Name;
barDisp2.ResultCase = loads[0].Loadcase.Number;
barDisp2.Position = 1;
barDisp2.Divisions = 2;
barForce2.FX = es[i];
barForce2.ObjectId = bars[i].Name;
barForce2.ResultCase = loads[0].Loadcase.Number;
barForce2.Position = 1;
barForce2.Divisions = 2;
outResult2.barDisplacement = barDisp2;
outResult2.barForce = barForce2;
outResults.Add(outResult1);
outResults.Add(outResult2);
}
return(outResults);
}
/***************************************************/
private IEnumerable <IResult> ExtractBarStrain(List <int> ids, List <int> loadcaseIds)
{
List <BarStrain> barStrains = new List <BarStrain>();
IFView view = m_LusasApplication.getCurrentView();
IFResultsContext resultsContext = m_LusasApplication.newResultsContext(view);
string entity = "Strain - Thick 3D Beam";
string location = "Feature extreme";
foreach (int loadcaseId in loadcaseIds)
{
IFLoadset loadset = d_LusasData.getLoadset(loadcaseId);
if (!loadset.needsAssociatedValues())
{
resultsContext.setActiveLoadset(loadset);
}
IFUnitSet unitSet = d_LusasData.getModelUnits();
double forceSIConversion = 1 / unitSet.getForceFactor();
double lengthSIConversion = 1 / unitSet.getLengthFactor();
List <string> components = new List <string>()
{
"Ex", "Ey", "Ez", "Bx", "By", "Bz"
};
d_LusasData.startUsingScriptedResults();
Dictionary <string, IFResultsComponentSet> resultsSets = GetResultsSets(entity, components, location, resultsContext);
foreach (int barId in ids)
{
Dictionary <string, double> featureResults = GetFeatureResults(components, resultsSets, unitSet, barId, "L", 6);
List <string> keys = featureResults.Keys.ToList();
double eX = 0; double eY = 0; double eZ = 0; double bX = 0; double bY = 0; double bZ = 0;
featureResults.TryGetValue("Ex", out eX); featureResults.TryGetValue("Ey", out eY); featureResults.TryGetValue("Ez", out eZ);
featureResults.TryGetValue("Bx", out bX); featureResults.TryGetValue("By", out bY); featureResults.TryGetValue("Bz", out bZ);
//TODO: resolve below identifiers extractable through the API
int mode = -1;
double timeStep = 0;
double position = 0;
int divisions = 0;
BarStrain barStrain = new BarStrain(
barId,
Adapters.Lusas.Convert.GetName(loadset.getName()),
mode,
timeStep,
position,
divisions,
eX,
eY,
eZ,
0,
0,
0,
0,
0,
0);
BH.Engine.Base.Compute.RecordWarning("Please note only axial and shear strains will be returned when pulling BarStrain results.");
barStrains.Add(barStrain);
}
d_LusasData.stopUsingScriptedResults();
d_LusasData.flushScriptedResults();
}
return(barStrains);
}
