/// <summary>
/// Case 17: Calculation of capacity state in symmetric section for fixed moment
/// </summary>
public void Case17()
{
// Case 17a
// geometry definition
Geometry geometry = new Geometry();
geometry.Add(0.0, 0.0);
geometry.Add(0.0, 0.6);
geometry.Add(0.3, 0.6);
geometry.Add(0.3, 0.0);
// rebars definition
List <Rebar> rebars = new List <Rebar>();
double rebarArea = 0.040 * 0.040 * Math.PI / 4.0;
rebars.Add(new Rebar(0.03, 0.03, rebarArea));
rebars.Add(new Rebar(0.03, 0.57, rebarArea));
rebars.Add(new Rebar(0.27, 0.57, rebarArea));
rebars.Add(new Rebar(0.27, 0.03, rebarArea));
// concrete parameters
Concrete concrete = new Concrete();
concrete.SetStrainStressModelRectangular(17.12e6, 0.0035, 32e9, 0.8);
// steel parameters
Steel steel = new Steel();
steel.SetModelIdealElastoPlastic(310e6, 0.025, 200e9);
// solver creation and parameterization
RCSolver solver = RCSolver.CreateNewSolver(geometry);
solver.SetRebars(rebars);
solver.SetConcrete(concrete);
solver.SetSteel(steel);
//calulation
solver.SolveResistanceF(184.9E3, Axis.x, true);
// result for rebars
SetOfForces forcesRebar = solver.GetInternalForces(ResultType.Rebars);
// result for concrete
SetOfForces forcesConcrete = solver.GetInternalForces(ResultType.Concrete);
Point2D Gcc = solver.GetStressGravityCenter(ResultType.Concrete);
double Acc = solver.GetConcreteStressArea();
// result for RC section
SetOfForces forces = solver.GetInternalForces(ResultType.Section);
double angle = solver.GetNeutralAxisAngle();
double dist = solver.GetNeutralAxisDistance();
// result presentation
sb.AppendLine();
sb.AppendLine(decoration);
sb.AppendLine("Case 17a: Calculation of capacity state in symmetric section for fixed moment");
sb.AppendLine(decoration);
sb.AppendLine(FormatOutput("Ns", forcesRebar.AxialForce, 6));
sb.AppendLine(FormatOutput("Mxs", forcesRebar.MomentX, 6));
sb.AppendLine(FormatOutput("Mys", forcesRebar.MomentY, 6));
sb.AppendLine(FormatOutput("Ac", Acc, 6));
sb.AppendLine(FormatOutput("Gcx", Gcc.X, 6));
sb.AppendLine(FormatOutput("Gcy", Gcc.Y, 6));
sb.AppendLine(FormatOutput("Nc", forcesConcrete.AxialForce, 6));
sb.AppendLine(FormatOutput("Mxc", forcesConcrete.MomentX, 6));
sb.AppendLine(FormatOutput("Myc", forcesConcrete.MomentY, 6));
sb.AppendLine(FormatOutput("N", forces.AxialForce, 6));
sb.AppendLine(FormatOutput("Mx", forces.MomentX, 6));
sb.AppendLine(FormatOutput("My", forces.MomentY, 6));
sb.AppendLine(FormatOutput("dist", dist, 6));
sb.AppendLine(FormatOutput("angle", angle, 6));
// Case 17b
// rebars definition
rebars.Clear();
rebarArea = 0.032 * 0.032 * Math.PI / 4.0;
rebars.Add(new Rebar(0.03, 0.03, rebarArea));
rebars.Add(new Rebar(0.03, 0.57, rebarArea));
rebars.Add(new Rebar(0.27, 0.57, rebarArea));
rebars.Add(new Rebar(0.27, 0.03, rebarArea));
solver.SetRebars(rebars);
//calulation
solver.SolveResistanceF(138.67E3, Axis.y, true);
// result for rebars
forcesRebar = solver.GetInternalForces(ResultType.Rebars);
// result for concrete
forcesConcrete = solver.GetInternalForces(ResultType.Concrete);
Gcc = solver.GetStressGravityCenter(ResultType.Concrete);
Acc = solver.GetConcreteStressArea();
// result for RC section
forces = solver.GetInternalForces(ResultType.Section);
angle = solver.GetNeutralAxisAngle();
dist = solver.GetNeutralAxisDistance();
// result presentation
sb.AppendLine();
sb.AppendLine(decoration);
sb.AppendLine("Case 17b: Calculation of capacity state in symmetric section for fixed moment");
sb.AppendLine(decoration);
sb.AppendLine(FormatOutput("Ns", forcesRebar.AxialForce, 6));
sb.AppendLine(FormatOutput("Mxs", forcesRebar.MomentX, 6));
sb.AppendLine(FormatOutput("Mys", forcesRebar.MomentY, 6));
sb.AppendLine(FormatOutput("Ac", Acc, 6));
sb.AppendLine(FormatOutput("Gcx", Gcc.X, 6));
sb.AppendLine(FormatOutput("Gcy", Gcc.Y, 6));
sb.AppendLine(FormatOutput("Nc", forcesConcrete.AxialForce, 6));
sb.AppendLine(FormatOutput("Mxc", forcesConcrete.MomentX, 6));
sb.AppendLine(FormatOutput("Myc", forcesConcrete.MomentY, 6));
sb.AppendLine(FormatOutput("N", forces.AxialForce, 6));
sb.AppendLine(FormatOutput("Mx", forces.MomentX, 6));
sb.AppendLine(FormatOutput("My", forces.MomentY, 6));
sb.AppendLine(FormatOutput("dist", dist, 6));
sb.AppendLine(FormatOutput("angle", angle, 6));
}
