CC BY
33
3. Для решения краевой задачи разработан метод, сочетающий модифицированный метод Ньютона и метод градиентного спуска. Метод реализован в программе на языке PASCAL. Ниже приведены результаты расчёта для случая, когда ai=0.2, аг=40, е=0.2. В начальный момент времени управляемый КА движется по круговой орбите Земли и начинает движение для мягкой встречи с КА, который движется по эллиптической орбите в плоскости несовпадающей с плоскостью орбиты Земли В таблице в безразмерных переменных приведены координаты положения и вектора скорости космических аппаратов
Таблица фазовых состояний КА (в безразмерных переменных)
*3 vi V2 V3
1.0 0.0 0.0 0.0 1.0 0.0
2.3100 0.7500 0.1750 -0 1006 06037 0 1409
2.3320 0 1830 0.0427 00557 0 6348 0.1481
В первой строке таблицы указано начальное состояние управляемого аппарата, во второй - неуправляемого, а в третьей - их состояние в момент мягкой встречи В безразмерных переменных длительность перелёта равна 21 6153 или 3.4402 земных года.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1 Челноков Ю Н. Применение кватернионов в теории орбитального движения искусственного спутника И//Космич исслед. 1993 Т. 31, вып 3 С. 3 - 15
2. Сапунков ЯГ Применение КЭ-переменных к задаче оптимального управления космическим аппаратом // Космич исслед 1996 Т 34, вып 4. С 428 - 433
УДК 533.6011
Я. Г. Сапунков, Г. П. Шиндяпнн, В. А. Норшнев, В. Н. Федорец
ВЫЧИСЛИТЕЛЬНЫЙ ЭКСПЕРИМЕНТ НО ОПРЕДЕЛЕНИЮ ГАЗОДИНАМИЧЕСКИХ ПАРАМЕТРОВ В КАМЕРЕ ДЕТОНАНИОННОГО ДВИГАТЕЛЯ'
В статье метод расчёта движения продуктов детонации в цилиндрической детонационной камере длиной 1\ (1 ] обобщается на случай, когда к камере присоединен диффузор длиной /о с углом полураствора а, заполненный нейтральным газом (воздухом).
' Работа выполнена при поддержке Министерства образования РФ, научно-технической программы, проект № 01 01 030
1. Постановка задачи. Ось Ох направляется вдоль оси камеры и диффузора Сечение х=0 соответствует закрытому концу камеры, сечение х -1\ разделяет в начальный момент взрывчатую смесь, находящуюся в детонационной камере левее этого сечения, и воздух в диффузоре Сечение х /1+/0 соответствует открытому концу диффузора. У закрытого конца камеры создается детонационная волна Чспмена-Жуге, которая движется по камере в положительном направлении оси Ох. Па рисунке изображена картина течения продуктов детонации и воздуха на плоскости Ох1 Указаны основные линии сильных и слабых разрывов
В момент когда детонационная волна О достигнет сечения х~1\, возникает произвольный разрыв, который распадается на ударную волну Л\ которая движется по воздуху, находящемуся в диффузоре, и приводит его в движение, контактный разрыв С,', который разделяет воздух и продукты детонации, и волну разрежения в продуктах детонации, фронт которой /'2.
На рисунке область 0 соответствует диффузору, который в начальный момент заполнен воздухом, область 1 - камере с взрывчатой смесью, область 2 - центрированной волне разрежения за фронтом детонации Линия /•'[ - задний фронт центрированной волны разрежения, он отделяет её от области покоя 3 в продуктах детонации. Область 4 соответствует возмущенному движению продуктов детонации, а область 5 - воздуху, прошедшему через ударную волну Л'.
2. Уравнения движения, начальные и граничные условия. Обозначим через Б(х) площадь сечения камеры и диффузора, через /?3, я3 давление, плотность и скорость звука в области покоя 3 соответственно Движение продуктов детонации в камере и диффузоре оказывается изэн-тропическим и описывается следующей системой уравнений:
В системе уравнений (1) через р, р, v, у2 обозначены плотность, давление, скорость и отношение теплоемкостей продуктов детонации соответственно
Движение воздуха в диффузоре описывается уравнениями
?P + v?r + w±(vS)=о (2)
dt дх dt дbe р дх dt дх S дх
В системе уравнений (2) р, р, v, у0 - плотность, давление, скорость и отношение теплоемкостей воздуха соответственно.
Движение детонационной волны и продуктов детонации на промежутке времени 0<t<t\ описаны в [1]. Состояние продуктов детонации в момент времени является начальным состоянием для расчета их дальнейшего движения, начальные условия для воздуха при t=t\ имеют вид
l^<x<J)+lo, Р=Ро, P=Po, v=0, (3)
где ро, Ро - атмосферное давление и плотность воздуха при атмосферном давлении соответственно.
Решение систем дифференциальных уравнений (1), (2) после момента времени t\ должно удовлетворять граничным условиям на закрытом конце камеры, на контактном разрыве, на ударной волне и на выходном сечении диффузора На дне камеры выполняется условие
при х=0, f>/|, v=0. (4)
Скорость движения контактного разрыва обозначим через Uc, через Vc4, Рсл - скорость и давление в области 4 в продуктах детонации непосредственно перед контактным разрывом, а через vC5, рс5 - скорость и давление воздуха в области 5 непосредственно за контактным разрывом. Тогда условия на контактном разрыве при t>t\, jr=jtc, где x=x^í) -закон движения контактного разрыва, запишутся в виде
V« = VC5 =t/c; Ры = Pd (5)
Скорость движения ударной волны S обозначим через Us, скорость звука в области 0 перед ударной волной обозначим через а0. Обозначив через ví5, ps5, pt5 скорость, давление и плотность в воздухе в области 5 непосредственно за ударной волной получим
Ví5 = Us-7 (l-^o). Ps s = Ра
У о + 1
Í 2Уо иг Yo-О
Уо + 1 Го + и
Prf-Po/ - 0 0 = 1 Yo-' ()
(Уо-^о+З a¿ V Ро
В выходном сечении диффузора выполняется условие для давления
t>0,x=l\+l0, р=ро- (7)
При этом проверяется выполнение ограничения на скорость: скорость течения газов в выходном сечении меньше или равна скорости звука Если это условие нарушается, то давление в выходном сечении диффузора определяется из условия, что скорость течения газов в выходном сечении равна скорости звука
Таким образом, параметры течения продуктов детонации и воздуха в детонационном двигателе определяются в результате решения краевой задачи для систем дифференциальных уравнений (1) и (2) для продуктов детонации и воздуха соответственно с начальными условиями (3) и граничными условиями (4) - (7).
3. Результаты расчёта Для решения краевой задачи с использованием ЭВМ составлена программа на языке PASCAL, в которой реализован метод Лакса-Вендроффа. Расчёты и проведенные эксперименты показывают, что наличие диффузора увеличивает суммарный импульс давления детонационного двигателя на 30%.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1 Сапунков Я.Г, Шиндяпин Г П., Поршнев НА, Федоре ц ОН. Математическая модель детонационного двигателя // Математика Механика Сб. науч тр Саратов Изд-во Сарат ун-та, 2001 Вып. 3. С. 178-181
УДК 232.5,232.135
М. И. Сафрончик
СВЕДЕНИЕ НЕСТАЦИОНАРНОЙ ЗАДАЧИ О ДВИЖЕНИИ ВЯЗКОПЛАСТИЧНОЙ СРЕДЫ К СИСТЕМЕ ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
Рассма1ривается течение вязкопластичной среды по наклонной плоскости под действием силы тяжести. При постановке задачи учитывается неодинаковое поведение материала при нагружении и разгрузке (реологическая модель содержит два предела текучести: статический тст и динамический та) Основной математической трудностью при решении задачи является наличие неизвестной, изменяющейся во времени границы области вязкопластичного течения
Ниже предлагается метод решения задачи, являющийся несколько измененным методом "мгновенных собственных функций", разработанным В.Г. Меламедом для решения задачи Стефана [1].
