View the mode shape of an eigen analysis
Moderators: silvia, selimgunay, Moderators
-
- Posts: 109
- Joined: Sat May 05, 2007 12:28 pm
- Location: Houston, TX
# Procedure to calculate the periods and eigen vectors of a bridge
#
# BGN, Mar 2004, Mostly taken from OSP.
#
proc ModalAnalysis { nModes name } {
#Wipe the analysis
wipeAnalysis
#Perform Eigen Value analysis
set eigenvalues [eigen $nModes]
#Set up recorders for eigen vectors
set jjj 0
for {set jjj 0} {$jjj < $nModes} {incr jjj 1} {
recorder Node -file tempFile[expr {$jjj}].tmp -node all -dof 1 2 3 "eigen [expr $jjj + 1]"
}
#Run a one step gravity load with no loading.
test EnergyIncr 1e-6 20 0
algorithm Newton
system SparseGeneral
integrator LoadControl 1 1 .0000001 1
numberer Plain
constraints Transformation
analysis Static
analyze 1
#Wipe the analysis
wipeAnalysis
#Remove the recorders
remove recorders
#Create .eig file
set _fID [open $name/$name.eig w]
#Open the temp files
for {set jjj 0} {$jjj < $nModes} {incr jjj 1} {
set _fIDt [open tempFile[expr {$jjj}].tmp]
set tempPeriod [expr 2*3.1416/sqrt([lindex $eigenvalues $jjj])]
set tempValues [read -nonewline $_fIDt]
puts $tempPeriod;# $tempValues
puts $_fID "$tempPeriod $tempValues"
close $_fIDt
file delete tempFile[expr {$jjj}].tmp
}
#Close the .eig File
close $_fID
}
I hope it's clear now how to create an ".eig" file now. Take care,
#
# BGN, Mar 2004, Mostly taken from OSP.
#
proc ModalAnalysis { nModes name } {
#Wipe the analysis
wipeAnalysis
#Perform Eigen Value analysis
set eigenvalues [eigen $nModes]
#Set up recorders for eigen vectors
set jjj 0
for {set jjj 0} {$jjj < $nModes} {incr jjj 1} {
recorder Node -file tempFile[expr {$jjj}].tmp -node all -dof 1 2 3 "eigen [expr $jjj + 1]"
}
#Run a one step gravity load with no loading.
test EnergyIncr 1e-6 20 0
algorithm Newton
system SparseGeneral
integrator LoadControl 1 1 .0000001 1
numberer Plain
constraints Transformation
analysis Static
analyze 1
#Wipe the analysis
wipeAnalysis
#Remove the recorders
remove recorders
#Create .eig file
set _fID [open $name/$name.eig w]
#Open the temp files
for {set jjj 0} {$jjj < $nModes} {incr jjj 1} {
set _fIDt [open tempFile[expr {$jjj}].tmp]
set tempPeriod [expr 2*3.1416/sqrt([lindex $eigenvalues $jjj])]
set tempValues [read -nonewline $_fIDt]
puts $tempPeriod;# $tempValues
puts $_fID "$tempPeriod $tempValues"
close $_fIDt
file delete tempFile[expr {$jjj}].tmp
}
#Close the .eig File
close $_fID
}
I hope it's clear now how to create an ".eig" file now. Take care,
Bayram Aygun
Graduate Student, Civil&Env. Eng.
Rice University
Graduate Student, Civil&Env. Eng.
Rice University
Hello Bayram Aygun!
I tried to view mode shapes through OSP for several times but I am not getting it. Can u send an example file that works on OSP to view mode shapes?
to
thanks
I tried to view mode shapes through OSP for several times but I am not getting it. Can u send an example file that works on OSP to view mode shapes?
to
thanks
Last edited by kishor on Mon Dec 15, 2008 9:36 am, edited 1 time in total.
Kishor
Research Assitant
McGill University
Canada
Research Assitant
McGill University
Canada
-
- Posts: 109
- Joined: Sat May 05, 2007 12:28 pm
- Location: Houston, TX
I was using OSP to visualize my bridge-foundation-soil model but I don't feel comfortable with the mode shapes I'm getting with OSP, so I turned to the viewing codes Silvia posted here.
I'm using the latest version (v2.0.0) and I'm able to run the procedures without getting an error, BUT when the figures start appearing my computer crashes down I see "NOT RESPONDING" signs on the figures. Every time I try to visualize the mode shapes OpenSees shuts down. Maybe some of the commands in the procedures are not compatible with the latest version, I don't know.
I'd really appreciate if you could help me out on this,
Thanks,
I'm using the latest version (v2.0.0) and I'm able to run the procedures without getting an error, BUT when the figures start appearing my computer crashes down I see "NOT RESPONDING" signs on the figures. Every time I try to visualize the mode shapes OpenSees shuts down. Maybe some of the commands in the procedures are not compatible with the latest version, I don't know.
I'd really appreciate if you could help me out on this,
Thanks,
Bayram Aygun
Graduate Student, Civil&Env. Eng.
Rice University
Graduate Student, Civil&Env. Eng.
Rice University
-
- Posts: 109
- Joined: Sat May 05, 2007 12:28 pm
- Location: Houston, TX
-
- Posts: 14
- Joined: Sun Nov 16, 2008 6:55 pm
- Location: KU
Code for visualization of modes
[quote="silvia"]yes, but only with the latest version of OS.
[put these procs into a file and then source it.]
Hello,
I did run my model as well as the code mentioned here by Silvia Mazzoni.
I got thee windows: deformed shape, node number and the third one seems to be for the modes. However, they appear in the monitor for a very short time (maybe a couple of dt) and I can't visualize them properly. It is possible to frozen the images or save them in a file?? maybe I'm doing something wrong...
[put these procs into a file and then source it.]
Hello,
I did run my model as well as the code mentioned here by Silvia Mazzoni.
I got thee windows: deformed shape, node number and the third one seems to be for the modes. However, they appear in the monitor for a very short time (maybe a couple of dt) and I can't visualize them properly. It is possible to frozen the images or save them in a file?? maybe I'm doing something wrong...
-
- Posts: 73
- Joined: Tue Feb 20, 2007 6:48 am
- Location: IIEES
- Contact:
-
- Posts: 73
- Joined: Tue Feb 20, 2007 6:48 am
- Location: IIEES
- Contact:
Re: View the mode shape of an eigen analysis
Hi everyone,
I am wondering to know if it is possible to change the color of either nodes or deformed shape for viewing the mode shapes?
I am wondering to know if it is possible to change the color of either nodes or deformed shape for viewing the mode shapes?