public void SetUp()
{
nodeFactory = new NodeFactory(ModelType.Truss1D);
start = nodeFactory.Create(0);
end = nodeFactory.Create(1);
elementFactory = new ElementFactory(ModelType.Truss1D);
material = new GenericElasticMaterial(0, 0.1, 0, 0);
section = new SolidRectangle(0.1, 1);
SUT = elementFactory.CreateLinearTruss(start, end, material, section);
}
public void HashCode_changes_with_Material_and_CrossSection()
{
int SUTOriginalHash = SUT.GetHashCode();
LinearTruss equal = elementFactory.CreateLinearTruss(start, end, material, section);
Assert.AreEqual(SUTOriginalHash, equal.GetHashCode());
GenericElasticMaterial material2 = new GenericElasticMaterial(0, 2, 0, 0);
LinearTruss unequalMaterial = elementFactory.CreateLinearTruss(start, end, material2, section);
Assert.AreNotEqual(SUTOriginalHash, unequalMaterial.GetHashCode());
SolidRectangle section2 = new SolidRectangle(2, 1);
LinearTruss unequalCrossSection = elementFactory.CreateLinearTruss(start, end, material, section2);
Assert.AreNotEqual(SUTOriginalHash, unequalCrossSection.GetHashCode());
}
public void Equality_depends_on_material_and_CrossSection()
{
Assert.IsTrue(SUT.Equals(SUT));
LinearTruss equal = elementFactory.CreateLinearTruss(start, end, material, section);
Assert.IsTrue(SUT.Equals(equal));
GenericElasticMaterial material2 = new GenericElasticMaterial(0, 2, 0, 0);
LinearTruss unequalMaterial = elementFactory.CreateLinearTruss(start, end, material2, section);
Assert.IsFalse(SUT.Equals(unequalMaterial));
SolidRectangle section2 = new SolidRectangle(2, 1);
LinearTruss unequalCrossSection = elementFactory.CreateLinearTruss(start, end, material, section2);
Assert.IsFalse(SUT.Equals(unequalCrossSection));
}
/// <summary>
///
/// </summary>
/// <returns></returns>
private double CalculateStiffnessConstant()
{
try
{
LinearConstantSpring spring = this.Element as LinearConstantSpring;
if (spring != null)
{
return(spring.SpringConstant);
}
}
catch (InvalidCastException)
{
// do nothing
}
LinearTruss truss = this.Element as LinearTruss;
return(truss.Material.YoungsModulus * truss.CrossSection.Area / this.Element.OriginalLength);
}
public void AnalyzeSimple2DFrame()
{
FiniteElementModel model = new FiniteElementModel(ModelType.Truss2D);
//build geometric model and constraints
FiniteElementNode node1 = model.NodeFactory.CreateFor2DTruss(0, 0);
model.ConstrainNode(node1, DegreeOfFreedom.X);
model.ConstrainNode(node1, DegreeOfFreedom.Z);
FiniteElementNode node2 = model.NodeFactory.CreateFor2DTruss(0, 3);
FiniteElementNode node3 = model.NodeFactory.CreateFor2DTruss(3, 0);
FiniteElementNode node4 = model.NodeFactory.CreateFor2DTruss(3, 3);
FiniteElementNode node5 = model.NodeFactory.CreateFor2DTruss(6, 0);
model.ConstrainNode(node5, DegreeOfFreedom.Z);
FiniteElementNode node6 = model.NodeFactory.CreateFor2DTruss(6, 3);
IMaterial material = new GenericElasticMaterial(0, 70000000, 0, 0);
ICrossSection section = new SolidRectangle(0.03, 0.01);
LinearTruss truss1 = model.ElementFactory.CreateLinearTruss(node1, node2, material, section);
LinearTruss truss2 = model.ElementFactory.CreateLinearTruss(node1, node3, material, section);
LinearTruss truss3 = model.ElementFactory.CreateLinearTruss(node2, node3, material, section);
LinearTruss truss4 = model.ElementFactory.CreateLinearTruss(node2, node4, material, section);
LinearTruss truss5 = model.ElementFactory.CreateLinearTruss(node1, node4, material, section);
LinearTruss truss6 = model.ElementFactory.CreateLinearTruss(node3, node4, material, section);
LinearTruss truss7 = model.ElementFactory.CreateLinearTruss(node3, node6, material, section);
LinearTruss truss8 = model.ElementFactory.CreateLinearTruss(node4, node5, material, section);
LinearTruss truss9 = model.ElementFactory.CreateLinearTruss(node4, node6, material, section);
LinearTruss truss10 = model.ElementFactory.CreateLinearTruss(node3, node5, material, section);
LinearTruss truss11 = model.ElementFactory.CreateLinearTruss(node5, node6, material, section);
//apply forces
ForceVector force50Z = model.ForceFactory.CreateForTruss(0, -50000);
model.ApplyForceToNode(force50Z, node2);
model.ApplyForceToNode(force50Z, node6);
ForceVector force100Z = model.ForceFactory.CreateForTruss(0, -100000);
model.ApplyForceToNode(force100Z, node4);
//solve model
IFiniteElementSolver solver = new MatrixInversionLinearSolver(model);
FiniteElementResults results = solver.Solve();
//assert results
ReactionVector reactionAtNode1 = results.GetReaction(node1);
Assert.AreEqual(0, reactionAtNode1.X, 1);
Assert.AreEqual(100000, reactionAtNode1.Z, 1);
ReactionVector reactionAtNode5 = results.GetReaction(node1);
Assert.AreEqual(0, reactionAtNode5.X, 1);
Assert.AreEqual(100000, reactionAtNode5.Z, 1);
DisplacementVector displacementAtNode2 = results.GetDisplacement(node2);
Assert.AreEqual(7.1429, displacementAtNode2.X, 0.001);
Assert.AreEqual(-9.0386, displacementAtNode2.Z, 0.001);
DisplacementVector displacementAtNode3 = results.GetDisplacement(node3);
Assert.AreEqual(5.2471, displacementAtNode3.X, 0.001);
Assert.AreEqual(-16.2965, displacementAtNode3.Z, 0.001);
DisplacementVector displacementAtNode4 = results.GetDisplacement(node4);
Assert.AreEqual(5.2471, displacementAtNode4.X, 0.001);
Assert.AreEqual(-20.0881, displacementAtNode4.Z, 0.001);
DisplacementVector displacementAtNode5 = results.GetDisplacement(node5);
Assert.AreEqual(10.4942, displacementAtNode5.X, 0.001);
Assert.AreEqual(0, displacementAtNode5.Z, 0.001);
DisplacementVector displacementAtNode6 = results.GetDisplacement(node6);
Assert.AreEqual(3.3513, displacementAtNode6.X, 0.001);
Assert.AreEqual(-9.0386, displacementAtNode6.Z, 0.001);
}
public void AnalyzeSimple2DBuilding()
{
FiniteElementModel model = new FiniteElementModel(ModelType.Truss2D);
//add nodes
FiniteElementNode node0 = model.NodeFactory.CreateFor2DTruss(0, 0);
FiniteElementNode node1 = model.NodeFactory.CreateFor2DTruss(3, 0);
FiniteElementNode node2 = model.NodeFactory.CreateFor2DTruss(6, 0);
FiniteElementNode node3 = model.NodeFactory.CreateFor2DTruss(9, 0);
FiniteElementNode node4 = model.NodeFactory.CreateFor2DTruss(0, 3);
FiniteElementNode node5 = model.NodeFactory.CreateFor2DTruss(3, 3);
FiniteElementNode node6 = model.NodeFactory.CreateFor2DTruss(6, 3);
FiniteElementNode node7 = model.NodeFactory.CreateFor2DTruss(9, 3);
FiniteElementNode node8 = model.NodeFactory.CreateFor2DTruss(0, 6);
FiniteElementNode node9 = model.NodeFactory.CreateFor2DTruss(3, 6);
FiniteElementNode node10 = model.NodeFactory.CreateFor2DTruss(6, 6);
FiniteElementNode node11 = model.NodeFactory.CreateFor2DTruss(9, 6);
//add constraints
model.ConstrainNode(node0, DegreeOfFreedom.X);
model.ConstrainNode(node0, DegreeOfFreedom.Z);
model.ConstrainNode(node1, DegreeOfFreedom.X);
model.ConstrainNode(node1, DegreeOfFreedom.Z);
model.ConstrainNode(node2, DegreeOfFreedom.X);
model.ConstrainNode(node2, DegreeOfFreedom.Z);
model.ConstrainNode(node3, DegreeOfFreedom.X);
model.ConstrainNode(node3, DegreeOfFreedom.Z);
//properties
IMaterial material = new GenericElasticMaterial(0, 70000000, 0, 0);
ICrossSection section = new SolidRectangle(0.03, 0.01);
//add horizontal elements
LinearTruss truss0 = model.ElementFactory.CreateLinearTruss(node0, node1, material, section);
LinearTruss truss1 = model.ElementFactory.CreateLinearTruss(node1, node2, material, section);
LinearTruss truss2 = model.ElementFactory.CreateLinearTruss(node2, node3, material, section);
LinearTruss truss3 = model.ElementFactory.CreateLinearTruss(node4, node5, material, section);
LinearTruss truss4 = model.ElementFactory.CreateLinearTruss(node5, node6, material, section);
LinearTruss truss5 = model.ElementFactory.CreateLinearTruss(node6, node7, material, section);
LinearTruss truss6 = model.ElementFactory.CreateLinearTruss(node8, node9, material, section);
LinearTruss truss7 = model.ElementFactory.CreateLinearTruss(node9, node10, material, section);
LinearTruss truss8 = model.ElementFactory.CreateLinearTruss(node10, node11, material, section);
//add vert elements
LinearTruss truss9 = model.ElementFactory.CreateLinearTruss(node0, node4, material, section);
LinearTruss truss10 = model.ElementFactory.CreateLinearTruss(node4, node8, material, section);
LinearTruss truss11 = model.ElementFactory.CreateLinearTruss(node1, node5, material, section);
LinearTruss truss12 = model.ElementFactory.CreateLinearTruss(node5, node9, material, section);
LinearTruss truss13 = model.ElementFactory.CreateLinearTruss(node2, node6, material, section);
LinearTruss truss14 = model.ElementFactory.CreateLinearTruss(node6, node10, material, section);
LinearTruss truss15 = model.ElementFactory.CreateLinearTruss(node3, node7, material, section);
LinearTruss truss16 = model.ElementFactory.CreateLinearTruss(node7, node11, material, section);
//add bracing
LinearTruss truss17 = model.ElementFactory.CreateLinearTruss(node1, node6, material, section);
LinearTruss truss18 = model.ElementFactory.CreateLinearTruss(node2, node5, material, section);
LinearTruss truss19 = model.ElementFactory.CreateLinearTruss(node5, node10, material, section);
LinearTruss truss20 = model.ElementFactory.CreateLinearTruss(node6, node9, material, section);
ForceVector force50Z = model.ForceFactory.CreateForTruss(0, 50000);
model.ApplyForceToNode(force50Z, node11);
//model.ApplyForceToNode(force50Z, node8);
IFiniteElementSolver solver = new MatrixInversionLinearSolver(model);
FiniteElementResults results = solver.Solve();
DisplacementVector displacementAtNode11 = results.GetDisplacement(node11);
double dispX = displacementAtNode11.X;
double dispY = displacementAtNode11.Y;
//Assert.AreEqual(7.1429, displacementAtNode11.X, 0.001);
//Assert.AreEqual(-9.0386, displacementAtNode11.Z, 0.001);
//FiniteElementNode node2 = model.NodeFactory.CreateFor2DTruss(0, 3);
}
