Ejemplo n.º 1
0
        public void Calculate_Program()
        {
            string design_file = Path.Combine(system_path, "DESIGN.FIL");
            string boq_file    = Path.Combine(system_path, "BoQ.FIL");

            StreamWriter sw        = new StreamWriter(new FileStream(rep_file_name, FileMode.Create));
            StreamWriter sw_design = new StreamWriter(new FileStream(design_file, FileMode.Create));
            StreamWriter sw_boq    = new StreamWriter(new FileStream(boq_file, FileMode.Create));

            try
            {
                #region Write Report
                sw.WriteLine("\t\t*******************************************************");
                sw.WriteLine("\t\t*                      ASTRA Pro                      *");
                sw.WriteLine("\t\t*       TechSOFT Engineering Services (I) Pvt. Ltd.   *");
                sw.WriteLine("\t\t*                                                     *");
                sw.WriteLine("\t\t*                  DESIGN OF SINGLE SPAN              *");
                sw.WriteLine("\t\t*  ONE WAY CONTINUOUS RCC SLAB BY LIMIT STATE METHOD  *");
                sw.WriteLine("\t\t*******************************************************");
                sw.WriteLine("\t\t-------------------------------------------------------");
                sw.WriteLine("\t\t    THIS RESULT CREATED ON " + System.DateTime.Now.ToString("dd.MM.yyyy  AT HH:mm:ss") + " ");
                sw.WriteLine("\t\t-------------------------------------------------------");
                //sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("USER'S DATA");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("Span Length = {0} m", LengthString);
                sw.WriteLine("Assume Slab Thickness = {0} mm", AssumeSlabThickness);
                sw.WriteLine("Slab Load = {0} kN/sq.m", SlabLoad);
                sw.WriteLine("Finish Thickness = {0} mm", FinishThickness);
                sw.WriteLine("Finish Load = {0} kN/sq.m", FinishLoad);
                sw.WriteLine("Permanent Load Factor = {0}", FixedLoadfactor);
                sw.WriteLine("Imposed Load = {0} kN/sq.m", ImposedLoad);
                sw.WriteLine("Imposed Load Factor = {0}", ImposedLoadFactor);
                sw.WriteLine("Dia of Main Reinforcement = {0} mm", DiaMainReinforcement);
                sw.WriteLine("Dia of Distribution/Temperature Reinforcement = {0} mm", DiaDistributionReinforcement);
                sw.WriteLine("Percentage Distribution Reinforcement = {0} %", PercentageDistReinforcement.ToString("0.00"));
                sw.WriteLine("Grade of Concrete = M {0} N/sq.mm", ConcreteGrade);
                sw.WriteLine("Grade of Steel = Fe {0} N/sq.mm", SteelGrade);
                sw.WriteLine("α = {0}", Alpha);
                sw.WriteLine("β = {0}", Beta);
                sw.WriteLine("γ = {0}", Gamma);
                sw.WriteLine("δ = {0}", Delta);
                sw.WriteLine("λ = {0}", Lamda);
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("DESIGN CALCULATIONS");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("STEP 1 ");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();

                double end_span = 0.0;
                double val      = 0.0;
                if (Lengths.Count > 1)
                {
                    if (Lengths[0] > Lengths[Lengths.Count - 1])
                    {
                        end_span = Lengths[0];
                    }
                    else
                    {
                        end_span = Lengths[Lengths.Count - 1];
                    }
                }

                double mid_span = 0.0d;
                for (int i = 1; i < Lengths.Count - 1; i++)
                {
                    if (Lengths[i] > mid_span)
                    {
                        mid_span = Lengths[i];
                    }
                }


                sw.WriteLine();
                sw.WriteLine("Consider a 1m wide strip and assume overall thickness ");
                sw.WriteLine("of Slab as {0} mm. Owing to symmetry only half the slab ", d_Do);
                sw.WriteLine("need be considered.");
                sw.WriteLine();

                double DL = (d_Do / 1000.0) * 1.0 * d_SL;
                sw.WriteLine("   Dead load of slab (DL) = {0:f3} * 1.0 * {1:f2}", (d_Do / 1000.0), d_SL);
                sw.WriteLine("                          = {0} kN/m", DL);

                double thick_dead_load = (d_Fo / 1000.0) * d_FL;
                sw.WriteLine();
                sw.WriteLine("   Dead load of {0:f2} cm thick finish = {1:f3} * 1.0 * {2}", (d_Fo / 10.0), (d_Fo / 1000.0), d_FL);
                sw.WriteLine("                                     = {0:f3} kN/m", thick_dead_load);
                sw.WriteLine();

                double total_dead_load = DL + thick_dead_load;
                sw.WriteLine("   Total dead load = {0:f2} + {1:f2}", DL, thick_dead_load);
                sw.WriteLine("                   = {0:f2} kN/m", total_dead_load);

                double IL = d_LL * 1.0;
                sw.WriteLine();
                sw.WriteLine("   Imposed Load (IL) = {0} * 1.0 ", d_LL);
                sw.WriteLine("                     = {0:f2} kN/m", IL);

                double factored_dead_load = d_DF * total_dead_load;
                sw.WriteLine();
                sw.WriteLine("   Factored dead load  = {0:f2} * {1:f2}", d_DF, total_dead_load);
                sw.WriteLine("                       = {0:f2} kN/m", factored_dead_load);

                double v1 = factored_dead_load * Lengths[0];
                sw.WriteLine("                       = {0:f2} * {1:f2}", factored_dead_load, Lengths[0]);
                sw.WriteLine("                       = {0:f2} kN/m", v1);

                double factored_imposed_load = ImposedLoadFactor * ImposedLoad;
                sw.WriteLine();
                sw.WriteLine("   Factored imposed load = {0:f2} * {1:f2} ", ImposedLoadFactor, ImposedLoad);
                sw.WriteLine("                         = {0:f2} kN/m", factored_imposed_load);

                double v2 = factored_imposed_load * Lengths[0];
                sw.WriteLine();
                sw.WriteLine("                         = {0:f2} * {1:f2}", factored_imposed_load, Lengths[0]);
                sw.WriteLine("                         = {0} kN", v2);

                sw.WriteLine();
                sw.WriteLine("    Moment at End Span = (DL * l)/12 + (LL * l)/10");
                sw.WriteLine("                       = ({0} * {1})/12 + ({2} * {1})/10", v1, end_span, v2);
                double E_end_moment = ((v1 * end_span) / 12) + ((v2 * end_span) / 10);
                sw.WriteLine("                       = {0} kN/m", E_end_moment.ToString("0.00"));


                sw.WriteLine();
                sw.WriteLine("    Moment at Interior Span = (DL * l)/24 + (LL * l)/12");
                sw.WriteLine("                            = ({0} * {1})/24 + ({2} * {1})/12", v1, mid_span, v2);

                double F_interior_moment = ((v1 * mid_span) / 16) + ((v2 * mid_span) / 12);
                sw.WriteLine("                            = {0} kN/m", F_interior_moment.ToString("0.00"));

                sw.WriteLine();
                sw.WriteLine("    Moment at support next to end support");
                sw.WriteLine("          = (DL * l)/10 + (LL * l)/9");
                sw.WriteLine("          = ({0} * {1})/10 + ({2} * {1})/9", v1, end_span, v2);

                double B_end_support_moment = ((v1 * end_span) / 10) + ((v2 * end_span) / 9);
                sw.WriteLine("          = {0} kN/m", B_end_support_moment.ToString("0.00"));


                sw.WriteLine();
                double max_moment = 0.0;

                if (E_end_moment > F_interior_moment && E_end_moment > B_end_support_moment)
                {
                    max_moment = E_end_moment;
                }
                else if (F_interior_moment > E_end_moment && F_interior_moment > B_end_support_moment)
                {
                    max_moment = F_interior_moment;
                }
                else if (B_end_support_moment > E_end_moment && B_end_support_moment > F_interior_moment)
                {
                    max_moment = B_end_support_moment;
                }

                //sw.WriteLine("Maximum bending moments at critical sections are computed using bending moment");
                ////Table ?
                //sw.WriteLine("coefficients given in Table 11.1 and are tabulated in Table 14.2");

                sw.WriteLine();
                sw.WriteLine("    For a balanced design, Mu = 0.138 * σ_ck * b * d * d");

                double d = Math.Sqrt((max_moment * 10E+5) / (0.138 * d_sigma_ck * 1000));

                sw.WriteLine();
                sw.WriteLine("    d =√(({0} * 10E+5)/(0.138 * {1} * 1000))", max_moment.ToString("0.00"), ConcreteGrade.ToString("0"));
                sw.WriteLine("      = {0:f2} mm", d);

                val  = (int)d / 10;
                val += 2;
                val  = val * 10;

                sw.WriteLine();
                if (val < 100)
                {
                    d = 100;
                    sw.WriteLine("This effective depth is too small, ");
                    sw.WriteLine();
                    sw.WriteLine("Let us adopt an over all depth D = 100 mm");
                    sw.WriteLine();
                    sw.WriteLine("Assume {0} mm bars with {1} mm clear cover.", DiaMainReinforcement, ClearCover);
                }
                else
                {
                    d = val;
                }

                double eff_depth = d - ClearCover - d_d1 / 2;
                sw.WriteLine();
                sw.WriteLine("   Effective depth = {0:f2} - {1:f2} - {2}/2 ", d, ClearCover, d_d1);
                sw.WriteLine("                   = {0:f2} mm", eff_depth);
                sw.WriteLine();


                // At, B


                double a, b, c, At1;

                a = (0.87 * d_sigma_y * d_sigma_y) / (d_sigma_ck * 1000);
                b = 0.87 * d_sigma_y * eff_depth;
                c = B_end_support_moment * 10E+5;
                double b_ac = (b * b) - 4 * a * c;
                At1 = (b) - Math.Sqrt(Math.Abs(b_ac));
                At1 = At1 / (2 * a);

                At1  = (int)At1;
                At1 += 1;
                sw.WriteLine("At B, the area of steel is given by ");
                sw.WriteLine();
                sw.WriteLine("          Mu = 0.87 * σ_y * At(d-((σ_y*At)/(σ_ck*b))");
                sw.WriteLine();
                sw.WriteLine("    {0:f3} * 10E+5 = 0.87 * {1:f2} * At({2:f2} - (({1:f2} * At )/({3:f2} * 1000))",
                             B_end_support_moment,
                             d_sigma_y,
                             eff_depth,
                             d_sigma_ck);
                sw.WriteLine();
                sw.WriteLine("          At = {0:f3} sq.mm/m", At1);
                sw.WriteLine();

                double area1        = Math.PI * d_d1 * d_d1 / 4;
                int    nos_main_rod = (int)(At1 / area1);
                nos_main_rod++;

                double main_spacing = (double)1000.0 / nos_main_rod;
                main_spacing = (int)main_spacing / 10;
                main_spacing = main_spacing * 10;

                double At1_ = area1 * 1000 / main_spacing;

                sw.WriteLine();
                sw.WriteLine("Use {0} mm bars @ {1} mm c/c.", d_d1, main_spacing);
                sw.WriteLine();

                if (At1_ > At1)
                {
                    sw.WriteLine(" Area of steel provided = {0:f0} sq.mm/m > {1:f0} sq.mm/m, OK", At1_, At1);
                }
                else
                {
                    sw.WriteLine(" Area of steel provided = {0:f0} sq.mm/m > {1:f0} sq.mm/m, NOT OK", At1_, At1);
                }


                double At2;

                a    = (0.87 * d_sigma_y * d_sigma_y) / (d_sigma_ck * 1000);
                b    = 0.87 * d_sigma_y * eff_depth;
                c    = E_end_moment * 10E+5;
                b_ac = (b * b) - 4 * a * c;
                At2  = (b) - Math.Sqrt(Math.Abs(b_ac));
                At2  = At2 / (2 * a);

                At2 = (int)At2;
                //At2 += 1;
                sw.WriteLine();
                sw.WriteLine(" At E, the area of steel is given by ");
                sw.WriteLine();
                sw.WriteLine("Or,  {0} * 10E+5 = 0.87 * {1} * At(({2} - ({3}*At)/({4} * 1000))",
                             E_end_moment.ToString("0.00"),
                             d_sigma_y,
                             eff_depth,
                             d_sigma_y,
                             d_sigma_ck);

                sw.WriteLine();
                sw.WriteLine("Or, As = {0:f2} sq.mm/m", At2);

                int nos_dist_rod = (int)(At2 / area1);
                nos_dist_rod++;

                double main_spacing_2 = (double)1000.0 / nos_dist_rod;
                main_spacing_2 = (int)main_spacing_2 / 10;
                main_spacing_2++;
                main_spacing_2 = main_spacing_2 * 10;
                double At2_ = area1 * 1000 / main_spacing_2;
                sw.WriteLine();
                sw.WriteLine("Use {0} mm bars @ {1} mm c/c, As = {2:f0} sq.mm/m",
                             d_d1,
                             main_spacing_2,
                             At2_);

                // At, F

                sw.WriteLine();
                sw.WriteLine("At F, The Area of steel is given by,");
                sw.WriteLine();
                sw.WriteLine("  {0} * 10E+5 = 0.87 * {1} * At({2} - (({3} * At)/({4} * 1000))",
                             F_interior_moment.ToString("0.000"),
                             d_sigma_y,
                             eff_depth,
                             d_sigma_y,
                             d_sigma_ck);


                double At3;

                a    = (0.87 * d_sigma_y * d_sigma_y) / (d_sigma_ck * 1000);
                b    = 0.87 * d_sigma_y * eff_depth;
                c    = F_interior_moment * 10E+5;
                b_ac = (b * b) - 4 * a * c;
                At3  = (b) - Math.Sqrt(Math.Abs(b_ac));
                At3  = At3 / (2 * a);

                At3 = (int)At3;

                int nos_main_rod_1 = (int)(At3 / area1);
                nos_main_rod_1++;

                double main_spacing_1 = (double)1000.0 / nos_main_rod_1;
                main_spacing_1 = (int)main_spacing_1 / 10;
                main_spacing_1 = (main_spacing_1 * 10);
                main_spacing_1--;

                double At3_ = (area1 * 1000) / main_spacing_1;

                sw.WriteLine();
                sw.WriteLine("Use {0:f0} mm bars @ {1:f0} mm c/c",
                             d_d1,
                             main_spacing_1);

                sw.WriteLine("As = {0:f2} sq.mm/m.", At3_);
                sw.WriteLine();
                sw.WriteLine("The Code requires that the maximum bar spacing should not exceed 3d");
                sw.WriteLine();
                sw.WriteLine("Minimum area of high strength main steel");
                sw.WriteLine("       = {0} % of the gross concrete", PercentageDistReinforcement);

                double strength_main_steel = (PercentageDistReinforcement / 100.0) * 1000 * d;

                sw.WriteLine("       = ({0:f2} /100) * 1000 * {1:f2}", PercentageDistReinforcement, d);
                sw.WriteLine("       = {0:f2} sq.mm.m", strength_main_steel);

                sw.WriteLine();
                if (At3_ > strength_main_steel)
                {
                    sw.WriteLine("Thus the actual area provided is OK");
                }
                else
                {
                    sw.WriteLine("Thus the actual area provided is NOT OK");
                }


                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("Temperature reinforcement equal to 0.15% of the gross concrete");
                sw.WriteLine("area will be provided in the longitudinal direction");
                sw.WriteLine();
                //
                double dist_area = (0.15 / 100) * 1000 * 100;
                sw.WriteLine("  (0.15/100) * 1000 * 100 = {0:f3} sq.mm/m", dist_area);
                //area2
                double area2 = Math.PI * d_d2 * d_d2 / 4.0;

                int no_dist_rod = (int)(dist_area / area2);
                no_dist_rod++;

                double dist_spacing = 1000 / no_dist_rod;
                if (dist_spacing > 120)
                {
                    dist_spacing = 120;
                }
                //
                sw.WriteLine();
                sw.WriteLine("Use {0} mm MS bars @ {1} mm c/c", d_d2, dist_spacing);
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("STEP 2 : Check for shear");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine(" The maximum shear force occurs at the first interior support B in the Span AB,");
                sw.WriteLine();

                double Vu = 0.6 * v1 + 0.6 * v2;
                sw.WriteLine(" SF = 0.6 * {0} + 0.6 * {1} = {2} kN",
                             v1, v2, Vu);

                sw.WriteLine();
                double tau_v = Vu / (1000 * eff_depth);
                sw.WriteLine("    Nominal shear stress τ_v = Vu/b*d");
                sw.WriteLine("                             = {0} * 1000 / 1000 * {1}", Vu, eff_depth);
                sw.WriteLine("                             = {0:f4}", tau_v);

                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("At B, percentage of tension steel,");

                double B_percent_tension_steel = (100 * At1_) / (1000 * eff_depth);

                sw.WriteLine("   100 As/b*d = (100 * {0})/(1000 * {1}) = {2}%",
                             At1_.ToString("0.00"), eff_depth, B_percent_tension_steel.ToString("0.0000"));


                ShearValue shV   = new ShearValue();
                double     tau_c = shV.Get_M15(B_percent_tension_steel);
                sw.WriteLine("So, τ_c = {0} N/sq.mm", tau_c);

                double k_value = kV.Get_KValue(d);

                sw.WriteLine("and    k = {0}", k_value);
                double tau_c_dash = k_value * tau_c;

                sw.WriteLine("τ_c` = {0} * {1} > τ_v    OK", k_value, tau_c);
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("STEP 3 : Check for development length");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine("The slab is checked for development length. Let us take the end support A.");

                //
                sw.WriteLine();
                sw.WriteLine("  φ < ((4 * τ_bd)/σ_s)*((1.3*(M1/V) + Lo))");
                //
                //

                double percent_curtain_bar = 50.0;

                sw.WriteLine(" Let us curtail {0}% bars at 0.1 l from the centre of support A.", percent_curtain_bar);
                sw.WriteLine();
                //sw.WriteLine("ταβγδλγσ");
                sw.WriteLine(" M1 = 0.87 * σ_y * At((d - (σ_y * At)/(σ_ck * b))");
                percent_curtain_bar = percent_curtain_bar / 100;

                double M1 = 0.87 * d_sigma_y * percent_curtain_bar * At2_ * ((eff_depth - (d_sigma_y * At2_ * percent_curtain_bar) / (d_sigma_ck * 1000)));

                sw.WriteLine("   = 0.87 * {0} * {1} * {2} * (({3} - ({4} * {5} * {6})/({7} * 1000))",
                             d_sigma_y,
                             percent_curtain_bar,
                             At2_.ToString("0.00"),
                             eff_depth,
                             d_sigma_y,
                             At2_.ToString("0.00"),
                             percent_curtain_bar,
                             d_sigma_ck);
                M1 = M1 / 10E+5;
                sw.WriteLine("    = {0:f2} * 10E+5 N-mm", M1);
                sw.WriteLine("    = {0:f2} kN-m", M1);

                sw.WriteLine();
                double V = 0.4 * v1 + 0.45 * v2;
                sw.WriteLine("   Shear force V = 0.4 * {0} + 0.45 * {1}", v1, v2);
                sw.WriteLine("                 = {0:f2} kN", V);
                sw.WriteLine();
                sw.WriteLine("Assuming Lo = 200 mm including anchor value of U-hook");

                double phi = ((4 * 1.6) / (0.87 * d_sigma_y)) * (((k_value * M1 * 1000) / V) + 200);

                //
                // CC
                sw.WriteLine();
                sw.WriteLine("  φ <= ((4 * 1.6)/(0.87*{0}))*(({1}*{2}*1000)/{3} + 200)",
                             d_sigma_y,
                             k_value,
                             M1.ToString("0.00"),
                             V.ToString("0.00"));

                sw.WriteLine("   = {0:f2} mm.", phi);

                //
                sw.WriteLine();
                sw.WriteLine(" Since the actual bar diameter is less than {0} mm, the slab is safe", phi.ToString("0.000"));
                sw.WriteLine("in development length. The bars must be embedded in the support by");
                sw.WriteLine("at least Ld/3 distance and a U-hook must be provided.");
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("STEP 4 : Check for Deflection.");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine(" The Slab is checked for deflection. Let us consider the end span.");
                sw.WriteLine();
                sw.WriteLine(" L/D <= α * β * γ * δ * λ");
                sw.WriteLine("   α  = {0} for continuous span", Alpha);
                sw.WriteLine("   β  = 1 since L < 10 m");
                sw.WriteLine();

                //double Pt = (100 * (50 * 1000 / 170)) / (1000 * 81);
                double Pt = (100 * (area1 * 1000 / main_spacing_1)) / (1000 * eff_depth);
                // CC
                sw.WriteLine("Pt = 100 * ({0:f0} * 1000/{1:f0})/(1000 * {2:f2})",
                             area1.ToString("0"),
                             main_spacing_1,
                             eff_depth);

                sw.WriteLine("   = {0:f2}%", Pt);
                sw.WriteLine();

                double gama = ModificationFactor.GetGamma(Pt);
                sw.WriteLine(" For a service stress of 240 Mpa in Fe 415 grade steel bar");
                sw.WriteLine(" the value of γ = {0}", gama.ToString("0.00"));
                sw.WriteLine("              δ = 1 since compressive reinforcement is zero.");
                sw.WriteLine("and           λ = 1");
                sw.WriteLine();

                double l_d = Alpha * gama;
                sw.WriteLine();
                sw.WriteLine(" Allowable L/d = {0} * {1:f3}", Alpha, gama.ToString("0.000"));
                sw.WriteLine("               = {0:f2}", l_d);
                sw.WriteLine();
                sw.WriteLine(" Actual    L/d = {0}/{1}", (end_span * 1000), eff_depth);

                val = (end_span * 1000.0 / eff_depth);
                if (val < l_d)
                {
                    sw.WriteLine("               = {0:f2} < {1:f2}, OK", val, l_d);
                }
                else
                {
                    sw.WriteLine("               = {0:f2} > {1:f2}, NOT OK", val, l_d);
                }
                sw.WriteLine();

                #endregion

                sw.WriteLine();
                sw.WriteLine("---------------------------------------------------------------------------");
                sw.WriteLine("---------------------       END OF REPORT        --------------------------");
                sw.WriteLine("---------------------------------------------------------------------------");


                #region Write DESIGN File
                sw_design.WriteLine("SLAB DESIGN 03");
                sw_design.WriteLine();
                sw_design.WriteLine("lengths = {0}", LengthString);
                //lengths = 3000, 3000, 3000, 3000, 3000
                sw_design.WriteLine("span_bredth = {0}", d_b);
                //span_bredth = 8000;
                sw_design.WriteLine("main_dia = {0}", d_d1);
                //main_dia = 8.0;
                sw_design.WriteLine("dist_dia = {0}", d_d2);
                //dist_dia = 6.0;
                sw_design.WriteLine("end_span_reinforcement = {0}", main_spacing_2 * 2);
                //end_span_reinforcement = 340;
                sw_design.WriteLine("support_reinforcement = {0}", main_spacing);
                //support_reinforcement = 140;
                sw_design.WriteLine("interior_reinforcement = {0}", main_spacing_1 * 2);
                //interior_reinforcement = 480;
                sw_design.WriteLine("start_centre_distance = {0}", 300);
                //start_centre_distance = 300;
                sw_design.WriteLine("start_wall_distance = {0}", 450);
                //start_wall_distance = 450;
                sw_design.WriteLine("clear_cover = {0}", ClearCover);
                //clear_cover = 15;
                sw_design.WriteLine("upper_distance = {0}", 900);
                //upper_distance = 900;
                sw_design.WriteLine("lower_distance = {0}", 600);
                //lower_distance = 600;
                sw_design.WriteLine("depth = {0}", d_Do);
                //depth = 120;
                sw_design.WriteLine("wall_thickness = {0}", WallThickness);
                //wall_thickness = 250;
                sw_design.WriteLine("bar_hill = {0}", 25.4);
                //bar_hill = 25.4;
                sw_design.WriteLine("pillar_length = {0}", 500);
                //pillar_length = 700;
                sw_design.WriteLine();
                sw_design.WriteLine("FINISH");
                #endregion

                string size1 = "9.1 x 1.0";
                string bgd   = "Fe 415";
                double d1    = 8;
                double d2    = 6;
                double dist1 = 2935;
                double dist2 = 2635;
                double dist3 = 1800;
                double dist4 = 1000;
                double bno1  = 9;
                double bno2  = 9;
                double bno3  = 16;
                double bno4  = 69;
                double bno5  = 38;


                double bwt1 = 0.00616 * d1 * d1 * dist1 * bno1;

                bwt1 = bwt1 / 10e5;


                double bwt2 = 0.00616 * d1 * d1 * dist2 * bno1;
                bwt2 = bwt2 / 10e5;
                double bwt3 = 0.00616 * d1 * d1 * dist3 * bno2;
                bwt3 = bwt3 / 10e5;
                double bwt4 = 0.00616 * d2 * d2 * dist4 * bno3;
                bwt4 = bwt4 / 10e5;
                double bwt5 = 0.00616 * d2 * d2 * dist4 * bno4;
                bwt5 = bwt5 / 10e5;

                double depth  = 90;
                double shape1 = 2935;
                double shape2 = 2635;
                double shape3 = 1800;
                double shape4 = 1000;


                #region BoQ File


                sw_boq.WriteLine("BoQ Text");
                sw_boq.WriteLine("size1={0}", size1);
                sw_boq.WriteLine("bgd={0}", bgd);
                sw_boq.WriteLine("d1={0}", d1);
                sw_boq.WriteLine("d2={0}", d2);
                sw_boq.WriteLine("dist1={0}", dist1);
                sw_boq.WriteLine("dist2={0}", dist2);
                sw_boq.WriteLine("dist3={0}", dist3);
                sw_boq.WriteLine("dist4={0}", dist4);
                sw_boq.WriteLine("bno1={0}", bno1);
                sw_boq.WriteLine("bno2={0}", bno2);
                sw_boq.WriteLine("bno3={0}", bno3);
                sw_boq.WriteLine("bno4={0}", bno4);
                sw_boq.WriteLine("bno5={0}", bno5);
                sw_boq.WriteLine("bwt1={0:f4}", bwt1);
                sw_boq.WriteLine("bwt2={0:f4}", bwt2);
                sw_boq.WriteLine("bwt3={0:f4}", bwt3);
                sw_boq.WriteLine("bwt4={0:f4}", bwt4);
                sw_boq.WriteLine("bwt5={0:f4}", bwt5);
                sw_boq.WriteLine("depth={0}", depth);
                sw_boq.WriteLine("shape1={0}", shape1);
                sw_boq.WriteLine("shape2={0}", shape2);
                sw_boq.WriteLine("shape3={0}", shape3);
                sw_boq.WriteLine("shape4={0}", shape4);
                sw_boq.WriteLine("FINISH");

                #endregion
            }
            catch (Exception ex)
            {
            }
            finally
            {
                sw.Flush();
                sw.Close();
                sw_design.Flush();
                sw_design.Close();
                sw_boq.Flush();
                sw_boq.Close();
            }
        }
Ejemplo n.º 2
0
        public void Calculate_Program()
        {
            //rep_file_name = Path.Combine(user_path, "DESIGN_OF_ONE_RCC_SLAB.TXT");
            string view_file = Path.Combine(system_path, "DESIGN.FIL");
            string boq_file  = Path.Combine(system_path, "BoQ.FIL");

            StreamWriter sw_view = new StreamWriter(new FileStream(view_file, FileMode.Create));
            StreamWriter sw_boq  = new StreamWriter(new FileStream(boq_file, FileMode.Create));
            StreamWriter sw      = new StreamWriter(new FileStream(rep_file_name, FileMode.Create));

            try
            {
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("\t\t***********************************************");
                sw.WriteLine("\t\t*            ASTRA Pro Release 21             *");
                sw.WriteLine("\t\t*        TechSOFT Engineering Services        *");
                sw.WriteLine("\t\t*                                             *");
                sw.WriteLine("\t\t*           DESIGN OF SINGLE SPAN             *");
                sw.WriteLine("\t\t*  ONE WAY RCC SLAB BY LIMIT STATE METHOD     *");
                sw.WriteLine("\t\t***********************************************");
                sw.WriteLine("\t\t-----------------------------------------------");
                sw.WriteLine("\t\tTHIS RESULT CREATED ON " + System.DateTime.Now.ToString("dd.MM.yyyy  AT HH:mm:ss") + " ");
                sw.WriteLine("\t\t----------------------------------------------");
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");

                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("USER'S DATA");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine(" Length Perpendicular to Span [L] = {0:f2} m", L);
                sw.WriteLine(" Design Span for Slab [B]= {0:f2} m", B);
                sw.WriteLine(" Super imposed / Live Load [LL]= {0:f2} kN/sq.m", LL);
                sw.WriteLine(" Slab Load = {0:f2} kN/sq.m", Slab_load);
                sw.WriteLine(" Concrete Grade [f_ck] = M {0} N/sq.mm", sigma_ck);
                sw.WriteLine(" steel Grade [f_y] = Fe {0} N/sq.mm", sigma_y);
                sw.WriteLine(" Diameter of Main Reinforcement [d1] = {0} mm", d1);
                sw.WriteLine(" Diameter of Distribution Reinforcement [d2] = {0} mm", d2);
                sw.WriteLine(" Clear Cover [h1] = {0} mm", h1);
                sw.WriteLine(" End Cover [h2] = {0} mm", h2);
                sw.WriteLine(" Provide Distribution Reinforcement [ads] = {0} %", ads);
                sw.WriteLine(" Shear Strength of Concrete as % of Steel = {0} %", tc);
                sw.WriteLine(" α = {0}", alpha);
                sw.WriteLine(" β = {0}", beta);
                sw.WriteLine(" γ = {0}", gamma);
                sw.WriteLine(" δ = {0}", delta);
                sw.WriteLine(" λ = {0}", lamda);
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("DESIGN CALCULATIONS");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("■ STEP 1: Calculations for Overall and effective depth");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine("  Since length of the slab is more than twice the width, it is a");
                sw.WriteLine("ONE WAY RCC SLAB. Load will be transferred to the supports along");
                sw.WriteLine("the shorter span.");
                sw.WriteLine();
                sw.WriteLine("Consider a 100 cm wide strip of the slab parallel to its shorter span.");
                sw.WriteLine();

                double d, lowest_Span;

                double val1, val2;
                lowest_Span = (L > B) ? B * 1000 : L * 1000;
                d           = (lowest_Span / (alpha * beta * gamma * delta * lamda));

                sw.WriteLine("Minimum depth of slab");
                sw.WriteLine("    d = L /(α * β * γ * δ * λ)");
                sw.WriteLine();
                sw.WriteLine("Let α = {0}, β = {1}, γ = {2}, δ = {3} and λ = {4}", alpha, beta, gamma, delta, lamda);
                sw.WriteLine();
                sw.WriteLine("So, d = {0}/{1} = {2} mm", lowest_Span.ToString("0.0"), (alpha * beta * gamma * delta * lamda).ToString("0.00"), d.ToString("0.00"));

                double D = (d + h1);
                sw.WriteLine();
                sw.WriteLine("Let us adopt overall depth   D = {0} mm.", D.ToString("0.00"));
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("■ STEP 2: Calculations for Design Load, Moment and Shear");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                double deadLoad_slab = (D / 1000) * 1.0 * Slab_load;
                sw.WriteLine("Dead Load of slab = {0} * {1} * {2} = {3} kN/m.", (D / 1000).ToString("0.00"), 1.0, Slab_load, deadLoad_slab.ToString("0.00"));
                sw.WriteLine();
                sw.WriteLine("Superimposed load = {0} * 1 = {1} kN/m.", ll.ToString("0.00"), ll.ToString("0.00"));
                double totalLoad = deadLoad_slab + ll;
                sw.WriteLine();
                sw.WriteLine("Total Load  = {0:f2} + {1:f2} = {2:f2} kN/m.", deadLoad_slab, LL, totalLoad);


                sw.WriteLine();
                double factoredLoad, loadFactor;
                loadFactor   = 1.5d;
                factoredLoad = totalLoad * loadFactor;
                sw.WriteLine("Factored load if the load factor is {0}", loadFactor.ToString("0.00"));
                sw.WriteLine();
                sw.WriteLine("                               = {0} * {1} = {2} kN/m", loadFactor.ToString("0.00"), totalLoad.ToString("0.00"), factoredLoad.ToString("0.00"));
                sw.WriteLine();
                sw.WriteLine("Maximum BM at centre of shorter span");
                sw.WriteLine("                       = (Wu * l * l) / 8");
                sw.WriteLine();
                sw.WriteLine("Assume steel consist of {0} mm bars with {1} mm clear cover.", d1, h1);

                double half_depth = d1 / 2;
                double eff_depth  = D - h1 - half_depth;
                sw.WriteLine();
                sw.WriteLine("Effective depth = {0} - {1} - {2} = {3} mm", D.ToString("0.0"), h1.ToString("0.0"), half_depth.ToString("0.00"), eff_depth.ToString("0.00"));

                lowest_Span = lowest_Span / 1000;

                double eff_span = (lowest_Span) + (eff_depth / 1000);
                sw.WriteLine();
                sw.WriteLine("Effective Span of Slab = {0} + d = {0} + {1} = {2} m", (lowest_Span).ToString("0.00"), (eff_depth / 1000).ToString("0.00"), eff_span.ToString("0.00"));
                double BM = factoredLoad * eff_span * eff_span / 8;
                //sw.WriteLine("So, BM = ({0} * {1} * {1}) / 8 = {2} kNm", factoredLoad.ToString("0.00"), eff_span.ToString("0.00"), BM.ToString("0.00"));
                sw.WriteLine();
                sw.WriteLine("So, BM = M = ({0} * {1} * {1}) / 8 = {2} kNm", factoredLoad.ToString("0.00"), eff_span.ToString("0.00"), BM.ToString("0.00"));
                sw.WriteLine();
                sw.WriteLine("Max shear force = Vu = (Wn * lc) / 2");
                double max_shear_force = (factoredLoad * lowest_Span / 2.0);
                sw.WriteLine("                = {0} * {1}/2 = {2} kN = {3} N", factoredLoad.ToString("0.00"), lowest_Span.ToString("0.0"), max_shear_force.ToString("0.0"), (max_shear_force * 1000).ToString("0.00"));

                max_shear_force *= 1000;


                sw.WriteLine();
                sw.WriteLine("Depth of the slab is given by");
                sw.WriteLine("              BM = 0.138 * σ_ck * b * d* d");

                double M = (BM * 10E+5);
                d = ((BM * 10E+5) / (0.138 * sigma_ck * 1000));
                d = Math.Sqrt(d);
                sw.WriteLine();
                sw.WriteLine("or      d = √(({0} * 10E+5)/(0.138 * {1} * 1000)) = {2} mm", BM.ToString("0.00"), sigma_ck.ToString("0.00"), d.ToString("0.00"));

                d  = (int)(d / 10);
                d += 2;
                d *= 10;
                sw.WriteLine();
                sw.WriteLine("Adopt effective depth d = {0} mm and over all depth", d.ToString("0.0"));
                sw.WriteLine("                      D = {0} mm", eff_depth);


                sw.WriteLine();
                sw.WriteLine("Adopt of tension steel is given by ");
                sw.WriteLine("             M   = 0.87 * σ_y * A_t( d - ((σ_y * A_t)/(σ_ck * b))");

                double a, b, c, At;

                a = (0.87 * sigma_y * sigma_y) / (sigma_ck * 1000);
                b = 0.87 * sigma_y * d;
                c = M;
                double b_ac = (b * b) - 4 * a * c;
                At = (b) - Math.Sqrt(Math.Abs(b_ac));
                At = At / (2 * a);

                At  = (int)At / 10;
                At += 1;
                At  = At * 10;
                sw.WriteLine();
                sw.WriteLine("   {0} * 10E+5 = 0.87 * {1} * At * ({2} - {3} * At / ({4} * 1000))",
                             BM.ToString("0.00"), sigma_y.ToString("0.00"), d.ToString("0.0"), sigma_y.ToString("0.00"), sigma_ck.ToString("0.00"));
                sw.WriteLine();
                sw.WriteLine("or          At   = {0} sq.mm", At);
                sw.WriteLine();
                sw.WriteLine("Use {0} mm bars @ {1} mm c/c giving total area ", d1, d.ToString("0.0"));

                double est_value = (Math.PI * (d1 * d1) / 4) * (1000 / d);

                val1 = est_value;
                val2 = At;
                sw.WriteLine();
                if (val1 > val2)
                {
                    sw.WriteLine("                         = {0} sq.mm. > {1} sq.mm      OK", est_value.ToString("0"), At.ToString("0"));
                }
                else
                {
                    sw.WriteLine("                         = {0} sq.mm. < {1} sq.mm      NOT OK", est_value.ToString("0"), At.ToString("0"));
                }
                double n_rod = At / ((Math.PI * (d1 * d1) / 4));

                n_rod  = (int)n_rod;
                n_rod += 1;

                sw.WriteLine();
                sw.WriteLine("    Bend alternate bars at L/{0} from the face of support where ", n_rod);
                sw.WriteLine("moment reduces to less than half its maximum value. Temperature ");
                sw.WriteLine("reinforcement equal to {0}% of the gross concrete area will be", ads);
                sw.WriteLine("provided in the longitudinal direction.");
                sw.WriteLine();

                double dirArea = (ads / 100) * 1000 * eff_depth;
                sw.WriteLine("       = {0:f4} * 1000 * {1:f2}", (Ads / 100), eff_depth);
                sw.WriteLine("       = {0:f3} sq.mm.", dirArea);
                sw.WriteLine();
                sw.WriteLine("Use {0} mm MS bars @ 100 mm c/c giving total area ", d2.ToString("0"));
                sw.WriteLine();

                double a_st = Math.PI * d2 * d2 / 4;

                val1 = a_st * 10;
                val2 = dirArea;
                if (val1 > val2)
                {
                    sw.WriteLine("   = {0} * (1000/100) = {1} sq.mm.  > {2} sq.mm            OK", a_st.ToString("0.00"), (a_st * 10).ToString("0.00"), dirArea.ToString("0.0"));
                }
                else
                {
                    sw.WriteLine("   = {0} * (1000/100) = {1} sq.mm.  < {2} sq.mm           NOT OK", a_st.ToString("0.00"), (a_st * 10).ToString("0.00"), dirArea.ToString("0.0"));
                }


                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("■ STEP 3: Check for Shear");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine("Percent tension steel = (100 * At)/ (b * d)");
                a_st = Math.PI * d1 * d1 / 4;



                double percent = (100 * (a_st * (1000 / 300))) / (1000 * d);
                sw.WriteLine("                      = (100 * ({0} * (1000/300)) / (1000 * {1}) = {2}%", a_st.ToString("0.0"), d.ToString("0.0"), percent.ToString("0.00"));
                sw.WriteLine();
                sw.WriteLine("Shear strength of concrete for {0}% steel", percent.ToString("0.00"));
                ShearValue sh    = new ShearValue();
                double     tau_c = 0.0;

                tau_c = sh.Get_M15(percent);

                sw.WriteLine("     τ_c = {0} N/sq.mm.", tau_c);


                double k       = kVal_Table.Get_KValue(eff_depth);
                double tc_dash = k * tau_c;



                //tau_c = tau_c;
                sw.WriteLine();
                sw.WriteLine("For {0} mm thick slab, k = {1}", eff_depth.ToString("0.00"), k.ToString("0.00"));
                sw.WriteLine();
                sw.WriteLine("         So, τ_c` = k * Tc = {0} * {1} = {2} N/sq.mm", k.ToString("0.00"), tau_c.ToString("0.00"), tc_dash.ToString("0.00"));
                double Vu  = max_shear_force;
                double t_v = Vu / (1000 * d);

                sw.WriteLine();
                sw.WriteLine("Nominal shear stress Tv = Vu / b * d = {0}/(1000 * {1}) = {2} N/sq.mm", Vu, d, t_v.ToString("0.00"));

                sw.WriteLine();
                sw.WriteLine("The Slab is safe in shear.");
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("■ STEP 4: Check for development length");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();
                sw.WriteLine("Moment of resistance offered by {0} mm bars @ {1} mm c/c", d1, d * 2);
                sw.WriteLine();
                sw.WriteLine("M1 = 0.87 * σ_y * At * (d - (σ_y * At / σ_ck * b))");
                sw.WriteLine("   = 0.87 * {0:f2} * {1:f2} * (1000/300) * ", sigma_y, a_st);
                sw.WriteLine("     ({0:f2} - ({1:f2} * {2:f2} * (1000 / 300)) / {3:f2} * 1000)",
                             d, sigma_y, a_st, Sigma_ck);

                double M1 = 0.87 * sigma_y * a_st * (1000.0 / 300.0) * (d - (sigma_y * a_st * (1000.0d / 300.0d) / (sigma_ck * 1000.0)));



                sw.WriteLine();
                sw.WriteLine("   = {0:F2} N mm", M1);
                sw.WriteLine("Vu = {0:F2} N", Vu);
                sw.WriteLine();
                sw.WriteLine("Let us assume anchorage length Lo = 0");
                sw.WriteLine("                  Ld <= 1.3 * (M1/Vu)");
                sw.WriteLine("                  56φ <= 1.3 * ({0}/{1})", M1.ToString("0.00"), Vu.ToString("0.00"));

                double phi = (1.3 * (M1 / Vu) / 56.0);


                val1 = phi;
                val2 = d1;
                sw.WriteLine();
                if (val1 > val2)
                {
                    sw.WriteLine("                  φ < {0}", phi.ToString("0.00"));
                    sw.WriteLine("We have provided φ = {0} mm, So OK.", d1.ToString("0.00"));
                }
                else
                {
                    sw.WriteLine("                  φ > {0}", phi.ToString("0.00"));
                    sw.WriteLine("We have provided φ = {0} mm, So NOT OK.", d1.ToString("0.00"));
                }

                sw.WriteLine();
                sw.WriteLine("  The Code requires that bars must be carried");
                sw.WriteLine("into the supports by atleast Ld / 3 = 190 mm");
                sw.WriteLine();
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine("■ STEP 5: Check for deflection");
                sw.WriteLine("------------------------------------------------------------");
                sw.WriteLine();

                sw.WriteLine("Percent tension steel at midspan");
                sw.WriteLine("               = (100 * As) / (b * d)");
                sw.WriteLine("               = (100 * {0:F2} * 1000 / 150) / (1000 * 150)", a_st);
                double pps = (100 * a_st * 1000 / 150) / (1000 * 150);
                sw.WriteLine("               = {0}%", pps.ToString("0.00"));

                double gama = ModificationFactor.GetGamma(pps);
                sw.WriteLine();
                sw.WriteLine("             γ = {0}%", gama.ToString("0.00"));
                //sw.WriteLine(" σσσγγβαδλ■  √");
                sw.WriteLine();
                sw.WriteLine("  β = {0}, δ = {1} and λ = {2}", beta, delta, lamda);

                // Constant 20
                sw.WriteLine();
                sw.WriteLine("Allowable L/d = 20 * {0} = {1}", gama.ToString("0.00"), (20 * gama).ToString("0.00"));

                val1 = (eff_span * 1000) / d;
                val2 = (20 * gama);

                sw.WriteLine();
                if (val1 < val2)
                {
                    sw.WriteLine("Actual L/d = {0} / {1} = {2} < {3}       OK", (eff_span * 1000).ToString("0.00"), d, (eff_span * 1000 / d).ToString("0.00"), (20 * gama).ToString("0.00"));
                }
                else
                {
                    sw.WriteLine("Actual L/d = {0} / {1} = {2} > {3}       NOT OK", (eff_span * 1000).ToString("0.00"), d, (eff_span * 1000 / d).ToString("0.00"), (20 * gama).ToString("0.00"));
                }
                sw.WriteLine();
                sw.WriteLine("---------------------------------------------------------------------------");
                sw.WriteLine("---------------------       END OF REPORT        --------------------------");
                sw.WriteLine("---------------------------------------------------------------------------");

                double fd = 500;

                sw_view.WriteLine("SLAB DESIGN 02");
                sw_view.WriteLine("L = {0}", lowest_Span * 1000);
                sw_view.WriteLine("D = {0}", eff_depth);
                sw_view.WriteLine("b1 = {0}", d);
                sw_view.WriteLine("b2 = 100");
                sw_view.WriteLine("h1 = {0}", h1);
                sw_view.WriteLine("h2 = {0}", h2);
                sw_view.WriteLine("d1 = {0}", d1);
                sw_view.WriteLine("d2 = {0}", d2);
                sw_view.WriteLine("fd = {0:f0}", fd);
                sw_view.WriteLine("w1 = {0}", w1);
                sw_view.WriteLine("w2 = {0}", w1);
                sw_view.WriteLine("END");

                //BoQ CODE
                //S_No 1.,2.,3.
                //Member 3.5 8.0
                //Bar_Mark 01,02,03
                //Bar_CODE T10_B1,T6_B2,T10_B1
                //Bar_Grade Fe415,Fe415,Fe415
                //Bar_Dia 10,6,10
                //Bar_Length 3470,7970,3470
                //Bar_Nos 56,35,56
                //Bar_Weight 0.330,0.220,0.330
                //Bar_Shape 3019,120,1209
                //Bar_Shape 1209,120,3019
                //END BoQ

                int    main_bar_nos, dist_bar_nos;
                double wt_main_bar, wt_dist_bar;


                sw_boq.WriteLine("BoQ Code");
                sw_boq.WriteLine("S_No {0}.,{1}.,{2}.", 1, 2, 3);
                sw_boq.WriteLine("Member {0} {1}", B.ToString("0.0"), L.ToString("0.0"));
                sw_boq.WriteLine("Bar_Spacing {0},{1}", d.ToString("0.0"), "100.0");
                sw_boq.WriteLine("Bar_Mark 01,02,03");
                sw_boq.WriteLine("Bar_Code T{0}_B1,T{1}_B2,T{0}_B1", d1, d2);
                sw_boq.WriteLine("BAR_Grade Fe{0},Fe{0},Fe{0}", sigma_y);
                sw_boq.WriteLine("Bar_Dia {0},{1},{0}", d1, d2);
                sw_boq.WriteLine("Bar_Length {0},{1},{0}", (B * 1000 - 2 * h1), (L * 1000 - 2 * h1));

                main_bar_nos = (int)(L * 1000 / d);
                dist_bar_nos = (int)(B * 1000 / 100.0d);

                wt_main_bar  = 0.00616 * d1 * d1 * (B)*main_bar_nos;
                wt_dist_bar  = 0.00616 * d2 * d2 * (L)*dist_bar_nos;
                wt_main_bar /= 1000.0d;
                wt_dist_bar /= 1000.0d;
                //sw_boq.WriteLine("Bar_Nos {0},{1},{0}", "56", "35", "56");
                sw_boq.WriteLine("Bar_Nos {0},{1},{0}", main_bar_nos, dist_bar_nos);
                sw_boq.WriteLine("Bar_Weight {0},{1},{0}", wt_main_bar.ToString("0.000"), wt_dist_bar.ToString("0.000"));
                //sw_boq.WriteLine("Bar_Weight {0},{1},{0}", "0.330", "0.220");

                double k_d = eff_depth - (2 * ((h1 + (d1 / 2))));

                double aaa = Math.Sqrt((k_d * k_d) + (k_d * k_d));

                double tot_len = (lowest_Span * 1000.0) + 2 * w1 - 2 * h1;

                double sh1 = tot_len - (fd + (w1 - h1));
                double sh2 = (fd + (w1 - h1)) - eff_depth;



                sw_boq.WriteLine("Bar_Shape {0},{1:f0},{2}", sh1, aaa, sh2);
                //sw_boq.WriteLine("Bar_Shape {0},{1:f0},{2}", 3140, aaa, 1290);
                sw_boq.WriteLine("Bar_Shape {0}", ((L * 1000.0) - (2 * h1)));

                sw_boq.WriteLine("Bar_Shape {0},{1},{2}", sh2, aaa, sh1);
                sw_boq.WriteLine("END BoQ");
            }
            catch (Exception exx)
            {
            }
            finally
            {
                GC.Collect();
                sw.Flush();
                sw.Close();
                sw_view.Flush();
                sw_view.Close();
                sw_boq.Flush();
                sw_boq.Close();
            }
        }