Bouc-wen Model

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

Moderators: silvia, selimgunay, Moderators

Post Reply
fanyijin
Posts: 1
Joined: Mon Sep 12, 2022 5:10 pm

Bouc-wen Model

Post by fanyijin » Mon Sep 12, 2022 5:19 pm

Hi
I first simulate reinforced concrete with concrete 01 and steel 01 constitutive models, FC = 30MPa, FY = 335mpa, esteel = 2.0e11pa, and the cross-section is a circular pier with a diameter of 1m.
If I want to use bouc Wen model instead of reinforcement and concrete constitutive model to establish the same pier column above.
I modeled a single degree of freedom column with the Bouc-Wen model ,wheteher to define fiber section and to assigin Bouc-Wen material to each fiber of the section?
Is the code I wrote below correct?If it is correct, is it regular to assign values to these parameters when simulating RC columns? If the code is wrong, how should I adjust it?
In addition, is it the same method to build a two-dimensional model and a three-dimensional model?
Thank you very much for your reply!

#units:KN,m
wipe
model basic -ndm 2 -ndf 3
if { [file exists output] ==0 } {
file mkdir output;
}
node 1 0 0
node 2 0 0.5
node 3 0 1.0
node 4 0 1.5
node 5 0 2.0
node 6 0 2.5
node 7 0 3.0
node 8 0 3.5
node 9 0 4.0

fix 1 1 1 1

geomTransf Linear 1
uniaxialMaterial BoucWen 1 0.06 250000 1 -400 520 1 0.2 0.8 0.4

section Fiber 1 {
patch circ 1 40 10 0 0 0 0.5 0 360
}

#柱
element dispBeamColumn 1 1 2 3 1 1
element dispBeamColumn 2 2 3 3 1 1
element dispBeamColumn 3 3 4 3 1 1
element dispBeamColumn 4 4 5 3 1 1
element dispBeamColumn 5 5 6 3 1 1
element dispBeamColumn 6 6 7 3 1 1
element dispBeamColumn 7 7 8 3 1 1
element dispBeamColumn 8 8 9 3 1 1


recorder Node -file output/disp_9.out -time -node 9 -dof 1 disp

pattern Plain 2 Linear {
load 9 1.0 0.0 0.0
}
constraints Plain
numberer RCM
system BandGeneral
test EnergyIncr 1.0e-6 200
algorithm Newton
analysis Static
foreach disp {50 100 200 400 800 1000 } {
integrator DisplacementControl 9 1 0.0001;
analyze [expr $disp]
integrator DisplacementControl 9 1 -0.0001;
analyze [expr $disp]
integrator DisplacementControl 9 1 -0.0001;
analyze [expr $disp]
integrator DisplacementControl 9 1 0.0001;
analyze [expr $disp]
integrator DisplacementControl 9 1 0.0001;
analyze [expr $disp]
integrator DisplacementControl 9 1 -0.0001;
analyze [expr $disp]
integrator DisplacementControl 9 1 -0.0001;
analyze [expr $disp]
integrator DisplacementControl 9 1 0.0001;
analyze [expr $disp]

}

Post Reply