spsw
Moderators: silvia, selimgunay, Moderators
Re: spsw
please check my example file and tell me why beamcolumn element not jointing with shell element.thanks
wipe
file mkdir data
model basic -ndm 3 -ndf 6
node 1 0 0 0
node 2 5 0 0
node 3 0 0 5
node 4 5 0 5
puts "node ok"
fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
equalDOF 3 4 1
equalDOF 3 5 1
equalDOF 3 15 1
equalDOF 3 115 1
equalDOF 3 125 1
geomTransf Linear 1 0 -1 0
geomTransf Linear 2 0 1 0
uniaxialMaterial Steel01 1 2.354e8 2e11 0.02
section Fiber 1 {
patch quad 1 6 2 0.0443 0.0275 0.0443 -0.0275 0.05 -0.0275 0.05 0.0275
patch quad 1 2 6 -0.0443 0.00205 -0.0443 -0.00205 0.0443 -0.00205 0.0443 0.00205
patch quad 1 6 2 -0.05 0.0275 -0.05 -0.0275 -0.0443 -0.0275 -0.0443 0.0275
}
set e 2.e11
set nu 0.3
set G [expr $e/(2.*(1+$nu))]
set j 6.920e-8
set Tmattag 2
uniaxialMaterial Elastic $Tmattag [expr $G*$j]
section Aggregator 4 $Tmattag T -section 1
## beam and column element
element nonlinearBeamColumn 1 1 3 10 4 2
element nonlinearBeamColumn 2 2 4 10 4 2
element nonlinearBeamColumn 3 3 4 10 4 2
puts "element ok"
####steel shear wall element
set t1 0.0003
nDMaterial ElasticIsotropic 5 2.1e7 0.3
nDMaterial PlateFiber 6 5
section PlateFiber 7 6 $t1
block2D 10 10 5 4 ShellMITC4 7 {
1 0. 0. 0.
2 5. 0. 0.
3 5. 0. 5.
4 0. 0. 5.
}
puts "section ok"
recorder display model 700 100 500 500 -wipe
prp 5 5 1
vup 0 1 0
vpn 0 0 1
display 1 2 10
recorder Element -file data/eleoutput2.txt -time -ele 2 globalForce
recorder Element -file data/eleoutput10.txt -time -ele 10 globalForce
recorder Node -file data/node4disp.out -time -node 4 -dof 1 2 3 disp
recorder Node -file data/node10disp.out -time -node 10 -dof 1 2 3 disp
recorder Node -file data/node1reac.out -time -node 1 -dof 1 2 3 reaction
recorder Node -file data/node2reac.out -time -node 2 -dof 1 2 3 reaction
puts "recorder ok"
set PushOverLoadPattern 2
pattern Plain $PushOverLoadPattern Linear {
eleLoad -ele 2 -type -beamUniform 0. 0. -4500.
load 3 6000. 0. 0. 0. 0. 0.
}
puts "load ok"
constraints Plain
numberer RCM
system BandGeneral
test NormUnbalance 1.e-6 100
algorithm Newton
integrator DisplacementControl 4 1 0.001
analysis Static
analyze 100
loadConst -time 0.0
wipe
file mkdir data
model basic -ndm 3 -ndf 6
node 1 0 0 0
node 2 5 0 0
node 3 0 0 5
node 4 5 0 5
puts "node ok"
fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 1 1
equalDOF 3 4 1
equalDOF 3 5 1
equalDOF 3 15 1
equalDOF 3 115 1
equalDOF 3 125 1
geomTransf Linear 1 0 -1 0
geomTransf Linear 2 0 1 0
uniaxialMaterial Steel01 1 2.354e8 2e11 0.02
section Fiber 1 {
patch quad 1 6 2 0.0443 0.0275 0.0443 -0.0275 0.05 -0.0275 0.05 0.0275
patch quad 1 2 6 -0.0443 0.00205 -0.0443 -0.00205 0.0443 -0.00205 0.0443 0.00205
patch quad 1 6 2 -0.05 0.0275 -0.05 -0.0275 -0.0443 -0.0275 -0.0443 0.0275
}
set e 2.e11
set nu 0.3
set G [expr $e/(2.*(1+$nu))]
set j 6.920e-8
set Tmattag 2
uniaxialMaterial Elastic $Tmattag [expr $G*$j]
section Aggregator 4 $Tmattag T -section 1
## beam and column element
element nonlinearBeamColumn 1 1 3 10 4 2
element nonlinearBeamColumn 2 2 4 10 4 2
element nonlinearBeamColumn 3 3 4 10 4 2
puts "element ok"
####steel shear wall element
set t1 0.0003
nDMaterial ElasticIsotropic 5 2.1e7 0.3
nDMaterial PlateFiber 6 5
section PlateFiber 7 6 $t1
block2D 10 10 5 4 ShellMITC4 7 {
1 0. 0. 0.
2 5. 0. 0.
3 5. 0. 5.
4 0. 0. 5.
}
puts "section ok"
recorder display model 700 100 500 500 -wipe
prp 5 5 1
vup 0 1 0
vpn 0 0 1
display 1 2 10
recorder Element -file data/eleoutput2.txt -time -ele 2 globalForce
recorder Element -file data/eleoutput10.txt -time -ele 10 globalForce
recorder Node -file data/node4disp.out -time -node 4 -dof 1 2 3 disp
recorder Node -file data/node10disp.out -time -node 10 -dof 1 2 3 disp
recorder Node -file data/node1reac.out -time -node 1 -dof 1 2 3 reaction
recorder Node -file data/node2reac.out -time -node 2 -dof 1 2 3 reaction
puts "recorder ok"
set PushOverLoadPattern 2
pattern Plain $PushOverLoadPattern Linear {
eleLoad -ele 2 -type -beamUniform 0. 0. -4500.
load 3 6000. 0. 0. 0. 0. 0.
}
puts "load ok"
constraints Plain
numberer RCM
system BandGeneral
test NormUnbalance 1.e-6 100
algorithm Newton
integrator DisplacementControl 4 1 0.001
analysis Static
analyze 100
loadConst -time 0.0
-
- Posts: 115
- Joined: Mon Feb 20, 2017 1:19 am
- Location: RMIT University, Melbourne, Australia
Re: spsw
I wish if I could help you, but I have never come across such a problem. You may post your question in researchgate.net, so you will have more people to view your modelling problem. Good luck.
-
- Posts: 2
- Joined: Sat Nov 11, 2017 11:15 am
Re: spsw
Hi did you find any solution for your problem?
I am looking for some information about SPSW too in 3D. I would be helpful to know if you found something already.
Thanks in advance
I am looking for some information about SPSW too in 3D. I would be helpful to know if you found something already.
Thanks in advance