$************************************************
$    модель силового взаимодействия точки тела
$    с ограничивающей ровной горизонтальной
$    поверхностью, вызывающего вертикальную
$    силу и горизональную силу трения.
$    Только для плоских моделей.
$    Используется гистерезисная модель трения, с начальным
$    упругим участком
$************************************************
$      формальные параметры шага МОДЕЛЬ
$ p1  - номер тела неподвижной поверхности
$ p2  - номер узла тела неподвижной поверхности
$ p3  - номер контактирующего тела
$ p4  - номер узла контактирующего тела
$
BL14 NPDS  NTIP   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   102    0   0   0   0   0   0
BL16 NMX NAPR  TUF  KPAR  TDF  KPAR
@  'параметры силы трения'
      nmx1   1    492  4    50   1
@ 'характеристика вертикальной силы'
      nmx2   3    2    4    50    1
@ 'характеристика горизонтальной силы'
      nmx3   1    1    1    50    1
BL17  NSE  NT1  NUZ1  NT2  NUZ2  NMX
@ 'описание вертикального контакта'
      nse1 p1  p2   p3   p4   nmx2
@ 'описание горизонтального контакта'
      nse2 p1  p2   p3   p4   nmx3
@ 'вертикальная сила - фик. звено'
      nse3  nse1  -3   npu4   1   3
$@ 'горизонтальное перем. - многополюсник'
$      nse4  p3    p4   npu1   2   1
$@ 'горизонтальная скорость - многополюсник'
$      nse5  p3    p4   npu1   3  -1
@ 'горизонтальное перем. - многополюсник'
      nse4  nse2  -1   npu1   2   1
@ 'горизонтальная скорость - многополюсник'
      nse5  nse2  -2   npu1   3   1
@ 'сила трения - фик. звено'
      nse6  npu1   1   npu2   1   1
@ 'фик. звено - фик. звено'
      nse7  npu2   1   npu3   1   1
@ 'фик. звено - горизонтальная сила'
      nse8  npu3   1   nse2  -1   1
@ 'вертикальная сила - фик. звено'
      nse9  npu4   1   npu1   1   1
@ 'вертикальная реакция в опорном контакте'
#wkt.mac  3  1  nse1  1   3   1
@ 'сила трения в опорном контакте'
#wkt.mac  3  1  nse6  1   3   1
#tspring.mac  nse2
$************************************************
$      формальные параметры шага РАСЧЕТ
$ p1  - жесткость  опоры
$ p2  - коэффициент вязкого сопротивления по вертикали
$ p3  - коэффициент трения
$ p4  - жесткость  упругого участка трения
$ p5  - величина вязкого трения по горизонтали
BL31 NPDS    MAS     JNX     JNY     JNZ
BL33    NXAR    PRSOST    PRSL      KSTROK
@  'параметры силы трения'
        nmx1     1        1          1
          p3     p4         p5      0.0
        nmx1     1        2          1
          0.     0.   0.   0.
@ 'характеристика вертикальной силы'
        nmx2     3        1          1
          p1    0.00000001   0.00000001   0.
        nmx2     3        2          1
          p2     0.   0.   0.
@ 'характеристика горизонтальной силы'
        nmx3     1        1          1
          0.0    0.0      0.0        0.0
        nmx3     1        2          1
          0.0    0.   0.   0.