-
Notifications
You must be signed in to change notification settings - Fork 22
/
FiberRHSS.tcl
42 lines (39 loc) · 1.56 KB
/
FiberRHSS.tcl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
##################################################################################################################
# FiberRHSS.tcl
#
# SubRoutine to construct a fiber section: Rectangular HSS section
#
##################################################################################################################
#
# Input Arguments:
#------------------
# secID Section ID
# matID Material ID
# d Section depth
# t Tube tickness
# nfdy Number of fibers along depth that goes along local y axis
# nfty Number of fibers along thickness that goes along local y axis
# nfdz Number of fibers along depth that goes along local z axis
# nftz Number of fibers along thickness that goes along local z axis
#
# Written by: Dr. Ahmed Elkady, University of Southampton, UK
#
##################################################################################################################
proc FiberRHSS {secID matID d t nfdy nfty nfdz nftz} {
set dw [expr $d - 2 * $t]
set y1 [expr -$d/2]
set y2 [expr -$dw/2]
set y3 [expr $dw/2]
set y4 [expr $d/2]
set z1 [expr -$d/2]
set z2 [expr -$dw/2]
set z3 [expr $dw/2]
set z4 [expr $d/2]
section fiberSec $secID -GJ 1.e10 {
# nfIJ nfJK yI zI yJ zJ yK zK yL zL
patch quadr $matID $nftz $nfdy $y2 $z4 $y2 $z3 $y3 $z3 $y3 $z4
patch quadr $matID $nftz $nfdy $y2 $z2 $y2 $z1 $y3 $z1 $y3 $z2
patch quadr $matID $nfdz $nfty $y1 $z4 $y1 $z1 $y2 $z1 $y2 $z4
patch quadr $matID $nfdz $nfty $y3 $z4 $y3 $z1 $y4 $z1 $y4 $z4
}
}