1# $Header: /var/cvs/mbdyn/mbdyn/mbdyn-1.0/modules/module-mds/mds,v 1.8 2017/01/12 14:54:57 masarati Exp $ 2# MBDyn (C) is a multibody analysis code. 3# http://www.mbdyn.org 4# 5# Copyright (C) 1996-2017 6# 7# Pierangelo Masarati <masarati@aero.polimi.it> 8# Paolo Mantegazza <mantegazza@aero.polimi.it> 9# 10# Dipartimento di Ingegneria Aerospaziale - Politecnico di Milano 11# via La Masa, 34 - 20156 Milano, Italy 12# http://www.aero.polimi.it 13# 14# Changing this copyright notice is forbidden. 15# 16# This program is free software; you can redistribute it and/or modify 17# it under the terms of the GNU General Public License as published by 18# the Free Software Foundation (version 2 of the License). 19# 20# 21# This program is distributed in the hope that it will be useful, 22# but WITHOUT ANY WARRANTY; without even the implied warranty of 23# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 24# GNU General Public License for more details. 25# 26# You should have received a copy of the GNU General Public License 27# along with this program; if not, write to the Free Software 28# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 29# 30 31begin: data; 32 problem: initial value; 33end: data; 34 35begin: initial value; 36 initial time: 0.; 37 final time: 10.; 38 time step: 1.e-2; 39 40 tolerance: 1.e-6; 41 max iterations: 10; 42 43 method: ms, .6; 44 45 # derivatives coefficient: 1.e-12; 46 derivatives tolerance: 1.e-3; 47 derivatives max iterations: 10; 48 49 # linear solver: umfpack, cc; 50 51 output: iterations; 52 #output: residual; 53 #output: solution; 54 #output: jacobian; 55end: initial value; 56 57begin: control data; 58 loadable elements: 1; 59end: control data; 60 61set: real M = 1.; 62set: real OMEGA = 2*pi; 63set: real XI = .1; 64set: real K = M*OMEGA^2; 65set: real D = 2*XI*M*OMEGA; 66 67begin: elements; 68 module load: "libmodule-mds"; 69 user defined: 1, mds, 70 M, D, K, 71 # sine, 0., .2*OMEGA, 1., forever, 0., 72 sine, 0., 10.*OMEGA, 1., one, 0., 73 output, yes; 74end: elements; 75 76# vim:ft=mbd 77