problem while run parameter study

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

Moderators: silvia, selimgunay, Moderators

Post Reply
alaa
Posts: 21
Joined: Mon Jul 18, 2005 3:29 am
Location: egypt , cairo

problem while run parameter study

Post by alaa »

hi all
i have amodel like example 5.1 3d rigid frame .
but with bidirection eccentricity ex and ey
iam intersiting in master node displacement
i can run the model and get resonable result for each ex an ey
but when i try to run a parameter study while varaying ex and ey there is a problem.
well
i use this script to run a parameter study

test.tcl
set node132 0 ;# node 13 disp in direction 2 will count from(1 to 2) no.
set node1 3 ;# node1 reaction from (4 to 5)
set node2 6 ;# node 2 reaction from (7 to 8)
set node3 9 ;# node 3 reaction from ( 10 to 11 )
set node4 12 ;# node 4 reaction from(13 to 14 )
foreach ex { 0.02 0.1 ......... } ey { 0.02 0.1 ....... } {
wipe
set node132 [ expr $node132+1]
set node1 [ expr $node1+1]
set node2 [ expr $node2+1]
set node3 [ expr $node3+1]
set node4 [ expr $node4+1]
source 3dframe.tcl
}


then i run the model
for first values of ex and ey ( 0.02 0.02)
i got the same result if i run a normal model include only one value
for ex and ey ( 0.02 0.02 )
but stasting from second value for ex and ey
there is a different in results
start softly with the earthquake result then increse sharbly at the
end of earthquake
plz any one can help me

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

Post by silvia »

i am confused, please mail me your scripts.
mazzoni@berkeley_NOSPAM_.edu (remove the _NOSPAM_)
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
silvia
Posts: 3909
Joined: Tue Jan 11, 2005 7:44 am
Location: Degenkolb Engineers
Contact:

Post by silvia »

ALAA,
I have looked over the input file.
The problem with your input is that the ground motion is too large and the element is unable to recover when it has been pushed beyond it’s crushing strength.
You need to scale the ground motion down to have reasonable results. I put this into your 3dfame file.

Also,
I cleaned up your parameter-study input file. I think it makes the filenames a little easier to understand and it allows you to run as many analyses as you like.
This required a modification to the 3dframe, too, in the recorder area.

So I am attaching the two files, and will post this to the board for records.

Let me know if you need further help.

here is what the main file looks like:

Code: Select all

set iex {0.01 0.02 0.03 0.04 0.05 0.1} 
set iey {0.01 0.02 0.03 0.04 0.05 0.1}

set Nanal [llength $iex]

for {set ianal 1} {$ianal <= $Nanal} {incr ianal 1} {
	wipe
	puts Case$ianal
	set counter [expr $ianal-1]
	set ex [lindex $iex $counter]
	set ey [lindex $iey $counter]
	set endword Case$ianal
	source 3dframe.tcl
}
and this is the modification in the frame file:

Code: Select all

# Set up the acceleration records for Tabas fault normal and fault parallel
set GMfact [expr 0.5*$g];		
set tabasFP "Path -filePath tabasFP.txt -dt 0.02 -factor $GMfact";
and:

Code: Select all

recorder Node -file Rnode1$endword.xls -time -node 1 -dof 1 2 3 reaction;	
recorder Node -file Rnode2$endword.xls -time -node 2 -dof 1 2 3 reaction;	
recorder Node -file Rnode3$endword.xls -time -node 3 -dof 1 2 3 reaction;	
recorder Node -file Rnode4$endword.xls -time -node 4 -dof 1 2 3 reaction; 	
recorder Node -file DYnode13$endword.xls -time -node 13 -dof 2 disp	;
Silvia Mazzoni, PhD
Structural Consultant
Degenkolb Engineers
235 Montgomery Street, Suite 500
San Francisco, CA. 94104
alaa
Posts: 21
Joined: Mon Jul 18, 2005 3:29 am
Location: egypt , cairo

Post by alaa »

wow
it is already work
thanx so much
silvia
Post Reply