Problem
Moderators: silvia, selimgunay, Moderators
-
- Posts: 93
- Joined: Mon Apr 14, 2008 10:53 pm
Problem
Hello
when i run my script, opensees.exe goes out. would you please tell me what is the problem?
tnx
--------------------------------
#Units: N,m
wipe
model BasicBuilder -ndm 3 -ndf 6
# nodal coordinates:
node 1 0 0 6.1
node 2 12.19 0 4.55
node 3 6.095 10.55684967 4.55
node 4 -6.095 10.55684967 4.55
node 5 -12.19 0 4.55
node 6 -6.095 -10.55684967 4.55
node 7 6.095 -10.55684967 4.55
node 8 24.38 0 0
node 9 12.19 21.11369934 0
node 10 -12.19 21.11369934 0
node 11 -24.38 0 0
node 12 -12.19 -21.11369934 0
node 13 12.19 -21.11369934 0
# ---------------------------------------------
# nodal restraint:
fix 8 1 1 1 1 1 1
fix 9 1 1 1 1 1 1
fix 10 1 1 1 1 1 1
fix 11 1 1 1 1 1 1
fix 12 1 1 1 1 1 1
fix 13 1 1 1 1 1 1
# ---------------------------------------------
nDMaterial ElasticIsotropic 1 2E11 0.3
set ys "-VM"
set ps "-VM"
set eps "-NOS 1 1"
nDMaterial Template3Dep 2 1 -YS $ys -PS $ps -EPS $eps
# ---------------------------------------------
geomTransf Corotational 1 1 1 0
# ---------------------------------------------
# Define section:
set Width 1
set Depth 1.279
# some variables derived from the parameters
set y1 [expr $Depth/2.0]
set z1 [expr $Width/2.0]
section Fiber 1 {
patch rect 1 10 10 [expr -$y1] [expr -$z1] $y1 $z1
}
# Add Torsion property to existing sections:
section Aggregator 2 1 T -section 1
# ---------------------------------------------
# element:
element nonlinearBeamColumn 1 1 2 5 2 2
element nonlinearBeamColumn 2 1 3 5 2 2
element nonlinearBeamColumn 3 1 4 5 2 2
element nonlinearBeamColumn 4 1 5 5 2 2
element nonlinearBeamColumn 5 1 6 5 2 2
element nonlinearBeamColumn 6 1 7 5 2 2
element nonlinearBeamColumn 7 2 3 5 2 2
element nonlinearBeamColumn 8 3 4 5 2 2
element nonlinearBeamColumn 9 4 5 5 2 2
element nonlinearBeamColumn 10 5 6 5 2 2
element nonlinearBeamColumn 11 6 7 5 2 2
element nonlinearBeamColumn 12 7 2 5 2 2
element nonlinearBeamColumn 13 2 8 5 2 2
element nonlinearBeamColumn 14 3 9 5 2 2
element nonlinearBeamColumn 15 4 10 5 2 2
element nonlinearBeamColumn 16 5 11 5 2 2
element nonlinearBeamColumn 17 6 12 5 2 2
element nonlinearBeamColumn 18 7 13 5 2 2
#-----------------------------------------------
#load
pattern Plain 1 Linear {
load 1 0 0 -21.4E6 0 0 0
load 2 0 0 -21.4E6 0 0 0
load 3 0 0 -21.4E6 0 0 0
load 4 0 0 -21.4E6 0 0 0
load 5 0 0 -21.4E6 0 0 0
load 6 0 0 -21.4E6 0 0 0
load 7 0 0 -21.4E6 0 0 0
}
#-----------------------------------------------
recorder Node -file node2_dome.out -node 2 -time -dof 3 disp
recorder Element -time -file force1.out -ele 1 force
#----------------------------------------------------------
constraints Transformation
numberer RCM
system BandGeneral
test NormDispIncr 1.0e-6 30
algorithm Newton
integrator LoadControl 0.01
analysis Static
analyze 100
when i run my script, opensees.exe goes out. would you please tell me what is the problem?
tnx
--------------------------------
#Units: N,m
wipe
model BasicBuilder -ndm 3 -ndf 6
# nodal coordinates:
node 1 0 0 6.1
node 2 12.19 0 4.55
node 3 6.095 10.55684967 4.55
node 4 -6.095 10.55684967 4.55
node 5 -12.19 0 4.55
node 6 -6.095 -10.55684967 4.55
node 7 6.095 -10.55684967 4.55
node 8 24.38 0 0
node 9 12.19 21.11369934 0
node 10 -12.19 21.11369934 0
node 11 -24.38 0 0
node 12 -12.19 -21.11369934 0
node 13 12.19 -21.11369934 0
# ---------------------------------------------
# nodal restraint:
fix 8 1 1 1 1 1 1
fix 9 1 1 1 1 1 1
fix 10 1 1 1 1 1 1
fix 11 1 1 1 1 1 1
fix 12 1 1 1 1 1 1
fix 13 1 1 1 1 1 1
# ---------------------------------------------
nDMaterial ElasticIsotropic 1 2E11 0.3
set ys "-VM"
set ps "-VM"
set eps "-NOS 1 1"
nDMaterial Template3Dep 2 1 -YS $ys -PS $ps -EPS $eps
# ---------------------------------------------
geomTransf Corotational 1 1 1 0
# ---------------------------------------------
# Define section:
set Width 1
set Depth 1.279
# some variables derived from the parameters
set y1 [expr $Depth/2.0]
set z1 [expr $Width/2.0]
section Fiber 1 {
patch rect 1 10 10 [expr -$y1] [expr -$z1] $y1 $z1
}
# Add Torsion property to existing sections:
section Aggregator 2 1 T -section 1
# ---------------------------------------------
# element:
element nonlinearBeamColumn 1 1 2 5 2 2
element nonlinearBeamColumn 2 1 3 5 2 2
element nonlinearBeamColumn 3 1 4 5 2 2
element nonlinearBeamColumn 4 1 5 5 2 2
element nonlinearBeamColumn 5 1 6 5 2 2
element nonlinearBeamColumn 6 1 7 5 2 2
element nonlinearBeamColumn 7 2 3 5 2 2
element nonlinearBeamColumn 8 3 4 5 2 2
element nonlinearBeamColumn 9 4 5 5 2 2
element nonlinearBeamColumn 10 5 6 5 2 2
element nonlinearBeamColumn 11 6 7 5 2 2
element nonlinearBeamColumn 12 7 2 5 2 2
element nonlinearBeamColumn 13 2 8 5 2 2
element nonlinearBeamColumn 14 3 9 5 2 2
element nonlinearBeamColumn 15 4 10 5 2 2
element nonlinearBeamColumn 16 5 11 5 2 2
element nonlinearBeamColumn 17 6 12 5 2 2
element nonlinearBeamColumn 18 7 13 5 2 2
#-----------------------------------------------
#load
pattern Plain 1 Linear {
load 1 0 0 -21.4E6 0 0 0
load 2 0 0 -21.4E6 0 0 0
load 3 0 0 -21.4E6 0 0 0
load 4 0 0 -21.4E6 0 0 0
load 5 0 0 -21.4E6 0 0 0
load 6 0 0 -21.4E6 0 0 0
load 7 0 0 -21.4E6 0 0 0
}
#-----------------------------------------------
recorder Node -file node2_dome.out -node 2 -time -dof 3 disp
recorder Element -time -file force1.out -ele 1 force
#----------------------------------------------------------
constraints Transformation
numberer RCM
system BandGeneral
test NormDispIncr 1.0e-6 30
algorithm Newton
integrator LoadControl 0.01
analysis Static
analyze 100
-
- Posts: 93
- Joined: Mon Apr 14, 2008 10:53 pm
-
- Posts: 93
- Joined: Mon Apr 14, 2008 10:53 pm