steel box section

If you have a script you think might be useful to others post it
here. Hopefully we will be able to get the most useful of these incorporated in the manuals.

Moderators: silvia, selimgunay, Moderators

Post Reply
basir
Posts: 32
Joined: Mon Apr 13, 2009 9:55 pm
Location: IRAN university of science an technology

steel box section

Post by basir » Mon Apr 20, 2009 7:55 am

hi
how can i create a steel box section
or is there any ready tcl file like Wsection

silvia
Posts: 3912
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia » Fri Apr 24, 2009 9:20 am

modify Wsection to what you need.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104

saeedwizard
Posts: 4
Joined: Sun Mar 14, 2010 9:16 am
Location: KNtu

Post by saeedwizard » Thu Mar 18, 2010 9:04 am

[quote="silvia"]modify Wsection to what you need.[/quote]

How can i change it ?

in w sections we have 1 web in center , but in box sections ... .

PLZ help

more info please :cry:

fmk
Site Admin
Posts: 5866
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Post by fmk » Thu Mar 18, 2010 12:04 pm

it's only 4 patch commands instead of 3!!

:evil: here:

[code]
proc BoxSection {secID matID D B T nfTBl nfTBt nfSl nfSt} {
# ###################################################################
# Wsection $secID $matID $d $bf $tf $tw $nfTBl $nfTBt $nfSl $nfSt
# ###################################################################
# create a simple rectangular section given the nominal section properties
# WARNING: does not take into account rounded corners.
#
# written: Frank McKenna
#
# input parameters:
# secID - section ID number
# matID - material ID number
# H = nominal height
# B = nominal breadth
# T = thickness
# nfTBl = number of fibers along Top/Bottom length
# nfTBt = number of fibers through Top/Bottom thickness
# nfSl = number of fibers along sides length
# nfSt = number of fibers through sides thickness

set y1 [expr $D/2.0]
set y2 [expr $y1 - $T]

set z1 [expr $B/2.0]
set z2 [expr $z1-$T]

section fiberSec $secID {
# nfIJ nfJK yI zI yJ zJ yK zK yL zL
# top and bottom
patch quadr $matID $nfTBl $nfTBt $y2 $z1 $y2 -$z1 $y1 -$z1 $y1 $z1
patch quadr $matID $nfTBl $nfTBt -$y1 $z1 -$y1 -$z1 -$y2 -$z1 -$y2 $z1

# sides
patch quadr $matID $nfSt $nfSl -$y2 $z1 -$y2 $z2 $y2 $z2 $y2 $z1
patch quadr $matID $nfSt $nfSl -$y2 -$z2 -$y2 -$z1 $y2 -$z1 $y2 -$z2
}
}

[/code]

and here is a test:

[code]
# Simple Elastic Column Test of Box Section
# Written: fmk

set P 10.0
set L 100.0
set E 29000.0

set D 2.0
set B 1.0
set T 0.1
set nfTl 1
set nfTt 10
set nfSl 20
set nfSt 1

set y1 [expr $D]
set y2 [expr $D-2.0*$T]
set z1 [expr $B]
set z2 [expr $B-2.0*$T]

# theoretical displacement calc
set I [expr $z1*$y1*$y1*$y1/12.0 - $z2*$y2*$y2*$y2/12.0] ;
set dispT [expr $P*$L*$L*$L/(3.0*$E*$I)];

# opensees model
model BasicBuilder -ndm 2 -ndf 3

node 1 0.0 0.0
node 2 0.0 $L

fix 1 1 1 1

uniaxialMaterial Elastic 1 $E

source BoxSection.tcl
BoxSection 1 1 $D $B $T $nfTl $nfTt $nfSl $nfSt

geomTransf Linear 1
element forceBeamColumn 1 1 2 5 1 1

pattern Plain 1 "Linear" {
load 2 $P 0.0 0.0
}

integrator LoadControl 1
test NormDispIncr 1.0e-12 3 1
algorithm Newton
system ProfileSPD
constraints Plain
numberer Plain
analysis Static
analyze 1

set dispM [nodeDisp 2 1]
puts "Theoretical: $dispT Model: $dispM"

[/code]

saeedwizard
Posts: 4
Joined: Sun Mar 14, 2010 9:16 am
Location: KNtu

Post by saeedwizard » Thu Mar 18, 2010 10:18 pm

wow

tnx, verrry much

Post Reply