Сборка и испытания IDA-71S1 (этап 1). Проектирование этапа 2

Сообщение

Листригон · Добавлено: 17-10-2005 12:51

--

Что это было?

Lt_Flash · Добавлено: 17-10-2005 12:53

Так...Что-то я не то делаю или не понимаю...
Смотри, вот у нас кусок

Код: Выделить всё

  Compartment[i]:=PbeginA[i]+(Pgas-PbeginA[i])*(1-Power(2,(-(Temin/ThtA[i]))));
  PambA[i]:=(Compartment[i]-ValA[i])*ValB[i];
  SDepth[i]:=(10*PambA[i])-10;

Отсюда приравниваем глубину к нулю, а Pamb к 1, получается что Pcmp:=(1/ValB)+ValA; Далее, берем форумулу Ежика и туда подставляем...Я взял 0.60 потому что считал ради теста для 0.40 кислорода, 0.60 азота...Разложил на коэффициенты А, Б и С (где С уже конкретно НДЛ) и получил...полную фигню Где я что делаю неверно?
Pcmp у меня - это максимальное Pcomp для Pamb = 1;

Код: Выделить всё

  Pcmp:=(1/ValB[i])+ValA[i];
  A:=Power(2,1/ThtA[i]);
  B:=1/-((Pcmp-0.60)/(Pgas-0.60))+1;
  C:=LogN(A,B);

Листригон · Добавлено: 17-10-2005 13:07

Разбираюсь

Листригон · Добавлено: 17-10-2005 14:02

Ошибка по моему в том, что вычитать надо текущее значение Pcomp(накопленное на данный момент)

Lt_Flash · Добавлено: 17-10-2005 14:03

Листригон писал(а):Ошибка по моему в том, что вычитать надо текущее значение Pcomp(накопленное на данный момент)

Куда вычитать и откуда???
У меня вообще вот такое получилось:

Код: Выделить всё

// PambA = 1 => 1 = (Compartment[i]-ValA[i])*ValB[i];
// => 1/ValB[i] = Compartment[i]-ValA[i];
// => (1/ValB[i])+ValA[i] = Compartment[i];
//  Pcomp.max[i]=(1/ValB[i])+ValA[i];

Листригон · Добавлено: 17-10-2005 14:34

Кажется опять придеться лезть к М числам.

Lt_Flash · Добавлено: 17-10-2005 14:43

Листригон писал(а):Кажется опять придеться лезть к М числам.

Зачем к ним лезть? Где изъяны в моей логике, явно же где-то я накосячил! Что-то не туда выразил или еще как-то...

Листригон · Добавлено: 17-10-2005 15:03

У нас какая-то принципиальная ошибка. Обрати внимание на то ,что в

PambA = 1 => 1 = (Compartment-ValA)*ValB;

Мы получим для одно и то же значение, вне зависимости от глубины, а это не верно - а М числа показывают, что зависит (там более понятно), я кажется понял что надо делать, но более стройно попробую объяснить позже-надо самому яснее осознать.

Lt_Flash · Добавлено: 17-10-2005 15:08

Листригон писал(а):У нас какая-то принципиальная ошибка. Обрати внимание на то ,что в
PambA = 1 => 1 = (Compartment-ValA)*ValB;

Мы получим для одно и то же значение, вне зависимости от глубины, а это не верно - а М числа показывают, что зависит (там более понятно), я кажется понял что надо делать, но более стройно попробую объяснить позже-надо самому яснее осознать.

Видимо дело в том, что в Compartment у нас ТЕКУЩЕЕ насыщение тканей...Только тогда как нам все это приравнять к нужной нам единице, которая соответствует глубине в 0 метров?
Думаю, может вместо Compartment надо подставить пересчитанное Pcomp для данной ткани для текущей глубины? Или как?

Листригон · Добавлено: 17-10-2005 15:28

Все понял.
Опять М числа помогли- формулу Рamb=(Pcomp-a)*b, можно переписать ввиде Рcomp=Pamb/b +a , где Pcomp максимальное давление в ткани при внешнем давлении Pamb, позволяющее подняться до поверхости без декомпрессии. Таким образом нам надо приравнивать не к 1 , а к текущему внешнему давлению, чтобы вычислить это максимальное значение. Далее подставляем в формулу Ежика вместо Рcomp, а в качестве Pbegin текущее значение Рсомр.
Вроде так.

Lt_Flash · Добавлено: 17-10-2005 15:32

Листригон писал(а):Все понял.
Опять М числа помогли- формулу Рamb=(Pcomp-a)*b, можно переписать ввиде Рcomp=Pamb/b +a , где Pcomp максимальное давление в ткани при внешнем давлении Pamb, позволяющее подняться до поверхости без декомпрессии. Таким образом нам надо приравнивать не к 1 , а к текущему внешнему давлению, чтобы вычислить это максимальное значение. Далее подставляем в формулу Ежика вместо Рcomp, а в качестве Pbegin текущее значение Рсомр.
Вроде так.

Щас попробуем...
Че-то все равно не то...Можешь прямо здесь поправить, что я делаю не так?

Код: Выделить всё

  Pcmp:=(((1+(StrToInt(DepthL.Text)*0.1)))/ValB[i])+ValA[i];
  A:=Power(2,1/ThtA[i]);
  B:=1/-((Pcmp-Compartment[i])/(Pgas-Compartment[i]))+1;
  C:=LogN(A,B);

А то что-то я сегодня не соображаю...

Получается такая задача...
У нас есть Pcomp, в котором есть Te - которое является временем нахождения на данной глубине. Далее, у нас есть Pamb - которое показывает давление, до которого мы можем всплыть без появления пузырей (грубо говоря). Из всего этого нам надо выразить Te причем такое, при котором Pamb = 1. Вот тут я и запутался

Листригон · Добавлено: 17-10-2005 15:50

Так и я на этом же попался.
В коде вроде не вижу ошибки пока.
А что при вычислениях получается?

Lt_Flash · Добавлено: 17-10-2005 15:52

Листригон писал(а):Так и на этом же попался.
В коде вроде не вижу ошибки пока.
А что при вычислениях получается.

Ну смотри на фтп bullmann.exe... Четвертая колонка - НДЛ который получился, пятая - это я просто проверял как он берет корень степени tht из 2х...

Листригон · Добавлено: 17-10-2005 16:00

Очень интересно, все минусом .

Lt_Flash · Добавлено: 17-10-2005 16:03

Листригон писал(а):Очень интересно, все минусом .

Угу...Мне тоже интересно.Да даже если минус отбросить - тут прямая последовательность - от 1й до 16й ткани все больше и больше НДЛ. Короче, где-то что-то не так мы считаем. У Ежика вроде верная формула, подставляю тоже верно...

Листригон · Добавлено: 17-10-2005 16:12

Попробую вечером поподробнее посмотреть, сейчас убегаю. Код скинь.
И еще, может попробовать всетаки подбором P comp для данной глубины известно, текущая тоже (Pbegin) и далее итерациями, скажем через 10 секунд ( 1/6 min) или даже менее точно по 0,5 min?.

Lt_Flash · Добавлено: 19-10-2005 21:07

Думал про НДЛ. Если взять Pcomp макс для данной глубины и вычесть Pcomp текущий, то нихрена хорошего все равно не получается. Это то же самое что взять Te для заданной глубины (которое как бы должно быть НДЛ для нее). А Те у нас нихрена и не получается. Смотри, возьмем Pamb=1, формулу Ежика. Что подставляем в Pbegin? Что в Pgas? Я так думаю сначала надо вывести это самое Те, которое НДЛ для строго фиксированной глубины, иначе нифига не понятно...

Листригон · Добавлено: 20-10-2005 09:57

Думаю так надо делать.
Имеем выражение:
Pcomp = Pbegin + [ Pgas - Pbegin ] x [ 1 - 2 -te/tht ], где Pbegin=Pcomp(текущее), Рcomp получаем из выражнения
Pcomp.max=(Pamb(текущее) /b)+a; т.е Рcomp=Pcomp.max,
все члены известны, кроме te, которое в данном случае и есть NDL,
чтобы пока не связываться с формулами, просто подобрать te=NDL итерациями в цикле, допустим через 0,5min.

Lt_Flash · Добавлено: 20-10-2005 10:22

Листригон писал(а):Думаю так надо делать.
Имеем выражение:
Pcomp = Pbegin + [ Pgas - Pbegin ] x [ 1 - 2 -te/tht ], где Pbegin=Pcomp(текущее), Рcomp получаем из выражнения
Pcomp.max=(Pamb(текущее) /b)+a; т.е Рcomp=Pcomp.max,
все члены известны, кроме te, которое в данном случае и есть NDL,
чтобы пока не связываться с формулами, просто подобрать te=NDL итерациями в цикле, допустим через 0,5min.

Не понял, как это Pcomp=Pcomp(текущее) и одновременно Pcomp=Pcomp.max???

Lt_Flash · Добавлено: 24-10-2005 08:27

Листригон, сегодня только начинаю заниматься НДЛ. Сейчас и приступаю. У тебя никакой новой инфы нету?