normmat_new.cpp

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

                                  k)*\[CurlyEpsilon]n) - 
                      h*(4*(-1 + k)*\[CurlyEpsilon]c + (6 - 
                                  9*k)*\[CurlyEpsilon]c1 + 
                            4*(-1 + 
                                  k)*\[CurlyEpsilon]n)*(\[CurlyEpsilon]c - (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*ch + As1*d)*E0c*k*
                                  Power(\[CurlyEpsilon]c1, 2)))))/(k*
                  Power(\[CurlyEpsilon]c1, 2)*
                  Power(\[CurlyEpsilon]c - 
                      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)))))))), 2))), mhr2 = 1000000000)

Null

ac1 = ((12*Sqrt(3)*
              E0c*((-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*
              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)) - 4*\[CurlyEpsilon]c1*((6*(-ch + d)*E0c)/
                ch + (3*
                    E0c*(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)))*