Fiber section quad compared fiber

Forum for OpenSees users to post questions, comments, etc. on the use of the OpenSees interpreter, OpenSees.exe

Moderators: silvia, selimgunay, Moderators

Post Reply
Swanrhcp
Posts: 5
Joined: Tue Aug 06, 2019 12:17 pm

Fiber section quad compared fiber

Post by Swanrhcp »

Hi,
I have modelled a 3d structure with STKO and after I have modified .tcl file manually. I have a question: I have modelled Ipe550 section with STKO triangular fiber and I got a different result compared to the case in which the same section has modelled with quad patch model (writing manually section.tcl for opensees). Why? Is It correct? I have copied the code of Wsection.tcl from Wiki for the patch quad. Thx

this is my script:
wipe
model basic -ndm 3 -ndf 6
# timeSeries Constant
timeSeries Constant 2
# timeSeries Linear
timeSeries Linear 3
uniaxialMaterial Steel02 1 500000.0 200000000.0 0.005 20.0 0.925 0.15 0.1 7.0 0.1 7.0 0.0
#ipe550
set b 0.21
set tw 0.0111
set tf 0.0172
set H 0.55
#ipe300
set b_3 0.15
set tw_3 0.0071
set tf_3 0.0107
set H_3 0.30
set nfdw 16
set nftw 4
set nfbf 16
set nftf 4

set dw [expr $H - 2 * $tf]
set y1 [expr -$H/2]
set y2 [expr -$dw/2]
set y3 [expr $dw/2]
set y4 [expr $H/2]

set z1 [expr -$b/2]
set z2 [expr -$tw/2]
set z3 [expr $tw/2]
set z4 [expr $b/2]

section fiberSec 2 {
# nfIJ nfJK yI zI yJ zJ yK zK yL zL
patch quadr 1 $nfbf $nftf $y1 $z4 $y1 $z1 $y2 $z1 $y2 $z4
patch quadr 1 $nftw $nfdw $y2 $z3 $y2 $z2 $y3 $z2 $y3 $z3
patch quadr 1 $nfbf $nftf $y3 $z4 $y3 $z1 $y4 $z1 $y4 $z4
}
#dimensioni ipe300
set dw_3 [expr $H_3 - 2 * $tf_3]
set y1_3 [expr -$H_3/2]
set y2_3 [expr -$dw_3/2]
set y3_3 [expr $dw_3/2]
set y4_3 [expr $H_3/2]

set z1_3 [expr -$b_3/2]
set z2_3 [expr -$tw_3/2]
set z3_3 [expr $tw_3/2]
set z4_3 [expr $b_3/2]

section fiberSec 3 {
# nfIJ nfJK yI zI yJ zJ yK zK yL zL
patch quadr 1 $nfbf $nftf $y1_3 $z4_3 $y1_3 $z1_3 $y2_3 $z1_3 $y2_3 $z4_3
patch quadr 1 $nftw $nfdw $y2_3 $z3_3 $y2_3 $z2_3 $y3_3 $z2_3 $y3_3 $z3_3
patch quadr 1 $nfbf $nftf $y3_3 $z4_3 $y3_3 $z1_3 $y4_3 $z1_3 $y4_3 $z4_3
}
# tag x y z
node 1 20 0 0
node 2 20 0 5
node 3 0 0 0
node 4 5 0 0
node 5 15 0 0
node 6 10 0 0
node 7 10 0 5
node 8 0 0 5
node 9 15 0 5
node 10 5 0 5
node 11 20 0 1.25
node 12 20 0 2.5
node 13 20 0 3.75
node 14 1.25 0 0
node 15 2.5 0 0
node 16 3.75 0 0
node 17 16.25 0 0
node 18 17.5 0 0
node 19 18.75 0 0
node 20 6.25 0 0
node 21 7.5 0 0
node 22 8.75 0 0
node 23 11.25 0 0
node 24 12.5 0 0
node 25 13.75 0 0
node 26 10 0 1.25
node 27 10 0 2.5
node 28 10 0 3.75
node 29 1.25 0 3.75
node 30 2.5 0 2.5
node 31 3.75 0 1.25
node 32 13.75 0 3.75
node 33 12.5 0 2.5
node 34 11.25 0 1.25
node 35 6.25 0 3.75
node 36 7.5 0 2.5
node 37 8.75 0 1.25
node 38 18.75 0 3.75
node 39 17.5 0 2.5
node 40 16.25 0 1.25
node 41 15 0 1.25
node 42 15 0 2.5
node 43 15 0 3.75
node 44 5 0 1.25
node 45 5 0 2.5
node 46 5 0 3.75
node 47 0 0 1.25
node 48 0 0 2.5
node 49 0 0 3.75
node 50 1.25 0 5
node 51 2.5 0 5
node 52 3.75 0 5
node 53 16.25 0 5
node 54 17.5 0 5
node 55 18.75 0 5
node 56 11.25 0 5
node 57 12.5 0 5
node 58 13.75 0 5
node 59 6.25 0 5
node 60 7.5 0 5
node 61 8.75 0 5

geomTransf Linear 1 1.0 0.0 -0.0
element forceBeamColumn 1 1 11 1 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 2 1.0 0.0 -0.0
element forceBeamColumn 2 11 12 2 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 3 1.0 0.0 -0.0
element forceBeamColumn 3 12 13 3 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 4 1.0 0.0 -0.0
element forceBeamColumn 4 13 2 4 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 5 0.0 0.0 1.0
element forceBeamColumn 5 3 14 5 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 6 0.0 0.0 1.0
element forceBeamColumn 6 14 15 6 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 7 0.0 0.0 1.0
element forceBeamColumn 7 15 16 7 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 8 0.0 0.0 1.0
element forceBeamColumn 8 16 4 8 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 9 0.0 0.0 1.0
element forceBeamColumn 9 5 17 9 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 10 0.0 0.0 1.0
element forceBeamColumn 10 17 18 10 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 11 0.0 0.0 1.0
element forceBeamColumn 11 18 19 11 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 12 0.0 0.0 1.0
element forceBeamColumn 12 19 1 12 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 13 0.0 0.0 1.0
element forceBeamColumn 13 4 20 13 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 14 0.0 0.0 1.0
element forceBeamColumn 14 20 21 14 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 15 0.0 0.0 1.0
element forceBeamColumn 15 21 22 15 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 16 0.0 0.0 1.0
element forceBeamColumn 16 22 6 16 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 17 0.0 0.0 1.0
element forceBeamColumn 17 6 23 17 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 18 0.0 0.0 1.0
element forceBeamColumn 18 23 24 18 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 19 0.0 0.0 1.0
element forceBeamColumn 19 24 25 19 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 20 0.0 0.0 1.0
element forceBeamColumn 20 25 5 20 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 21 1.0 0.0 -0.0
element forceBeamColumn 21 6 26 21 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 22 1.0 0.0 -0.0
element forceBeamColumn 22 26 27 22 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 23 1.0 0.0 -0.0
element forceBeamColumn 23 27 28 23 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 24 1.0 0.0 -0.0
element forceBeamColumn 24 28 7 24 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 25 0.7071067811865475 0.0 0.7071067811865476
element forceBeamColumn 25 8 29 25 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 26 0.7071067811865476 0.0 0.7071067811865475
element forceBeamColumn 26 29 30 26 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 27 0.7071067811865476 0.0 0.7071067811865476
element forceBeamColumn 27 30 31 27 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 28 0.7071067811865476 0.0 0.7071067811865476
element forceBeamColumn 28 31 4 28 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 29 -0.7071067811865476 0.0 0.7071067811865476
element forceBeamColumn 29 9 32 29 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 30 -0.7071067811865477 0.0 0.7071067811865475
element forceBeamColumn 30 32 33 30 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 31 -0.7071067811865476 0.0 0.7071067811865476
element forceBeamColumn 31 33 34 31 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 32 -0.7071067811865475 0.0 0.7071067811865477
element forceBeamColumn 32 34 6 32 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 33 0.7071067811865476 0.0 0.7071067811865476
element forceBeamColumn 33 10 35 33 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 34 0.7071067811865477 0.0 0.7071067811865475
element forceBeamColumn 34 35 36 34 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 35 0.7071067811865476 0.0 0.7071067811865476
element forceBeamColumn 35 36 37 35 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 36 0.7071067811865475 0.0 0.7071067811865477
element forceBeamColumn 36 37 6 36 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 37 -0.7071067811865476 0.0 0.7071067811865476
element forceBeamColumn 37 2 38 37 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 38 -0.7071067811865477 0.0 0.7071067811865475
element forceBeamColumn 38 38 39 38 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 39 -0.7071067811865476 0.0 0.7071067811865476
element forceBeamColumn 39 39 40 39 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 40 -0.7071067811865475 0.0 0.7071067811865477
element forceBeamColumn 40 40 5 40 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 41 1.0 0.0 -0.0
element forceBeamColumn 41 5 41 41 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 42 1.0 0.0 -0.0
element forceBeamColumn 42 41 42 42 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 43 1.0 0.0 -0.0
element forceBeamColumn 43 42 43 43 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 44 1.0 0.0 -0.0
element forceBeamColumn 44 43 9 44 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 45 1.0 0.0 -0.0
element forceBeamColumn 45 4 44 45 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 46 1.0 0.0 -0.0
element forceBeamColumn 46 44 45 46 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 47 1.0 0.0 -0.0
element forceBeamColumn 47 45 46 47 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 48 1.0 0.0 -0.0
element forceBeamColumn 48 46 10 48 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 49 1.0 0.0 -0.0
element forceBeamColumn 49 3 47 49 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 50 1.0 0.0 -0.0
element forceBeamColumn 50 47 48 50 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 51 1.0 0.0 -0.0
element forceBeamColumn 51 48 49 51 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 52 1.0 0.0 -0.0
element forceBeamColumn 52 49 8 52 Lobatto 3 5
# Geometric transformation command
geomTransf Linear 53 0.0 0.0 1.0
element forceBeamColumn 53 8 50 53 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 54 0.0 0.0 1.0
element forceBeamColumn 54 50 51 54 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 55 0.0 0.0 1.0
element forceBeamColumn 55 51 52 55 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 56 0.0 0.0 1.0
element forceBeamColumn 56 52 10 56 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 57 0.0 0.0 1.0
element forceBeamColumn 57 9 53 57 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 58 0.0 0.0 1.0
element forceBeamColumn 58 53 54 58 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 59 0.0 0.0 1.0
element forceBeamColumn 59 54 55 59 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 60 0.0 0.0 1.0
element forceBeamColumn 60 55 2 60 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 61 0.0 0.0 1.0
element forceBeamColumn 61 7 56 61 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 62 0.0 0.0 1.0
element forceBeamColumn 62 56 57 62 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 63 0.0 0.0 1.0
element forceBeamColumn 63 57 58 63 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 64 0.0 0.0 1.0
element forceBeamColumn 64 58 9 64 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 65 0.0 0.0 1.0
element forceBeamColumn 65 10 59 65 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 66 0.0 0.0 1.0
element forceBeamColumn 66 59 60 66 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 67 0.0 0.0 1.0
element forceBeamColumn 67 60 61 67 Lobatto 2 5
# Geometric transformation command
geomTransf Linear 68 0.0 0.0 1.0
element forceBeamColumn 68 61 7 68 Lobatto 2 5

# Constraints.sp fix
fix 3 1 1 1 1 0 1
fix 1 0 1 1 1 0 1

# Patterns.addPattern loadPattern
pattern Plain 2 2 {

# Loads.Force NodeForce
load 2 0.0 0.0 -11.5 0.0 0.0 0.0
load 8 0.0 0.0 -11.5 0.0 0.0 0.0
load 7 0.0 0.0 -22.3 0.0 0.0 0.0
load 9 0.0 0.0 -22.3 0.0 0.0 0.0
load 10 0.0 0.0 -22.3 0.0 0.0 0.0
}


file mkdir Risultati

recorder Node -file Risultati/abbassamento_nodo_cont.txt -time -node 7 -dof 3 disp;
recorder Node -file Risultati/Reazio_cern.txt -time -node 3 -dof 3 reaction
recorder Element -file Risultati/ele_globalForce.txt -time -ele 49 globalForce
recorder Element -file Risultati/ele_sectionForce.txt -time -ele 49 section 1 force
recorder Element -file Risultati/ele_sectiondeformation.txt -time -ele 49 section 1 deformation;


# analyses command
constraints Plain
numberer RCM
system BandGeneral
test NormDispIncr 0.0001 40
algorithm Newton
integrator LoadControl 1.0
analysis Static
analyze 1

loadConst -time 0

# Patterns.addPattern loadPattern
pattern Plain 6 3 {

# Loads.Force NodeForce
load 7 0.0 0.0 -10.0 0.0 0.0 0.0
load 9 0.0 0.0 -10.0 0.0 0.0 0.0
load 10 0.0 0.0 -10.0 0.0 0.0 0.0
}

# analyses command
set Maxdisp 0.5
set disp_incr 0.01
set n_step [expr int($Maxdisp/$disp_incr)]

constraints Plain
numberer RCM
system BandGeneral
test NormDispIncr 0.0001 40
algorithm Newton
integrator DisplacementControl 7 3 [expr -$disp_incr]
analysis Static
analyze $n_step

puts "fatto"
Top
Post Reply

Return to “OpenSees.exe Users”