BoucWen model 's parameter of the relationship problems

[Wujianbin]

I want to simulate Steel01 material in Hysteresis behavior by the BoucWen model (No stiffness or strength degradation). The stiffness ratio is 0, and the initial stiffness and yield displacement are known. Juan David Morales once proposed the relationship between yield displacement Uy and γ and β, but I could not fit the actual hysteresis curve by using this relationship. Therefore, I want to know whether there is something wrong with my model or the problem with my parameter setting.
wipe
puts "System"
model basic -ndm 3 -ndf 6
puts "restraint"
node 1 1.000E+003 0.000E+000 0.000E+000
node 2 0.000E+000 0.000E+000 0.000E+000
puts "rigidDiaphragm"
puts "mass"
##mass 1 5.197E-002 5.197E-002 0.000E+000 0.000E+000 0.000E+000 0.000E+000
##mass 2 5.197E-002 5.197E-002 0.000E+000 0.000E+000 0.000E+000 0.000E+000
puts "node"
fix 1 0 0 0 0 0 0;
fix 2 1 1 1 1 1 1;
puts "Equal DOF"
puts "material"
uniaxialMaterial BoucWen 1 0 2735680 1 0.35 0.35 1 0 0.000001 0.0000001
##uniaxialMaterial Steel01 1 295 206000 0.001
##uniaxialMaterial Elastic 1 2.060E+005
uniaxialMaterial Elastic 2 2.482E+004
uniaxialMaterial Elastic 3 1.999E+005
uniaxialMaterial Elastic 201 2.641E+008
uniaxialMaterial Elastic 301 6.338E+008
uniaxialMaterial Elastic 401 7.547E+010
##HW400X300X10X16
section Fiber 1 -GJ 4.832E+013 {
fiber -1.200E+002 -1.920E+002 9.600E+002 1
fiber -6.000E+001 -1.920E+002 9.600E+002 1
fiber 0.000E+000 -1.920E+002 9.600E+002 1
fiber 6.000E+001 -1.920E+002 9.600E+002 1
fiber 1.200E+002 -1.920E+002 9.600E+002 1
fiber -1.200E+002 1.920E+002 9.600E+002 1
fiber -6.000E+001 1.920E+002 9.600E+002 1
fiber 0.000E+000 1.920E+002 9.600E+002 1
fiber 6.000E+001 1.920E+002 9.600E+002 1
fiber 1.200E+002 1.920E+002 9.600E+002 1
fiber 0.000E+000 -1.472E+002 7.360E+002 1
fiber 0.000E+000 -7.360E+001 7.360E+002 1
fiber 0.000E+000 0.000E+000 7.360E+002 1
fiber 0.000E+000 7.360E+001 7.360E+002 1
fiber 0.000E+000 1.472E+002 7.360E+002 1
}
section Aggregator 1001 201 Vy 301 Vz 401 T -section 1

puts "transformation"
geomTransf PDelta 1 0.000 0.000 1.000
puts "element"
element nonlinearBeamColumn 1 2 1 3 1001 1
puts "shell element"
puts "SOLID element"
puts "recorder"
recorder Node -file node0.out -time -nodeRange 1 2 -dof 1 2 3 disp
recorder Node -file node2.out -time -node 1 -dof 1 2 3 disp
recorder Node -file reaction2.out -time -node 2 -dof 1 2 3 reaction
puts "loading"
## Load Case = PUSH
pattern Plain 3 Linear {
load 1 1.000E+000 0.000E+000 0.000E+000 0.000E+000 0.000E+000 0.000E+000
}
puts "analysis"
constraints Plain
numberer Plain
system BandGeneral
test NormDispIncr 1.0e-1 2000
algorithm Newton
analysis Static
integrator DisplacementControl 1 1 0.1
analyze 100
integrator DisplacementControl 1 1 -0.2
analyze 100
integrator DisplacementControl 1 1 0.3
analyze 100
integrator DisplacementControl 1 1 -0.4
analyze 100
integrator DisplacementControl 1 1 0.5
analyze 100
integrator DisplacementControl 1 1 -0.6
analyze 100
integrator DisplacementControl 1 1 0.7
analyze 100
integrator DisplacementControl 1 1 -0.8
analyze 100

[whoppingtiara]

You have also the problem ?
geometry dash subzero

[Wujianbin]

Thanks for your reply, I have already solve my questiopns, this method can achieve γ and β