Hi dear Silvia
please check my model and tell me why committed displacement is zero!
# 1STORY FRAME MODELING
# Written: Seyyed Abbas Hoseyni
# Date : 2009
#################################
model BasicBuilder -ndm 2 -ndf 3;
# Units: kg,m
#################################
#nodal coordinates
#node ntag x y z
node 1 0 0
node 2 4 0
node 3 0 3.2
node 4 4 3.2
#----------------------------
node 5 0 0.2
node 6 0 0.2
node 7 0 3.0
node 8 0 3.0
#----------------------------
node 9 4 0.2
node 10 4 0.2
node 11 4 3.0
node 12 4 3.0
#----------------------------
node 13 0.2 3.2
node 14 0.2 3.2
node 15 3.8 3.2
node 16 3.8 3.2
#################################
#boundary conditions coordinates
#fix ntag TX TY TZ RX RY RZ
fix 1 1 1 1
fix 2 1 1 1
###################################
#nodal masses:
#mass nTag MX MY MZ IX IY IZ
mass 3 14000 0 0
mass 4 14000 0 0
###################################
#transformation:
geomTransf PDelta 1
#geomTransf PDelta 2
###################################
#MATERIALS
set fc 280
#set povason 0.2
set E [expr 15100*pow($fc,0.5)]
#set G [expr $E/(2*(1+$Povason))]
set col1 1
set beam1 2
######################################
#uniaxialMaterial Hysteretic $matTag $s1p $e1p $s2p $e2p $s3p $e3p $s1n $e1n $s2n $e2n $s3n $e3n $pinchX $pinchY $damage1 $damage2 $beta
uniaxialMaterial Hysteretic $col1 102902.0545 0.009545603 122416.8083 0.078248345 0.01 0.178248345 -102902.0545 -0.009545603 -122416.8083 -0.078248345 0.01 -0.178248345 0.1 0.1 0.01 0.01 0
uniaxialMaterial Hysteretic $beam1 17559.01157 0.010467317 21376.74949 0.082591519 0.01 0.182591519 -17559.01157 -0.010467317 -21376.74949 -0.082591519 0.01 -0.182591519 0.1 0.1 0.01 0.01 0
###################################
###################################
set WidthCol1 0.5
set DepthCol1 0.5
###################################
set WidthBeam1 0.5
set DepthBeam1 0.4
###################################
set Acol1 [expr $WidthCol1*$DepthCol1]
#-----------------------
set Abeam1 [expr $WidthBeam1*$DepthBeam1]
###################################
set IzCol1 [expr $WidthCol1*pow($DepthCol1,3)/12]
#-----------------------
set IzBeam1 [expr $WidthBeam1*pow($DepthBeam1,3)/12]
###################################
#Column Elastic Elements
######################################
#element elasticBeamColumn $eleTag $iNode $jNode $A $E $Iz $transfTag
#C1------------------
element elasticBeamColumn 1 1 5 $Acol1 $E $IzCol1 1
element elasticBeamColumn 2 6 7 $Acol1 $E $IzCol1 1
element elasticBeamColumn 3 8 3 $Acol1 $E $IzCol1 1
#C2------------------
element elasticBeamColumn 4 2 9 $Acol1 $E $IzCol1 1
element elasticBeamColumn 5 10 11 $Acol1 $E $IzCol1 1
element elasticBeamColumn 6 12 4 $Acol1 $E $IzCol1 1
############################################################################################
source rotSpring2D.tcl
#rotSpring2D $eleID $nodeR $nodeC $matID
#########################################
#Column RotSpring Elements
#########################################
#C1------------------
#rotSpring2D $eleID $nodeR $nodeC $matID
rotSpring2D 10 5 6 $col1
rotSpring2D 11 7 8 $col1
#C2------------------
rotSpring2D 12 9 10 $col1
rotSpring2D 13 11 12 $col1
########################################################################################################
#Beam Elastic Elements
#########################################################################
#element elasticBeamColumn $eleTag $iNode $jNode $A $E $Iz $transfTag
element elasticBeamColumn 7 3 13 $Abeam1 $E $IzBeam1 1
element elasticBeamColumn 8 13 15 $Abeam1 $E $IzBeam1 1
element elasticBeamColumn 9 16 4 $Abeam1 $E $IzBeam1 1
###################################################################################################
# Beam RotSpring Elements
########################################
#rotSpring2D $eleID $nodeR $nodeC $matID
rotSpring2D 14 13 14 $beam1
rotSpring2D 15 15 16 $beam1
###########################################################
source DisplayPlane.tcl;
source DisplayModel2D.tcl;
set ViewScale 1; # amplify display of deformed shape
DisplayModel2D DeformedShape $ViewScale ; # display deformed shape, the scaling factor needs to be adjusted for each model
puts "Model built!!!"
puts "Congratulations!"
#source Eigen.tcl
#########################################################
# Define gravity loads
# --------------------
# Set Parameters for the Gravity load
#Nodal Loads
set OutNodeLoad 8635
#Distributed Loads
set BotElementLoad 1925
# Create a Plain load pattern with a Linear TimeSeries
pattern Plain 1 "Linear" {
# Create nodal loads
# nd FX FY MZ ; # nd FX FY MZ
# Out Nodes Loads
load 3 0.00 [expr -$OutNodeLoad] 0.00 ; load 4 0.0 [expr -$OutNodeLoad] 0.0
# Create Distributed loads
# BotElementLoad
#eleLoad -ele $eleTag1 <$eleTag2 ....> -type -beamUniform $Wz
eleLoad -ele 7 8 9 -type -beamUniform $BotElementLoad
}
puts "Congratulations!"
constraints Transformation
##numberer RCM
#system BandGeneral
#test NormDispIncr 1.0e-5 100
#algorithm Newton
#integrator LoadControl 0.1
#integrator MinUnbalDispNorm .1
#analysis Static
#analyze 10
#wipeAnalysis
set DmaxPush 1
set Nsteps 100
set DxPush [expr int($DmaxPush/$Nsteps)]
#set PUSHOVER "DispControl"; # run displacement-controlled static pushover analysis
pattern Plain 2 "Linear" {
#load $nodeTag (ndf $LoadValues)
load 4 100.0 0.0 0.0
}
recorder Node -file node4-disp.out -node 4 -dof 1 disp
constraints Plain;
#numberer RCM;
numberer Plain
system BandGeneral
#system BandGeneral
test NormDispIncr 1.0e-2 100 0;
algorithm Newton;
#integrator DisplacementControl $nodeTag $dofTag $dU1 <$Jd $minDu $maxDu>
integrator DisplacementControl 4 1 $DxPush 1 $DxPush $DxPush
analysis Static;
analyze $Nsteps
print node 4
--------------------------------
results
Node: 4
Coordinates : 4 3.2
commitDisps: 0 0 0
unbalanced Load: 0 0 0
Mass :
14000 0 0
0 0 0
0 0 0
ID : 3 4 5
pushover analysis doesn't effect
Moderators: silvia, selimgunay, Moderators
-
- Posts: 5
- Joined: Thu Apr 23, 2009 9:37 am
- Location: Tehran University