{"comment":{"url":"https://api.github.com/repos/Rosilopes/Naca/comments/13126967","html_url":"https://github.com/Rosilopes/Naca/commit/2f2cdfff63381460afb328f2739db7b501b6c001#commitcomment-13126967","id":13126967,"user":{"login":"Rosilopes","id":14096483,"avatar_url":"https://avatars.githubusercontent.com/u/14096483?v=3","gravatar_id":"","url":"https://api.github.com/users/Rosilopes","html_url":"https://github.com/Rosilopes","followers_url":"https://api.github.com/users/Rosilopes/followers","following_url":"https://api.github.com/users/Rosilopes/following{/other_user}","gists_url":"https://api.github.com/users/Rosilopes/gists{/gist_id}","starred_url":"https://api.github.com/users/Rosilopes/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/Rosilopes/subscriptions","organizations_url":"https://api.github.com/users/Rosilopes/orgs","repos_url":"https://api.github.com/users/Rosilopes/repos","events_url":"https://api.github.com/users/Rosilopes/events{/privacy}","received_events_url":"https://api.github.com/users/Rosilopes/received_events","type":"User","site_admin":false},"position":2,"line":2,"path":"README.md","commit_id":"2f2cdfff63381460afb328f2739db7b501b6c001","created_at":"2015-09-09T03:27:29Z","updated_at":"2015-09-09T03:27:29Z","body":"% PERFIL NACA 2412\r\n% 34 PAINÉIS\r\n% ELEMENTOS DE SINGULARIDADES UTILIZADOS: FONTES E VÓRTICES\r\n\r\n% ESTRUTURA DO CÓDIGO\r\n    % 1.  VARIÁVEIS DE ENTRADA\r\n    % 2.  PERFIL\r\n    % 3.  GEOMETRIA DOS PAINÉIS\r\n    % 4.  DETERMINAÇÃO DOS VETORES (r_i), (r_i+1) E DOS ÂNGULOS BETAS (B)\r\n    % 5.  COMPONENTES DE VELOCIDADES INDUZIDAS PELAS SINGULARIDADES NO SISTEMA LOCAL\r\n    % 6.  COMPONENTES DE VELOCIDADES INDUZIDAS PELAS SINGULARIDADES NO SISTEMA GLOBAL\r\n    % 7.  CÁLCULO DA MATRIZ A PARA DETERMINAR AS SINGULARIDADES\r\n    % 8.  CÁLCULO DA MATRIZ b PARA DETERMINAR AS SINGULARIDADES\r\n    % 9.  CÁLCULO DAS SINGULARIDADES\r\n    % 10. VELOCIDADE TANGECIAL \r\n    % 11. COEFICIENTE DE PRESSÃO\r\n    % 12. COEFICIENTES AERODINÂMICOS\r\n    % 13. GRÁFICOS\r\n\r\n%__________________________________________________________________________\r\n\r\n% 1 VARIÁVEIS DE ENTRADA\r\n\r\n% 1.1 COORDENADAS CARTESIANAS DO PERFIL (34 PAINÉIS, 35 PONTOS DE FRONTEIRA)\r\nP=[1.0000\t-0.00\r\n0.9500\t-0.0048\r\n0.9000\t-0.0082\r\n0.8000\t-0.015\r\n0.7000\t-0.0214\r\n0.6000\t-0.0276\r\n0.5000\t-0.0334\r\n0.4000\t-0.038\r\n0.3000\t-0.0412\r\n0.2500\t-0.0422\r\n0.2000\t-0.0423\r\n0.1500\t-0.041\r\n0.1000\t-0.0375\r\n0.0750\t-0.0346\r\n0.0500\t-0.0301\r\n0.0250\t-0.0227\r\n0.0125\t-0.0165\r\n0.0000\t0\r\n0.0125\t0.0215\r\n0.0250\t0.0299\r\n0.0500\t0.0413\r\n0.0750\t0.0496\r\n0.1000\t0.0563\r\n0.1500\t0.0661\r\n0.2000\t0.0726\r\n0.2500\t0.0767\r\n0.3000\t0.0788\r\n0.4000\t0.078\r\n0.5000\t0.0724\r\n0.6000\t0.0636\r\n0.7000\t0.0518\r\n0.8000\t0.0375\r\n0.9000\t0.0208\r\n0.9500\t0.0114\r\n1.0000\t0.00\r\n];\r\n\r\n% 1.2 ÂNGULO THETA ENTRE OS PAINEIS E O EIXO DAS ABSCISSAS (radianos)\r\ntheta=[3.237299362;\r\n3.209488133;\r\n3.209488133;\r\n3.205505487;\r\n3.203513394;\r\n3.199527748;\r\n3.18756025;\r\n3.173581738;\r\n3.161589988;\r\n3.143592651;\r\n3.11559851;\r\n3.071706652;\r\n3.026108792;\r\n2.963499716;\r\n2.853809612;\r\n2.681150162;\r\n2.219128316;\r\n1.044169055;\r\n0.591685828;\r\n0.427832312;\r\n0.320550075;\r\n0.261846791;\r\n0.193546464;\r\n0.129275004;\r\n0.081816949;\r\n0.04197533;\r\n6.275185479;\r\n6.227243737;\r\n6.195411416;\r\n6.165728455;\r\n6.141148256;\r\n6.117712324;\r\n6.097354383;\r\n6.059017264];\r\n\r\n% 1.3 VELOCIDADE DO ESCOAMENTO\r\nV=10                        % Velocidade de corrente livre (m/s)\r\n\r\n% 1.4 ÂNGULOS DE ATAQUE (VARIAÇÃO DE -10° À 15°)\r\na=-10:1:15;                         % Ângulo de ataque (graus)\r\n \r\nalpha1=(a(1,1)/180)*pi;             % Ângulo de ataque (radianos)  -10°\r\nalpha2=(a(1,2)/180)*pi;             % Ângulo de ataque (radianos)   -9°\r\nalpha3=(a(1,3)/180)*pi;             % Ângulo de ataque (radianos)   -8°\r\nalpha4=(a(1,4)/180)*pi;             % Ângulo de ataque (radianos)   -7°\r\nalpha5=(a(1,5)/180)*pi;             % Ângulo de ataque (radianos)   -6°\r\nalpha6=(a(1,6)/180)*pi;             % Ângulo de ataque (radianos)   -5°\r\nalpha7=(a(1,7)/180)*pi;             % Ângulo de ataque (radianos)   -4°\r\nalpha8=(a(1,8)/180)*pi;             % Ângulo de ataque (radianos)   -3°\r\nalpha9=(a(1,9)/180)*pi;             % Ângulo de ataque (radianos)   -2°\r\nalpha10=(a(1,10)/180)*pi;           % Ângulo de ataque (radianos)   -1°\r\nalpha11=(a(1,11)/180)*pi;           % Ângulo de ataque (radianos)    0°\r\nalpha12=(a(1,12)/180)*pi;           % Ângulo de ataque (radianos)    1°\r\nalpha13=(a(1,13)/180)*pi;           % Ângulo de ataque (radianos)    2°\r\nalpha14=(a(1,14)/180)*pi;           % Ângulo de ataque (radianos)    3°\r\nalpha15=(a(1,15)/180)*pi;           % Ângulo de ataque (radianos)    4°\r\nalpha16=(a(1,16)/180)*pi;           % Ângulo de ataque (radianos)    5°\r\nalpha17=(a(1,17)/180)*pi;           % Ângulo de ataque (radianos)    6°\r\nalpha18=(a(1,18)/180)*pi;           % Ângulo de ataque (radianos)    7°\r\nalpha19=(a(1,19)/180)*pi;           % Ângulo de ataque (radianos)    8°\r\nalpha20=(a(1,20)/180)*pi;           % Ângulo de ataque (radianos)    9°\r\nalpha21=(a(1,21)/180)*pi;           % Ângulo de ataque (radianos)   10°\r\nalpha22=(a(1,22)/180)*pi;           % Ângulo de ataque (radianos)   11°\r\nalpha23=(a(1,23)/180)*pi;           % Ângulo de ataque (radianos)   12°\r\nalpha24=(a(1,24)/180)*pi;           % Ângulo de ataque (radianos)   13°\r\nalpha25=(a(1,25)/180)*pi;           % Ângulo de ataque (radianos)   14°\r\nalpha26=(a(1,26)/180)*pi;           % Ângulo de ataque (radianos)   15°\r\n\r\nalpha=(a/180)*pi;\r\n\r\n%__________________________________________________________________________\r\n\r\n% 2 PERFIL\r\n\r\n% 2.1 COORDENADAS DOS PAINÉIS \r\nX_in=P(1:18,1);      % Coordenada x do intradorso\r\nY_in=P(1:18,2);      % Cooerdenada y do intradorso\r\nX_ex=P(18:35,1);     % Coordenada x do extradorso\r\nY_ex=P(18:35,2);     % Cooerdenada y do extradorso\r\n\r\n% 2.2 DESENHO DO PERFIL \r\nfigure;plot(X_ex,Y_ex,'b','linewidth',2);hold on; plot(X_in,Y_in,'r','linewidth',2), grid\r\ntitle('PERFIL NACA 2412')\r\naxis([0 1 -.5 .5])\r\nxlabel('x')\r\nylabel('y')\r\n\r\n%__________________________________________________________________________\r\n\r\n% 3. GEOMETRIA DOS PAINÉIS\r\n\r\n% 3.1 PONTOS DE FRONTEIRA\r\nX_f=P(:,1);          % Coordenada x dos Pontos de fronteiras \r\nY_f=P(:,2);          % Coordenada y dos Pontos de fronteiras\r\nX_i=P(1:34,1);      % Coordenada x do ponto i\r\nX_i1=P(2:35,1);     % Coordenada x do ponto i+1\r\nY_i=P(1:34,2);      % Coordenada y do ponto i\r\nY_i1=P(2:35,2);     % Coordenada y do ponto i+1\r\n\r\n% 3.2 PONTOS DE CONTROLE (Ponto médio dos painéis)\r\nX_c=(X_i+X_i1)./2;   %Coordenada x dos pontos de controle\r\nY_c=(Y_i+Y_i1)./2;   %Coordenada y dos pontos de controle\r\n\r\n% 3.3 COMPRIMENTO DOS PAINÉIS\r\nL=((X_i1-X_i).^2 + (Y_i1-Y_i).^2).^.5;\r\n\r\n%__________________________________________________________________________\r\n\r\n% 4. DETERMINAÇÃO DOS VETORES (r_i), (r_i+1) E DOS ÂNGULOS BETAS (B)\r\n\r\n% 4.1 MATRIZ DOS VETORES r_i\r\nr_i(:,1)=sqrt((X_i(1,1)-X_c).^2+(Y_i(1,1)-Y_c).^2);\r\nr_i(:,2)=sqrt((X_i(2,1)-X_c).^2+(Y_i(2,1)-Y_c).^2);\r\nr_i(:,3)=sqrt((X_i(3,1)-X_c).^2+(Y_i(3,1)-Y_c).^2);\r\nr_i(:,4)=sqrt((X_i(4,1)-X_c).^2+(Y_i(4,1)-Y_c).^2);\r\nr_i(:,5)=sqrt((X_i(5,1)-X_c).^2+(Y_i(5,1)-Y_c).^2);\r\nr_i(:,6)=sqrt((X_i(6,1)-X_c).^2+(Y_i(6,1)-Y_c).^2);\r\nr_i(:,7)=sqrt((X_i(7,1)-X_c).^2+(Y_i(7,1)-Y_c).^2);\r\nr_i(:,8)=sqrt((X_i(8,1)-X_c).^2+(Y_i(8,1)-Y_c).^2);\r\nr_i(:,9)=sqrt((X_i(9,1)-X_c).^2+(Y_i(9,1)-Y_c).^2);\r\nr_i(:,10)=sqrt((X_i(10,1)-X_c).^2+(Y_i(10,1)-Y_c).^2);\r\nr_i(:,11)=sqrt((X_i(11,1)-X_c).^2+(Y_i(11,1)-Y_c).^2);\r\nr_i(:,12)=sqrt((X_i(12,1)-X_c).^2+(Y_i(12,1)-Y_c).^2);\r\nr_i(:,13)=sqrt((X_i(13,1)-X_c).^2+(Y_i(13,1)-Y_c).^2);\r\nr_i(:,14)=sqrt((X_i(14,1)-X_c).^2+(Y_i(14,1)-Y_c).^2);\r\nr_i(:,15)=sqrt((X_i(15,1)-X_c).^2+(Y_i(15,1)-Y_c).^2);\r\nr_i(:,16)=sqrt((X_i(16,1)-X_c).^2+(Y_i(16,1)-Y_c).^2);\r\nr_i(:,17)=sqrt((X_i(17,1)-X_c).^2+(Y_i(17,1)-Y_c).^2);\r\nr_i(:,18)=sqrt((X_i(18,1)-X_c).^2+(Y_i(18,1)-Y_c).^2);\r\nr_i(:,19)=sqrt((X_i(19,1)-X_c).^2+(Y_i(19,1)-Y_c).^2);\r\nr_i(:,20)=sqrt((X_i(20,1)-X_c).^2+(Y_i(20,1)-Y_c).^2);\r\n"}}