/// <summary>
/// Case 9: Calculation of internal forces in symmetric section for given state of strain with inclined neutral axis
/// </summary>
public void Case9()
{
// 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.012 * 0.012 * Math.PI / 4.0;
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 concrete = new Concrete();
concrete.SetStrainStressModelBiLinear(25e6, 0.003, 28e9, 0.002);
// steel parameters
Steel steel = new Steel();
steel.SetModelIdealElastoPlastic(550e6, 0.1, 200e9);
// solver creation and parameterization
RCSolver solver = RCSolver.CreateNewSolver(geometry);
solver.SetRebars(rebars);
solver.SetConcrete(concrete);
solver.SetSteel(steel);
//calulation
solver.SolveForces(ResultType.Section, 0.0025, -0.00383423, 0.34906585);
// 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);
// result presentation
sb.AppendLine();
sb.AppendLine(decoration);
sb.AppendLine("Case 9: 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));
}
/// <summary>
/// Case 4: Calculation of internal forces for given state of strain in the section and bilinear model of concrete
/// Case 4a, Case 4b and Case 4c - common geometry, concrete and steel parameters different rebars and calculation
/// </summary>
public void Case4()
{
// Case 4a:
// 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.016 * 0.016 * Math.PI / 4.0;
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 concrete = new Concrete();
concrete.SetStrainStressModelBiLinear(30e6, 0.0035, 32e9, 0.002);
// steel parameters
Steel steel = new Steel();
steel.SetModelIdealElastoPlastic(400e6, 0.1, 205e9);
// solver creation and parameterization
RCSolver solver = RCSolver.CreateNewSolver(geometry);
solver.SetRebars(rebars);
solver.SetConcrete(concrete);
solver.SetSteel(steel);
//calulation
solver.SolveForces(ResultType.Section, 0.0035, 0.0005);
// 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);
// result presentation
sb.AppendLine();
sb.AppendLine(decoration);
sb.AppendLine("Case 4a: Calculation of internal forces for given state of strain in the section and bilinear model of concrete ");
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));
output.Text = output.Text + sb.AppendLine();
//
// Case 4b
rebars.Clear();
rebarArea = 0.032 * 0.032 * Math.PI / 4.0;
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));
solver.SetRebars(rebars);
//calulation
solver.SolveForces(ResultType.Section, 0.0035, -0.02936558);
// 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 4b: Calculation of internal forces for given state of strain in the section and bilinear model of concrete ");
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));
output.Text = output.Text + sb.AppendLine();
//
// Case 4c
rebars.Clear();
rebarArea = 0.032 * 0.032 * Math.PI / 4.0;
rebars.Add(new Rebar(0.05, 0.05, rebarArea));
rebars.Add(new Rebar(0.25, 0.05, rebarArea));
solver.SetRebars(rebars);
//calulation
solver.SolveForces(ResultType.Section, 0.0015, -0.002);
// 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 4c: Calculation of internal forces for given state of strain in the section and bilinear model of concrete ");
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));
output.Text = output.Text + sb.AppendLine();
}
/// <summary>
/// Case 7 Calculation of internal forces for given state of strain in the section and inclined branch model of steel
/// </summary>
public void Case7()
{
// 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>();
rebars.Add(new Rebar(0.05, 0.05, 0.032 * 0.032 * Math.PI / 4.0));
rebars.Add(new Rebar(0.15, 0.05, 0.032 * 0.032 * Math.PI / 4.0));
rebars.Add(new Rebar(0.25, 0.05, 0.032 * 0.032 * Math.PI / 4.0));
rebars.Add(new Rebar(0.05, 0.15, 0.020 * 0.020 * Math.PI / 4.0));
rebars.Add(new Rebar(0.25, 0.15, 0.020 * 0.020 * Math.PI / 4.0));
rebars.Add(new Rebar(0.05, 0.55, 0.012 * 0.012 * Math.PI / 4.0));
rebars.Add(new Rebar(0.25, 0.55, 0.012 * 0.012 * Math.PI / 4.0));
// concrete parameters
Concrete concrete = new Concrete();
concrete.SetStrainStressModelLinear(20e6, 0.0035, 30e9);
// steel parameters
Steel steel = new Steel();
steel.SetModelWithHardening(500e6, 0.075, 200e9, 1.05);
// solver creation and parameterization
RCSolver solver = RCSolver.CreateNewSolver(geometry);
solver.SetRebars(rebars);
solver.SetConcrete(concrete);
solver.SetSteel(steel);
//calulation
solver.SolveForces(ResultType.Section, 0.0035, -0.0070);
// 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);
// result presentation
sb.AppendLine();
sb.AppendLine(decoration);
sb.AppendLine("Case 7 Calculation of internal forces for given state of strain in the section and inclined branch model of steel ");
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));
}