$************************************************
$    Модель соединения двух балочных элементов
$************************************************
$      формальные параметры шага МОДЕЛЬ
$ p1  - номер тела 1 балки
$ p2  - номер узла 1 балки
$ p3  - номер тела 2 балки
$ p4  - номер узла 2 балки
$
$BL14 NPDS  NTIP   LX  LY  LZ  UX  UY  UZ
$@  'вспомогательная масса контрактильного элемента'
$        npt1   npt1+100  1   0   0   0   0   0
BL16 NMX NAPR  TUF  KPAR  TDF  KPAR
@  'Характеристика соединения'
      nmx1   1    1    1    50    1
      nmx1   2    1    1    50    1
      nmx1   3    1    1    50    1
      nmx1   4    1    1    50    1
      nmx1   5    1    1    50    1
      nmx1   6    1    1    50    1
BL17  NSE  NT1  NUZ1  NT2  NUZ2  NMX
@  'Соединение балок'
      nse1 p1  p2   p3   p4   nmx1
$************************************************
$      формальные параметры шага РАСЧЕТ
$ p1  - длина 1 балки
$ p2  - площадь  поперечного  сечения  1  балки
$ p3  - Jx 1  балки
$ p4  - Jy 1  балки
$ p5  - Jz 1  балки
$ p6  - признак направления 2 балки (1-3)
$ p7  - длина 2 балки
$ p8  - площадь  поперечного  сечения  2  балки
$ p9  - Jx 2  балки
$ p10 - Jy 2  балки
$ p11 - Jz 2  балки
$ p12 - признак направления 2 балки (1-3)
$
BL31 NPDS    MAS     JNX     JNY     JNZ
d_factor=0.01
e=e_modul
g=g_modul
BL33    NXAR    PRSOST    PRSL      KSTROK
if p6 == 1 then
  cx1=2*e*p2/p1
endif
if p6 == 2 then
  cx1=12*8*e*p5/(p1*p1*p1)
endif
if p6 == 3 then
$?cx1=12*8*e*p4/(p1*p1*p1)
endif
if p12 == 1 then
  cx2=2*e*p8/p7
endif
if p12 == 2 then
  cx2=12*8*e*p11/(p7*p7*p7)
endif
if p12 == 3 then
$?cx2=12*8*e*p10/(p7*p7*p7)
endif
cx=cx1*cx2/(cx1+cx2)
$************************************
if p6 == 1 then
  cy1=e*12*8*p5/(p1*p1*p1)
endif
if p6 == 2 then
  cy1=2*e*p2/p1
endif
if p6 == 3 then
$?cy1=e*12*8*p3/(p1*p1*p1)
endif
if p12 == 1 then
  cy2=e*12*8*p11/(p7*p7*p7)
endif
if p12 == 2 then
  cy2=2*e*p2/p7
endif
if p12 == 3 then
$?  cy2=12*e*p10/(p7*p7*p7)
endif
cy=cy1*cy2/(cy1+cy2)
$************************************
if p6 == 1 then
  cz1=e*12*8*p4/(p1*p1*p1)
endif
if p6 == 2 then
  cz1=e*12*8*p4/(p1*p1*p1)
endif
if p6 == 3 then
$?cz1=e*12*8*p3/(p1*p1*p1)
endif
if p12 == 1 then
  cz2=e*12*8*p10/(p7*p7*p7)
endif
if p12 == 2 then
  cz2=e*12*8*p10/(p7*p7*p7)
endif
if p12 == 3 then
$?  cz2=12*e*p10/(p7*p7*p7)
endif
cz=cz1*cz2/(cz1+cz2)
$************************************
if p6 == 1 then
  cux1=2*g*p3/p1
endif
if p6 == 2 then
  cux1=2*4*e*p4/p1
endif
if p6 == 3 then
$?cux1=e*12*8*p3/(p1*p1*p1)
endif
if p12 == 1 then
  cux2=2*g*p9/p7
endif
if p12 == 2 then
  cux2=2*4*e*p4/p7
endif
if p12 == 3 then
$?  cux2=12*e*p10/(p7*p7*p7)
endif
cux=cux1*cux2/(cux1+cux2)
$************************************
if p6 == 1 then
  cuy1=2*4*e*p4/p1
endif
if p6 == 2 then
  cuy1=2*g*p3/p1
endif
if p6 == 3 then
$?cuy1=e*12*8*p3/(p1*p1*p1)
endif
if p12 == 1 then
  cuy2=2*4*e*p10/p7
endif
if p12 == 2 then
  cuy2=2*g*p3/p7
endif
if p12 == 3 then
$?  cuy2=12*e*p10/(p7*p7*p7)
endif
cuy=cuy1*cuy2/(cuy1+cuy2)
$************************************
if p6 == 1 then
  cuz1=2*4*e*p5/p1
endif
if p6 == 2 then
  cuz1=2*4*e*p5/p1
endif
if p6 == 3 then
$?cuz1=e*12*8*p3/(p1*p1*p1)
endif
if p12 == 1 then
  cuz2=2*4*e*p11/p7
endif
if p12 == 2 then
  cuz2=2*4*e*p11/p7
endif
if p12 == 3 then
$?  cuz2=12*e*p10/(p7*p7*p7)
endif
cuz=cuz1*cuz2/(cuz1+cuz2)
$  'Характеристика соединения'
        nmx1     1        1          1
          cx    0.    0.    0.
        nmx1     1        2          1
          sqrt(cx)*d_factor   0.   0.   0.
        nmx1     2        1          1
          cy    0.    0.    0.
        nmx1     2        2          1
          sqrt(cy)*d_factor   0.   0.   0.
        nmx1     3        1          1
          cz    0.    0.    0.
        nmx1     3        2          1
          sqrt(cz)*d_factor   0.   0.   0.
        nmx1     4        1          1
          cux    0.    0.    0.
        nmx1     4        2          1
          sqrt(cux)*d_factor   0.   0.   0.
        nmx1     5        1          1
          cuy    0.    0.    0.
        nmx1     5        2          1
          sqrt(cuy)*d_factor   0.   0.   0.
        nmx1     6        1          1
          cuz  0.    0.    0.
        nmx1     6        2          1
          sqrt(cuz)*d_factor   0.   0.   0.