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