bug in time outputs

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

Moderators: silvia, selimgunay, Moderators

Post Reply
dgale
Posts: 55
Joined: Fri Jan 15, 2016 4:01 am
Location: UPM

bug in time outputs

Post by dgale » Thu May 24, 2018 8:27 am

Hello,
Could you help to understand what happen with time in my outputs?

I post the code below, just run it and see the output disp2.out at the first column.

My idea with this code is to give an input of sp displacement in a node and to obtain the forces in the spring. If you think I need a special advice, please don't be shameful.

Thanks in advance

#---------------------------------------------------------------------------------------------------------------------------------------
# MODELO MUELLE
wipe;
model basic -ndm 3 -ndf 6;
# NUDOS
node 1 0 0 0
node 2 0 0 0
# CONDICIONES DE CONTORNO EN NODOS:
fix 1 1 1 1 1 1 1
# MATERIALES:
# Hormigon
set Hormigon 1;
set Ec 30000;
set fc0 -43;
set epsc0 -0.002;
set fu0 -35.0;
set epsu0 -0.006;
set ft 2.8;

# $matTag $fc $ec $ecu $Ec $fct $et $beta
uniaxialMaterial Concrete04 $Hormigon $fc0 $epsc0 $epsu0 $Ec $ft [expr $ft/$Ec] 0.1

# Hormigon elastico
set HormigonElastico 2;
uniaxialMaterial Elastic $HormigonElastico $Ec


# ELEMENTOS:
equalDOF 1 2 2 3 4 5 6
# element zeroLength 1 1 2 -mat $HormigonElastico -dir 1
element zeroLength 1 1 2 -mat $Hormigon -dir 1

# OUTPUTS
set dtR 0.005; # Delta t del recording
set nSD 8

recorder Node -file disp2.out -precision $nSD -time -closeOnWrite -dT [expr $dtR-0.00000001] -node 2 -dof 1 disp;
recorder Node -file node1reac.out -precision $nSD -time -closeOnWrite -dT [expr $dtR-0.00000001] -node 1 -dof 1 reaction
recorder Element -file elem1force.out -precision $nSD -time -closeOnWrite -dT [expr $dtR-0.00000001] -ele 1 force;
# ANALISIS
constraints Penalty 1e16 1e16
system BandGeneral; # Tipo de matriz de rigidez a resolver
numberer Plain; #Numeración del modelo

# Displacement control--> Pushover
set maxU -0.006; # Desplaziento final
set dU -0.0001; # Incremento de desplazamiento


# timeSeries Linear 1; #Aplicación del pseudotiempo <------------------------Correct outputs
# pattern Plain 1 1 {
# sp 2 1 $dU
# }


set filePath Dfile2.txt;
timeSeries Path 1 -dt 1.0 -filePath $filePath <------------------------weird outputs in the time column
pattern Plain 1 1 -fact 1.0 {
sp 2 1 $dU
}

integrator DisplacementControl 2 1 $dU 1 $dU $dU;
test NormDispIncr 1.0e-10 50 2
algorithm KrylovNewton
analysis Static
analyze 100
#------------------------------------------------------------------------------------------------------------------------------------------------------
Dfile2.txt

0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
#-------------------------------------------------------------------------------------------------------------------------------------------------------
Results disp2.out LOOK AT 1º COLUMN, that should be time

7 -0.0001
22 -0.0002
42 -0.0003
67 -0.0004
97 -0.0005
132 -0.0006
172 -0.0007
217 -0.0008
267 -0.0009
322 -0.001
382 -0.0011
447 -0.0012
517 -0.0013
592 -0.0014
672 -0.0015
757 -0.0016
847 -0.0017
942 -0.0018
1042 -0.0019
.....

fmk
Site Admin
Posts: 5883
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley
Contact:

Re: bug in time outputs

Post by fmk » Fri Jun 08, 2018 8:11 am

what exactly are you expecting .. the time, given your model, is the load factor applied to enforce the constraint.

Post Reply