/// <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)); }