Sinusoidal displacement loading protocol

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

Moderators: silvia, selimgunay, Moderators

Post Reply
giovannimilan
Posts: 71
Joined: Wed Apr 05, 2017 12:49 am
Location: ARUP, amsterdam

Sinusoidal displacement loading protocol

Post by giovannimilan » Sun Sep 10, 2017 11:59 pm

Hi all,

is there an easy-way to prescribe a displacement history at specific nodes of a structure, referring to a .tcl file containing the protocol?

I'd like use a sinusoidal loading protocol used during the lab test.

Cheers,

Jhno
Posts: 214
Joined: Sat May 05, 2012 2:55 pm
Location: Laval

Re: Sinusoidal displacement loading protocol

Post by Jhno » Mon Sep 11, 2017 4:21 am


giovannimilan
Posts: 71
Joined: Wed Apr 05, 2017 12:49 am
Location: ARUP, amsterdam

Re: Sinusoidal displacement loading protocol

Post by giovannimilan » Thu Sep 14, 2017 4:59 am

Jhno wrote:
> sp
> http://opensees.berkeley.edu/wiki/index.php/Sp_Command

Thanks for the reply.

Does it require a non-Plain constrain handler?

I ve got a model with zero length element and EqualDOF and i get a warning that I Cannot run the analysis.

the thing is that i switched to lagrange handler but the deform shape does not make sense. it is messing with the others constrains I have in the model ( i reckon EqualDOF in aprticular).

see code below:
if {$analysis == "protocol"} {
puts "Running the sinusoidal loading protocol...\n"

timeSeries Path 1 -fileTime "protocol_time.tcl" -filePath "protocol_value.tcl";

set controlnode 5; # top node belonging to a rotational zero length spring

pattern Plain 1 1 {
sp $controlnode 1 1.0
}

set numberOfsteps 1343; # corresponding to 300 sec

set incrOutput 0.01

# Cyclic loading to the model by Disp. control
constraints Lagrange
numberer Plain
system BandGeneral
set TolStatic 1.e-7; # Convergence Test: tolerance
set maxNumIterStatic 500; # Convergence Test: maximum number of iterations
set printFlagStatic 0; # Convergence Test: flag used to print information on convergence (optional) # 1: print information on each step;
test EnergyIncr $TolStatic $maxNumIterStatic $printFlagStatic;
algorithm KrylovNewton;
integrator LoadControl $incrOutput
analysis Static
analyze $numberOfsteps; #steps

set currentTime [getTime];

puts "current time is $currentTime in sec\n"

puts "the loading protocol has reached the end...\n"
}

Jhno
Posts: 214
Joined: Sat May 05, 2012 2:55 pm
Location: Laval

Re: Sinusoidal displacement loading protocol

Post by Jhno » Thu Sep 14, 2017 6:43 am

I am sorry I wouldn't know exactly the solution to your problems.

If you have alot of constraint in the model, I don't think the plain handler will work, especially adding more SP_Constraint with the SP command.
Maybe you can try using the penalty method (large values, say 1.e18) and see what the deform. shapes look likes? Gl.

giovannimilan
Posts: 71
Joined: Wed Apr 05, 2017 12:49 am
Location: ARUP, amsterdam

Re: Sinusoidal displacement loading protocol

Post by giovannimilan » Thu Sep 14, 2017 7:04 am

Jhno wrote:
> I am sorry I wouldn't know exactly the solution to your problems.
>
> If you have alot of constraint in the model, I don't think the plain
> handler will work, especially adding more SP_Constraint with the SP
> command.
> Maybe you can try using the penalty method (large values, say 1.e18) and
> see what the deform. shapes look likes? Gl.

I made it run with constraints Transformation

Cheers,

Post Reply