recorder Element command problem

[camphilli]

I am using UCSD soil model and elements to model soil liquefaction. I need to verify the stress-strain relationship of the material, but, the recorder Element command is not working, therefore, I am not able to obtain as an output neither the stress nor the strain time histories. I already have the acceleration, velocity and displacement time histories but as I mentioned before I can not construct the backbone curve (this command is not working also) or the hysteresis loops for a given element.

I am using the following command line (the input file is too long to post it here):

recorder Element 1 -time -file stress1 -dT $dT material 1 stress

Thanks much for the help

Camilo Phillips

[fmk]

just post the element and material command you are using, so i can see ele type and material type.

[amin_asareh]

you have to run this recorder in opensees version 1.7.1. it will not give an answer in new versions. I don't know why!

[fmk]

if someone would post the element & material ommand i can either fix it or tell you what's wrong with the command.

[zhenxinli1999]

recorder EnvelopeNode -file $dataDir/Nenvy.out -time -node 6 -dof 1 2 3 disp
recorder EnvelopeElement -file $dataDir/FEnvy175.out -time -ele 175 localForce
Is there something wrong?
When I use the Example07 I can get the output for Node recorder and Element recorder.

[silvia]

frank has asked that you post the command where you define the elements and materials, please.
also, please remove the -time in the recorder command -- follow the manual format.

[amin_asareh]

the following recorders have some problems when opensees 1.7.3 is used. the output for these recorders does not have any information and it is 0Kbs. these recorders only work in 1.7.1 version.

recorder Element $ele -time -file $name11 material 1 stress
recorder Element $ele -time -file $name22 material 1 strain

[rjaeger]

Which element are you using?

[amin_asareh]

We are using Pressure dependant multi yield materials. I got the stress strain recorders in version 1.7.1 with no problem. but they do not work with the new versions.

[fmk]

POST a 1 element 1 material example. there
are too many options for me to try and figure out what you are doing. Again, PLEASE POST a simple 1 element 1 node example.

[bayram_aygun]

The code below has only one soil element and one P-y spring. It looks long but it's fairly simple. My integrator doesn't work properly but I post this code to help Dr. McKenna figure out why this "recorder" command doesn't record!

Thanks,



#One Pile-Soil Interaction Model By Bayram Aygun, NOV2007
#Units (m,kN,ton)

wipe

set NumSteps 2000
set StepSize [expr 1.0/$NumSteps]
#
##Build Soil Model

model basic -ndm 2 -ndf 2


############################################

#define loose sand stratum parameters

############################################

set rho 1.7000;
set refShearModul 55000.;
set refBulkModul 150000.;
set frictionAng 29.;
set peakShearStra 0.1;
set refPress 80. ;
set pressDependCoe 0.5;
set PTAng 29.;
set contrac 0.21;
set dilat1 0.;
set dilat2 0.;
set liquefac1 10.;
set liquefac2 0.02;
set liquefac3 1.0;
set e 0.85;

# $tag $nd

nDMaterial PressureDependMultiYield 3 2 $rho $refShearModul $refBulkModul $frictionAng $peakShearStra $refPress $pressDependCoe $PTAng 0.17 $dilat1 $dilat2 $liquefac1 $liquefac2 $liquefac3 $e

updateMaterialStage -material 3 -stage 0 #PDMY

################################################

# Define Fluid-Solid Coupling for Loose Sand

################################################

set combinedBulkModul_1 2.2D6; #Combined undrained bulk modulus for weak soil

#nDMaterial FluidSolidPorousMaterial $tag $nd $soilMatTag $combinedBulkModul

nDMaterial FluidSolidPorous 5 2 3 $combinedBulkModul_1

updateMaterialStage -material 5 -stage 0 #FSP

#############################################

# Define P-y Spring Parameters

#############################################

# Define PyLiq1

set SoilType 2;
set P_ult_1 65.;
set Y_50_1 0.085;
set C_d_1 0.07;
set C_1 0.0;
set pRes_1 6.5;

# uniaxialMaterial PyLiq1 $matTag $soilType $pult $Y50 $Cd $c $pRes $solidElem1 $solidElem2

uniaxialMaterial PyLiq1 7 $SoilType $P_ult_1 $Y_50_1 $C_d_1 $C_1 $pRes_1 1 1


updateMaterialStage -material 7 -stage 0

#Define Soil Nodes
# x y

node 1 0. 0.
node 2 10. 0.
node 3 10. 10.
node 4 0. 10.


################################################

#define the actual soil element

################################################

set inclination 0.
set massDen 2.0 ;# solid mass density
set fluidDen 1.0 ;# fluid mass density
set pi 3.1415926535

set unitWeightX [expr ($massDen-$fluidDen)*9.81*sin($inclination/180.0*$pi)] ;# unit weight in X direction
set unitWeightY [expr -($massDen-$fluidDen)*9.81*cos($inclination/180.0*$pi)] ;# unit weight in Y direction


# This command is used to construct a FourNodeQuad element object which uses a bilinear isoparametric formulation.

#element quad $eleTag $iNode $jNode $kNode $lNode $thick $type $matTag <$pressure $rho $b1 $b2>



element quad 1 1 2 3 4 1.0 "PlaneStrain" 5 0.0 0 $unitWeightX [expr $unitWeightY*1.]



# fix the soil base

fix 1 1 1
fix 2 1 1

#tie nodes 3 and 4

equalDOF 3 4 1 2

##dummy node
node 33 10 10



# element zeroLength $eleTag $iNode $jNode -mat $matTag1 $matTag2 ... -dir $dir1 $dir2

element zeroLength 47 33 3 -mat 7 -dir 1

#fix the dummy node in vertical direction

fix 33 0 1


################################
# Gravity Application
###############################

system ProfileSPD
test NormUnbalance 1.0e-5 10 0
constraints Transformation
integrator LoadControl 1 1 1 1

algorithm Newton
numberer RCM

analysis Static

analyze 5

recorder Node -file SoilDisp.dat -time -node 3 -dof 1 2 disp
print -node 3 4


updateMaterialStage -material 5 -stage 1 #FSPM
updateMaterialStage -material 3 -stage 1 #PDMY

updateMaterials -material 3 bulkModulus 100000
updateMaterialStage -material 7 -stage 1 #PyLiq1

model basic -ndm 2 -ndf 3

node 5 10 0
node 6 10 10
node 7 10 20

#fix the pile tip in every possible DOF
fix 5 1 1 1

# fixing the pile nodes in the vertical direction

fix 6 0 1 0
fix 7 0 1 0
#
#create geometric transformation for pile elements

geomTransf Linear 1

# define materials for pile (for constitutive nonlinear beam-column elements)

#uniaxialMaterial Concrete02 $matTag $fpc $epsc0 $fpcu $epsU $lambda $ft $Ets
uniaxialMaterial Concrete02 1 -27580. -0.003 -5516. -0.01 0.1 3681. 1930544
#build cover concrete (unconfined)


#uniaxialMaterial Steel02 $matTag $Fy $E $b $R0 $cR1 $cR2 $a1 $a2 $a3 \
#$a4
#use Steel01 next time.!!!!!!!!!!!!!
uniaxialMaterial Steel02 2 460572.64 199949200. 0.01 18 0.925 0.15
# build rein. matr


#section Fiber $secTag {
section Fiber 109 {

# patch quad $matTag $numSubdivIJ $numSubdivJK $yI $zI $yJ $zJ $yK $zK $yL $zL
patch quad 1 4 16 -0.30 0.3 -0.3 -0.3 0.3 -0.3 0.3 -0.3

# layer straight $matTag $numBars $areaBar $yStart $zStart $yEnd $zEnd
layer straight 2 16 0.00145 0.25 0.25 0.25 -0.25; # top layer rein
layer straight 2 16 0.00145 -0.25 0.25 -0.25 -0.25; # bottom layer rein

}; # end of fiber section definition



#element nonlinearBeamColumn $eleTag $iNode $jNode $numIntgrPts $secTag $transfTag <-mass $massDens>\
element nonlinearBeamColumn 17 5 6 5 109 1 #<-iter $maxIters $tol>
element nonlinearBeamColumn 39 6 7 5 109 1

#Slave the pile node in 3DOF to the dummy node in 2DOF in the x & y direction
#equalDOF $rNodeTag $cNodeTag $dof1 $dof2

equalDOF 33 6 1 2


mass 7 1. 0. 0. 0. 0. 0.

###########################
#Now DEFINE TRANSIENT lOADS
#############################

set ground_motion "Path -filePath EQ_Rix41.acn -dt 0.005 -factor [expr 1.00*384.6]"
pattern UniformExcitation 2 1 -accel $ground_motion



#rezero time
setTime 0.0
wipe analysis #This command does not destroy the elements, nodes, materials, etc. It does...
#destroy the solution strategies: the algorithm, analysis, eq'n solver, constraint handler, etc.


constraints Transformation
test NormDispIncr 1.e-2 1000 0
algorithm Newton
numberer RCM
system ProfileSPD

set NewmarkGamma 0.5; # Newmark-integrator gamma parameter (also HHT)
set NewmarkBeta 0.3025; # Newmark-integrator beta parameter
integrator Newmark $NewmarkGamma $NewmarkBeta
analysis Transient
analyze $NumSteps 0.0005
#0.010 0.00010 0.010 10

recorder Node -file PileTopDisp.out -node 7 -dof 1 2 -dT 0.01 disp
print -node 3 4 7 6