I've designed reinforced concrete beam and I intend to apply the damaged , remove the mesh. Of course I would enter this damage only part of the beam، not across the concrete section.
Thank you for your help
damege of beam
Moderators: silvia, selimgunay, Moderators
Re: damege of beam
would you please explain your question better?
sorry I can not understand what is your problem?
sorry I can not understand what is your problem?
Re: damege of beam
I really grateful for your response dear nahid2011
I have designed a reinforced concrete beam in the Opensees software. By removing a part of the mesh in beam, I would to apply the damage. I did it by Using the fiber, but removing the mesh would happen throughout the beam.
Please help me، how can I apply damage in the part of concrete beam ?
I have designed a reinforced concrete beam in the Opensees software. By removing a part of the mesh in beam, I would to apply the damage. I did it by Using the fiber, but removing the mesh would happen throughout the beam.
Please help me، how can I apply damage in the part of concrete beam ?
Re: damege of beam
if you have applied a fiber section .. then it seems that you are using a frame element
IMO you'd better divide your beam and use multi elements ... then modify the fiber section in the damage zone ... (I mean define a new modified fiber section for the element in the damage zone)
I hope I've got your question right
IMO you'd better divide your beam and use multi elements ... then modify the fiber section in the damage zone ... (I mean define a new modified fiber section for the element in the damage zone)
I hope I've got your question right
Re: damege of beam
Dear nahid2011 thank you a lot for your Guidance.
I have designed a Cantilever beam.
my program is this. but I do not know how apply damage in the beam such as the picture.
How can I put photos of the beam in here?
wipe
source beam_recorder.txt
file mkdir eigen
set txt [open eigen/eigen.out w+]
file mkdir data3
# Define model builde
model basic -ndm 2 -ndf 3
# Define nodes
node 1 0. 0.
node 2 1. 0.
node 3 1.1 0.
node 4 3. 0.
# Fix all degrees of freedom at node 1
fix 1 1 1 1
mass 2 1.e-8 2000 1.e-8
mass 3 1.e-8 2000 1.e-8
# Define geometric transformation
geomTransf Linear 1
section Elastic 1 2e11 1e10 5.625e-4
# Define Materials for nonlinear Beam
# Concrete tag fpc epsc0 fpcu epsU
# Core Concrete(confined)
uniaxialMaterial Concrete01 1 -28.e6 -0.0024 -5.6e6 -0.015
# Cover Concrete(unconfined)
uniaxialMaterial Concrete01 2 -24.e6 -0.002 -4.8e6 -0.005
# Steel tag Fy E0 b
uniaxialMaterial Steel01 4 400.e6 2.e11 0.01
# Define Cross-Section
set hbf 0.4
set hbfi [expr $hbf-2*0.05]
set bbf 0.3
set bbfi [expr $bbf-2*0.05]
section Fiber 2 {
# Create the concrete core fibers
patch quad 1 20 20 [expr $hbfi/2] [expr $bbfi/2] [expr -$hbfi/2] [expr $bbfi/2] [expr -$hbfi/2] [expr -$bbfi/2] [expr $hbfi/2] [expr -$bbfi/2]
# Create the concrete cover fibers[top,bottom,left,righ]
# Top
patch quad 2 20 6 [expr $hbfi/2] [expr $bbf/2] [expr $hbfi/2] [expr -$bbf/2] [expr $hbf/2] [expr -$bbf/2] [expr $hbf/2] [expr $bbf/2]
# Right
patch quad 2 6 20 [expr $hbfi/2] [expr $bbfi/2] [expr $hbfi/2] [expr $bbf/2] [expr -$hbfi/2] [expr $bbf/2] [expr -$hbfi/2] [expr $bbfi/2]
# Bottom
patch quad 2 6 20 [expr -$hbfi/2] [expr $bbf/2] [expr -$hbf/2] [expr $bbf/2] [expr -$hbf/2] [expr -$bbfi/2] [expr -$hbfi/2] [expr -$bbfi/2]
# Left
patch quad 2 6 20 [expr -$hbfi/2] [expr -$bbfi/2] [expr -$hbfi/2] [expr -$bbf/2] [expr $hbfi/2] [expr -$bbf/2] [expr $hbfi/2] [expr -$bbfi/2]
# Create the reinforcing fibers(top,bottom)
# Top
layer straight 4 2 [expr pow(10.e-3,2)*2*asin(1)] [expr $hbfi/2] [expr $bbfi/2] [expr $hbfi/2] [expr -$bbfi/2]
# Bottom
layer straight 4 5 [expr pow(10.e-3,2)*2*asin(1)] [expr -$hbfi/2] [expr $bbfi/2] [expr -$hbfi/2] [expr -$bbfi/2]
}
# Define element
# tag inode jnode numIntgrPts secTag tranfTag
element nonlinearBeamColumn 1 1 2 10 2 1
element nonlinearBeamColumn 2 2 3 10 1 1
element nonlinearBeamColumn 3 3 4 10 2 1
set eig [eigen 2]
puts $txt "Eigen Values are: $eig"
puts $txt "First Mode period is:[expr 2*3.1415/pow([lindex $eig 0],0.5)]"
puts $txt "Second Mode period is:[expr 2*3.1415/pow([lindex $eig 1],0.5)]"
eigen generalized 2
ModalAnalysis 2
# Define load
pattern Plain 1 Linear {
load 4 0. -100000. 0.
}
# Define recorder
recorder Node -file data3/node2disp.out -time -node 2 -dof 1 2 3 disp
recorder Node -file data3/node1reac.out -time -node 1 -dof 1 2 3 reaction
recorder Element -file ele2sec1Defo.out -time -ele 2 section 1 deformation
# Define Analysis objects
constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.e-8 6
algorithm Newton
integrator MinUnbalDispNorm 0.01
analysis Static
analyze 1
I have designed a Cantilever beam.
my program is this. but I do not know how apply damage in the beam such as the picture.
How can I put photos of the beam in here?
wipe
source beam_recorder.txt
file mkdir eigen
set txt [open eigen/eigen.out w+]
file mkdir data3
# Define model builde
model basic -ndm 2 -ndf 3
# Define nodes
node 1 0. 0.
node 2 1. 0.
node 3 1.1 0.
node 4 3. 0.
# Fix all degrees of freedom at node 1
fix 1 1 1 1
mass 2 1.e-8 2000 1.e-8
mass 3 1.e-8 2000 1.e-8
# Define geometric transformation
geomTransf Linear 1
section Elastic 1 2e11 1e10 5.625e-4
# Define Materials for nonlinear Beam
# Concrete tag fpc epsc0 fpcu epsU
# Core Concrete(confined)
uniaxialMaterial Concrete01 1 -28.e6 -0.0024 -5.6e6 -0.015
# Cover Concrete(unconfined)
uniaxialMaterial Concrete01 2 -24.e6 -0.002 -4.8e6 -0.005
# Steel tag Fy E0 b
uniaxialMaterial Steel01 4 400.e6 2.e11 0.01
# Define Cross-Section
set hbf 0.4
set hbfi [expr $hbf-2*0.05]
set bbf 0.3
set bbfi [expr $bbf-2*0.05]
section Fiber 2 {
# Create the concrete core fibers
patch quad 1 20 20 [expr $hbfi/2] [expr $bbfi/2] [expr -$hbfi/2] [expr $bbfi/2] [expr -$hbfi/2] [expr -$bbfi/2] [expr $hbfi/2] [expr -$bbfi/2]
# Create the concrete cover fibers[top,bottom,left,righ]
# Top
patch quad 2 20 6 [expr $hbfi/2] [expr $bbf/2] [expr $hbfi/2] [expr -$bbf/2] [expr $hbf/2] [expr -$bbf/2] [expr $hbf/2] [expr $bbf/2]
# Right
patch quad 2 6 20 [expr $hbfi/2] [expr $bbfi/2] [expr $hbfi/2] [expr $bbf/2] [expr -$hbfi/2] [expr $bbf/2] [expr -$hbfi/2] [expr $bbfi/2]
# Bottom
patch quad 2 6 20 [expr -$hbfi/2] [expr $bbf/2] [expr -$hbf/2] [expr $bbf/2] [expr -$hbf/2] [expr -$bbfi/2] [expr -$hbfi/2] [expr -$bbfi/2]
# Left
patch quad 2 6 20 [expr -$hbfi/2] [expr -$bbfi/2] [expr -$hbfi/2] [expr -$bbf/2] [expr $hbfi/2] [expr -$bbf/2] [expr $hbfi/2] [expr -$bbfi/2]
# Create the reinforcing fibers(top,bottom)
# Top
layer straight 4 2 [expr pow(10.e-3,2)*2*asin(1)] [expr $hbfi/2] [expr $bbfi/2] [expr $hbfi/2] [expr -$bbfi/2]
# Bottom
layer straight 4 5 [expr pow(10.e-3,2)*2*asin(1)] [expr -$hbfi/2] [expr $bbfi/2] [expr -$hbfi/2] [expr -$bbfi/2]
}
# Define element
# tag inode jnode numIntgrPts secTag tranfTag
element nonlinearBeamColumn 1 1 2 10 2 1
element nonlinearBeamColumn 2 2 3 10 1 1
element nonlinearBeamColumn 3 3 4 10 2 1
set eig [eigen 2]
puts $txt "Eigen Values are: $eig"
puts $txt "First Mode period is:[expr 2*3.1415/pow([lindex $eig 0],0.5)]"
puts $txt "Second Mode period is:[expr 2*3.1415/pow([lindex $eig 1],0.5)]"
eigen generalized 2
ModalAnalysis 2
# Define load
pattern Plain 1 Linear {
load 4 0. -100000. 0.
}
# Define recorder
recorder Node -file data3/node2disp.out -time -node 2 -dof 1 2 3 disp
recorder Node -file data3/node1reac.out -time -node 1 -dof 1 2 3 reaction
recorder Element -file ele2sec1Defo.out -time -ele 2 section 1 deformation
# Define Analysis objects
constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.e-8 6
algorithm Newton
integrator MinUnbalDispNorm 0.01
analysis Static
analyze 1