fiber section cant creat.why?

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

Moderators: silvia, selimgunay, Moderators

Post Reply
amceuma
Posts: 23
Joined: Sat May 26, 2007 4:48 am
Location: iran

fiber section cant creat.why?

Post by amceuma »

dear friends, why this proc for RC fiber section cant create?

#Define sections----------------------------------------------------------------------------------------------------------------------------
proc buildrcrectsection {ID H B cover $concrete $steel numbartop1 numbartop2 numbarbot1 numbarbot2 numbarint barareatop1 barareatop2 barareabot1 barareabot2 barareaint nfy nfz} {
set covery [expr $H/2];
set coverz [expr $B/2];
set corey [expr $covery-$cover];
set corez [expr $coverz-$cover];
set locaddedbar [expr $corez/2];
section fiber $ID {
patch quadr $concrete 2 20 -$corey $corez -$corey -$corez $corey -$corez $corey $corez;
patch quadr $concrete 20 20 -$covery $coverz -$corey $corez $corey $corez $covery $coverz;
patch quadr $concrete 2 20 -$corey -$corez -$covery -$coverz $covery -$coverz $corey $corez;
patch quadr $concrete 20 2 -$covery $coverz -$covery -$coverz -$corey -$corez -$corey $corez;
patch quadr $concrete 20 2 $corey $corez $corey -$corez $covery -$coverz $covery $coverz;

layer straight $steel $numbarint $barareaint -$corey $corez $corey $corez;
layer straight $steel $numbarint $barareaint -$corey -$corez $corey -$corez;
layer straight $steel $numbartop1 $barareatop1 $corey $corez $corey -$corez;
layer straight $steel $numbartop2 $barareatop2 $corey $locaddedbar $corey -$locaddedbar;
layer straight $steel $numbarbot1 $barareabot1 -$corey $corez -$corey -$corez;
layer straight $steel $numbarbot2 $barareabot2 -$corey $locaddedbar -$corey -$locaddedbar;
};
};

----------------------------------------------------------------------------------
#Build model-----------------------------------------------------------------------------------------------------------------------------
wipe
model BasicBuilder -ndm 2 -ndf 3;
source libunits.tcl;
source DisplayPlane.tcl; # procedure for displaying a plane in model
source DisplayModel2D.tcl; # procedure for displaying 2D perspective of model
source buildrcrectsection.tcl;

#Define Geometry;
set Lcol [expr 3*$m];
set Lbeam [expr 5*$m];
set Nstory 4;
set Nbay 3;

#Define nodal coordinates
for {set level 1} {$level <= [expr $Nstory+1]} {incr level 1} {
set y [expr ($level-1)*$Lcol];
for {set pier 1} {$pier <= [expr $Nbay+1]} {incr pier 1} {
set x [expr ($pier-1)*$Lbeam];
set nodeID [expr $level*10+$pier];
node $nodeID $x $y;
}
};
fix 11 1 1 1;
fix 12 1 1 1;
fix 13 1 1 1;
fix 14 1 1 1;

#Define sections
source libmaterials.tcl;
buildrcrectsection 44 [expr 35*$cm] [expr 25*$cm] [expr 5.3*$cm] 1 2 3 0 2 0 0 [expr 2*$cm2] 0 [expr 2*$cm2] 0 0 20 20;
buildrcrectsection 33 [expr 35*$cm] [expr 25*$cm] [expr 5.3*$cm] 1 2 4 0 3 0 0 [expr 2*$cm2] 0 [expr 2*$cm2] 0 0 20 20;
buildrcrectsection 22 [expr 35*$cm] [expr 25*$cm] [expr 5.3*$cm] 1 2 2 2 3 0 0 [expr 2*$cm2] [expr 3.14*$cm2] [expr 2*$cm2] 0 0 20 20;
buildrcrectsection 11 [expr 35*$cm] [expr 25*$cm] [expr 5.3*$cm] 1 2 2 2 3 0 0 [expr 2*$cm2] [expr 3.14*$cm2] [expr 2*$cm2] 0 0 20 20;
buildrcrectsection 00 [expr 35*$cm] [expr 35*$cm] [expr 5.4*$cm] 1 2 3 0 3 0 3 [expr 2.54*$cm2] 0 [expr 2.54*$cm2] 0 [expr 2.54*$cm2] 20 20;
.
.
.
.
.
.
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

try and remove the dollar signs from the input arguments:
cover $concrete $steel numbartop1
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
amceuma
Posts: 23
Joined: Sat May 26, 2007 4:48 am
Location: iran

Post by amceuma »

Dear silvia,thanks for reply
but it doesnt work yet!
Last edited by amceuma on Tue Apr 29, 2008 8:50 am, edited 1 time in total.
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

you need to debug it carefully, then, please.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
amceuma
Posts: 23
Joined: Sat May 26, 2007 4:48 am
Location: iran

Post by amceuma »

Dear silvia
I change the fiber to fiberSec and then warning "subcommand `patch` is only valid inside a `section` command" and warning "error reading information in {}" appeared. I dont know what is my mistake! and what is the problem? :(
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

how about using what the manual says: Fiber, with the upper case.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
amceuma
Posts: 23
Joined: Sat May 26, 2007 4:48 am
Location: iran

Post by amceuma »

dear silvia
the previous warning appeared,again!
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

try to simplify the model and build it up.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
amceuma
Posts: 23
Joined: Sat May 26, 2007 4:48 am
Location: iran

Post by amceuma »

when i use the buildrcrectsection.tcl that exist in example files the script run and analyze correctly but when i use the section procedure written by me the program gives wrong error. can i send my files foy you to check them? I cant find the reason.
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

i am sorry, but i can't do that kind of checking.
as i said, please rebuild it carefully.
even use my file as a template.
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
amceuma
Posts: 23
Joined: Sat May 26, 2007 4:48 am
Location: iran

Post by amceuma »

i found the reason of this problem. the section id 00 is incorrect!!!
buildrcrectsection 00 [expr 35*$cm] [expr 35*$cm] [expr 5.4*$cm] 1 2 3 0 3 0 3 [expr 2.54*$cm2] 0 [expr 2.54*$cm2] 0 [expr 2.54*$cm2] 20 20
silviaTemp
Posts: 5
Joined: Thu May 01, 2008 6:21 am
Location: UC Berkeley

Post by silviaTemp »

good job.
OpenSees User-Support Manager
UC Berkeley
Post Reply