Regards , please someone help
the recorder element shell does not work ? files are empty always.
I'm wrong
PD: displacement is accurate
version = Opensees 2.2.1
wipe
# unidades kg , mts , sec.
# coordenadas y grados de libertad
model basic -ndm 3 -ndf 6
# parametros geometricos, ancho mts y largo mts
set x 9
set y 3
set largo $x.
set alto $y.
# segmentos a cada 0.1 mts o sea que habran 10 segmentos en un metro
set segmentosX [expr {10*$x}]
set segmentosY [expr {10*$y}]
set nx [expr $segmentosX]
set ny [expr $segmentosY]
set nodostotal [expr (($nx)+1)*(($ny)+1)];
set area_nodo [expr (($alto/$ny)*($largo/$nx))]
# generacion de nodos y elementos
# argumento para bloque2d , en este caso el tipo de material 1.
set eleArgs "1"
# Definicion de nodos del bloque2D
set blockData "
1 0. 0. 0.
2 $x 0. 0.
3 $x $y 0.
4 0. $y 0.
"
# material
# $tag_de_seccion $Modulo_de_elasticidad $Relacion_de_poisson $Espesor_de_la_seccion $densidad_relativa
# $secTag $E $nu $h $rho
section ElasticMembranePlateSection 1 2.535e9 0.2 0.25 2.4
# Definicion de bloque 2D para generar mesh
# block2d $nx $ny $e1 $n1 elemento (argumentos del elemento)
# $nx $ny numeros de elementos en x y y local respectivamente
# $e1 $n1 elemento de inicio y nodo inicial para la generacion respectivamente
block2D $nx $ny 1 1 ShellMITC4 $eleArgs $blockData
#condiciones de frontera empotrado en los lados y al fondo del muro
fixX 0.0 1 1 1 1 1 1
fixX $x. 1 1 1 1 1 1
fixY 0.0 1 1 1 1 1 1
# CARGAS
# Patron de cargas para presion de agua
pattern Plain 1 Linear {
for {set i 1} {$i <= [expr $ny+1]} {incr i} {
for {set j 1} {$j <= [expr $nx+1]} {incr j} {
set nodo [expr ($j+($i-1)*($nx+1))]
set presion [expr (($ny+1-$i)*($alto/$ny)*($area_nodo)*1000)]
load $nodo 0.0 0.0 $presion 0.0 0.0 0.0
}
}
}
# Archivos de salida
file mkdir data
set FileName0 "data/"
set FileName1 NodeDisp
set FileNameNODOS nodostotal
set FileNameEle elementos_fuerza
recorder Node -file $FileName0$FileName1.out -nodeRange 1 $nodostotal -dof 1 2 3 4 5 6 disp;
recorder Element -file $FileName0$FileNameEle.out -eleRange 1 2700 section 1 forces;
recorder Element -file data/stresses_elemento.out -ele all section 1 stresses;
recorder Element -file dato/ele_force.out -eleRange 1 2700 section 1 force;
recorder Element -file data/stresses_elemento.out -ele all section 1 stresse;
#analysis
integrator LoadControl 1.0 1 1.0 1 10
test EnergyIncr 1.0e-15 20 1
algorithm Newton
numberer RCM
constraints Plain
system BandGeneral
analysis Static
analyze 1
#resultados
# for {set ki 46} {$ki <= [expr $nodostotal]} {incr ki 91} {
# print node $ki
# }
print $FileName0$FileNameNODOS.ops
exit
recorder shell PROBLEM ?
Moderators: silvia, selimgunay, Moderators
-
- Posts: 2
- Joined: Thu May 28, 2009 9:43 am
- Location: nicaragua
- Contact:
recorder shell PROBLEM ?
________________
ing. Marcel ToruNo
ing. Marcel ToruNo