normmat_new.cpp

Finite element program for mechanical problem. It can solve various problem in solid problem

                      Power(\[CurlyEpsilon]c1, 
                        2)*(fy - 
                          E0c*\[CurlyEpsilon]n)*(ch*fy + (-ch + d)*
                            E0c*\[CurlyEpsilon]n) + 
                    d*(6*fy*k*N0*Power(\[CurlyEpsilon]c1, 2) - 
                          6*As1*fy*k*
                            Power(\[CurlyEpsilon]c1, 
                              2)*(fy - E0c*\[CurlyEpsilon]n) + 
                          E0c*(-6*k*N0*
                                  Power(\[CurlyEpsilon]c1, 
                                    2)*\[CurlyEpsilon]n + 
                                b*d*fc*(-4*(-1 + k)*
                                        Power(\[CurlyEpsilon]c1, 3) + 
                                      3*(-2 + 3*k)*\[CurlyEpsilon]c1*
                                        Power(\[CurlyEpsilon]n, 2) - 
                                      4*(-1 + k)*
                                        Power(\[CurlyEpsilon]n, 3)))))) < 
        0, \[CurlyEpsilon]ct = -1, \[CurlyEpsilon]ct = (-3*E0c*k*
              Power(\[CurlyEpsilon]c1, 
                2)*(-2*As2*ch*fy - As1*d*fy + As2*d*fy + d*N0 + 
                  2*As2*(ch - d)*E0c*\[CurlyEpsilon]n) + 
            3*b*Power(d, 2)*E0c*
              fc*(2*(-1 + k)*
                    Power(\[CurlyEpsilon]c1, 
                      2) + (-2 + 3*k)*\[CurlyEpsilon]c1*\[CurlyEpsilon]n - 
                  2*(-1 + k)*Power(\[CurlyEpsilon]n, 2)) + 
            Sqrt(3)*Sqrt(E0c*(3*E0c*
                        Power(k*
                              Power(\[CurlyEpsilon]c1, 
                                2)*(-2*As2*ch*fy - As1*d*fy + As2*d*fy + 
                                  d*N0 + 
                                  2*As2*(ch - d)*E0c*\[CurlyEpsilon]n) + 
                            b*Power(d, 2)*
                              fc*(-2*(-1 + k)*
                                    Power(\[CurlyEpsilon]c1, 
                                      2) + (2 - 
                                        3*k)*\[CurlyEpsilon]c1*\[CurlyEpsilon]\
n + 2*(-1 + k)*Power(\[CurlyEpsilon]n, 2)), 
                          2) - (2*As2*(-ch + d)*E0c*k*
                              Power(\[CurlyEpsilon]c1, 2) - 
                            b*Power(d, 2)*
                              fc*((-2 + k)*\[CurlyEpsilon]c1 + 
                                  4*(-1 + k)*\[CurlyEpsilon]n))*(-6*As2*k*
                              Power(\[CurlyEpsilon]c1, 
                                2)*(fy - 
                                  E0c*\[CurlyEpsilon]n)*(ch*fy + (-ch + d)*
                                    E0c*\[CurlyEpsilon]n) + 
                            d*(6*fy*k*N0*Power(\[CurlyEpsilon]c1, 2) - 
                                  6*As1*fy*k*
                                    Power(\[CurlyEpsilon]c1, 
                                      2)*(fy - E0c*\[CurlyEpsilon]n) + 
                                  E0c*(-6*k*N0*
                                        Power(\[CurlyEpsilon]c1, 
                                        2)*\[CurlyEpsilon]n + 
                                        b*d*
                                        fc*(-4*(-1 + k)*
                                        Power(\[CurlyEpsilon]c1, 3) + 
                                        3*(-2 + 3*k)*\[CurlyEpsilon]c1*
                                        Power(\[CurlyEpsilon]n, 2) - 
                                        4*(-1 + k)*
                                        Power(\[CurlyEpsilon]n, 3))))))))/(3.*
            E0c*(2*As2*(ch - d)*E0c*k*Power(\[CurlyEpsilon]c1, 2) + 
                b*Power(d, 2)*
                  fc*((-2 + k)*\[CurlyEpsilon]c1 + 
                      4*(-1 + k)*\[CurlyEpsilon]n))))

If(mhr1 == 1000000000 || 
      E0c*(3*E0c*
                Power(As2*k*
                      Power(\[CurlyEpsilon]c1, 
                        2)*(-2*ch*fy + d*fy + 2*ch*E0c*\[CurlyEpsilon]n - 
                          2*d*E0c*\[CurlyEpsilon]n) + 
                    d*(-(As1*fy*k*Power(\[CurlyEpsilon]c1, 2)) + 
                          2*b*d*fc*(Power(\[CurlyEpsilon]ct, 
                                    2) + (\[CurlyEpsilon]c1 - \
\[CurlyEpsilon]n)*\[CurlyEpsilon]n) + 
                          k*(N0*Power(\[CurlyEpsilon]c1, 2) + 
                                b*d*fc*(-2*Power(\[CurlyEpsilon]ct, 2) - 
                                      3*\[CurlyEpsilon]c1*\[CurlyEpsilon]n + 
                                      2*Power(\[CurlyEpsilon]n, 2)))), 
                  2) - (2*As2*(-ch + d)*E0c*k*Power(\[CurlyEpsilon]c1, 2) + 
                    b*Power(d, 2)*
                      fc*((-2 + 3*k)*\[CurlyEpsilon]c1 - 
                          4*(-1 + 
                                k)*(\[CurlyEpsilon]ct + \
\[CurlyEpsilon]n)))*(-6*As2*k*
                      Power(\[CurlyEpsilon]c1, 
                        2)*(fy - 
                          E0c*\[CurlyEpsilon]n)*(ch*fy + (-ch + d)*
                            E0c*\[CurlyEpsilon]n) + 
                    d*(6*fy*k*N0*Power(\[CurlyEpsilon]c1, 2) - 
                          6*As1*fy*k*
                            Power(\[CurlyEpsilon]c1, 
                              2)*(fy - E0c*\[CurlyEpsilon]n) + 
                          E0c*(-6*k*N0*
                                  Power(\[CurlyEpsilon]c1, 
                                    2)*\[CurlyEpsilon]n + 
                                b*d*fc*(-4*(-1 + k)*
                                        Power(\[CurlyEpsilon]ct, 3) + 
                                      3*(-2 + 3*k)*\[CurlyEpsilon]c1*
                                        Power(\[CurlyEpsilon]n, 2) - 
                                      4*(-1 + k)*
                                        Power(\[CurlyEpsilon]n, 3)))))) < 
        0 || \[CurlyEpsilon]ct == -1, \[CurlyEpsilon]c = -1, \[CurlyEpsilon]c \
= (-3*As2*E0c*k*
              Power(\[CurlyEpsilon]c1, 
                2)*(-2*ch*fy + d*fy + 2*ch*E0c*\[CurlyEpsilon]n - 
                  2*d*E0c*\[CurlyEpsilon]n) + 
            3*d*E0c*(k*(As1*fy - N0)*Power(\[CurlyEpsilon]c1, 2) - 
                  2*b*d*fc*Power(\[CurlyEpsilon]ct, 2) + 
                  2*b*d*fc*k*Power(\[CurlyEpsilon]ct, 2) + 
                  b*d*fc*(-2 + 3*k)*\[CurlyEpsilon]c1*\[CurlyEpsilon]n - 
                  2*b*d*fc*(-1 + k)*Power(\[CurlyEpsilon]n, 2)) + 
            Sqrt(3)*Sqrt(E0c*(3*E0c*
                        Power(As2*k*
                              Power(\[CurlyEpsilon]c1, 
                                2)*(-2*ch*fy + d*fy + 
                                  2*ch*E0c*\[CurlyEpsilon]n - 
                                  2*d*E0c*\[CurlyEpsilon]n) + 
                            d*(-As1*fy*k*Power(\[CurlyEpsilon]c1, 2) + 
                                  2*b*d*
                                    fc*(Power(\[CurlyEpsilon]ct, 
                                        2) + (\[CurlyEpsilon]c1 - \
\[CurlyEpsilon]n)*\[CurlyEpsilon]n) + 
                                  k*(N0*Power(\[CurlyEpsilon]c1, 2) + 
                                        b*d*fc*(-2*
                                        Power(\[CurlyEpsilon]ct, 2) - 
                                        3*\[CurlyEpsilon]c1*\[CurlyEpsilon]n \
+ 2*Power(\[CurlyEpsilon]n, 2)))), 
                          2) - (2*As2*(-ch + d)*E0c*k*
                              Power(\[CurlyEpsilon]c1, 2) + 
                            b*Power(d, 2)*
                              fc*((-2 + 3*k)*\[CurlyEpsilon]c1 - 
                                  4*(-1 + 
                                        k)*(\[CurlyEpsilon]ct + \
\[CurlyEpsilon]n)))*(-6*As2*k*
                              Power(\[CurlyEpsilon]c1, 
                                2)*(fy - 
                                  E0c*\[CurlyEpsilon]n)*(ch*fy + (-ch + d)*
                                    E0c*\[CurlyEpsilon]n) + 
                            d*(6*fy*k*N0*Power(\[CurlyEpsilon]c1, 2) - 
                                  6*As1*fy*k*
                                    Power(\[CurlyEpsilon]c1, 
                                      2)*(fy - E0c*\[CurlyEpsilon]n) + 
                                  E0c*(-6*k*N0*
                                        Power(\[CurlyEpsilon]c1, 
                                        2)*\[CurlyEpsilon]n + 
                                        b*d*fc*(-4*(-1 + k)*
                                        Power(\[CurlyEpsilon]ct, 3) + 
                                        3*(-2 + 3*k)*\[CurlyEpsilon]c1*
                                        Power(\[CurlyEpsilon]n, 2) - 
                                        4*(-1 + k)*
                                        Power(\[CurlyEpsilon]n, 3))))))))/(3.*
            E0c*(2*As2*(ch - d)*E0c*k*Power(\[CurlyEpsilon]c1, 2) + 
                b*Power(d, 2)*
                  fc*((2 - 3*k)*\[CurlyEpsilon]c1 + 
                      4*(-1 + k)*(\[CurlyEpsilon]ct + \[CurlyEpsilon]n)))))

\[CurlyEpsilon]chr2 = \[CurlyEpsilon]c

If(\[CurlyEpsilon]c < 0 && \[CurlyEpsilon]c != -1, 
    mhr2 = 0.08333333333333333*(6*As1*
              E0c*(-2*cd + h)*(\[CurlyEpsilon]n + 
                  1/(6.*(As2*ch + As1*d)*E0c*k*
                          Power(\[CurlyEpsilon]c1, 2))*(6*As2*ch*E0c*
                          k*\[CurlyEpsilon]c*Power(\[CurlyEpsilon]c1, 2) + 
                        3*d*k*Power(\[CurlyEpsilon]c1, 
                            2)*(N0 + 
                              As1*E0c*(\[CurlyEpsilon]c - \[CurlyEpsilon]n) - 
                              As2*E0c*(\[CurlyEpsilon]c + \[CurlyEpsilon]n)) +
                         Sqrt(3)*
                          Sqrt(k*
                              Power(\[CurlyEpsilon]c1, 
                                2)*(3*k*Power(\[CurlyEpsilon]c1, 2)*
                                    Power(-2*As2*ch*E0c*\[CurlyEpsilon]c + 
                                        d*(-N0 + 
                                        E0c*(As1*(-\[CurlyEpsilon]c + \
\[CurlyEpsilon]n) + As2*(\[CurlyEpsilon]c + \[CurlyEpsilon]n))), 2) - 
                                  2*(As2*ch + As1*d)*
                                    E0c*(b*Power(d, 2)*fc*
                                        Power(\[CurlyEpsilon]c + \
\[CurlyEpsilon]n, 2)*(4*(-1 + k)*\[CurlyEpsilon]c + (6 - 
                                        9*k)*\[CurlyEpsilon]c1 + 
                                        4*(-1 + k)*\[CurlyEpsilon]n) + 
                                        6*k*\[CurlyEpsilon]c*
                                        Power(\[CurlyEpsilon]c1, 
                                        2)*(As2*ch*E0c*\[CurlyEpsilon]c + 
                                        d*(N0 - E0c*(As1*\[CurlyEpsilon]n + 
                                        As2*(\[CurlyEpsilon]c + \
\[CurlyEpsilon]n))))))))) - 
            6*As2*E0c*(-2*ch + 
                  h)*(\[CurlyEpsilon]c - (ch*\[CurlyEpsilon]c)/
                    d + \[CurlyEpsilon]n + (ch*(6*As2*ch*E0c*
                              k*\[CurlyEpsilon]c*
                              Power(\[CurlyEpsilon]c1, 2) + 
                            3*d*k*Power(\[CurlyEpsilon]c1, 
                                2)*(N0 + 
                                  As1*E0c*(\[CurlyEpsilon]c - \
\[CurlyEpsilon]n) - As2*E0c*(\[CurlyEpsilon]c + \[CurlyEpsilon]n)) + 
                            Sqrt(3)*
                              Sqrt(k*
                                  Power(\[CurlyEpsilon]c1, 
                                    2)*(3*k*Power(\[CurlyEpsilon]c1, 2)*
                                        Power(-2*As2*ch*E0c*\[CurlyEpsilon]c +
                                         d*(-N0 + 
                                        E0c*(As1*(-\[CurlyEpsilon]c + \
\[CurlyEpsilon]n) + As2*(\[CurlyEpsilon]c + \[CurlyEpsilon]n))), 2) - 
                                      2*(As2*ch + As1*d)*
                                        E0c*(b*Power(d, 2)*fc*
                                        Power(\[CurlyEpsilon]c + \
\[CurlyEpsilon]n, 2)*(4*(-1 + k)*\[CurlyEpsilon]c + (6 - 
                                        9*k)*\[CurlyEpsilon]c1 + 
                                        4*(-1 + k)*\[CurlyEpsilon]n) + 
                                        6*k*\[CurlyEpsilon]c*
                                        Power(\[CurlyEpsilon]c1, 
                                        2)*(As2*ch*E0c*\[CurlyEpsilon]c + 
                                        d*(N0 - E0c*(As1*\[CurlyEpsilon]n + 
                                        As2*(\[CurlyEpsilon]c + \
\[CurlyEpsilon]n)))))))))/(6.*d*(As2*ch + As1*d)*E0c*k*
                        Power(\[CurlyEpsilon]c1, 2))) - (b*d*fc*
                  Power(\[CurlyEpsilon]c + \[CurlyEpsilon]n, 
                    2)*(2*d*(\[CurlyEpsilon]c + \[CurlyEpsilon]n)*((-1 + 
                                  k)*\[CurlyEpsilon]c + (2 - 
                                  3*k)*\[CurlyEpsilon]c1 + (-1 +