Вычислительная математика - ИСЭМ СО РАН
Вычислительная математика - ИСЭМ СО РАН
Вычислительная математика - ИСЭМ СО РАН
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Рис. 1: Изменение скорости движения снаряда<br />
Рис. 2: Изменение давления пороховых газов<br />
Для сокращения объема вычислений, связанных c обращением матрицы, можно воспользоваться<br />
доказанной леммой и предложить эквивалентную схеме (6) формулу:<br />
k∑<br />
β j B(x i+1−j , t i+1−j )<br />
j=1<br />
3. Численные раcчеты<br />
k∑<br />
α j x i+1−j = h<br />
j=0<br />
k∑<br />
β j g(x i+1−j , t i+1−j ). (8)<br />
Рассмотрим численное решение системы (3), редуцированной к виду (5), многошаговыми<br />
методами. Постоянные коэффициенты задачи описывают характеристики авиационной<br />
пушки калибра 30 мм; порохового заряда, состоящего из семиканальных зерен пороха<br />
марки 4/7 Цгр (целлюлозно-гранулированный) и снаряда ОФЗ (осколочно-фугасный зажигательный)<br />
и принимают следующие значения (единицы измерения СИ):<br />
s = 7, 3 · 10 −4 ; W 0 = 1, 09 · 10 −4 ; m = 0, 394; ω = 0, 095; f = 9 · 10 6 ;<br />
α = 9, 49 · 10 −4 ; ρ = 1640; u 1 = 10 −9 ; k = 1, 2; e 1 = 5 · 10 −4 .<br />
j=1<br />
115