Текст программы
dimension a(11), b(11), c(11), d(11), e(11), f(11), t(11)
open(2,file='kurs4.res')
n=11
fom=0.6
dfo=0.06
dx=0.1
Bi1=0.215054
Bi2=0.021505
fo=0
t0=0.04
tcp1=1
tcp2=0.2
do i=1,n
t(i)=t0
end do
a(1)=-1/dx
b(1)=1/dx+Bi1
c(1)=0
do i=2,n-1
a(i)=-1/(dx**2)
b(i)=2/(dx**2)+1/dfo
c(i)=-1/(dx**2)
end do
a(n)=0
b(n)=1/dx+Bi2
c(n)=-1/dx
1 fo=fo+dfo
d(1)=-Bi1*tcp1
d(n)=-Bi2*tcp2
do i=2,n-1
d(i)=-t(i)/dfo
end do
e(1)=-a(1)/b(1)
f(1)=-d(1)/b(1)
do
i=2,n-1
e(i)=-a(i)/(b(i)+e(i-1)*c(i))
f(i)=-(d(i)+f(i-1)*c(i))/(b(i)+e(i-1)*c(i))
end do
t(n)=-(d(n)+f(n-1)*c(n))/(b(n)+c(n)*e(n-1))
do i=n-1,1,-1
t(i)=e(i)*t(i+1)+f(i)
end do
write(2,5) fo
write(2,5) t
5 format(2x,11f7.3)
if (fo.lt.fom) go to 1
stop
end
Комментарии к программе
Определяющую
температуру t*=tср1=500
0С,
тогда
,
а
,
.
Предварительно вычисляем числа Био и Фурье:
,
.
Шаг по координате выбираем h=0.1. При этом интервал [0;1]разбиваем на 10 участков, то есть число узлов N=11.
Шаг по времени выбираем равным fi=0.06. При этом время нагрева делится на М=10 точек.
Результаты
X/Fo |
0,06 |
0,12 |
0,18 |
0,24 |
0,3 |
0,36 |
0,42 |
0,48 |
0,54 |
0,6 |
0 |
0,082 |
0,098 |
0,111 |
0,123 |
0,133 |
0,143 |
0,153 |
0,162 |
0,171 |
0,18 |
0,1 |
0,068 |
0,085 |
0,098 |
0,11 |
0,12 |
0,13 |
0,14 |
0,149 |
0,158 |
0,167 |
0,2 |
0,059 |
0,074 |
0,087 |
0,098 |
0,109 |
0,119 |
0,129 |
0,138 |
0,147 |
0,156 |
0,3 |
0,053 |
0,066 |
0,078 |
0,089 |
0,1 |
0,11 |
0,119 |
0,129 |
0,138 |
0,147 |
0,4 |
0,049 |
0,06 |
0,071 |
0,082 |
0,092 |
0,102 |
0,111 |
0,121 |
0,13 |
0,139 |
0,5 |
0,047 |
0,056 |
0,066 |
0,076 |
0,086 |
0,096 |
0,105 |
0,115 |
0,124 |
0,133 |
0,6 |
0,045 |
0,053 |
0,062 |
0,072 |
0,082 |
0,091 |
0,101 |
0,11 |
0,119 |
0,128 |
0,7 |
0,045 |
0,052 |
0,06 |
0,069 |
0,079 |
0,088 |
0,098 |
0,107 |
0,116 |
0,125 |
0,8 |
0,045 |
0,052 |
0,06 |
0,069 |
0,078 |
0,087 |
0,096 |
0,105 |
0,114 |
0,123 |
0,9 |
0,046 |
0,053 |
0,061 |
0,069 |
0,078 |
0,087 |
0,096 |
0,105 |
0,114 |
0,123 |
1 |
0,049 |
0,055 |
0,063 |
0,071 |
0,08 |
0,089 |
0,098 |
0,106 |
0,115 |
0,124 |
