Outline
OpenSees mini-Workshop tutorial
by Patxi Uriz
University of California, Berkeley
Introduction ! Example brace model ! Post earthquake lateral resisting capacity ! Large displacement analysis of frame
!
Brace Modeling
!
Brace Modeling
!
Menegotto-Pinto Steel Material Model Steel02
W Section
Box Section
Pipe Section
Gusset plate Section
Brace Modeling
!
Brace Modeling
!
Bracing Elements
Gusset plate elements
Corotational Transformation { geomTransf Corotational $mtag } j
OpenSees Model
Initial deformation - from allowable tolerances i 2t Rigid element
Brace Modeling
!
Modeling with limited local buckling
Fracture Loss of compression strength from local buckling
Section X-X
Predicted Fracture
Brace Modeling
!
Effect of modeling
Example
!
Hazard Assessment
Goal: Determine post-EQK capacity
! Lateral ! Vertical ! Etc.
Example ([Link])
model BasicBuilder -ndm 2 -ndf 3 ## Nodes node 1 0.0 0.0 ## Mass mass 7 1.0 0.0 0.0 #Boundary Conditions fix 1 1 1 1 ## Materials uniaxialMaterial Steel02 1 40.0 29800.0 0.003 20 \ 0.925 0.15 0.0005 0.01 0.0005 0.01 ## Pipe Section (4" diameter, 1/2" wall thickness) section fiberSec 1 { patch circ 1 20 4 0.0 0.0 2.0 1.5 360.0 0.0 } ## Gusset Plate Section (1" thick, 20" wide)
Example
## Transformation geomTransf Corotational 1 ## Define Model element elasticBeamColumn 1 1 2 1e3 1e8 1e6 1 ;# 'Rigid' offset element nonlinearBeamColumn 2 2 3 3 2 1 ;# Gusset plate element nonlinearBeamColumn 3 3 4 3 1 1 ;# Brace element nonlinearBeamColumn 4 4 5 3 1 1 ;# Brace element nonlinearBeamColumn 5 5 6 3 2 1 ;# Gusset plate element elasticBeamColumn 6 6 7 1e3 1e8 1e6 1 ;# 'Rigid' offset
Example
## Source in the model source [Link] ## Apply the nodal Load pattern Plain 1 Linear { load 7 1.0 0.0 0.0 } ## Recorder recorder Node -file [Link] -time -node 7 -dof 1 disp ## Static Analysis parameters test EnergyIncr 1.0e-8 30 0 algorithm Newton system UmfPack numberer RCM constraints Plain analysis Static set peaks [ list 0.25 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 ] for {set i 1 } { $i <= 9 } {incr i } { set dU [expr -1.0*pow((-1.0),$i)*[lindex $peaks [expr $i-1] ]/100.0 ] integrator DisplacementControl 7 1 $dU 1 $dU $dU analyze 100 }
Example
## Source in the model source [Link] ## Load profile pattern Plain 1 Linear { load 7 1.0 0.0 0.0 } ## Recorder recorder Node -file [Link] -time -node 7 -dof 1 disp ## Static Analysis parameters test EnergyIncr 1.0e-8 30 0 algorithm Newton system UmfPack numberer RCM constraints Plain analysis Static ## Push past buckling set dU -0.01 integrator DisplacementControl 7 1 $dU 1 $dU $dU analyze 300 ## close out model and recorders wipe
Example (cont)
## Find peak compressive capacity set Capacity 0.0 set fid [open [Link] r] foreach ln [split [read $fid] \n] { set tmpMin [lindex $ln 0] if {[llength $ln] == 0} { continue } elseif { $tmpMin < $Capacity } { set Capacity $tmpMin } } close $fid
Exaple (cont)
## Run dynamic analysis source [Link] recorder Node -file [Link] -time -node 7 -dof 1 disp ## Transient Analysis parameters set SineSeries "Sine 0.0 2.0 0.80 -factor 200 " pattern UniformExcitation 1 1 -accel $SineSeries integrator Newmark 0.5 0.25 0.05 0.05 0.0 0.0 test EnergyIncr 1.0e-8 300 0 algorithm KrylovNewton system UmfPack numberer RCM constraints Plain analysis Transient analyze 40000 0.0001 ## Keep structure data and rename the file wipeAnalysis remove loadPattern 1
Example (cont_
## Load profile pattern Plain 1 Linear { load 7 1.0 0.0 0.0 } ## Recorder recorder Node -file [Link] -time -node 7 -dof 1 disp ## Static Analysis parameters test EnergyIncr 1.0e-8 30 0 algorithm Newton system UmfPack numberer RCM constraints Plain analysis Static ## Push past buckling set dU -0.01 integrator DisplacementControl 7 1 $dU 1 $dU $dU analyze 700 ## close out model and recorders wipe
Example (cont)
## Find peak compressive capacity set PostEQKCap 0.0 set fid [open [Link] r] foreach ln [split [read $fid] \n] { set tmpMin [lindex $ln 0] if {[llength $ln] == 0} { continue } elseif { $tmpMin < $PostEQKCap} { set PostEQKCap $tmpMin } } close $fid puts [ format " \n\n\n\n \ Prior to loading the buckling capacity was: %3.4f kips \n \ After dynamic time-history analysis, \n \ the buckling capacity is now : %3.4f kips \n " $Capacity $PostEQKCap ]
Example
!
Hazard Assessment
Lateral pushover, before and after earthquake
Hazard Assessment
Hazard Assessment
P[ DM < dm | IM = im] =
edp
! P[ DM < dm | EDP = edp] dP[ EDP < edp | IM = im]dedp
Reference
!
OpenSees Examples
Uriz, P., and Mahin, S. A. (2004). "Seismic Vulnerability Assessment of Concentrically Braced Steel Frames." International Journal of Steel Structures, 4(4), 239-248.
![Example (cont)
Example (cont)
## Find peak compressive capacity
set Capacity 0.0
set fid [open StaticPush.dat r]
foreach ln [](/p?url=https%3A%2F%2Fscreenshots.scribd.com%2FScribd%2F252_100_85%2F194%2F61782648%2F5.jpeg&__src=https%3A%2F%2Fwww.scribd.com%2Fdocument%2F61782648%2FOpen-Sees-Example&__type=image)
