normmat_new.cpp

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

                                      Power(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n, 2)))))/((As2*ch + 
                        As1*d)*k*Power(\[CurlyEpsilon]c1, 2))) + (4*(3*As2*ch*
                    d*k*Power(\[CurlyEpsilon]c1, 
                      2)*(4*fy + 
                        E0c*(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n)) + 
                  3*As1*d*k*
                    Power(\[CurlyEpsilon]c1, 
                      2)*(-(ch*
                              E0c*(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n)) + 
                        2*d*(2*fy + E0c*\[CurlyEpsilon]c1 + 
                              2*E0c*\[CurlyEpsilon]n)) + 
                  2*ch*(3*d*k*N0*Power(\[CurlyEpsilon]c1, 2) + 
                        Sqrt(3)*
                          Sqrt(k*
                              Power(\[CurlyEpsilon]c1, 
                                2)*(3*k*Power(\[CurlyEpsilon]c1, 2)*
                                    Power(As2*ch*
                                        E0c*\[CurlyEpsilon]c1 + (d*(2*N0 + 
                                        E0c*(As1*(\[CurlyEpsilon]c1 - 
                                        2*\[CurlyEpsilon]n) - 
                                        As2*(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n))))/2., 2) - 
                                  2*(As2*ch + As1*d)*
                                    E0c*(-(b*Power(d, 2)*fc*
                                        Power(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n, 
                                        2)*(-4*\[CurlyEpsilon]c1 + 
                                        7*k*\[CurlyEpsilon]c1 + 
                                        4*\[CurlyEpsilon]n - 
                                        4*k*\[CurlyEpsilon]n))/4. + 
                                        3*k*Power(\[CurlyEpsilon]c1, 
                                        3)*((As2*ch*E0c*\[CurlyEpsilon]c1)/
                                        2. + 
                                        d*(N0 - (E0c*(As2*\[CurlyEpsilon]c1 + 
                                        2*As1*\[CurlyEpsilon]n + 
                                        2*As2*\[CurlyEpsilon]n))/
                                        2.))))))))/(ch*(As2*ch + As1*d)*k*
              Power(\[CurlyEpsilon]c1, 2)))/(48.*E0c)

ac2 = (-2 + (2*d)/
          ch - (2*Sqrt(3)*((-8*b*(As2*ch + As1*d)*E0c*
                          fc*(\[CurlyEpsilon]c1 - \[CurlyEpsilon]n)*
                          Power(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n, 2) + 
                        k*(12*Power(N0, 2)*Power(\[CurlyEpsilon]c1, 2) - 
                              12*(As1 + As2)*E0c*N0*
                                Power(\[CurlyEpsilon]c1, 
                                  2)*(\[CurlyEpsilon]c1 + 
                                    2*\[CurlyEpsilon]n) + 
                              E0c*(3*Power(As2, 2)*E0c*
                                      Power(\[CurlyEpsilon]c1, 2) + 
                                    2*As2*(3*As1*E0c*
                                        Power(\[CurlyEpsilon]c1, 2) + 
                                        b*ch*
                                        fc*(7*\[CurlyEpsilon]c1 - 
                                        4*\[CurlyEpsilon]n)) + 
                                    As1*(3*As1*E0c*
                                        Power(\[CurlyEpsilon]c1, 2) + 
                                        2*b*d*
                                        fc*(7*\[CurlyEpsilon]c1 - 
                                        4*\[CurlyEpsilon]n)))*
                                Power(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n, 
                                  2)))*(Power(As2, 2)*E0c*k*
                          Power(\[CurlyEpsilon]c1, 2) + 
                        2*As2*(As1*E0c*k*Power(\[CurlyEpsilon]c1, 2) + 
                              b*ch*fc*(k*\[CurlyEpsilon]c1 + 
                                    4*\[CurlyEpsilon]n - 
                                    4*k*\[CurlyEpsilon]n)) + 
                        As1*(As1*E0c*k*Power(\[CurlyEpsilon]c1, 2) + 
                              2*b*d*fc*(k*\[CurlyEpsilon]c1 + 
                                    4*\[CurlyEpsilon]n - 
                                    4*k*\[CurlyEpsilon]n))) - 
                  3*E0c*Power(Power(As2, 2)*E0c*k*
                          Power(\[CurlyEpsilon]c1, 
                            2)*(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n) + 
                        As1*(2*b*d*
                                fc*(\[CurlyEpsilon]c1 + 
                                    2*\[CurlyEpsilon]n)*((-1 + 
                                        2*k)*\[CurlyEpsilon]c1 - 
                                    2*(-1 + k)*\[CurlyEpsilon]n) + 
                              k*Power(\[CurlyEpsilon]c1, 
                                  2)*(-2*N0 + 
                                    As1*E0c*(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n))) + 
                        2*As2*(b*ch*
                                fc*(\[CurlyEpsilon]c1 + 
                                    2*\[CurlyEpsilon]n)*((-1 + 
                                        2*k)*\[CurlyEpsilon]c1 - 
                                    2*(-1 + k)*\[CurlyEpsilon]n) + 
                              k*Power(\[CurlyEpsilon]c1, 
                                  2)*(-N0 + 
                                    As1*E0c*(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n))), 2))*
              Sqrt(k*Power(\[CurlyEpsilon]c1, 
                    2)*(3*k*Power(\[CurlyEpsilon]c1, 2)*
                        Power(As2*ch*
                              E0c*\[CurlyEpsilon]c1 + (d*(2*N0 + 
                                      E0c*(As1*(\[CurlyEpsilon]c1 - 
                                        2*\[CurlyEpsilon]n) - 
                                        As2*(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n))))/2., 2) - 
                      2*(As2*ch + As1*d)*
                        E0c*(-(b*Power(d, 2)*fc*
                                    Power(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n, 
                                      2)*((-4 + 7*k)*\[CurlyEpsilon]c1 - 
                                        4*(-1 + k)*\[CurlyEpsilon]n))/4. + 
                            3*k*Power(\[CurlyEpsilon]c1, 
                                3)*((As2*ch*E0c*\[CurlyEpsilon]c1)/2. + 
                                  d*(N0 - (E0c*(As2*\[CurlyEpsilon]c1 + 
                                        2*As1*\[CurlyEpsilon]n + 
                                        2*As2*\[CurlyEpsilon]n))/
                                        2.))))))/((As2*ch + As1*d)*
              k*\[CurlyEpsilon]c1*
              Power(-8*b*(As2*ch + As1*d)*E0c*
                    fc*(\[CurlyEpsilon]c1 - \[CurlyEpsilon]n)*
                    Power(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n, 2) + 
                  k*(12*Power(N0, 2)*Power(\[CurlyEpsilon]c1, 2) - 
                        12*(As1 + As2)*E0c*N0*
                          Power(\[CurlyEpsilon]c1, 
                            2)*(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n) + 
                        E0c*(3*Power(As2, 2)*E0c*
                                Power(\[CurlyEpsilon]c1, 2) + 
                              2*As2*(3*As1*E0c*Power(\[CurlyEpsilon]c1, 2) + 
                                    b*ch*fc*(7*\[CurlyEpsilon]c1 - 
                                        4*\[CurlyEpsilon]n)) + 
                              As1*(3*As1*E0c*Power(\[CurlyEpsilon]c1, 2) + 
                                    2*b*d*
                                      fc*(7*\[CurlyEpsilon]c1 - 
                                        4*\[CurlyEpsilon]n)))*
                          Power(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n, 2)), 
                2)) + (2*As2*ch*k*Power(\[CurlyEpsilon]c1, 2) + (As1 - As2)*d*
                k*Power(\[CurlyEpsilon]c1, 
                  2) + (2*
                    Sqrt(3)*(Power(As2, 2)*E0c*k*
                          Power(\[CurlyEpsilon]c1, 
                            2)*(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n) + 
                        As1*(2*b*d*
                                fc*(\[CurlyEpsilon]c1 + 
                                    2*\[CurlyEpsilon]n)*((-1 + 
                                        2*k)*\[CurlyEpsilon]c1 - 
                                    2*(-1 + k)*\[CurlyEpsilon]n) + 
                              k*Power(\[CurlyEpsilon]c1, 
                                  2)*(-2*N0 + 
                                    As1*E0c*(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n))) + 
                        2*As2*(b*ch*
                                fc*(\[CurlyEpsilon]c1 + 
                                    2*\[CurlyEpsilon]n)*((-1 + 
                                        2*k)*\[CurlyEpsilon]c1 - 
                                    2*(-1 + k)*\[CurlyEpsilon]n) + 
                              k*Power(\[CurlyEpsilon]c1, 
                                  2)*(-N0 + 
                                    As1*E0c*(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n))))*
                    Sqrt(k*Power(\[CurlyEpsilon]c1, 
                          2)*(3*k*Power(\[CurlyEpsilon]c1, 2)*
                              Power(As2*ch*
                                    E0c*\[CurlyEpsilon]c1 + (d*(2*N0 + 
                                        E0c*(As1*(\[CurlyEpsilon]c1 - 
                                        2*\[CurlyEpsilon]n) - 
                                        As2*(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n))))/2., 2) - 
                            2*(As2*ch + As1*d)*
                              E0c*(-(b*Power(d, 2)*fc*
                                        Power(\[CurlyEpsilon]c1 + 
                                        2*\[CurlyEpsilon]n, 
                                        2)*((-4 + 7*k)*\[CurlyEpsilon]c1 - 
                                        4*(-1 + k)*\[CurlyEpsilon]n))/4. + 
                                  3*k*Power(\[CurlyEpsilon]c1, 
                                      3)*((As2*ch*E0c*\[CurlyEpsilon]c1)/2. + 
                                        d*(N0 - (E0c*(As2*\[CurlyEpsilon]c1 + 
                                        2*As1*\[CurlyEpsilon]n + 
                                        2*As2*\[CurlyEpsilon]n))/2.))))))/(-8*
                  
                    b*(As2*ch + As1*d)*E0c*
                      fc*(\[CurlyEpsilon]c1 - \[CurlyEpsilon]n)*
                      Power(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n, 2) + 
                    k*(12*Power(N0, 2)*Power(\[CurlyEpsilon]c1, 2) - 
                          12*(As1 + As2)*E0c*N0*
                            Power(\[CurlyEpsilon]c1, 
                              2)*(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n) + 
                          E0c*(3*Power(As2, 2)*E0c*
                                  Power(\[CurlyEpsilon]c1, 2) + 
                                2*As2*(3*As1*E0c*Power(\[CurlyEpsilon]c1, 2) +
                                       b*ch*
                                        fc*(7*\[CurlyEpsilon]c1 - 
                                        4*\[CurlyEpsilon]n)) + 
                                As1*(3*As1*E0c*Power(\[CurlyEpsilon]c1, 2) + 
                                      2*b*d*
                                        fc*(7*\[CurlyEpsilon]c1 - 
                                        4*\[CurlyEpsilon]n)))*
                            Power(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n, 
                              2))))/((As2*ch + As1*d)*k*
              Power(\[CurlyEpsilon]c1, 2)))/2.

ac3 = (Sqrt(3)*((-8*b*(As2*ch + As1*d)*E0c*
                    fc*(\[CurlyEpsilon]c1 - \[CurlyEpsilon]n)*
                    Power(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n, 2) + 
                  k*(12*Power(N0, 2)*Power(\[CurlyEpsilon]c1, 2) - 
                        12*(As1 + As2)*E0c*N0*
                          Power(\[CurlyEpsilon]c1, 
                            2)*(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n) + 
                        E0c*(3*Power(As2, 2)*E0c*
                                Power(\[CurlyEpsilon]c1, 2) + 
                              2*As2*(3*As1*E0c*Power(\[CurlyEpsilon]c1, 2) + 
                                    b*ch*fc*(7*\[CurlyEpsilon]c1 - 
                                        4*\[CurlyEpsilon]n)) + 
                              As1*(3*As1*E0c*Power(\[CurlyEpsilon]c1, 2) + 
                                    2*b*d*
                                      fc*(7*\[CurlyEpsilon]c1 - 
                                        4*\[CurlyEpsilon]n)))*
                          Power(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n, 
                            2)))*(Power(As2, 2)*E0c*k*
                    Power(\[CurlyEpsilon]c1, 2) + 
                  2*As2*(As1*E0c*k*Power(\[CurlyEpsilon]c1, 2) + 
                        b*ch*fc*(k*\[CurlyEpsilon]c1 + 4*\[CurlyEpsilon]n - 
                              4*k*\[CurlyEpsilon]n)) + 
                  As1*(As1*E0c*k*Power(\[CurlyEpsilon]c1, 2) + 
                        2*b*d*fc*(k*\[CurlyEpsilon]c1 + 4*\[CurlyEpsilon]n - 
                              4*k*\[CurlyEpsilon]n))) - 
            3*E0c*Power(Power(As2, 2)*E0c*k*
                    Power(\[CurlyEpsilon]c1, 
                      2)*(\[CurlyEpsilon]c1 + 2*\[CurlyEpsilon]n) +