$************************************************
$    модель силового взаимодействия точки тела
$    с ограничивающей неподвижной ровной горизонтальной
$    поверхностью, вызывающего вертикальную
$    силу и горизональную силу трения.
$    Для плоских и пространственных моделей.
$    Используется гистерезисная модель трения.
$    вязким участком. Ц.м. неподвижного тела должен иметь нулевые координаты
$    лежать в плоскости ограничивающей поверхности.
$    Нормальная реакция вычисляется с помощью пружины в шарнирах 
$    "бесконечной" длины
$************************************************
$      формальные параметры шага МОДЕЛЬ
$ 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   111    0   0   0   0   0   0
BL16 NMX NAPR  TUF  KPAR  TDF  KPAR
@  'параметры силы трения'
      nmx1   1    494  4    50   1
$@ 'характеристика вертикальной силы'
$      nmx2   1    2    4    50    1
@ 'характеристика вертикальной силы'
      nmx2   1    7    3    50    1
BL17  NSE  NT1  NUZ1  NT2  NUZ2  NMX
@ 'описание вертикального контакта fzspat'
      nse1 p1  p2   p3   p4   nmx2
@ 'вертикальная сила - вычислитель силы 1'
      nse2  nse1  -5   npu1   1   1
@ 'перем. по X - многополюсник'
      nse3  p3    p4   npu1   2   1
@ 'перем. по Y - многополюсник'
      nse4  p3    p4   npu1   3   2
@ 'скорость по X - многополюсник'
      nse5  p3    p4   npu1   4  -1
@ 'скорость по Y - многополюсник'
      nse6  p3    p4   npu1   5  -2
@ 'сила трения по X - фик. звено'
      nse7  npu1   1   npu2   1   1
@ 'сила трения по Y - фик. звено'
      nse8  npu1   2   npu3   1   1
@ 'модуль силы трения - фик. звено'
      nse9  npu1   3   npu4   1   1
@ 'сила трения по X - тело 2'
      nse10 npu2   1   p3    p4   1
@ 'сила трения по Y - тело 2'
      nse11 npu3   1   p3    p4   2
@ 'вертикальная реакция в опорном контакте'
#wkt.mac  3  1  nse1  1   1   1
@ 'сила трения в опорном контакте'
#wkt.mac  3  1  nse9  1   1   1
@ 'описание вертикального контакта fzspat'
#tspring.mac  nse1
$************************************************
$      формальные параметры шага РАСЧЕТ
$ p1  - жесткость  опоры
$ p2  - коэффициент вязкого сопротивления по вертикали
$ p3  - коэффициент трения
$ p4  - не используется
$ p5  - номер тела неподвижной поверхности
$ p6  - номер узла тела неподвижной поверхности
$ p7  - X координата узла тела неподвижной поверхности
$ p8  - Y координата узла тела неподвижной поверхности
$ p9  - Z координата узла тела неподвижной поверхности
$
BL32 NPDS    NUZ      X     Y    Z
r=500
      p5      p6      p7    p8     -r
BL33    NXAR    PRSOST    PRSL      KSTROK
@  'параметры силы трения'
        nmx1     1        1          1
          p3     p1        p2      0.0
        nmx1     1        2          1
          0.     0.   0.   0.
@ 'характеристика вертикальной силы'
        nmx2     1        1          1
         p1     -0.00000001   0.    0.
        nmx2     1        2          1
          p2     0.   0.   0.
$       nmx2     1        1          1
$        0.00000001   p1  0.00000001   0.
$       nmx2     1        2          1
$         p2     0.   0.   0.
@ 'описание вертикального контакта fzspat'
#lspring.mac  nse1 r+p9