Talk:Corotational Truss Element

From OpenSeesWiki
Revision as of 22:24, 5 July 2011 by Dborello (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

I removed the reference to the doRayleigh flag. That does not seem to be implemented in the source code.

corotTruss appears to implement Rayleigh damping by default. Truss does not. Someone should verify these findings.

I used the following (MATLAB) script for verification:


<nowiki>


clear all; clc %% Damp = 1; Corot = 1;

fout = fopen('Tester\Driver.tcl','w');

fprintf(fout,'wipe\n'); fprintf(fout,'model BasicBuilder -ndm 2 -ndf 3\n'); fprintf(fout,'file mkdir Data/\n');

fprintf(fout,'node 1 0.0 0.0\n'); fprintf(fout,'node 2 1.0 0.0\n'); fprintf(fout,'node 3 2.0 0.0\n'); fprintf(fout,'node 4 3.0 0.0\n'); fprintf(fout,'fix 1 1 1 1 \n'); fprintf(fout,'fix 2 0 1 1\n'); fprintf(fout,'fix 3 0 1 1\n'); fprintf(fout,'fix 4 0 1 1\n'); fprintf(fout,'mass 2 1.0 0.0 0.0\n'); fprintf(fout,'mass 3 1.0 0.0 0.0\n'); fprintf(fout,'mass 4 1.0 0.0 0.0\n');

fprintf(fout,'uniaxialMaterial Elastic 1 100.0\n');

if Corot

   fprintf(fout,'element corotTruss 1 1 2 1.0 1\n');
   fprintf(fout,'element corotTruss 2 2 3 1.0 1\n');
   fprintf(fout,'element corotTruss 3 3 4 1.0 1\n');

else

   fprintf(fout,'element Truss 1 1 2 1.0 1 -doRayleigh 1 \n');
   fprintf(fout,'element Truss 2 2 3 1.0 1 -doRayleigh 1 \n');
   fprintf(fout,'element Truss 3 3 4 1.0 1 -doRayleigh 1 \n');

end

if Damp

   fprintf(fout,'region 1 -ele 1 2 3 -rayleigh 0.0 0.1 0.0 0.0\n');

end

fprintf(fout,'timeSeries Rectangular 1 0.0 1.0 -factor 1.0e-4\n'); fprintf(fout,'pattern UniformExcitation 1 1 -accel 1\n');

fprintf(fout,'recorder Node -file Data/TestD.txt -time -nodeRange 1 4 -dof 1 disp\n'); fprintf(fout,'recorder Element -file Data/TestF.txt -time -eleRange 1 3 localForce\n');

fprintf(fout,'constraints Penalty 1.0e8 1.0e8\n'); fprintf(fout,'numberer RCM\n'); fprintf(fout,'system ProfileSPD\n'); fprintf(fout,'test EnergyIncr 1.0e-5 30 1\n'); fprintf(fout,'algorithm NewtonLineSearch\n');

fprintf(fout,'analysis Transient\n'); fprintf(fout,'integrator Newmark 0.5 0.25\n'); fprintf(fout,'analyze 500 0.02\n');

%% tic; cd('Tester'); !OpenSees_2_2_2 Driver.tcl toc cd ..

%% Read results raw = dlmread('Tester\Data\TestD.txt',' '); time = raw(:,1); x = raw(:,2:end);

raw = dlmread('Tester\Data\TestD.txt',' '); N = raw(:,2:end);

%% Plot figure(1); %clf; plot(time,x(:,end),'LineWidth',2); hold all <\nowiki>