Error parsing viscous body, linear viscous generic constitutive law
I am testing the viscous body. I assumed the linear viscous generic constitutive law would be a good fit for this, but when I run the smple model below (one mass, gravity and a viscous body joint). It fails with the following error:
warning: "naive" solver should be used with "colamd"
Creating scalar solver with Naive linear solver
Reading Structural(1)
Reading Body(2)
Reading Joint(3)
[/home/rcrozier/src/mbdyn-official-git/libraries/libmbutil/mathp.cc:5383,func=ExpressionElement* MathParser::stmt()] (unable to find variable "default::viscous") at line 57
HighParser::GetValue(): error return from MathParser at line 57, file </home/rcrozier/src/rnfoundry-hg/common/multibody/MBDyn/test/Test_viscousBody/input.mbd>
An error occurred during the execution of MBDyn; aborting...
input file:
begin: data;
problem: initial value;
end: data;
# initial value problem
begin: initial value;
initial time: 0.0;
final time: 5.0;
time step: 0.01;
tolerance: 0.000000001;
max iterations: 20;
linear solver: naive ;
end: initial value;
begin: control data;
structural nodes: 1;
rigid bodies: 1;
joints: 1;
gravity;
default orientation: orientation matrix;
output results: netcdf, no text;
end: control data;
begin: nodes;
# 6 DOF structural node
structural : 1, dynamic, # label, type
position, 0.0, 0.0, 0.0, # absolute position
orientation,
matr,
1.0, 0.0, 0.0,
0.0, 1.0, 0.0,
0.0, 0.0, 1.0, # absolute orientation
velocity, 0.0, 0.0, 0.0, # absolute velocity
angular velocity, 0.0, 0.0, 0.0 # absolute angular velocity
; # end structural node
end: nodes;
begin: elements;
# body
body : 2, 1, # label, node label
1.0, # mass
0.0, 0.0, 0.0, # relative centre of mass
matr,
0.004000000000000001, 0.0, 0.0,
0.0, 0.004000000000000001, 0.0,
0.0, 0.0, 0.004000000000000001 # inertia matrix
; # end one-mass body
# viscous body
joint : 3, viscous body,
1, # node label
linear viscous generic,
matr,
1.0, 0.0, 0.0, 0.0, 0.0, 0.0,
0.0, 1.0, 0.0, 0.0, 0.0, 0.0,
0.0, 0.0, 1.0, 0.0, 0.0, 0.0,
0.0, 0.0, 0.0, 1.0, 0.0, 0.0,
0.0, 0.0, 0.0, 0.0, 1.0, 0.0,
0.0, 0.0, 0.0, 0.0, 0.0, 1.0
; # end viscous body
gravity : uniform, single, 0.0, 0.0, -1.0, const, 9.810000000000000497;
end: elements;