hi dr fmk
i am sorry, i have a question again
i try to model tmd(tuned mass damper) at top of a building with 3 storey.its including of mass,damping and stiffness.in one node at top of building, i modeled stiffness and dampeing with zerolength.my problem is this that with changeing stiffness(K) and damper(C).i didnt see any change in recorde of zerolength, until i change mass they changed.
please help me
thank you.
may i post the scripts?
set m 1;
set kg 1;
set sec 1;
set sec2 [expr $sec*$sec];
set m2 [expr $m*$m];
set g [expr 9.81*$m/$sec2];
set dataDir Data;
file mkdir $dataDir;
set GMdir "../GMfiles/";
# #########################################################
set Fy1 [expr 35107570.2*$kg/$m2];
set Fy2 [expr 25277450.5*$kg/$m2];
set E [expr 20600000000*$kg/$m2];
set b 0.02;
set steel 1;
set ssteel 2;
set sssteel 3;
set K 47338.2;
set C 10648.8;
set alpha 1;
set Tol 1.0e-8;
constraints Plain;
numberer RCM;
system SparseGeneral;
test NormDispIncr $Tol 10;
algorithm Newton;
integrator LoadControl 0.1;
analysis Static;
analyze 10;
loadConst -time 0.0;
puts "$$$$$$"
wipeAnalysis
constraints Plain;
numberer RCM;
system BandGeneral;
test EnergyIncr 1.0e-8 10 0;
algorithm Newton;
integrator Newmark 0.5 0.25;
analysis Transient;
proc ReadSMDFile {inFilename outFilename dt} {
upvar $dt DT;
if [catch {open $inFilename r} inFileID] {
puts stderr "Cannot open $inFilename for reading"
} else {
set outFileID [open $outFilename w]
set flag 0
foreach line [split [read $inFileID] \n] {
if {[llength $line] == 0} {
continue
} elseif {$flag == 1} {
puts $outFileID $line
} else {
foreach word [split $line] {
if {$flag == 1} {
set DT $word
break
}
if {[string match $word "DT="] == 1} {set flag 1}
}
}
}
close $outFileID;
close $inFileID;
}
}
set GMdirection 1;
set GMfile "tabas.acc" ;
set GMfact 1;
set DtAnalysis [expr 0.02*$sec];
set TmaxAnalysis [expr 20*$sec];
set MpropSwitch 1.0;
set KcurrSwitch 0.0;
set KcommSwitch 1.0;
set KinitSwitch 0.0;
set nEigenI 1;
set nEigenJ 3;
set lambdaN [eigen [expr $nEigenJ]];
set lambdaI [lindex $lambdaN [expr $nEigenI-1]];
set lambdaJ [lindex $lambdaN [expr $nEigenJ-1]];
set omegaI [expr pow($lambdaI,0.5)];
set omegaJ [expr pow($lambdaJ,0.5)];
set alphaM [expr $MpropSwitch*$xDamp*(2*$omegaI*$omegaJ)/($omegaI+$omegaJ)];
set betaKcurr [expr $KcurrSwitch*2.*$xDamp/($omegaI+$omegaJ)];
set betaKcomm [expr $KcommSwitch*2.*$xDamp/($omegaI+$omegaJ)];
set betaKinit [expr $KinitSwitch*2.*$xDamp/($omegaI+$omegaJ)];
rayleigh $alphaM $betaKcurr $betaKinit $betaKcomm;
set IDloadTag 400;
set inFile tabas.acc;
set outFile tabas.th;
ReadSMDFile $inFile $outFile dt;
set GMfatt [expr $g*$GMfact];
set AccelSeries "Series -dt $dt -filePath $outFile -factor $GMfatt";
