# MBDyn (C) is a multibody analysis code. 
# http://www.mbdyn.org
# 
# Copyright (C) 1996-2022
# 
# Pierangelo Masarati	<masarati@aero.polimi.it>
# Paolo Mantegazza	<mantegazza@aero.polimi.it>
# 
# Dipartimento di Scienze e Tecnologie Aerospaziali - Politecnico di Milano
# via La Masa, 34 - 20156 Milano, Italy
# http://www.aero.polimi.it
# 
# Changing this copyright notice is forbidden.
# 
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation (version 2 of the License).
# 
# 
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
# 
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA

begin: data;
	problem: initial value; # the default
end: data;

begin: initial value;
	initial time: 0.;
	final time: 1.;
	time step: 1.e-3;

	max iterations: 10;
	tolerance: 1.e-6;

	derivatives coefficient: 1.e-6;
end: initial value;

begin: control data;
	structural nodes: 1;
	rigid bodies: 1;
	forces: 1;

	output results: netcdf; #, no text;
end: control data;

begin: nodes;
	# in zero, with horizontal speed
	structural: 
		1,		# the label of the node
		dynamic,  	# the type of node: dynamic bc we want to add inertia props
		null, 		# initial position (3x1 vector)
		eye, 		# initial orientation (orientation matrix). Alternatives:
		# matr, 1.,0.,0., 0.,1.,0., 0.,0.,1.,
		# 1, 1.,0.,0.,  2,0.,1.,0.,
		# 3, 0.,0.,1.,  1,1.,0.,0.,
		0.,1.,0., 	# initial velocity (3x1 vector): component in one of the two horizontal directions, 0 in the vertical 
		null; 		# initial angular velocity (3x1 vector)
end: nodes;

begin: elements;
	body: 
		1, 		# the label of body
		1, 		# the label of the node
		1., 		# the mass
		null, 		# the offset of its center of mass with respect to the node. Alernatives:
		# reference, node, null 
		eye;		# inertia tensor: for the moment i want it to be positive definite so that the problem is solvable, but I don't care about the coefficient as the body doesn't rotate for now
	
	force: 
		1, 		# the label of force	
		absolute, 	# the force type (absolute: direction of force remains cte even if node rotates)
		1,		# the label of the node it is attached to
		position, null,	# the offset from the node of the point where the force is applied
		0.,0.,1.,	# the direction of the force
		const, -9.81; 	# the magnitude
	/*
	 * we have better ways to define gravity, though
	 */
end: elements;

# vim:ft=mbd

/*
* to run use
mbdyn rigidbody.mbd
ls -lrt (to show the desired output)
*/