Commit c0f5691a authored by Andreas Zilian's avatar Andreas Zilian
Browse files

Minor modifcations of parameters.

parent c8b47e26
......@@ -15,7 +15,7 @@ s3.m = 0.00; % [kg]
s3.k = 1000; % [N]
%% Definition of time function for time-dependent forces/base motions
tF = mafe.TimeFunction( 2.5, 1.0, 0.0 ); % force: Omega = 2.5 [rad/s], a_c = 1.0, a_s = 0.0
tU = mafe.TimeFunction( 5.0, 0.0, 1.0 ); % base : Omega = 5.0 [rad/s], a_c = 0.0, a_s = 1.0
tU = mafe.TimeFunction( 5.86, 0.0, 1.0 ); % base : Omega = 5.0 [rad/s], a_c = 0.0, a_s = 1.0
tfuns = [tF, tU];
%% Definition of nodes in the system
% Node 1
......@@ -63,7 +63,7 @@ ana = mafe.FeAnalysisDynamicFD( fep, tfuns );
ana.analyse();
%% apply initial conditions
modeshape = real(ana.X(:,1:2:fep.ndofs)); % get some selected (or all) modes
u0 = 1.00 * modeshape(:,1) + 0.00 * modeshape(:,2);
u0 = 0.00 * modeshape(:,1) + 0.00 * modeshape(:,2);
v0 = 0.00 * modeshape(:,1) + 0.00 * modeshape(:,2);
ana.applyInitialConditions( u0, v0 );
......@@ -88,7 +88,7 @@ cla; clf;
fig = figure(1); subplot('Position',[0.05 0.05 0.90 0.90]); axis equal; hold on;
% evaluate and plot for given time instants
% first undamped eigenfrequency is omega = ~5.9 [rad/s] -> T = 2pi/omega = 1.1 [s]
for tt = 0:0.02:1.1*2
for tt = 0:0.02:1.1*20
ana.putSolToDof(tt);
clf; subplot('Position',[0.05 0.05 0.90 0.90]); axis equal; hold on;
xlim([-1 5])
......
......@@ -75,7 +75,7 @@ fep.plotSystem('deformed');
ana.putModeToDof(3, 'real');
fep.plotSystem('deformed');
fep.plotSystem('reference');
%
%% -----------------------------------------------------------------------------
disp('press key to proceed'); pause;
%% -----------------------------------------------------------------------------
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment