Hello everyone,
I am attempting to use a 2D cantilever beam that spans 40 units that stops displacing downward after reaching 1 unit downward. I want to use the elastic gapPP material and assign it to a zero length element. I have node 2 and 3 at 40.0 and 0.0 and fixed node 3 but the gap element does not work. What can I do to make the gap material work? Should I use it with a different element type? My code is pasted below. I would really appreciate any help with this matter.
-------------------------------------------------------------------------------------
wipe
model basic -ndm 2 -ndf 3
set dataDir TEST;
file mkdir $dataDir;
# DEFINE NODES
# node # X Y
node 1 0.0 0.0
node 2 40.0 0.0
node 3 40.0 0.0
puts "Nodes have been defined"
# DEFINE BOUNDARY CONDITIONS
fix 1 1 1 1
fix 3 1 1 1
puts "Boundary conditions have been defined"
# DEFINE MATERIAL PROPERTIES
##Beam Concrete (cover) Tag f'c ec0 fcu ecu
uniaxialMaterial Concrete01 1 4.8 0.002 0.205 0.0045
##Beam Concrete (confined core) Tag f'c ec0 fcu ecu
uniaxialMaterial Concrete01 2 5.3 0.00220 1.06 0.0085
uniaxialMaterial Elastic 999 1.0E-10; # Filler material
set E 28758.62;
#uniaxialMaterial Steel01 $matTag $Fy $E0 $b
uniaxialMaterial Steel01 3 83.4 $E 0.006
uniaxialMaterial ElasticPPGap 4 $E -50000000 -1.0
puts "Material properties have been defined"
# DEFINE SECTION PROPERTIES
set b [expr 3.94/2.0]; # b=half section width
set h1 [expr 6.302/2.0]; # h1=half core concrete height
set h1l 3.743;
set b1 [expr 3.152/2.0]; # b1=half core concrete width
set h [expr 7.09/2.0]; # h=half section height
set hl [expr 8.27/2.0];
set h2 2.7176; # To consider actual location of longitudinal bars
section Fiber 1 {
# MTag NIJ NJK $yI $zI $yJ $zJ $yK $zK $yL $zL
patch quad 1 20 5 -$b $h1 $b $h1 $b $h -$b $h; # concrete top cover
patch quad 2 20 40 -$b -$h $b -$h $b $h -$b $h; # concrete core
patch quad 1 20 5 -$b -$h $b -$h $b -$h1 -$b -$h1; # concrete bottom cover
patch quad 1 5 40 -$b -$h1 -$b1 -$h1 -$b1 $h1 -$b $h1; # concrete left cover
patch quad 1 5 40 $b1 -$h1 $b -$h1 $b $h1 $b1 $h1; # concrete right cover
# layer straight $matTag $numFiber $areaFiber $yStart $zStart $yEnd $zEnd
layer straight 3 2 0.122 -$b $h2 $b $h2; # top reinforcement
layer straight 3 2 0.122 -$b -$h2 $b -$h2; # bottom reinforcement
};
puts "Section has been defined"
# DEFINE ELEMENTS
geomTransf Linear 1
puts "Geometric transformation has been defined"
element dispBeamColumn 1 1 2 5 1 1
element zeroLength 2 2 3 -mat 4 -dir 2
puts "Elements have been defined"
# Create recorders
recorder Node -file $dataDir/Node2Tests.out -time -node 2 -dof 1 2 3 disp
recorder Element -file $dataDir/Ele1Force.out -time -ele 1 force
puts "Output files have been defined"
#----------------------------------------------------------------------------
# DEFINE ANALYSIS PARAMETERS
#----------------------------------------------------------------------------
pattern Plain 1 "Linear" {
load 2 0.0 1.0 0.0;
}
# Define analysis parameters
constraints Plain;
numberer Plain;
system BandGeneral;
test RelativeNormDispIncr 1.0 200000;
algorithm Newton;
integrator LoadControl -0.1;
analysis Static;
puts "Analysis parameters have been defined"
puts "Analysis starts. Please wait for results"
analyze 100;
puts "Analysis has been performed"
Elastic GapPP material to zero length element
Moderators: silvia, selimgunay, Moderators