\(\(\(\ \)\( (*\(\(\(\(\(\(\(\(\(\(\(\(\(\(\(\(\(\(\(\(\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--
																	PROGRAMA\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\ DE\ SOLUCOES\ ANALITICAS\ DE\ PLACAS\ 
	 (*\ LAMINADO\ RETANGULAR\ SIMETRICO\ ORTOTROPICO\ SIMPLESMENTE\ APOIADO\ 
			\ \ SOB\ CARGA\ UNIFORME\ qo\ OU\ CARGA\ SENOIDAL . \ 
		\ \ \ TEORIA\ DE\ PRIMEIRA\ \(\(ORDEM\)\(.\)\)\ *) 
	

                                        	
                                        Paulo\ de\ Tarso\ R . \ 
                                        Mendonca, \ \(21/01\)/
                                        1999. 
GRANTE\  - \ \(Grupo\ de\ Analise\ e\ Projeto\ Mecânici\)
\(Departamento\ de\ Engenharia\ Mecanica\  - \ UFSC\)
\(Florianopolis, \ SC\  - \ 
                                        CP\ 476\)
\(www . 
                                        grante . ufsc . br, \ 
                                        mendonca@grante . ufsc . 
                                        br\)
\(\((48)\)\ 3331\) - 9899
	
	INSTRUCOES\ DE\ \(USO : 
	
                                        A\)\)\()\)\)\ Entre\ com\ os\ dados\ na\ celula\ abaixo\ e\ execute\ o\ frontend\ 
		\ \ \((Kernel -> \(Evaluation -> Evaluate\ Notebook\))\) . 
			
	
                                        B\)\()\)\)\ Role\ para\ a\ ultima\ celula\ e\ escolha\ os\ resultados\ que\ deseja\ listados\ ou\ plotados . 
			\ Os\ comandos\ de\ listar\ e\ plotar\ ja\ estao\ prontos . \ 
                                        Basta\ altera - 
                                        los\ para
	\ \ \ \(o\ dado\ desejado\ e\ executar\ a\ celula\ usando\ Shift + 
                                        Enter . 
				
	
                                        C\)\)\()\)\)\ As\ saidas\ disponiveis\ \(sao : \ deslocamento\ wp\ no\ ponto\), \ wmax\ maximo, \ 
                                      momentos\ fletores\ e\ 
		\ \ \(torcores\ Mx\), My, \ Mxy\ e\ tensoes\ em\ cada\ interface, \ 
                                      nas\ direcoes\ x - y\ e\ 1 - 2, \ 
                                      de\ forma
	\ \ \ \(nao\ normalizada\ e\ normalizada . \ 
                                        Os\ nomes\ sao\ os\ \(seguintes : 
			\ \ wp\)\), \ \ \ wmax, \ \ mx, \ \ my, \ \ mxy, \ \ \ \ sigx, \ \ sigy, \ \ tauxy, \ 
                                      tauxz\ \ \ \ tauxy, \ \ \ \ \(\(->\)\(\ \)\(valores\ nao\ normalizados\)\), \ e
	\ \ \ \ wpN, \ \ wmaxN, \ 
                                      mx, \ \ my, \ \ mxy, \ \ sigxN, \ 
                                      sigyN, \ tauxyN, \ tauyzN, \ 
                                      tauxzN\ \ \ \ \ \  -> \ 
                                        valores\ normalizados . 
					Os\ valores\ acima\ de\ tensoes\ cisalhantes\ são\ obtidos\ pelas\ relacoes
	\ \ \ constitutivas . \ 
                                        Valores\ obtidos\ pela\ integracao\ das\ equacoes\ 
	\ \ \ \(diferenciais\ de\ equilibrio\ \(sao : 
			\ \ \ tauyzi\)\), \ tauxzi, \ tauyziN, \ 
                                      tauxziN . 
				\((E'\ ainda\ feita\ analise\ de\ seguranca, \ conforme\ itens\ H, \ 
                                        I\ e\ J\ abaixo)\)
				

                                        D\)\(	\)\()\)\)\ As\ normalizacoes\ sao\ aquelas\ usadas\ por\ \(Pagano : 
			\ \ \((\ sigxN, \ sigyN, \ 
                                        tauxyN)\)\)\  = \ \((sigx, \ \ sigy, \ \ tauxy)\) H\^2/\((qo\ a\^2)\), 
	\ \ \ \ wpN\  = \ 
                                  100\ E2\ w\ H\^3/\((qo\ a\^4)\), 
	\ \ \ \ \((mxN, \ myN, \ mxyN)\)\  = \ \((mxN, \ myN, \ 
                                      mxyN)\)\ H/\((qo\ a\^3)\); 
	\ \ \ \ \(onde\ E2, \ H, \ a\ e\ qo\ sao\ modulo, \ espessura\ do\ laminado, \ 
                                lado\ da\ placa\ e\ carga, \ 
	\ \ \ \ como\ definidos\ na\ celula\ abaixo . 
				
	
                                    E\)\)\()\)\)\ Os\ valores\ sao\ listados\ e\ plotados\ em\ termos\ das\ coordenadas\ z\ de\ cada\ interface . \ 
	\ \ \ As\ laminas\ sao\ contadas\ no\ sentido\ positivo\ de\ z . \ 
                                A\ origem\ e'\ na\ superficie\ media
	\ \ \ \(do\ laminado . \ 
                                  Os\ valores\ normalizados\ sao\ listados\ e\ plotados\ em\ termos\ de\ zN\), \ 
	\ \ \ definido\ como\ zN\  = \ z/H, \ 
                            variando\ de\  - 
                              1/2\ a\ 1/2. 

	
                                F\)\()\)\)\ Observe\ que\ o\ programa\ so\ funciona\ para\ laminados\ de\ laminas\ ortotropicas
			\ \(alinhadas\ aos\ eixos\ x - y\ ou\ laminas\ isotropicas . 
					
	G\)\)\()\)\) 
                      Pode - se\ entrar\ com\ laminas\ obliquas . \ 
                        As\ matrizes\ de\ material\ 
			serao\ calculadas\ corretamente\ mas\ nos\ calculos\ subsequentes\ os\ termos\ de\ 
		\ \(acoplamento\ sao\ ignorados . 
			
	
                          H\)\)\()\)\)\ A\ analise\ de\ resistencia\ e'\ feita\ nas\ superficies\ inferior\ e\ superior\ de\ 
\ 		\(cada\ lamina\), \ pelos\ metodos\ da\ Tensao\ Maxima\ \((TM)\)\ e\ Tsai - 
                Hill\ \((TH)\) . 
			Caso\ nao\ se\ deseje\ esta\ analise, \ 
              ou\ nao\ se\ disponha\ dos\ \ 
	\ 	\(valores\ de\ resistencia\), \ 
              basta\ colocar\ valores\ bastantes\ altos\ para\ 
		\ \(os\ dados\ de\ resistencia . 
			
	
                    I\)\)\()\)\)\ O\ resultado\ da\ analise\ de\ seguranca\ e'\ o\ FATOR\ DE\ SOLICITACAO
	\ \ \ \(f\ definido\), \ 
          dependendo\ do\ criterio\ \(como : 
			\ Tensao\ \(Maxima : \ 
                    f\)\) = \(Maximo {\((sig1\/Xt\ se\ sig1 > 
                        0\ ou\ sig1\/\(\(\ \)\(Xc\)\)\ se\ sig1 < 
                        0)\), \((sig2\/Yt\ se\ sig2 > 
                        0\ ou\ sig2\/Yc\ se\ sig2 < 0)\), 
                    sig12\/S}
		\ \ Tsai - \(Hill : \ \((sig1\/Xt)\)\^2 - \(\(sig1\ sig2\)\/\(Xt\ Yt\)\) \((sig2\/Yt)\)\^2 + \((sig12\/S)\)\^2\) = 
              f\^2\ caso\ sig1 > 0\ e\ sig2 > 0\), \ 
          etc . 
		
	
            J\)\()\)\)\ O\ coeficiente\ de\ seguranca\ no\ ponto\ \((xc, \ 
          yc)\), \ para\ cada\ criterio, \ 
		\ \ cstm, \ csth, \ 
      e'\ obtido\ por\ 1/
          Max[f]\ dentre\ todas\ as\ interfaces\ de\ laminas . 
			
	
	\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(---\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)*) \)\)\)

\(\( (*\ \ DADOS\ *) 
n\  = \ 6\ ; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ Numero\ de\ laminas\ *) 
teta = {0, 0, 90, 90, 0, 0} Pi/180; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ Angulos\ das\ laminas.\ A\ lamina\ 1\ e' inferior, \ 
      com\ cota\ z\ \(negativa.\)\ *) 

    e1 = {\ \ \ \ 172.25, \ \ 172.25, 172.25, \ 172.25, \ \ 172.25, \ \ 
          172.25} 10\^3; 
     (*\ Modulo\ na\ direcao\ 1\ de\ cada\ lamina.\ MPa\ *) 

    e2 = {\ \ \ 6890, \ 6890, 6890, 6890, \ 6890, \ 6890}; \ \ \ \ \ \ \ \ 
     (*\ Modulo\ na\ direcao\ 2\ de\ cada\ lamina\ *) 

    g12 = {3445, 3445, 3445, \ 3445, 3445, 3445}; \ \ \ \ \ \ 
     (*\ Modulo\ na\ direcao\ 12\ de\ cada\ lamina\ *) 

    nu = {\ \ 0.25, \ 0.25, 0.25, \ 0.25, \ 0.25, \ 0.25}; \ \ \ \ \ \ 
     (*\ Coef.Poisson\ na\ direcao\ 12\ de\ cada\ lamina\ *) 

    h = {10, 10, 10, 10, 10, 10\ }/2; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ Espessuras\ das\ laminas.\ mm.\ *) 

    g13 = {3445, 3445, 3445, \ 3445, 3445, 3445}; \ \ \ \ \ \ \ 
     (*\ Modulos\ \(transversais.\)\ \ \ \ \ \ \ \ *) 

    g23 = \ {\ \ 1378, \ 1378, 1378, \ 1378, \ 1378, \ 1378}; 

    e3 = {\ \ \ 6890, \ 6890, 6890, 6890, \ 6890, \ 6890}; 

    kapa\  = \ \ 5/6\ ; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  (*\ fator\ de\ cisalhamento\ 5/6\ *) 
mh\  = \ 1; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ numero\ de\ termos\ na\ serie\ em\ x\ *) 
nh\  = \ 1; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ numero\ de\ termos\ na\ serie\ em\ y\ *) 
qo\  = \ 0.001; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ Carga\ distribuida, \ F/Area.\ Se\ mh = \(nh = 1\), \ q0\ e'

      					\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
      o\ modulo\ de\ carga\ senoidal.\ Se\ mh, nh\ diferentes\ de

      	\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 1, 
      \ carga\ uniforme\ de\ valor\ qo\ *) 
a = 3000; 
b = 3000; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ Lados\ da\ placa.\ a\ na\ direcao\ x\ *) 
xc = \ a/2; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ Ponto\ para\ o\ calculo\ das\ tensoes\ coplanares, \ 
      deslocamento\ e\ \(momentos.\)*) 
yc = 0/2; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ A\ rigem\ dos\ eixos\ e' no\ vertice\ do\ retangulo\ *) 

    Xt = {\ 1000, 1000, 1000, 1000, \ 1000, 1000\ }; \ \ \ \ \ 
     (*\ Xt, \ Xc, \ Yt, \ Yc, \ S\ sao\ os\ valores\ de\ *) 

    Xc = {600, 600\ , 600, 600, 600, 600}; \ \ \ \ \ \ \ \ \ 
     (*\ \ resistencia\ da\ lamina, \ nas\ direcoes\ 1, 2, *) 

    Yt = {\ 30, \ 30, \ 30, 30, 30, \ 30}; \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ a\ tracao\ e\ compressao.\ S\ e'\ cisalhante\ *) 

    Yc = {\ 120, \ 120\ , 120, 120, 120, \ 120}; \ 

    S = {\ 70, \ 70, \ 70, 70, 70, \ 70}; 

     (*\  ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  === 
         ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  == 
		 === 
         ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  === 
         ===  ===  ===  ===  ===  ===  ===  ===  ===  ===  == *) \)\)

\(\( (*\ \ DADOS\ \ \ \ \ \ \ \ reserva\ \ *) 

     (*\ \ n\  = \ 4\ ; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  (*\ Numero\ de\ laminas\ *) 
teta = {0, 90, 90, 0} Pi/180; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ Angulos\ das\ laminas.\ A\ lamina\ 1\ e' inferior, \ 
        com\ cota\ z\ \(negativa.\)\ *) 

      e1 = {172.25, \ 172.25, \ \ 172.25, \ \ 172.25} 10\^3; 
       (*\ Modulo\ na\ direcao\ 1\ de\ cada\ lamina.\ MPa\ *) 

      e2 = {\ 6890, 6890, \ 6890, \ 6890}; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ Modulo\ na\ direcao\ 2\ de\ cada\ lamina\ *) 

      g12 = {3445, \ 3445, 3445, 3445}; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ Modulo\ na\ direcao\ 12\ de\ cada\ lamina\ *) 

      nu = {0.25, \ 0.25, \ 0.25, \ 0.25}; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ Coef.Poisson\ na\ direcao\ 12\ de\ cada\ lamina\ *) 

      h = {10, 10, 10, 10\ }; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  (*\ Espessuras\ das\ laminas.\ mm.\ *) 

      g13 = {3445, 3445, \ 3445, \ 3445\ }; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ Modulos\ \(transversais.\)\ \ \ \ \ \ \ \ *) 

      g23 = \ {1378, \ 1378, \ 1378, \ 1378}; 

      e3 = {\ 6890, 6890, \ 6890, \ 6890}; 
kapa\  = \ 5/6\ \ ; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \  (*\ fator\ de\ cisalhamento\ *) 
mh\  = \ 1; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ numero\ de\ termos\ na\ serie\ em\ x\ *) 
nh\  = \ 1; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ numero\ de\ termos\ na\ serie\ em\ y\ *) 
qo\  = \ 0.001; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ Carga\ distribuida, \ F/Area.\ Se\ mh = \(nh = 1\), \ q0\ e'

        					\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
        o\ modulo\ de\ carga\ senoidal.\ Se\ mh, nh\ diferentes\ de

        	\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
        1, \ carga\ uniforme\ de\ valor\ qo\ *) 
a = 400; 
b = 400; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ Lados\ da\ placa.\ a\ na\ direcao\ x\ *) 
xc = \ a/2; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ Ponto\ para\ o\ calculo\ das\ tensoes\ coplanares, \ 
        deslocamento\ e\ \(momentos.\)*) 
yc = 0/2; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ A\ rigem\ dos\ eixos\ e' no\ vertice\ do\ retangulo\ *) 

      Xt = {1000, 1000, \ 1000, 1000\ }; \ \ \ \ \ \ \ \ \ \ 
       (*\ Xt, \ Xc, \ Yt, \ Yc, \ S\ sao\ os\ valores\ de\ *) 

      Xc = {600, 600, 600, 600}; \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ \ resistencia\ da\ lamina, \ nas\ direcoes\ 1, 2, *) 

      Yt = {30, 30, 30, \ 30}; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
       (*\ a\ tracao\ e\ compressao.\ S\ e'\ cisalhante\ *) 

      Yc = {120, 120, 120, \ 120}; \ 
S = {70, 70, 70, \ 70}; 
      \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ *) \)\)

\(\( (*\ CALCULOS\ 

      	\ \ \ \ \ \ \ \ \ \ \ \ \ COTAS\ *) 
z = Table[0, {j, n + 1}]; 

    z0 = Table[0, {j, n}]; 
\ \ \ \ \ \ \ \ \ 
     (*\ Cotas\ da\ superficie\ inferior\ de\ cada\ lamina\ e\ de\ sua\ 
        sup.media\ *) 
H = Sum[h[\([j]\)], {j, 1, n}]; 
    \ \ \ \ \ \ \ \ \ \ \ \  (*\ espessura\ total\ do\ laminado\ *) 

    z[\([1]\)] = \(-H\)/2; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ cota\ sup.inf.\ do\ laminado\ *) 

    Do[z[\([j + 1]\)] = z[\([j]\)] + h[\([j]\)]; 
	\ \ \ 
      z0[\([j]\)] = z[\([j]\)] + h[\([j]\)]/2, {j, 1, n}]; \)\)

H

\(\( (*\ Inicializar\ areas\ *) 
q = Table[1, {i, 4}, {j, n}]; 
    \ \ \ \ \ \  (*\ Cada\ coluna\ contem\ matriz\ Q\ *) 

    qb = Table[1, {i, 6}, {j, n}]; \ \ \ \ \ 
     (*\ e\ Qbarra\ de\ cada\ lamina\ *) 
cx = Table[1, {i, 3}, {j, n}]; 
    \ \ \ \ \  (*\ rigidez\ transversal, \ C44, C45, C55\ \ *) 

    ma = {0.0, 0.0, 0.0, 0.0, 0.0, 0.0}; \ 

    mb = {0.0, 0.0, 0.0, 0.0, 0.0, 0.0}; 

    md = {0.0, 0.0, 0.0, 0.0, 0.0, 0.0}; 
mac = {0.0, 0.0, 0.0}; 

    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ Matriz\ D.\ Valores\ em\ coluna\ da\ parte\ 
        triang.sup.
			\ \ \ \ \ \ \ \ \ \ \ \ \ d11\ d12\ d22\ d16\ d26\ 
        d66.\ O\ mesmo\ com\ qb.
			\ \ \ \ \ \ \ \ \ \ \ \ \ q11\ q12\ 
        q22\ q66\ \ \ \ \ \ \ \ *) \)\)

\(\( (*\ Calculo\ de\ q\ em\ cada\ lamina\ *) 

    Do[q[\([1, j]\)] = 
        e1[\([j]\)]\^2\/\(e1[\([j]\)] - nu[\([j]\)]\^2\ e2[\([j]\)]\); 
	
      q[\([2, j]\)] = 
        \(nu[\([j]\)]\ e1[\([j]\)]\ e2[\([j]\)]\)\/\(e1[\([j]\)] - 
            nu[\([j]\)]\^2\ e2[\([j]\)]\); 
	
      q[\([3, j]\)] = 
        \(\ e1[\([j]\)]\ e2[\([j]\)]\)\/\(e1[\([j]\)] - 
            nu[\([j]\)]\^2\ e2[\([j]\)]\); 
		
      q[\([4, j]\)] = g12[\([j]\)], {j, 1, n}]\)\)

\(\( (*\ Calculo\ de\ qbarra\ em\ cada\ lamina\ *) 

    Do[c4 = Cos[teta[\([j]\)]\ ]^4; \ \ \ \ \ \ \ \ \ 
		\ \ 
      s4 = Sin[teta[\([j]\)]\ ]^4; 
	\ \ \ \ c2 = Cos[teta[\([j]\)]\ ]^2; 
      \ \ \ \ \ \ \ \ \ \ 
	\ 	\ s2 = Sin[teta[\([j]\)]\ ]^2; 
	\ \ \ \ 
      c3 = Cos[teta[\([j]\)]\ ]^3; \ \ \ \ \ \ \ \ \ \ 
\ \ \ \ 	
      s3 = Sin[teta[\([j]\)]\ ]^3; 
	\ \ \ \ co = Cos[teta[\([j]\)]\ ]; 
      \ \ \ \ \ \ \ \ \ \ 
	\ 	\ se = Sin[teta[\([j]\)]\ ]; 
	
      qb[\([1, j]\)] = 
        q[\([1, j]\)]\ c4 + 2 \((q[\([2, j]\)] + 2  q[\([4, j]\)])\) s2\ c2\ 
           + q[\([3, j]\)] s4; 
	
      qb[\([3, j]\)] = \ 
        q[\([1, j]\)] s4 + 2 \((q[\([2, j]\)] + 2  q[\([4, j]\)])\) s2\ c2 + 
          q[\([3, j]\)] c4\ \ \ ; 
	
      qb[\([6, j]\)] = 
        \((q[\([1, j]\)] + q[\([3, j]\)] - 
                2 \((q[\([2, j]\)] + q[\([4, j]\)])\))\) s2\ c2 + 
          q[\([4, j]\)] \((s4 + c4)\)\ \ ; 
	
      qb[\([2, j]\)] = 
        \((q[\([1, j]\)] + q[\([3, j]\)] - 4  q[\([4, j]\)])\) s2\ c2 + 
          q[\([2, j]\)] \((s4 + c4)\); 
	
      qb[\([4, j]\)] = 
        \((q[\([1, j]\)] - q[\([2, j]\)] - 2  q[\([4, j]\)])\) se\ c3 + 
          \((q[\([2, j]\)] - q[\([3, j]\)] + 2  q[\([4, j]\)])\) s3\ co; 
	
      qb[\([5, j]\)] = 
        \((q[\([1, j]\)] - q[\([2, j]\)] - 2  q[\([4, j]\)])\) s3\ co + 
          \((q[\([2, j]\)] - q[\([3, j]\)] + 2  q[\([4, j]\)])\) se\ c3; 

      \ \ \ \ \ \ \ cx[\([1, j]\)] = g23[\([j]\)]\ c2\  + \ g13[\([j]\)]\ s2; 
      		 (*\ Calculo\ de\ cx\ em\ cada\ lamina\ *) 
\ \ \ \ \ \ \ 
      cx[\([2, j]\)] = \((g13[\([j]\)]\  - \ g23[\([j]\)])\)\ se\ co; 

      \ \ \ \ \ \ \ cx[\([3, j]\)] = g23[\([j]\)]\ s2\  + \ g13[\([j]\)]\ c2, 
      \ \ \ \ \ \ \ \ \ \ \ \ {j, 1, n}]\)\)

TableForm[q]TableForm[qb]TableForm[cx]\(Dimensions[qb]
\)

\(\( (*\ CALCULO\ DAS\ MATRIZES\ DE\ RIGIDEZ\ A, \ B, \ 
      D\ e\ Ac\ de\ cisalhamento*) 

    Do[
	Do[
		ma[\([i]\)] = ma[\([i]\)] + qb[\([i, j]\)]\ h[\([j]\)]; 

        		mb[\([i]\)] = 
          mb[\([i]\)] + qb[\([i, j]\)]\ h[\([j]\)]\ z0[\([j]\)]; 
		
        md[\([i]\)] = 
          md[\([i]\)] + 
            qb[\([i, j]\)] 
              \((h[\([j]\)]\ z0[\([j]\)]\^2\  + h[\([j]\)]\^3/12.0)\)\ , 
        \ {j, 1, n}], {i, 1, 6}]; 
\ 
    Do[
\ 	Do[\ 
        mac[\([i]\)] = mac[\([i]\)] + kapa\ cx[\([i, j]\)]\ h[\([j]\)], 
        \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {j, 1, n}], {i, 1, 3}]; \)\)

TableForm[ma]TableForm[mb]TableForm[md]TableForm[mac]

\(\( (*\ 
      LAMINADO\ RETANGULAR\ SIMPLESMENTE\ APOIADO\ SOB\ CARGA\ UNIFORME\ 
        q0.\ TEORIA\ DE\ PRIM.ORDEM\ *) 
 (*\ inicializacoes\ *) 

    qmn = Table[0, {i, mh}, {j, nh}]; 
wmn = Table[0, {i, mh}, {j, nh}]; 
    \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ coeficientes\ do\ deslocam.\ w, \ psix\ e\ psiy\ *) 

    psixmn = Table[0, {i, mh}, {j, nh}]; 

    psiymn = Table[0, {i, mh}, {j, nh}]; \)\)

\(\( (*\ CARGA\ UNIFORME\ q0\ *) 

     (*\ \ Do[
\ \ \ \ 	
        Do[qmn[\([i, j]\)] = \(16\ qo\)\/\(Pi\^2\ i\ j\), {i, 1, mh, 2}], {j, 
          1, nh, 2}]; \ \ \ \ \ \ \ *) 

     (*\ CARGA\ SENOIDAL\ qo\ \ USADO\ SO\ NO\ CASO\ mh\  = \ 
        \(nh\  = \ 1\)\ \ \ \ \ *) 
\ 
    Do[
\ \ \ \ 	Do[qmn[\([i, j]\)] = qo, {i, 1, mh, 2}], {j, 1, nh, 2}]; 
    \ \ 
TableForm[qmn]\)\)

\(\( (*\ 
      Matriz\ mk\ e\ vetor\ vf\ para\ calculo\ dos\ coeficientes\ de\ 
        deslocamento\ *) 
\ mk = Table[0, {i, 3}, {j, 3}]; \ 

    vf = Table[0, {i, 3}]; 
vsol = Table[0, {i, 3}]; \ \ \ \ \ \ \ \ \)\)

\(Do[
	
      Do[\ pm = \(i\ Pi\)\/a; \ \ pn = \(j\ Pi\)\/b; 
		
        mk[\([1, 1]\)] = 
          md[\([1]\)] pm\^2\  + \ md[\([6]\)] pn\^2\  + \ mac[\([3]\)]; 
		
        mk[\([1, 2]\)] = \((md[\([2]\)]\  + \ md[\([6]\)])\)\ pm\ pn; 
	\ \ 
        mk[\([1, 3]\)] = mac[\([3]\)]\ pm; 
\ 	\ 
        mk[\([2, 2]\)] = 
          md[\([6]\)]\ pm\^2\  + \ md[\([3]\)] pn\^2\  + \ mac[\([1]\)]; 

        \ \ 	mk[\([2, 3]\)] = mac[\([1]\)]\ pn; 
\ \ 	
        mk[\([3, 3]\)] = mac[\([3]\)]\ pm\^2\  + mac[\([1]\)] pn\^2\ ; 

        	\ \ mk[\([2, 1]\)] = mk[\([1, 2]\)]; 
	\ \ 
        mk[\([3, 1]\)] = \ mk[\([1, 3]\)]; 
\ \ \ \ 
        mk[\([3, 2]\)] = mk[\([2, 3]\)]; 
		\ \ \ \ \ \ \ \ \ \ \ 
        vsol = LinearSolve[mk, {0, 0, qmn[\([i, j]\)]}\ ]; 
			
        wmn[\([i, j]\)]\ \ \ \ \ \ \  = vsol[\([3]\)]; 
		\ \ 
        psixmn[\([i, j]\)] = vsol[\([1]\)]; 
		\ \ 
        psiymn[\([i, j]\)] = vsol[\([2]\)], {i, 1, mh, 2}], {j, 1, nh, 2}]\)

TableForm[wmn]

\( (*\ Deslocamento\ no\ ponto\ e\ deslocamento\ maximo\ \ *) 

    wp = Sum[wmn[\([i, j]\)] Sin[\(i\ Pi\ xc\ \)\/a]\ 
          Sin[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {j, 1, nh, 2}]\)wpN = ((wp (H^3) e2[[1]] 100)/(qo (a^4)))wmax = 
      Sum[wmn[[i, j]] Sin[((i Pi )/2)] Sin[((j Pi )/2)], {i, 1, 
          mh, 2}, {j, 1, nh, 2}]wmaxN = ((wmax (H^3) e2[[1]] 100)/(qo (a^4)))

\(\( (*\ \ TESTE : \ DEVERIA\ DAR\ 1\ \ \ *) 

     (*\ Sum[Sum[
            \(16\ \ \ Sin[\(i\ Pi\ \)\/2]\ Sin[\(j\ Pi\ \)\/2]\)\/\(Pi\ Pi\ i
                \ j\), {i, 1, 99, 2}], {j, 1, 99, 2}]\  // 
        N\ \ \(TESTE : \ DEVERIA\ DAR\ 1\)\ *) \)\)

\(\( (*\ 
      \(--\(--\(--
                \(--\(--
                    \(--\(--
                        \(--\(--
                            \(--\(--
                                \(--\(--
                                    \(--\(--
                                        \(--\(--
                                        \(--\(--
                                        \(--\(--\(--\(--
                                        \(--\(--\(--\(--\(--\(--\(--\(--
                                        \(--\(--\(--\(--\(--\(--\(--\(-

                                        													CALCULO
                                        \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\ DOS\ DESLOCAMENTOS\ E\ ROTACOES\ NO
          \ PONTO\ \((x, y)\) = \((xc, yc)\)\ \ *) 

    wpx = Sum[
        wmn[\([i, j]\)] \((\(i\ Pi\ \)\/a)\) Cos[\(i\ Pi\ xc\ \)\/a]\ 
          Sin[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {j, 1, nh, 2}]; 

    wpxx = \(-
        Sum[wmn[\([i, j]\)] \(\((\(i\ Pi\ \)\/a)\)\^2\) 
            Sin[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    wpxxx = \(-
        Sum[wmn[\([i, j]\)] \(\((\(i\ Pi\ \)\/a)\)\^3\) 
            Cos[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    wpy = Sum[
        wmn[\([i, j]\)] \((\(j\ Pi\ \)\/b)\) Sin[\(i\ Pi\ xc\ \)\/a]\ 
          Cos[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {j, 1, nh, 2}]; 

    wpyy = \(-
        Sum[wmn[\([i, j]\)] \(\((\(j\ Pi\ \)\/b)\)\^2\) 
            Sin[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    wpyyy = \(-
        Sum[wmn[\([i, j]\)] \(\((\(j\ Pi\ \)\/b)\)\^3\) 
            Sin[\(i\ Pi\ xc\ \)\/a]\ Cos[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    wpxy = Sum[
        wmn[\([i, j]\)] \((\(i\ Pi\ \)\/a)\) \((\(j\ Pi\ \)\/b)\) 
          Cos[\(i\ Pi\ xc\ \)\/a] Cos[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {j, 
          1, nh, 2}]; 

    wpxxy = \(-
        Sum[wmn[\([i, j]\)] \(\((\(i\ Pi\ \)\/a)\)\^2\) \((\(j\ Pi\ \)\/b)\) 
            Sin[\(i\ Pi\ xc\ \)\/a] Cos[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    wpxyy = \(-
        Sum[wmn[\([i, j]\)] \((\(i\ Pi\ \)\/a)\) \(\((\(j\ Pi\ \)\/b)\)\^2\) 
            Cos[\(i\ Pi\ xc\ \)\/a] Sin[\(j\ Pi\ yc\ \)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 


    psx = Sum[
        psixmn[\([i, j]\)]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
          Cos[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {j, 
          1, nh, 2}]; 

    psxx = \(-
        Sum[
          psixmn[\([i, j]\)] \((\(i\ Pi\)\/a)\)\ \ \ \ 
            Sin[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    psxxx = \(-
        Sum[psixmn[\([i, j]\)] \(\((\(i\ Pi\)\/a)\)\^2\) 
            Cos[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    psxxy = \(-
        Sum[psixmn[\([i, j]\)] \((\(i\ Pi\)\/a)\) \((\(j\ Pi\)\/b)\) 
            Sin[\(i\ Pi\ xc\ \)\/a]\ Cos[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    psxy = Sum[
        psixmn[\([i, j]\)] \((\(j\ Pi\)\/b)\)\ \ \ \ \ \ \ \ \ 
          Cos[\(i\ Pi\ xc\ \)\/a]\ Cos[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {j, 
          1, nh, 2}]; 

    psxyy = \(-
        Sum[psixmn[\([i, j]\)] \((\(j\ Pi\)\/b)\)\^2\ \ \ \ 
            Cos[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    psxxxx = Sum[
        psixmn[\([i, j]\)] \(\((\(i\ Pi\)\/a)\)\^3\) Sin[\(i\ Pi\ xc\ \)\/a]\ 
          Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {j, 1, nh, 2}]; 

    psxxxy = \(-
        Sum[psixmn[\([i, j]\)] \(\((\(i\ Pi\)\/a)\)\^2\) \((\(j\ Pi\)\/b)\) 
            Cos[\(i\ Pi\ xc\ \)\/a]\ Cos[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    psxxyy = Sum[
        psixmn[\([i, j]\)] \(\((\(j\ Pi\)\/b)\)\^2\) \((\(i\ Pi\)\/a)\)\ \ 
          Sin[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {j, 
          1, nh, 2}]; \)\)

\(psy = 
      Sum[psiymn[\([i, j]\)]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
          Sin[\(i\ Pi\ xc\ \)\/a]\ Cos[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {j, 
          1, nh, 2}]; 

    psyx = Sum[
        psiymn[\([i, j]\)] \((\(i\ Pi\)\/a)\)\ \ \ \ Cos[\(i\ Pi\ xc\ \)\/a]\ 
          Cos[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {j, 1, nh, 2}]; 

    psyxx = \(-
        Sum[psiymn[\([i, j]\)] \(\((\(i\ Pi\)\/a)\)\^2\) 
            Sin[\(i\ Pi\ xc\ \)\/a]\ Cos[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    psyxy = \(-
        Sum[psiymn[\([i, j]\)] \((\(i\ Pi\)\/a)\) \((\(j\ Pi\)\/b)\) 
            Cos[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    psyy = \(-
        Sum[psiymn[\([i, j]\)] \((\(j\ Pi\)\/b)\)\ \ \ \ \ \ \ \ \ 
            Sin[\(i\ Pi\ xc\ \)\/a] Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {j, 
            1, nh, 2}]\); 

    psyyy = \(-
        Sum[psiymn[\([i, j]\)] \((\(j\ Pi\)\/b)\)\^2\ \ \ \ 
            Sin[\(i\ Pi\ xc\ \)\/a]\ Cos[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    psyxxy = Sum[
        psiymn[\([i, j]\)] \(\((\(i\ Pi\)\/a)\)\^2\) \((\(j\ Pi\)\/b)\) 
          Sin[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {j, 
          1, nh, 2}]; 

    psyxyy = \(-
        Sum[psiymn[\([i, j]\)] \((\(i\ Pi\)\/a)\)\ \((\(j\ Pi\)\/b)\)\^2\ \ \ 
            Cos[\(i\ Pi\ xc\ \)\/a]\ Cos[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {
            j, 1, nh, 2}]\); 

    psyyyy = Sum[
        psiymn[\([i, j]\)] \((\(j\ Pi\)\/b)\)\^3\ \ \ \ 
          Sin[\(i\ Pi\ xc\ \)\/a]\ Sin[\(j\ Pi\ yc\)\/b], {i, 1, mh, 2}, {j, 
          1, nh, 2}]; \)

\(TableForm[{{psxx, psxy, psxxx, psxxy}, 

        	\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {psxyy, psxxx, 
          psxxxy, psxxyy}, 

        		\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {psyx, psyy, psyxx, 
          psyxy}, 

        		\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ {psyyy, psyxxy, 
          psyxyy, psyyyy}}]\ \ \ \ \ \ \ 
 (*\ Curvaturas\ kapa\ \ *) \)\(kxy\  = \ psxy + psyx; 
kxyx\  = \ psxxy + psyxx; 

    kxyy\  = \ psxyy + psyxy; 

psyy\ \ ; 
TableForm[psiymn]\ ; \)

\(\( (*\ TENSOES\ *) 
sigx = Table[0, {i, 2  n}, {j, 2}]; \ \ 
     (*\ sigx\ na\ sup\ inf\ e\ sup\ de\ cada\ uma\ das\ n\ \(lam.\)*) 

    sigy = Table[0, {i, 2  n}, {j, 2}]; 

    sigxy = Table[0, {i, 2  n}, {j, 2}]; 

    tauxz = Table[0, {i, 2  n}, {j, 2}]; 

    tauyz = Table[0, {i, 2  n}, {j, 2}]; 

    tauxzi = Table[0, {i, 2  n}, {j, 2}]; \ \ 
     (*\ tensoes\ integradas\ \ *) 
tauyzi = Table[0, {i, 2  n}, {j, 2}]; 
    \)\)

\(\(Do[
      sigx[\([2  k - 1, 2]\)] = 
        z[\([k]\)] 
          \((qb[\([1, k]\)] psxx + qb[\([2, k]\)] psyy + 
              qb[\([4, k]\)] \((psxy + psyx)\))\); 
\ \ 
      sigx[\([2  k, 2]\)] = 
        z[\([k + 1]\)] 
          \((qb[\([1, k]\)] psxx + qb[\([2, k]\)] psyy + 
              qb[\([4, k]\)] \((psxy + psyx)\))\); 
	
\ 
      sigy[\([2  k - 1, 2]\)] = 
        z[\([k]\)] 
          \((qb[\([2, k]\)] psxx + qb[\([3, k]\)] psyy + 
              qb[\([5, k]\)] \((psxy + psyx)\))\); 
\ 
      sigy[\([2  k, 2]\)] = 
        z[\([k + 1]\)] 
          \((qb[\([2, k]\)] psxx + qb[\([3, k]\)] psyy + 
              qb[\([5, k]\)] \((psxy + psyx)\))\); 
	
\ 
      sigxy[\([2  k - 1, 2]\)] = 
        z[\([k]\)] 
          \((qb[\([4, k]\)] psxx + qb[\([5, k]\)] psyy + 
              qb[\([6, k]\)] \((psxy + psyx)\))\); 

      sigxy[\([2  k, 2]\)] = 
        z[\([k + 1]\)] 
          \((qb[\([4, k]\)] psxx + qb[\([5, k]\)] psyy + 
              qb[\([6, k]\)] \((psxy + psyx)\))\), 
      \ \ \ \ \ \ \ \ \ \ \ \ \ {k, n}]; \)\)

\(\( (*\ Tensoes\ cisalhantes\ obtidas\ pelas\ rel.constitutivas\ *) 


    \(Do[\ tauyz[\([2  k - 1, 2]\)] = 
        \((cx[\([1, k]\)] \((psy + wpy)\) + cx[\([2, k]\)] \((psx + wpx)\))
          \); 
\ tauyz[\([2  k, 2]\)] = tauyz[\([2  k - 1, 2]\)]; 
\ 
      tauxz[\([2  k - 1, 2]\)] = 
        \((cx[\([2, k]\)] \((psy + wpy)\) + cx[\([3, k]\)] \((psx + wpx)\))
          \); 
\ tauxz[\([2  k, 2]\)] = \ tauxz[\([2  k - 1, 2]\)], 
      \ \ \ \ \ \ \ \ \ \ \ \ \ {k, n}]; \)						\)\)

\(\( (*\ 
      Tensoes\ cisalhantes\ obtidas\ por\ integracao\ das\ eqs.equilibrio\ *) 
      

tauxzi[\([1, 2]\)] = 0; 
tauyzi[\([1, 2]\)] = 0; 

    Do[zk2 = z[\([k]\)]\^2/2; \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
      zk12 = z[\([k + 1]\)]\^2/2; \ \ \ \ \ 
	
      tauxzi[\([2  k, 2]\)] = \ 
        tauxzi[\([2  k - 1, 2]\)] - 
	\ \ \ \ \ 
          \((\ \ qb[\([1, k]\)] zk12\ psxxx + qb[\([2, k]\)] zk12\ psyxy + \ 
              qb[\([4, k]\)]\ zk12\ kxyx + 
		\ \ \ \ 
              qb[\([4, k]\)] zk12\ psxxy + qb[\([5, k]\)] zk12\ psyyy + \ 
              qb[\([6, k]\)]\ zk12\ kxyy\ \ )\) + 
		\ \ 
          \((\ \ qb[\([1, k]\)] zk2\ psxxx + qb[\([2, k]\)] zk2\ psyxy + 
              qb[\([4, k]\)]\ zk2\ kxyx + 
		\ \ \ \ 
              qb[\([4, k]\)] zk2\ psxxy + qb[\([5, k]\)] zk2\ psyyy + 
              qb[\([6, k]\)]\ zk2\ kxyy\ \ )\); 
	
       (*\ superf.inf.da\ prox.lam.*) 
		
      tauxzi[\([2  k + 1, 2]\)] = tauxzi[\([2  k, 2]\)]\ , {k, n}]; \ 


    Do[zk2 = z[\([k]\)]\^2/2; \ \ \ \ \ \ \ \ \ \ \ \ \ 
      zk12 = z[\([k + 1]\)]\^2/2; \ \ 
	
      tauyzi[\([2  k, 2]\)] = \ 
        tauyzi[\([2  k - 1, 2]\)] - 
	\ \ \ \ \ 
          \((\ \ qb[\([2, k]\)] zk12\ psxxy + qb[\([3, k]\)] zk12\ psyyy + \ 
              qb[\([5, k]\)]\ zk12\ kxyy + 
		\ \ \ \ 
              qb[\([4, k]\)] zk12\ psxxx + qb[\([5, k]\)] zk12\ psyxy + 
              qb[\([6, k]\)]\ zk12\ kxyx\ \ )\) + 
		\ \ 
          \((\ \ qb[\([2, k]\)] zk2\ psxxy + qb[\([3, k]\)] zk2\ psyyy + 
              qb[\([5, k]\)]\ zk2\ kxyy + 
		\ \ \ \ 
              qb[\([4, k]\)] zk2\ psxxx + qb[\([5, k]\)] zk2\ psyxy + 
              qb[\([6, k]\)]\ zk2\ kxyx\ \ )\); 
		
      tauyzi[\([2  k + 1, 2]\)] = 
        tauyzi[\([2  k, 2]\)]\ \ \ \ \ \ \ \ \ \ \ \ , {k, n}]; \)\)

\(sigxN = Table[0, {i, 2  n}, {j, 2}]; \ \ 
     (*\ sigx\ na\ sup\ inf\ e\ sup\ de\ cada\ uma\ das\ n\ \(lam.\)*) 

    sigyN = Table[0, {i, 2  n}, {j, 2}]; \ \ \ \ 
     (*\ valores\ normalizados\ *) 
sigxyN = Table[0, {i, 2  n}, {j, 2}]; 

    tauxzN = Table[0, {i, 2  n}, {j, 2}]; 

    tauyzN = Table[0, {i, 2  n}, {j, 2}]; 

    tauxziN = Table[0, {i, 2  n}, {j, 2}]; 

    tauyziN = Table[0, {i, 2  n}, {j, 2}]; 

    sigxN = sigx  H\^2\/\(qo\ a\^2\); 
sigyN = sigy  H\^2\/\(qo\ a\^2\); 

    sigxyN = sigxy  H\^2\/\(qo\ a\^2\); 
tauyzN = tauyz\ H\/\(qo\ a\); 
    \ \ \ \ \ \ \ \ \ tauxzN = tauxz\ H\/\(qo\ a\); 

    tauyziN = tauyzi\ H\/\(qo\ a\); \ \ \ \ \ \ \ \ \ 
    tauxziN = tauxzi\ H\/\(qo\ a\); 

    Do[tauyz[\([2  k - 1, 1]\)] = z[\([k]\)]; \ \ \ \ \ \ \ \ \ \ 
      tauyz[\([2  k, 1]\)]\ \ \  = z[\([k + 1]\)]; 
		\ \ 
      tauxz[\([2  k - 1, 1]\)] = z[\([k]\)]; \ \ \ \ \ \ \ \ \ \ 
      tauxz[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]; 
		\ \ 
      tauyzN[\([2  k - 1, 1]\)] = z[\([k]\)]/H; \ \ 
      tauyzN[\([2  k, 1]\)]\  = z[\([k + 1]\)]/H; 
		\ \ 
      tauxzN[\([2  k - 1, 1]\)] = z[\([k]\)]/H; \ \ 
      tauxzN[\([2  k, 1]\)]\ \  = z[\([k + 1]\)]/H; 
		
\ \ \ \ 	
      tauyzi[\([2  k - 1, 1]\)] = z[\([k]\)]; \ \ \ \ \ \ \ \ \ \ 
      tauyzi[\([2  k, 1]\)]\ \ \  = z[\([k + 1]\)]; 
		\ \ 
      tauxzi[\([2  k - 1, 1]\)] = z[\([k]\)]; \ \ \ \ \ \ \ \ \ 
      tauxzi[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]; 
		\ \ 
      tauyziN[\([2  k - 1, 1]\)] = z[\([k]\)]/H; \ 
      tauyziN[\([2  k, 1]\)]\  = z[\([k + 1]\)]/H; 
		\ \ 
      tauxziN[\([2  k - 1, 1]\)] = z[\([k]\)]/H; \ 
      tauxziN[\([2  k, 1]\)]\ \  = z[\([k + 1]\)]/H, 

      	\ \ \ \ \ \ \ \ {k, n}]; \)

\(\( (*\ TENSOES\ NAS\ DIRECOES\ PRINCIPAIS\ *) 

    sig1 = Table[0, {i, 2  n}, {j, 2}]; \ \ 
     (*\ sigx\ na\ sup\ inf\ e\ sup\ de\ cada\ uma\ das\ n\ \(lam.\)*) 

    sig2 = Table[0, {i, 2  n}, {j, 2}]; 

    sig12 = Table[0, {i, 2  n}, {j, 2}]; 

    Do[co = Cos[teta[\([k]\)]]; \ se = Sin[teta[\([k]\)]]; 
	\ \ \ 
      c2 = co\ co; s2 = se\ se; 
	\ \ 
      sig1[\([2  k - 1, 2]\)] = 
        c2\ sigx[\([2  k - 1, 2]\)] + s2\ sigy[\([2  k - 1, 2]\)] + 
          2\ se\ co\ sigxy[\([2  k - 1, 2]\)]; 
\ \ 	
      sig1[\([2  k, 2]\)] = 
        c2\ sigx[\([2  k, 2]\)] + s2\ sigy[\([2  k, 2]\)] + 
          2\ se\ co\ sigxy[\([2  k, 2]\)]; 
	
	
      sig2[\([2  k - 1, 2]\)] = 
        s2\ sigx[\([2  k - 1, 2]\)] + c2\ sigy[\([2  k - 1, 2]\)] - 
          2\ se\ co\ sigxy[\([2  k - 1, 2]\)]; 
	
      sig2[\([2  k, 2]\)] = 
        s2\ sigx[\([2  k, 2]\)] + c2\ sigy[\([2  k, 2]\)] - 
          2\ se\ co\ sigxy[\([2  k, 2]\)]; 
	
	
      sig12[\([2  k - 1, 2]\)] = 
        \(-se\)\ co\ sigx[\([2  k - 1, 2]\)] + 
          se\ co\ sigy[\([2  k - 1, 2]\)] + 
          \((c2 - s2)\)\ sigxy[\([2  k - 1, 2]\)]; 
		
      sig12[\([2  k, 2]\)] = 
        \(-se\)\ co\ sigx[\([2  k, 2]\)] + se\ co\ sigy[\([2  k, 2]\)] + 
          \((c2 - s2)\)\ sigxy[\([2  k, 2]\)], {k, n}]; \)\)

\(\(sig1N = Table[0, {i, 2  n}, {j, 2}]; \ \ 
     (*\ sig1\ na\ sup\ inf\ e\ sup\ de\ cada\ uma\ das\ n\ \(lam.\)*) 

    sig2N = Table[0, {i, 2  n}, {j, 2}]; \ \ \ \ \ \ 
     (*\ valores\ normalizados\ *) 
sig12N = Table[0, {i, 2  n}, {j, 2}]; 

    sig1N = sig1  H\^2\/\(qo\ a\^2\); 
sig2N = sig2  H\^2\/\(qo\ a\^2\); 

    sig12N = sig12  H\^2\/\(qo\ a\^2\); 
\)\)

\(\(Do[sigx[\([2  k - 1, 1]\)] = z[\([k]\)]; 
	\ \ \ \ 
      sigx[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]; 
		\ \ 
      sigy[\([2  k - 1, 1]\)] = z[\([k]\)]; 
	\ \ \ \ 
      sigy[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]; 
		\ \ 
      sigxy[\([2  k - 1, 1]\)] = z[\([k]\)]; 
	\ \ \ \ 
      sigxy[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]; 

      		\ \ \ \ \ \ \ \ \ \ \ sig1[\([2  k - 1, 1]\)] = z[\([k]\)]; 

      	\ \ \ \ sig1[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]; 
		\ \ 
      sig2[\([2  k - 1, 1]\)] = z[\([k]\)]; 
	\ \ \ \ 
      sig2[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]; 
		\ \ 
      sig12[\([2  k - 1, 1]\)] = z[\([k]\)]; 
	\ \ \ \ 
      sig12[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]; 

      		\ \ \ \ \ \ \ \ \ \ \ sigxN[\([2  k - 1, 1]\)] = z[\([k]\)]/H; 

      	\ \ \ \ sigxN[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]/H; 

      		\ \ sigyN[\([2  k - 1, 1]\)] = z[\([k]\)]/H; 
	\ \ \ \ 
      sigyN[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]/H; 
		\ \ 
      sigxyN[\([2  k - 1, 1]\)] = z[\([k]\)]/H; 
	\ \ \ \ 
      sigxyN[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]/H; 

      		\ \ \ \ \ \ \ \ \ \ \ sig1N[\([2  k - 1, 1]\)] = z[\([k]\)]/H; 

      	\ \ \ \ sig1N[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]/H; 

      		\ \ sig2N[\([2  k - 1, 1]\)] = z[\([k]\)]/H; 
	\ \ \ \ 
      sig2N[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]/H; 
		\ \ 
      sig12N[\([2  k - 1, 1]\)] = z[\([k]\)]/H; 
	\ \ \ \ 
      sig12N[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]/H, 
      \ \ \ \ \ \ \ \ {k, n}]; \)\)

\(\( (*\ MOMENTOS\ FLETORES\ E\ TORCORES, \ Mx, \ My, \ Mxy\ \ *) 

    mx = md[\([1]\)] psxx + md[\([2]\)] psyy + md[\([4]\)] \((psxy + psyx)\); 
    
\ my = md[\([2]\)] psxx + md[\([3]\)] psyy + 
        md[\([5]\)] \((psxy + psyx)\); 
\ 
    mxy = md[\([4]\)] psxx + md[\([5]\)] psyy + 
        md[\([6]\)] \((psxy + psyx)\); 
		mxN = \(mx\ H\)\/\(qo\ a\^3\); \ 
     (*\ Valores\ normalizados\ *) 
	myN = \(mx\ H\)\/\(qo\ a\^3\); 
	
    mxyN = \(mxy\ H\)\/\(qo\ a\^3\); \)\)

\( (*\ FATORES\ DE\ SOLICITACAO\ *) 


    ftm = Table[0, {i, 2  n}, {j, 2}]; \ \ 
     (*\ fatores\ em\ cada\ interface\ para\ tensao\ maxima\ *) 

    fth = Table[0, {i, 2  n}, {j, 2}]; \ \ \ \ \ 
     (*\ fator\ para\ Tsai - Hill\ \ *) 
ftw = Table[0, {i, 2  n}, {j, 2}]; 
    \ \ \ \ \  (*\ fator\ para\ Tsai - Wu\ \ *) 

    Do[
	 (*\ CRITERIO\ DA\ TENSAO\ MAXIMA\ *) 

      ftm[\([2  k - 1, 2]\)] = 
        Max[{If[\ sig1[\([2  k - 1, 2]\)] >= 0, 
              sig1[\([2  k - 1, 2]\)]\/Xt[\([k]\)], \ 
              \(-sig1[\([2  k - 1, 2]\)]\)\/Xc[\([k]\)]], 
            If[\ sig2[\([2  k - 1, 2]\)] >= 0, 
              sig2[\([2  k - 1, 2]\)]\/Yt[\([k]\)], \ 
              \(-sig2[\([2  k - 1, 2]\)]\)\/Yc[\([k]\)]]\ , 
            Abs[sig12[\([2  k - 1, 2]\)]]\/S[\([k]\)]}]; 
		
      ftm[\([2  k, 2]\)] = 
        Max[{If[\ sig1[\([2  k, 2]\)] >= 0, sig1[\([2  k, 2]\)]\/Xt[\([k]\)], 
              \ \(-sig1[\([2  k, 2]\)]\)\/Xc[\([k]\)]], 
            If[\ sig2[\([2  k, 2]\)] >= 0, sig2[\([2  k, 2]\)]\/Yt[\([k]\)], 
              \ \(-sig2[\([2  k, 2]\)]\)\/Yc[\([k]\)]]\ , 
            Abs[sig12[\([2  k, 2]\)]]\/S[\([k]\)]}]; 
	
       (*\ CRITERIO\ DE\ TSAI - HILL\ *) 
	
	
      fth[\([2  k, 2]\)] = 	
        If[\ sig1[\([2  k, 2]\)] >= 0, sig1[\([2  k, 2]\)]\/Xt[\([k]\)], \ 
              \(-sig1[\([2  k, 2]\)]\)\/Xc[\([k]\)]]^2 + 
			
          If[\ sig2[\([2  k, 2]\)] >= 0, sig2[\([2  k, 2]\)]\/Yt[\([k]\)], \ 
              \(-sig2[\([2  k, 2]\)]\)\/Yc[\([k]\)]]^2 + 
			
          \((sig12[\([2  k, 2]\)]\/S[\([k]\)])\)^2 - 
			
          If[\ sig1[\([2  k, 2]\)] >= 0, sig1[\([2  k, 2]\)]\/Xt[\([k]\)], \ 
              \(-sig1[\([2  k, 2]\)]\)\/Xc[\([k]\)]] 
            If[\ sig2[\([2  k, 2]\)] >= 0, sig2[\([2  k, 2]\)]\/Yt[\([k]\)], 
              \ \(-sig2[\([2  k, 2]\)]\)\/Yc[\([k]\)]]; 
	
		
      fth[\([2  k - 1, 2]\)] = 	
        If[\ sig1[\([2  k - 1, 2]\)] >= 0, 
              sig1[\([2  k - 1, 2]\)]\/Xt[\([k]\)], \ 
              \(-sig1[\([2  k - 1, 2]\)]\)\/Xc[\([k]\)]]^2 + 
			
          If[\ sig2[\([2  k - 1, 2]\)] >= 0, 
              sig2[\([2  k - 1, 2]\)]\/Yt[\([k]\)], \ 
              \(-sig2[\([2  k - 1, 2]\)]\)\/Yc[\([k]\)]]^2 + 
			
          \((sig12[\([2  k - 1, 2]\)]\/S[\([k]\)])\)^2 - 
			
          If[\ sig1[\([2  k - 1, 2]\)] >= 0, 
              sig1[\([2  k - 1, 2]\)]\/Xt[\([k]\)], \ 
              \(-sig1[\([2  k - 1, 2]\)]\)\/Xc[\([k]\)]] 
            If[\ sig2[\([2  k - 1, 2]\)] >= 0, 
              sig2[\([2  k - 1, 2]\)]\/Yt[\([k]\)], \ 
              \(-sig2[\([2  k - 1, 2]\)]\)\/Yc[\([k]\)]]; 
		
       (*\ CRITERIO\ TSAI - WU\ *) 

      ftw[\([2  k - 1, 2]\)] = 
        \((1\/Xt[\([k]\)] - 1\/Xc[\([k]\)])\) sig1[\([2  k - 1, 2]\)] + 
          \((1\/Yt[\([k]\)] - 1\/Yc[\([k]\)])\) sig2[\([2  k - 1, 2]\)] + 

          			
	sig1[\([2  k - 1, 2]\)]^2\/\(Xt[\([k]\)] Xc[\([k]\)]\) + 
          sig2[\([2  k - 1, 2]\)]^2\/\(Yt[\([k]\)] Yc[\([k]\)]\) + 
          \((sig12[\([2  k - 1, 2]\)]\/S[\([k]\)])\)^2 - 
			
          \(sig1[\([2  k - 1, 2]\)] 
              sig2[\([2  k - 1, 2]\)]\)\/\@\(Xt[\([k]\)] Xc[\([k]\)] 
                Yt[\([k]\)] Yc[\([k]\)]\); 
	
      ftw[\([2  k, 2]\)] = 
        \((1\/Xt[\([k]\)] - 1\/Xc[\([k]\)])\) sig1[\([2  k, 2]\)] + 
          \((1\/Yt[\([k]\)] - 1\/Yc[\([k]\)])\) sig2[\([2  k, 2]\)] + 
			
          
	sig1[\([2  k, 2]\)]^2\/\(Xt[\([k]\)] Xc[\([k]\)]\) + 
          sig2[\([2  k, 2]\)]^2\/\(Yt[\([k]\)] Yc[\([k]\)]\) + 
          \((sig12[\([2  k, 2]\)]\/S[\([k]\)])\)^2 - 
			
          \(sig1[\([2  k, 2]\)] 
              sig2[\([2  k, 2]\)]\)\/\@\(Xt[\([k]\)] Xc[\([k]\)] Yt[\([k]\)] 
                Yc[\([k]\)]\); 
	
	\ \ 
      fth[\([2  k - 1, 2]\)] = \@fth[\([2  k - 1, 2]\)]; 
		
      fth[\([2  k, 2]\)] = \@fth[\([2  k, 2]\)]; 
	\ 
      ftw[\([2  k - 1, 2]\)] = \@Abs[ftw[\([2  k - 1, 2]\)]]; 
		
      ftw[\([2  k, 2]\)] = \@Abs[ftw[\([2  k, 2]\)]]
	, 
			{k, n}]; 

    \ TableForm[ftm]\)\(TableForm[fth]
 (*\ TableForm[ftw]*) \)

\(\( (*\ Coeficiente\ de\ seguranca\ *) 
maxtm = Max[ftm]; 

    maxth = Max[fth]; 
If[maxtm == 0, cstm = 0, cstm = 1/maxtm]; 

    If[maxth == 0, csth = 0, csth = 1/maxth]; \)\)

\(\(Do[ftm[\([2  k - 1, 1]\)] = z[\([k]\)]; 
	\ \ \ \ 
      ftm[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)]; 
	
      fth[\([2  k - 1, 1]\)] = z[\([k]\)]; 
	\ \ \ \ 
      fth[\([2  k, 1]\)]\ \ \ \ \ \ \  = z[\([k + 1]\)], {k, n}]; \)\)

wmax
                                                  
     (* valor maximo do deslocamento,  no centro. *) wmaxNwp                                                   
     (* deslocamento no ponto (xc, yc) dado  *) wpNPrint["coef.seg tensao maxima ", cstm]Print["coef.seg Tsai-Hill     ", csth]

Print[" Momentos: Mx  = ", mx,   ", My =  ", my,   
      ", Mxy =  ", mxy]Print[" Momentos: MxN = ", mxN, ", MyN = ", myN, 
      ", MxyN = ", mxyN]

TableForm[sig1N]                     
      TableForm[sig2N]                   
      TableForm[sig1]

ListPlot[sigxN,PlotJoined -> True,GridLines -> Automatic,Frame -> True,AxesLabel -> {"sigxN", "zN"}]ListPlot[sig1N,PlotJoined -> True,GridLines -> Automatic,Frame -> True,AxesLabel -> {"sig1N", "zN"}]

TableForm[sigxN]                     
      \(TableForm[sigyN]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 
     (*\ tensoes\ \ na\ direcao\ \ y\ *) 
\)

TableForm[fth]	ListPlot[fth,PlotJoined -> True,GridLines -> Automatic,Frame -> True,AxesLabel -> {"Fator Tsai-Hill.", "z"}]

TableForm[tauyzN]ListPlot[tauyzN,PlotJoined -> True]

TableForm[tauyziN]ListPlot[tauyziN,PlotJoined -> True]

tauyzNtauyziN