$ lanbp1, modified from H. Chang's nlrgd10 file
sol 105
diag 8
include 'lanbpa.v2004' $ fix error in Lanczos
include 'lanbchka.v2004' $ check eigensolutions
cend
title = QA check for nonlinear rigid element, RBAR   lanbp1
subti = nonlinear static analysis - Lagrange multipler method
label = bars in series with rbar, baseline, no lanba alter
temp(init) = 10
disp = all
spc  = 10
rigid = lagr
MPCF=ALL $ SPCF=ALL
ELFORCE=ALL
subcase 1
  param, usetprt, 0
  param, usetsel, -1 $ print them all
  temp(load) = 20
  load = 100
subcase 2
  method = 10
begin bulk
spc1,	10,	123456,	1
eigrl,	10,	.1 	1000. $ two known roots in this range
force,100,5,0,-100.0,1.0,0.0,0.0
$force,100,41,0,-100.0,1.0,0.0,0.0
cord2r,10,,0.0,0.0,0.0,0.0,0.0,1.0,+co1
+co1,0.0,1.0,0.0
$celas2,101,100.0, 3,5
grid,	1,	,	1.0
grid, 3,, 2.0,0.0,0.0
grid, 4,, 4.0,0.0,0.0
grid,	5,	,	5.0
cbar,	13,	13,	1,	3,	0.	0.	1.
Cbar,	35,	13,	3,	5,	== $
Pbar,	13,	13,	1.	2.	3.	4.
Mat1,	13,	100.,	,	.3	
rbar, 4,3,4,,,,,2.0-6
tempd,10,45.0
temp,10,3,50.0
temp,10,4,50.0
tempd,20,49.0
temp,20,3,60.0
enddata