$***********************************************************
$   макро задает силу давления в цилиндре четырехтактного двс
$   вход регулирующего воздействия через фиктивное звено
$***********************************************************
$       формальные параметры шага МОДЕЛЬ
$ p1  - номер тела цилиндр
$ p2  - номер узла центра крышки цилиндра
$ p3  - номер тела поршень
$ p4  - номер узла центра верхней поверхности цилиндра
$ p5  - номер ф.з. на которое подается регулирующее воздействие
$ p6  - номер тела коленвал
$ p7  - направление оси вращения коленвала (4 - 6)
BL14  NPDS  TPDS  LX  LY  LZ  UX  UY  UZ
@  'вычисление индикаторных сил'
       npu1   nmx1   0   0   0   0   0   0
@  'фиктивное звено'
       npu2   111    0   0   0   0   0   0
@  'фиктивное звено'
       npu3   111    0   0   0   0   0   0
@  'фиктивное звено'
       npu4   111    0   0   0   0   0   0
@  'фиктивное звено'
       npu5   111    0   0   0   0   0   0
@  'инвертор'
       npu6   102    0   0   0   0   0   0
BL16 NMX NAPR  TUF  KPAR  TDF  KPAR
@   'многополюсник индикаторных сил двс'
       nmx1   1     470    10    50  1
@   'измеритель перемещения поршня'
       nmx2   1     1      1     50  1
BL17   NSE NPDS1  NUZ1 NPDS2  NUZ2  NMX
@     'измеритель перемещения поршня'
       nse1   p1  p2  p3   p4  nmx2
@     'Вход - перемещение поршня'
       nse2   nse1  -1   npu1 1   1
@     'Вход - скорость поршня'
       nse3   nse1  -2   npu1 2   1
@     'Вход - угол поворота коленвала'
       nse4   p6    1    npu1 3   p7
@     'Вход - регулирующее воздействие'
       nse5   p5    1    npu1 4   1
@     'Выход -  сила давления на поршень'
       nse6   npu1 1   npu2   1   1
@     'Выход -  сила давления на поршень'
       nse7   npu2 1   npu6   1   1
@     'передача силы на цилиндр и поршень'
       nse8   npu6 1   npu3   1    1
@     'передача смещения'
       nse9   nse1  -1  npu4   1   1
@     'передача смещения'
       nse10  npu4   1  npu5   1   1
@     'передача силы на цилиндр и поршень'
       nse11   npu3 1   nse1  -1    1
@     'индикаторная сила на поршень'
#wkt.mac  3  1  nse6  1  1  1
@     'индикаторная сила на поршень nse'
#wkt.mac  3  1  nse1  1  1  1
@     'смещение поршня'
#wkt.mac  3  1  nse10  1  1  1
@     'смещение поршня nse'
#wkt.mac  2  1  nse1  1  1  1
@  'пружина измеритель'
#tspring.mac nse1
$************************************************
$      формальные параметры шага РАСЧЕТ
$     p1  - площадь поршня
$     p2  - ход поршня
$     P3  - показатель адиабаты цикла сжатия
$     P4  - показатель адиабаты рабочего цикла
$     P5  - угол впрыска/зажигания (градусы)
$     P6  - время сгорания топлива - время скачка давления
$     P7  - величина скачка давления
$     P8  - степень сжатия
$     P9  - номер такта рабочего цикла ( 1 или 3)
$     P10 - давление всасывания
$     P11 - минимальное длина камеры сгорания
BL33    NXAR    PRSOST    PRSL      KSTROK
@   'многополюсник индикаторных сил двс'
         nmx1    1    1    3
        p1     p2      p3    p4
        p5     p6      p7    p8
        p9     p10     0.0   0.0
         nmx1  1    2    1
         0.0    0.       0.      0.
@   'измеритель перемещения поршня'
         nmx2    1    1    1
         0.0   0.0   0.0   0.0
         nmx2    1    2    1
         0.0    0.       0.      0.
@  'пружина измеритель'
#lspring.mac nse1 p11