mirror of
https://github.com/jeremytammik/RevitSdkSamples.git
synced 2026-08-13 16:27:42 +00:00
Jeremy Tammik
171 lines
8.1 KiB
C#
171 lines
8.1 KiB
C#
//
|
|
// (C) Copyright 2003-2013 by Autodesk, Inc.
|
|
//
|
|
// Permission to use, copy, modify, and distribute this software in
|
|
// object code form for any purpose and without fee is hereby granted,
|
|
// provided that the above copyright notice appears in all copies and
|
|
// that both that copyright notice and the limited warranty and
|
|
// restricted rights notice below appear in all supporting
|
|
// documentation.
|
|
//
|
|
// AUTODESK PROVIDES THIS PROGRAM "AS IS" AND WITH ALL FAULTS.
|
|
// AUTODESK SPECIFICALLY DISCLAIMS ANY IMPLIED WARRANTY OF
|
|
// MERCHANTABILITY OR FITNESS FOR A PARTICULAR USE. AUTODESK, INC.
|
|
// DOES NOT WARRANT THAT THE OPERATION OF THE PROGRAM WILL BE
|
|
// UNINTERRUPTED OR ERROR FREE.
|
|
//
|
|
// Use, duplication, or disclosure by the U.S. Government is subject to
|
|
// restrictions set forth in FAR 52.227-19 (Commercial Computer
|
|
// Software - Restricted Rights) and DFAR 252.227-7013(c)(1)(ii)
|
|
// (Rights in Technical Data and Computer Software), as applicable.
|
|
//
|
|
|
|
using System;
|
|
using System.Collections.Generic;
|
|
using System.Linq;
|
|
using System.Text;
|
|
using Autodesk.CodeChecking.Concrete;
|
|
|
|
namespace ConcreteCalculationsExample
|
|
{
|
|
partial class Example
|
|
{
|
|
|
|
/// <summary>
|
|
/// Case 6: Calculation of internal forces for given state of strain in the section and power-rectangular model of concrete
|
|
/// Case 6a, Case 6b and Case 6c - common geometry and steel parameters different rebars, concrete exponent and calculation
|
|
/// </summary>
|
|
public void Case6()
|
|
{
|
|
// Case 6a:
|
|
// 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.SetStrainStressModelPowerRectangular(30e6, 0.0035, 32e9, 0.002, 1.4);
|
|
// steel parameters
|
|
Steel steel = new Steel();
|
|
steel.SetModelIdealElastoPlastic(400e6, 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.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 6a: Calculation of internal forces for given state of strain in the section and power-rectangular 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));
|
|
//
|
|
// Case 6b
|
|
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);
|
|
// concrete parameters
|
|
concrete.Power = 1.8;
|
|
solver.SetConcrete(concrete);
|
|
//calulation
|
|
solver.SolveForces(ResultType.Section, 0.0035, -0.03034666);
|
|
// 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 6b: Calculation of internal forces for given state of strain in the section and power-rectangular 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));
|
|
//
|
|
// Case 6c
|
|
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);
|
|
// concrete parameters
|
|
concrete.Power = 1.2;
|
|
solver.SetConcrete(concrete);
|
|
//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 6c: Calculation of internal forces for given state of strain in the section and power-rectangular 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));
|
|
}
|
|
}
|
|
}
|