private static void SolveQuadCantileverDecompositionTest3()
{
#region Quad Cantilever Model
double youngModulus = 3.0e07;
double poissonRatio = 0.3;
double nodalLoad = 1000;
// Create a new elastic 2D material
var material = new ElasticMaterial2D(StressState2D.PlaneStress)
{
YoungModulus = youngModulus,
PoissonRatio = poissonRatio
};
// Model creation
Model model = new Model();
// Add a single subdomain to the model
model.SubdomainsDictionary.Add(0, new Subdomain(0));
// Add nodes to the nodes dictonary of the model
int indexNode = 0;
for (int i = 0; i < 6; i++)
{
for (int j = 0; j < 3; j++)
{
model.NodesDictionary.Add(indexNode, new Node(id: indexNode++, x: i, y: j, z: 0.0));
}
}
// Constrain left nodes of the model
for (int i = 0; i < 3; i++)
{
model.NodesDictionary[i].Constraints.Add(new Constraint {
DOF = StructuralDof.TranslationX
});
model.NodesDictionary[i].Constraints.Add(new Constraint {
DOF = StructuralDof.TranslationY
});
}
int indexElement = 0;
for (int i = 0; i < 5; i++)
{
for (int j = 0; j < 2; j++)
{
var element = new Element()
{
ID = indexElement,
ElementType = new Quad4(material)
};
element.AddNode(model.NodesDictionary[i * 3 + j]);
element.AddNode(model.NodesDictionary[(i + 1) * 3 + j]);
element.AddNode(model.NodesDictionary[(i + 1) * 3 + j + 1]);
element.AddNode(model.NodesDictionary[i * 3 + j + 1]);
model.ElementsDictionary.Add(indexElement, element);
model.SubdomainsDictionary[0].Elements.Add(element);
indexElement++;
}
}
// Add nodal load values to node 3
model.Loads.Add(new Load()
{
Amount = nodalLoad, Node = model.NodesDictionary[17], DOF = StructuralDof.TranslationY
});
// Needed in order to make all the required data structures
model.ConnectDataStructures();
#endregion
var domainDecomposer = new AutomaticDomainDecomposer(model, 3);
domainDecomposer.UpdateModel();
Dictionary <int, int[]> expectedSubdomains = new Dictionary <int, int[]>()
{
{ 0, new int[] { 0, 2, 1, 3 } },
{ 1, new int[] { 4, 6, 5, 7 } },
{ 2, new int[] { 8, 9 } }
};
for (int i = 0; i < expectedSubdomains.Count; i++)
{
var subdomainElements = model.SubdomainsDictionary[i].Elements;
Assert.Equal(expectedSubdomains[i].Length, model.SubdomainsDictionary[i].Elements.Count);
for (int j = 0; j < expectedSubdomains[i].Length; j++)
{
Assert.Equal(expectedSubdomains[i][j], subdomainElements[j].ID);
}
}
}
private static void SolveQuadCantileverDecompositionTest2()
{
#region Quad Cantilever Model
double youngModulus = 3.0e07;
double poissonRatio = 0.3;
double nodalLoad = 1000;
// Create a new elastic 2D material
var material = new ElasticMaterial2D(StressState2D.PlaneStress)
{
YoungModulus = youngModulus,
PoissonRatio = poissonRatio
};
// Model creation
Model model = new Model();
// Add a single subdomain to the model
model.SubdomainsDictionary.Add(0, new Subdomain(0));
// Add nodes to the nodes dictonary of the model
int indexNode = 0;
for (int i = 0; i < 25; i++)
{
for (int j = 0; j < 5; j++)
{
model.NodesDictionary.Add(indexNode, new Node(id: indexNode++, x: i, y: j, z: 0.0));
}
}
// Constrain left nodes of the model
for (int i = 0; i < 5; i++)
{
model.NodesDictionary[i].Constraints.Add(new Constraint {
DOF = StructuralDof.TranslationX
});
model.NodesDictionary[i].Constraints.Add(new Constraint {
DOF = StructuralDof.TranslationY
});
}
int indexElement = 0;
for (int i = 0; i < 24; i++)
{
for (int j = 0; j < 4; j++)
{
var element = new Element()
{
ID = indexElement,
ElementType = new Quad4(material)
};
element.AddNode(model.NodesDictionary[i * 5 + j]);
element.AddNode(model.NodesDictionary[(i + 1) * 5 + j]);
element.AddNode(model.NodesDictionary[(i + 1) * 5 + j + 1]);
element.AddNode(model.NodesDictionary[i * 5 + j + 1]);
model.ElementsDictionary.Add(indexElement, element);
model.SubdomainsDictionary[0].Elements.Add(element);
indexElement++;
}
}
// Add nodal load values to node 3
model.Loads.Add(new Load()
{
Amount = nodalLoad, Node = model.NodesDictionary[124], DOF = StructuralDof.TranslationY
});
// Needed in order to make all the required data structures
model.ConnectDataStructures();
#endregion
var domainDecomposer = new AutomaticDomainDecomposer(model, 8);
domainDecomposer.UpdateModel();
Dictionary <int, int[]> expectedSubdomains = new Dictionary <int, int[]>()
{
{ 0, new int[] { 0, 4, 1, 5, 8, 9, 2, 6, 10, 12, 13, 14 } },
{ 1, new int[] { 3, 7, 11, 15, 18, 19, 17, 21, 22, 23, 16, 20 } },
{ 2, new int[] { 24, 28, 25, 29, 32, 33, 26, 30, 34, 36, 37, 38 } },
{ 3, new int[] { 27, 31, 35, 39, 42, 43, 41, 45, 46, 47, 40, 44 } },
{ 4, new int[] { 48, 52, 49, 53, 56, 57, 50, 54, 58, 60, 61, 62 } },
{ 5, new int[] { 51, 55, 59, 63, 66, 67, 65, 69, 70, 71, 64, 68 } },
{ 6, new int[] { 72, 76, 73, 77, 80, 81, 74, 78, 82, 84, 85, 86 } },
{ 7, new int[] { 75, 79, 83, 87, 90, 91, 89, 93, 94, 95, 88, 92 } }
};
for (int i = 0; i < expectedSubdomains.Count; i++)
{
var subdomainElements = model.SubdomainsDictionary[i].Elements;
Assert.Equal(expectedSubdomains[i].Length, model.SubdomainsDictionary[i].Elements.Count);
for (int j = 0; j < expectedSubdomains[i].Length; j++)
{
Assert.Equal(expectedSubdomains[i][j], subdomainElements[j].ID);
}
}
}
