public void DGExample5_1ReturnsWarpingStressAtMidspan()
        {

            SetAiscDG9Example5_1Parameters();
            z = 0.5 * L;
            TorsionalFunctionFactory tf = new TorsionalFunctionFactory();
            ITorsionalFunction function = tf.GetTorsionalFunction(TorsionalFunctionCase.Case3, E, G, J, L, z, T_u, C_w, t, alpha);
            double theta_2der = function.Get_theta_2();
            SectionStressAnalysis st = new SectionStressAnalysis();
            double sigma = st.GetNormalStressDueToWarpingOpenSection(E, W_no, theta_2der);
            double refValue = 28.0;
            double actualTolerance = EvaluateActualTolerance(sigma, refValue);
            Assert.LessOrEqual(actualTolerance, tolerance);

        }
        public static Dictionary<string, object> ShearStressDueToWarping(double E, double S_ws,double t_el, double theta_3der)
        {
            //Default values
            double tau_w = 0;


            //Calculation logic:
            SectionStressAnalysis analysis = new SectionStressAnalysis();
            tau_w = analysis.GetShearStressDueToWarpingOpenSecton(E, S_ws, t_el, theta_3der);

            return new Dictionary<string, object>
            {
                { "tau_w", tau_w }
 
            };
        }
        public static Dictionary<string, object> NormalStressDueToAxialLoad(double P,double A)
        {
            //Default values
            double sigma_a = 0;


            //Calculation logic:
            SectionStressAnalysis analysis = new SectionStressAnalysis();
            sigma_a = analysis.GetNormalStressDueToAxialLoad(P, A);

            return new Dictionary<string, object>
            {
                { "sigma_a", sigma_a }
 
            };
        }
        public static Dictionary<string, object> NormalStressDueToWarping(double E, double W_ns,double theta_2der)
        {
            //Default values
            double sigma_ws = 0;


            //Calculation logic:
            SectionStressAnalysis analysis = new SectionStressAnalysis();
            sigma_ws = analysis.GetNormalStressDueToWarpingOpenSection(E, W_ns, theta_2der);
            
            
            return new Dictionary<string, object>
            {
                { "sigma_ws", sigma_ws }
 
            };
        }
        public static Dictionary<string, object> NormalStressDueToBending(double M,double y,double I)
        {
            //Default values
            double sigma_b = 0;


            //Calculation logic:
            SectionStressAnalysis analysis = new SectionStressAnalysis();
            sigma_b = analysis.GetNormalStressDueToBending(M, y, I);


            return new Dictionary<string, object>
            {
                { "sigma_b", sigma_b }
 
            };
        }