/// <summary>
/// Case 16: Calculation of capacity state in symmetric section for fixed axial force
/// </summary>
public void Case16()
{
// Case 16a
// 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.SolveResistanceM(678E3, Axis.x, false);
// 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 16a: Calculation of capacity state in symmetric section for fixed axial force");
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 16b
// geometry definition
geometry = new Geometry();
geometry.Add(0.0, 0.0);
geometry.Add(0.0, 0.3);
geometry.Add(0.6, 0.3);
geometry.Add(0.6, 0.0);
// rebars definition
rebars.Clear();
rebarArea = 0.036 * 0.036 * Math.PI / 4.0;
rebars.Add(new Rebar(0.09, 0.06, rebarArea));
rebars.Add(new Rebar(0.09, 0.12, rebarArea));
rebars.Add(new Rebar(0.09, 0.18, rebarArea));
rebars.Add(new Rebar(0.09, 0.24, rebarArea));
rebars.Add(new Rebar(0.51, 0.06, rebarArea));
rebars.Add(new Rebar(0.51, 0.12, rebarArea));
rebars.Add(new Rebar(0.51, 0.18, rebarArea));
rebars.Add(new Rebar(0.51, 0.24, rebarArea));
// steel parameters
steel = new Steel();
steel.DesignStrength = 420e6;
steel.HardeningFactor = 1.0;
steel.ModulusOfElasticity = 200e9;
steel.StrainUltimateLimit = 0.025;
// solver creation and parameterization
solver = RCSolver.CreateNewSolver(geometry);
solver.SetRebars(rebars);
solver.SetConcrete(concrete);
solver.SetSteel(steel);
//calulation
solver.SolveResistanceM(1700E3, 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 16b: Calculation of capacity state in symmetric section for fixed axial force");
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));
}
/// <summary>
/// Case 10: Calculation of capacity state in symmetric section for bending moment Mx
/// </summary>
public void Case10()
{
// Case 10a
// 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.020 * 0.020 * Math.PI / 4.0;
rebars.Add(new Rebar(0.05, 0.05, rebarArea));
rebars.Add(new Rebar(0.25, 0.05, rebarArea));
// concrete parameters
Concrete concrete = new Concrete();
concrete.SetStrainStressModelRectangular(20e6, 0.0035, 30e9, 0.9);
// steel parameters
Steel steel = new Steel();
steel.SetModelIdealElastoPlastic(500e6, 0.075, 200e9);
// solver creation and parameterization
RCSolver solver = RCSolver.CreateNewSolver(geometry);
solver.SetRebars(rebars);
solver.SetConcrete(concrete);
solver.SetSteel(steel);
//calulation
solver.SolveResistanceM(0, Axis.x, false);
// 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 10a: Calculation of internal forces in symmetric section for given state of strain with inclined neutral axis");
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 10b
// rebars definition
rebarArea = 0.032 * 0.032 * Math.PI / 4.0;
rebars.Clear();
rebars.Add(new Rebar(0.05, 0.05, rebarArea));
rebars.Add(new Rebar(0.05, 0.55, rebarArea));
rebars.Add(new Rebar(0.25, 0.55, rebarArea));
rebars.Add(new Rebar(0.25, 0.05, rebarArea));
// concrete parameters
concrete.SetStrainStressModelParabolicRectangular(30e6, 0.0035, 32e9, 0.0020);
// steel parameters
steel.DesignStrength = 400e6;
steel.StrainUltimateLimit = 0.1;
// solver parameterization
solver.SetRebars(rebars);
solver.SetConcrete(concrete);
solver.SetSteel(steel);
//calulation
solver.SolveResistanceM(0, Axis.x, false);
// 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);
// result presentation
sb.AppendLine();
sb.AppendLine(decoration);
sb.AppendLine("Case 10b: Calculation of internal forces in symmetric section for given state of strain with inclined neutral axis ");
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));
}