public static void EmbeddedExample(Model model)
{
HostElementsBuilder(model);
EmbeddedElementsBuilder(model);
//var embeddedGrouping = new EmbeddedGrouping(model, model.ElementsDictionary.Where(x => x.Key <= 2500000).Select(kv => kv.Value), model.ElementsDictionary.Where(x => x.Key > 2500000).Select(kv => kv.Value), true);
//var embeddedGrouping = new EmbeddedGrouping(model, model.ElementsDictionary.Where(x => x.Key <= 2500).Select(kv => kv.Value), model.ElementsDictionary.Where(x => x.Key > 2500).Select(kv => kv.Value), true);
//var embeddedGrouping = new EmbeddedGrouping(model, model.ElementsDictionary.Where(x => x.Key <= 312500).Select(kv => kv.Value), model.ElementsDictionary.Where(x => x.Key > 312500).Select(kv => kv.Value), true);
//var embeddedGrouping = new EmbeddedGrouping(model, model.ElementsDictionary.Where(x => x.Key <= 1250).Select(kv => kv.Value), model.ElementsDictionary.Where(x => x.Key > 1250).Select(kv => kv.Value), true);
//var embeddedGrouping = new EmbeddedGrouping(model, model.ElementsDictionary.Where(x => x.Key <= 90).Select(kv => kv.Value), model.ElementsDictionary.Where(x => x.Key > 90).Select(kv => kv.Value), true);
//var embeddedGrouping = new EmbeddedGrouping(model, model.ElementsDictionary.Where(x => x.Key <= 90).Select(kv => kv.Value), model.ElementsDictionary.Where(x => x.Key > 90).Select(kv => kv.Value), true);
var embeddedGrouping = new EmbeddedGrouping(model, model.ElementsDictionary.Where(x => x.Key <= 8100).Select(kv => kv.Value), model.ElementsDictionary.Where(x => x.Key > 8100).Select(kv => kv.Value), true);
}
private static Model GetModelFromNEUModel(NEUModel neuModel, bool includeRotationsWhenEmbedding)
{
var model = new Model();
foreach (var node in neuModel.Nodes)
{
model.NodesDictionary.Add(node.Key, node.Value);
}
foreach (var constraint in neuModel.Constraints)
{
model.NodesDictionary[constraint.Item1].Constraints.Add(constraint.Item2);
}
foreach (var neuElement in neuModel.Elements)
{
var neuNodes = neuElement.Value.Item1;
int propertyId = neuElement.Value.Item2;
NEUElementType neuElementType = neuElement.Value.Item3;
NEUTopology neuTopology = neuElement.Value.Item4;
var property = neuModel.Properties[propertyId];
var material = neuModel.Materials[property.MaterialID].GetMaterial(neuTopology);
IFiniteElement elementType;
switch (neuElementType)
{
case NEUElementType.Solid1:
elementType = new Hexa8((IFiniteElementMaterial3D)material);
break;
case NEUElementType.Solid2:
elementType = new Hexa8((IFiniteElementMaterial3D)material);
//elementType = new Hexa20((IFiniteElementMaterial3D)material);
break;
case NEUElementType.Bar:
case NEUElementType.Beam:
case NEUElementType.Beam82:
var beamProperty = (NEUBeamProperty)property;
elementType = new Beam3D(material)
{
SectionArea = beamProperty.Area,
MomentOfInertiaZ = beamProperty.I1,
MomentOfInertiaY = beamProperty.I2,
MomentOfInertiaPolar = beamProperty.J
};
break;
default:
throw new ArgumentException("Unsupported element type.");
}
var element = new Element()
{
ElementType = elementType, ID = neuElement.Key
};
for (int i = 0; i < neuNodes.Count; i++)
{
element.AddNode(model.NodesDictionary[neuNodes[i]]);
}
//element.AddNodes(model.NodesDictionary.Where(x =>
// {
// return neuNodes.IndexOf(x.Key) >= 0;
// }).Select(x => x.Value).ToList<Node>());
model.ElementsDictionary.Add(neuElement.Key, element);
}
var dofTypeTranslations = new Dictionary <DOFType, DOFType>()
{
{ DOFType.X, DOFType.RotX },
{ DOFType.Y, DOFType.RotY },
{ DOFType.Z, DOFType.RotZ }
};
foreach (var load in neuModel.Loads)
{
switch (load.Item2)
{
case NEULoadType.NodalForce:
model.Loads.Add(new Load()
{
Node = model.NodesDictionary[load.Item3], Amount = load.Item5, DOF = load.Item4
});
break;
case NEULoadType.NodalMoment:
model.Loads.Add(new Load()
{
Node = model.NodesDictionary[load.Item3], Amount = load.Item5, DOF = dofTypeTranslations[load.Item4]
});
break;
default:
throw new NotImplementedException("Load type not implemented.");
}
}
var hostGroups = neuModel.Groups.Where(x => x.Value.Item1.ToUpper().StartsWith(NEUMeshGenerator.hostGroupPrefix)).Select(x => x.Value).OrderBy(x => x.Item1).ToArray();
var embeddedGroups = neuModel.Groups.Where(x => x.Value.Item1.ToUpper().StartsWith(NEUMeshGenerator.embeddedGroupPrefix)).Select(x => x.Value).OrderBy(x => x.Item1).ToArray();
var groupSize = hostGroups.Length;
if (groupSize != embeddedGroups.Length)
{
throw new ArgumentException("Number of host groups is not equalt to number of embedded groups. Check group names.");
}
for (int i = 0; i < groupSize; i++)
{
var hostName = hostGroups[i].Item1.Substring(NEUMeshGenerator.hostGroupPrefix.Length);
var embeddedName = embeddedGroups[i].Item1.Substring(NEUMeshGenerator.embeddedGroupPrefix.Length);
if (hostName != embeddedName)
{
throw new ArgumentException(String.Format("Current host/embedded pair name does not match ({0} - {1}).", hostName, embeddedName));
}
var hostElements = model.ElementsDictionary.Where(x => hostGroups[i].Item2.Any(id => id == x.Key)).Select(x => x.Value);
var embeddedElements = model.ElementsDictionary.Where(x => embeddedGroups[i].Item2.Any(id => id == x.Key)).Select(x => x.Value);
var embeddedGrouping = new EmbeddedGrouping(model, hostElements, embeddedElements, includeRotationsWhenEmbedding);
}
var subdomain = new Subdomain()
{
ID = 1
};
foreach (var element in model.Elements)
{
subdomain.ElementsDictionary.Add(element.ID, element);
}
model.SubdomainsDictionary.Add(1, subdomain);
return(model);
}
public static void ExampleWithEmbedded(Model model)
{
HostElementsBuilder(model);
EmbeddedElementsBuilder(model);
var embeddedGrouping = new EmbeddedGrouping(model, model.ElementsDictionary.Where(x => x.Key == 1).Select(kv => kv.Value), model.ElementsDictionary.Where(x => x.Key == 2).Select(kv => kv.Value), true);
}
public static void EmbeddedExampleBuilder(Model model)
{
MatrixModelBuilder(model);
EBEModelBuilder(model);
var embeddedGrouping = new EmbeddedGrouping(model, model.ElementsDictionary.Where(x => x.Key <= 10).Select(kv => kv.Value), model.ElementsDictionary.Where(x => x.Key > 10).Select(kv => kv.Value), true);
}
