problems of modal analysis
Moderators: silvia, selimgunay, Moderators
-
- Posts: 9
- Joined: Wed Jun 17, 2020 6:54 pm
problems of modal analysis
Hello,I am building a model of four span continuous beam bridge recently.When I run the analysis of eigen,some errors has appeared. I have put my tcl.file here,can you give me some suggestions of this model, tank you !
OpenSees > # modify2.tcl
OpenSees > logFile "modify2.log"
OpenSees > set systemTime [clock seconds]
1600222769
OpenSees >
OpenSees > puts "==============================================================================="
===============================================================================
OpenSees > puts "OpenSees Navigator: Version 2.5.7"
OpenSees Navigator: Version 2.5.7
OpenSees > puts "File sourced: E:/µ²¿é/È«ÇÅÄ£Äâ/1/modify/modify2.tcl"
File sourced: E:/碌虏驴茅/脠芦脟脜脛拢脛芒/1/modify/modify2.tcl
OpenSees > puts "Starting OpenSees analysis: [clock format $systemTime -format "%d-%b-%Y %H:%M:%S"]"
Starting OpenSees analysis: 16-Sep-2020 10:19:29
OpenSees > puts "==============================================================================="
===============================================================================
OpenSees >
OpenSees > puts "\n########################"
########################
OpenSees > puts "# Analysis-Sequence 1 #"
# Analysis-Sequence 1 #
OpenSees > puts "########################"
########################
OpenSees >
OpenSees > # Start timing of this analysis sequence
OpenSees > # --------------------------------------
OpenSees > set tStart [clock clicks -milliseconds]
1600222769996
OpenSees >
OpenSees > # Start of model generation
OpenSees > # =========================
OpenSees >
OpenSees > # Create ModelBuilder
OpenSees > # -------------------
OpenSees > model BasicBuilder -ndm 3 -ndf 6
OpenSees >
OpenSees > # Define geometry
OpenSees > # ---------------
OpenSees > source NodeCoord.tcl
OpenSees >
OpenSees > # Define Single Point Constraints
OpenSees > # -------------------------------
OpenSees > source SPConstraint.tcl
42 43 44 45 46 47
OpenSees >
OpenSees > # Define nodal masses
OpenSees > # -------------------
OpenSees > source NodeMass.tcl
OpenSees >
OpenSees > # Define Multi Point Constraints
OpenSees > # ------------------------------
OpenSees > source MPConstraint.tcl
OpenSees >
OpenSees > # Define material(s)
OpenSees > # ------------------
OpenSees > source Materials.tcl
OpenSees >
OpenSees > # Define friction model(s)
OpenSees > # ------------------------
OpenSees > source FrictionModels.tcl
OpenSees >
OpenSees > # Define section(s)
OpenSees > # -----------------
OpenSees > source Sections.tcl
OpenSees >
OpenSees > # Define geometric transformation(s)
OpenSees > # ----------------------------------
OpenSees > source GeoTran.tcl
OpenSees >
OpenSees > # Define element(s)
OpenSees > # -----------------
OpenSees > source Elements.tcl
expected floating-point number but got "-shearDist"
OpenSees >
OpenSees > # Define damping parameters
OpenSees > # -------------------------
OpenSees >
OpenSees > # Define time series
OpenSees > # ------------------
OpenSees > source TimeSeries.tcl
OpenSees >
OpenSees > # Start of analysis generation
OpenSees > # ============================
OpenSees >
OpenSees > # Get Initial Stiffness
OpenSees > # ---------------------
OpenSees > initialize
OpenSees >
OpenSees > puts "o Analysis: EigenDefaultCase"
o Analysis: EigenDefaultCase
OpenSees > # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
OpenSees >
OpenSees > # Define load pattern
OpenSees > # -------------------
OpenSees > source LoadPattern_1.tcl
OpenSees >
OpenSees > # Define recorder(s)
OpenSees > # ------------------
OpenSees > source Recorder_1.tcl
2
OpenSees >
OpenSees > # Define analysis options
OpenSees > # -----------------------
OpenSees > source AnalysisOptn_1.tcl
OpenSees >
OpenSees > # Analyze model
OpenSees > # -------------
OpenSees > set pi [expr acos(-1.0)]
3.141592653589793
OpenSees > set eigFID [open EigenDefaultCase_Node_EigenVectors_EigenVal.out w]
file91a820
OpenSees > set lambda [eigen -fullGenLapack 3]
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 72
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 93
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 103
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 105
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 110
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 112
non-varyng assumed
-479898.10073334811022505164 -678.40598248962487559766 5.30501512295851362921
OpenSees > puts $eigFID " lambda omega period frequency"
OpenSees > foreach lambda $lambda {
set omega [expr sqrt($lambda)]
set period [expr 2.0*$pi/$omega]
set frequ [expr 1.0/$period]
puts $eigFID [format " %+1.6e %+1.6e %+1.6e %+1.6e" $lambda $omega $period $frequ]
}
domain error: argument not in valid range
OpenSees > close $eigFID
OpenSees >
OpenSees > # Record eigenvectors
OpenSees > # -------------------
OpenSees > record
OpenSees >
OpenSees > # Stop timing of this analysis sequence
OpenSees > # -------------------------------------
OpenSees > set tStop [clock clicks -milliseconds]
1600222770588
OpenSees > puts "o Time taken: [expr ($tStop-$tStart)/1000.0] sec"
o Time taken: 0.592 sec
OpenSees >
OpenSees > # Clean up
OpenSees > # --------
OpenSees > wipe
OpenSees > set systemTime [clock seconds]
1600222770
OpenSees > puts "\n==============================================================================="
===============================================================================
OpenSees > puts "Finished OpenSees analysis: [clock format $systemTime -format "%d-%b-%Y %H:%M:%S"]"
Finished OpenSees analysis: 16-Sep-2020 10:19:30
OpenSees > # modify2.tcl
OpenSees > logFile "modify2.log"
OpenSees > set systemTime [clock seconds]
1600222769
OpenSees >
OpenSees > puts "==============================================================================="
===============================================================================
OpenSees > puts "OpenSees Navigator: Version 2.5.7"
OpenSees Navigator: Version 2.5.7
OpenSees > puts "File sourced: E:/µ²¿é/È«ÇÅÄ£Äâ/1/modify/modify2.tcl"
File sourced: E:/碌虏驴茅/脠芦脟脜脛拢脛芒/1/modify/modify2.tcl
OpenSees > puts "Starting OpenSees analysis: [clock format $systemTime -format "%d-%b-%Y %H:%M:%S"]"
Starting OpenSees analysis: 16-Sep-2020 10:19:29
OpenSees > puts "==============================================================================="
===============================================================================
OpenSees >
OpenSees > puts "\n########################"
########################
OpenSees > puts "# Analysis-Sequence 1 #"
# Analysis-Sequence 1 #
OpenSees > puts "########################"
########################
OpenSees >
OpenSees > # Start timing of this analysis sequence
OpenSees > # --------------------------------------
OpenSees > set tStart [clock clicks -milliseconds]
1600222769996
OpenSees >
OpenSees > # Start of model generation
OpenSees > # =========================
OpenSees >
OpenSees > # Create ModelBuilder
OpenSees > # -------------------
OpenSees > model BasicBuilder -ndm 3 -ndf 6
OpenSees >
OpenSees > # Define geometry
OpenSees > # ---------------
OpenSees > source NodeCoord.tcl
OpenSees >
OpenSees > # Define Single Point Constraints
OpenSees > # -------------------------------
OpenSees > source SPConstraint.tcl
42 43 44 45 46 47
OpenSees >
OpenSees > # Define nodal masses
OpenSees > # -------------------
OpenSees > source NodeMass.tcl
OpenSees >
OpenSees > # Define Multi Point Constraints
OpenSees > # ------------------------------
OpenSees > source MPConstraint.tcl
OpenSees >
OpenSees > # Define material(s)
OpenSees > # ------------------
OpenSees > source Materials.tcl
OpenSees >
OpenSees > # Define friction model(s)
OpenSees > # ------------------------
OpenSees > source FrictionModels.tcl
OpenSees >
OpenSees > # Define section(s)
OpenSees > # -----------------
OpenSees > source Sections.tcl
OpenSees >
OpenSees > # Define geometric transformation(s)
OpenSees > # ----------------------------------
OpenSees > source GeoTran.tcl
OpenSees >
OpenSees > # Define element(s)
OpenSees > # -----------------
OpenSees > source Elements.tcl
expected floating-point number but got "-shearDist"
OpenSees >
OpenSees > # Define damping parameters
OpenSees > # -------------------------
OpenSees >
OpenSees > # Define time series
OpenSees > # ------------------
OpenSees > source TimeSeries.tcl
OpenSees >
OpenSees > # Start of analysis generation
OpenSees > # ============================
OpenSees >
OpenSees > # Get Initial Stiffness
OpenSees > # ---------------------
OpenSees > initialize
OpenSees >
OpenSees > puts "o Analysis: EigenDefaultCase"
o Analysis: EigenDefaultCase
OpenSees > # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
OpenSees >
OpenSees > # Define load pattern
OpenSees > # -------------------
OpenSees > source LoadPattern_1.tcl
OpenSees >
OpenSees > # Define recorder(s)
OpenSees > # ------------------
OpenSees > source Recorder_1.tcl
2
OpenSees >
OpenSees > # Define analysis options
OpenSees > # -----------------------
OpenSees > source AnalysisOptn_1.tcl
OpenSees >
OpenSees > # Analyze model
OpenSees > # -------------
OpenSees > set pi [expr acos(-1.0)]
3.141592653589793
OpenSees > set eigFID [open EigenDefaultCase_Node_EigenVectors_EigenVal.out w]
file91a820
OpenSees > set lambda [eigen -fullGenLapack 3]
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 72
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 93
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 103
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 105
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 110
non-varyng assumed
WARNING PlainHandler::handle() - constraint matrix not identity, ignoring constraint for node 112
non-varyng assumed
-479898.10073334811022505164 -678.40598248962487559766 5.30501512295851362921
OpenSees > puts $eigFID " lambda omega period frequency"
OpenSees > foreach lambda $lambda {
set omega [expr sqrt($lambda)]
set period [expr 2.0*$pi/$omega]
set frequ [expr 1.0/$period]
puts $eigFID [format " %+1.6e %+1.6e %+1.6e %+1.6e" $lambda $omega $period $frequ]
}
domain error: argument not in valid range
OpenSees > close $eigFID
OpenSees >
OpenSees > # Record eigenvectors
OpenSees > # -------------------
OpenSees > record
OpenSees >
OpenSees > # Stop timing of this analysis sequence
OpenSees > # -------------------------------------
OpenSees > set tStop [clock clicks -milliseconds]
1600222770588
OpenSees > puts "o Time taken: [expr ($tStop-$tStart)/1000.0] sec"
o Time taken: 0.592 sec
OpenSees >
OpenSees > # Clean up
OpenSees > # --------
OpenSees > wipe
OpenSees > set systemTime [clock seconds]
1600222770
OpenSees > puts "\n==============================================================================="
===============================================================================
OpenSees > puts "Finished OpenSees analysis: [clock format $systemTime -format "%d-%b-%Y %H:%M:%S"]"
Finished OpenSees analysis: 16-Sep-2020 10:19:30
-
- Posts: 917
- Joined: Mon Sep 09, 2013 8:50 pm
- Location: University of California, Berkeley
Re: problems of modal analysis
You may need to use constraints Transformation.
Are you using equaldof, rigid diaphragm or similar multi-dof constraints?
Are you using equaldof, rigid diaphragm or similar multi-dof constraints?
-
- Posts: 9
- Joined: Wed Jun 17, 2020 6:54 pm
Re: problems of modal analysis
thanks for your suggestion,actually I am using rigid links to simulate the connection between supscture and bearing.
-
- Posts: 9
- Joined: Wed Jun 17, 2020 6:54 pm
Re: problems of modal analysis
thanks for your suggestion,actually I am using rigid links to simulate the connection between supscture and bearing.selimgunay wrote: ↑Tue Sep 15, 2020 10:16 pm You may need to use constraints Transformation.
Are you using equaldof, rigid diaphragm or similar multi-dof constraints?
Re: problems of modal analysis
Regardless, use constraints Transformation as Selim suggested.
-
- Posts: 9
- Joined: Wed Jun 17, 2020 6:54 pm