항공대 열전달 평판온도측정 수치해석 Matlab Code 및 레포트
- 최초 등록일
- 2015.09.23
- 최종 저작일
- 2011.10
- 24페이지/ 한컴오피스
- 가격 2,500원
목차
없음
본문내용
k=237;h=10;Tinf=20; %% 열전도도,대류열전달계수,온도경계층밖온도
p=숫자지정; %% p행 p열 node로 나눔
delx=0.1/p; %% delx=등분길이
C=zeros(p^2,1); %% p^2 만큼의 방정식의 상수행렬
A=zeros(p^2,p^2); %% p^2 만큼의 방정식의 계수행렬
A(1,p*(p-1)+2)=2; A(1,p*(p-2)+1)=1; A(1,p*(p-1)+1)=-4; C(1,1)=-300; %% 9번 방정식 계수 맞춤
A(2,p*(p-1)+p-1)=1; A(2,p*(p-2)+p)=1; A(2,p*p)=-4; C(2,1)=-700; %% 8번 방정식 계수 맞춤
A(3,p-1)=1; A(3,p+p)=2; A(3,p)=-4-(2*h*delx/k); C(3,1)=-400-(2*h*delx*Tinf/k); %% 7번 방정식 계수 맞춤
A(4,p+1)=1; A(4,2)=1; A(4,1)=-2-(h*delx/k); C(4,1)=-h*delx*Tinf/k; %% 6번 방정식 계수 맞춤
for m=2:1:p-1 %% 5번 방정식 계수 맞춤
A(5+m-2,p*(m-2)+1)=1;
A(5+m-2,p*m+1)=1;
A(5+m-2,p*(m-1)+2)=2;
A(5+m-2,p*(m-1)+1)=-4;
C(5+m-2,1)=0;
<중 략>
p=5로 지정 하였을 때
349.0549 351.2169 357.7077 368.4419 382.9497 400.0000
347.1705 349.3320 355.8710 366.8492 381.9847 400.0000
340.9633 343.0696 349.5951 361.0992 378.1400 400.0000
330.5434 332.3878 338.3408 349.8123 369.4762 400.0000
316.4346 317.5976 321.5679 330.3331 349.9523 400.0000
300.0000 300.0000 300.0000 300.0000 300.0000 350.0000
참고 자료
없음