Search found 8 matches

by luiscelorrio
Mon Dec 01, 2014 5:58 am
Forum: OpenSees.exe Users
Topic: Trying nonlinear methods with integrator LoadControl
Replies: 1
Views: 2245

Trying nonlinear methods with integrator LoadControl

How I can implement a code like this for the case of the integrator command is
Integrator LoadControl ?. That is, If I use integrator LoadControl, is possible to change the algorithm through the integration process?
I have some convergence problems and I would like to try several nonlinear methods.
Thanks on advance
Luis

set algorithmType Newton
algorithm $algorithmType;
integrator DisplacementControl $IDctrlNode $IDctrlDOF $Dincr
analysis Static

# --------------------------------- perform Static Pushover Analysis
set Nsteps [expr int($Dmax/$Dincr)]; # number of pushover analysis steps
set ok [analyze $Nsteps]; # this will return zero if no convergence problems were encountered

# ---------------------------------- in case of convergence problems
if {$ok != 0} {
# change some analysis parameters to achieve convergence
# performance is slower inside this loop
set ok 0;
set controlDisp 0.0; # start from zero
set D0 0.0; # start from zero
set Dstep [expr ($controlDisp-$D0)/($Dmax-$D0)]
while {$Dstep < 1.0 && $ok == 0} {
set controlDisp [nodeDisp $IDctrlNode $IDctrlDOF ]
set Dstep [expr ($controlDisp-$D0)/($Dmax-$D0)]
set ok [analyze 1 ]
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 2000 0
algorithm Newton -initial
set ok [analyze 1 ]
test $TestType $Tol $maxNumIter 0
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 ]
algorithm $algorithmType
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch .8
set ok [analyze 1 ]
algorithm $algorithmType
}
}
}; # end if ok !0

puts "DonePushover”
by luiscelorrio
Wed Nov 12, 2014 6:00 am
Forum: Reliability Analysis
Topic: Assign Random Variable to Reinforced Steel
Replies: 1
Views: 4399

Re: Assign Random Variable to Reinforced Steel

Hello, all
I think that Steel02 material is not valid, as ii is implemented for structural reliability analysis.
Using uniaxialMaterial Hardening the error no objects were able to identify parameter
disappears, however convergence problems appear.

Thanks
Luis
by luiscelorrio
Sun Nov 02, 2014 4:18 am
Forum: Reliability Analysis
Topic: Assign Random Variable to Reinforced Steel
Replies: 1
Views: 4399

Assign Random Variable to Reinforced Steel

Hello, I have a question about assigning random variables to material properties of layer in a fiber section.
How I can consider the material properties of the reinforced Steel bars in a RC section as random variables? I am working to compute structural reliability of a RC frame. I have defined RC section as a fiber section.
Reinforced steel are defined with:
# bottom layer reinforcement
layer straight $IDSteel $numBarsColExt $barAreaColExt -$coreY $coreZ -$coreY -$coreZ;
# top layer reinforcement
layer straight $IDSteel $numBarsColExt $barAreaColExt $coreY $coreZ $coreY -$coreZ;

I have considered the material Steel02 for the reinforced bars.
When I assign random variables to material properties of reinforced steel I use:

addToParameter $tagECE element $e section $ColExtSecTag material $IDSteel E

However it doesn’t work and I obtain this error:
Parameters addComponent 121 -- no objects were able to identify parameter section 2 material 3 E

How can solve it?
Thanks in advance
Luis
by luiscelorrio
Mon May 19, 2014 4:20 am
Forum: Reliability Analysis
Topic: computeGradients
Replies: 2
Views: 8889

computeGradients

Hello,
In my code I have used the command computeGradients
This is a pierce of my code:

computeGradients
foreach p [lsort -integer [getParamTags]] {
puts $gradDESfile "[sensNodeDisp 2 2 $p]"
puts "[sensNodeDisp 2 2 $p] "
}

My question is: Are the gradients obtained, normalized gradients or real gradients?
I need these gradientes for optimization.
Thanks on advance
by luiscelorrio
Thu Mar 13, 2014 12:07 am
Forum: Reliability Analysis
Topic: Gradients of Reliability Index
Replies: 1
Views: 3918

Gradients of Reliability Index

Hello,
I would like to know:
Is there a function or command to compute the gradients of the reliability index (betaFORM) with relation to means of random variables?
I have also some doubts about some functions starting with “get..” like

getNodeTags
getLSFTags
getParamTags
getRVTags

Could you explain them or/and list all of them?
Thanks on advance

Luis
by luiscelorrio
Sat Jan 25, 2014 2:07 am
Forum: Reliability Analysis
Topic: Reliability
Replies: 11
Views: 42070

Re: Reliability

Hello, I have read the paper: "Shape sensitivities in the reliability analysis of nonlinear frame structures". In this paper two examples are explained. I thanks that Michael had uploaded the first example involving a steel structure. I have run it and it works very well. Can you please upload the second example of a reinforced concrete structure?

I have some questions about section command.
Is possible to use the command section WFSection2d to compute DDM based sensitivity?

I have found that the information about the command Section Aggregator in the user manual and in the Sensitivity Command Manual is different. What is the correct command?

Thanks in advance.
Luis
by luiscelorrio
Fri Jan 03, 2014 4:16 pm
Forum: OpenSees.exe Users
Topic: Double T section in Steel
Replies: 0
Views: 1593

Double T section in Steel

Hi, I would like to computer and simulate a steel frame that consist in doule T sections. I will consider the width and height of column secctions and beam section as parameters or design variables in an optimization problem. I am modelling these beams and columns as nonlinearBeamColumn elements and sections using a fiber type section with 3 patch comands (upper flange, botton flange and web). Is this correct?.
How could i compute the gradients with respect to these parameter (b = width and h = height)?
Thanks in advance.
by luiscelorrio
Wed Aug 07, 2013 7:08 am
Forum: OpenSees.exe Users
Topic: DAKOTA - OPENSEES interface
Replies: 0
Views: 1576

DAKOTA - OPENSEES interface

Hi, Opensees Users.
My name is Luis Celorrio, I came from La Rioja (Spain).
I would like to know if someone has been working in Optimization Uncer Uncertainty using DAKOTA fron Sandia National Laboratory as an Optimization package and OpenSees as a Simulation Code. It is posible to connect these two packages?. What type of interface could be used?
Thanks on advance