Static Reversed Cyclic Analysis

[zhangchuanchao]

Hi Silvia
I am analysing a column ，the model has been completed 。But I could not know clearly about Static Reversed Cyclic Analysis。When I study example 4 about the part of Static Reversed Cyclic Analysis ，I think it is useful 。
1 so I want to consult that LibAnalysisStaticParameters.tcl and
LibGeneratePeaks.tcl which can be used at first hand ？
2 In the Ex4.Portal2D.analyze.Static.Cycle.tcl which parameters should bemodify ？
3 set iDmax "0.005 0.01 0.035 0.1" in example 4 ，could you tell me the meaning ？
If you can help me solve these problems I would be very grateful for you。Thank you very much

[silvia]

you should follow the main file of the example step by step, as well as the procedures.

[zhangchuanchao]

Hi Silvia
Now I I encounter a trouble, spent a very long time did not solve the problem ，could you help me ?
This is my procedure and the showing of system ,I do not know how to
solve Thank you very much
Oniy the Static Reversed-Cyclic Analysis .It is a column
set IDctrlNode 2; # node where displacement is read for displacement control
set IDctrlDOF 1; # degree of freedom of displacement read for displacement contro
# characteristics of cyclic analysis
set iDmax "0.005 0.01 0.035 0.1"; # vector of displacement-cycle peaks, in terms of storey drift ratio
set Dincr [expr 0.001*$LCol]; # displacement increment for pushover. you want this to be very small, but not too small to slow down the analysis
set Fact $LCol; # scale drift ratio by storey height for displacement cycles
set CycleType Full; # you can do Full / Push / Half cycles with the proc
set Ncycles 1; # specify the number of cycles at each peak

# create load pattern for lateral pushover load
set Hload [expr $Weight/2]; # define the lateral load as a proportion of the weight so that the pseudo time equals the lateral-load coefficient when using linear load pattern
set iPushNode 2; # define nodes where lateral load is applied in static lateral analysis
pattern Plain 200 Linear {; # define load pattern -- generalized
foreach PushNode $iPushNode {
load $PushNode $Hload 0.0 0.0 0.0 0.0 0.0
}
}

# ----------- set up analysis parameters
source LibAnalysisStaticParameters.tcl; # constraintsHandler,DOFnumberer,system-ofequations,convergenceTest,solutionAlgorithm,integrator

# --------------------------------- perform Static Cyclic Displacements Analysis
source LibGeneratePeaks.tcl
set fmt1 "%s Cyclic analysis: CtrlNode %.3i, dof %.1i, Disp=%.4f %s"; # format for screen/file output of DONE/PROBLEM analysis
foreach Dmax $iDmax {
set iDstep [GeneratePeaks $Dmax $Dincr $CycleType $Fact]; # this proc is defined above
for {set i 1} {$i <= $Ncycles} {incr i 1} {
set zeroD 0
set D0 0.0
foreach Dstep $iDstep {
set D1 $Dstep
set Dincr [expr $D1 - $D0]
integrator DisplacementControl $IDctrlNode $IDctrlDOF $Dincr
analysis Static
# ----------------------------------------------first analyze command------------------------
set ok [analyze 1]
# ----------------------------------------------if convergence failure-------------------------
if {$ok != 0} {
# if analysis fails, we try some other stuff
# performance is slower inside this loop global maxNumIterStatic; # max no. of iterations performed before "failure to converge" is ret'd
if {$ok != 0} {
puts "Trying Newton with Initial Tangent .."
test NormDispIncr $Tol 2000 0
algorithm Newton -initial
set ok [analyze 1]
test $testTypeStatic $TolStatic $maxNumIterStatic 0
algorithm $algorithmTypeStatic
}
if {$ok != 0} {
puts "Trying Broyden .."
algorithm Broyden 8
set ok [analyze 1 ]
algorithm $algorithmTypeStatic
}
if {$ok != 0} {
puts "Trying NewtonWithLineSearch .."
algorithm NewtonLineSearch 0.8
set ok [analyze 1]
algorithm $algorithmTypeStatic
}
if {$ok != 0} {
set putout [format $fmt1 "PROBLEM" $IDctrlNode $IDctrlDOF [nodeDisp $IDctrlNode $IDctrlDOF] $LunitTXT]
puts $putout
return -1
}; # end if
}; # end if
# -----------------------------------------------------------------------------------------------------
set D0 $D1; # move to next step
}; # end Dstep
}; # end i
}; # end of iDmaxCycl
# -----------------------------------------------------------------------------------------------------
if {$ok != 0 } {
puts [format $fmt1 "PROBLEM" $IDctrlNode $IDctrlDOF [nodeDisp $IDctrlNode $IDctrlDOF]$LunitTXT]
} else {
puts [format $fmt1 "DONE" $IDctrlNode $IDctrlDOF [nodeDisp $IDctrlNode $IDctrlDOF]$LunitTXT]
}

---------------Uniaxial Inelastic Material, Fiber Section, Nonlinear Model ----
----------------
---------------Static Reversed-Cyclic Analysis -----------------
Model Built
WARNING: CTestEnergyIncr::test() - failed to converge
after: 6 iterations
NewtnRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor -0.020803
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
WARNING: CTestEnergyIncr::test() - failed to converge
after: 6 iterations
NewtnRaphson::solveCurrentStep() -the ConvergenceTest object failed in test()
StaticAnalysis::analyze() - the Algorithm failed at iteration: 0 with domain at
load factor 0.00235002
OpenSees > analyze failed, returned: -3 error flag
Trying Newton with Initial Tangent ..
expected floating-point number but got " -0.00000000000034106051mm"
while executing
"format $fmt1 "DONE" $IDctrlNode $IDctrlDOF [nodeDisp $IDctrlNode $IDctrlDOF]$L
unitTXT"
invoked from within
"if {$ok != 0 } {
puts [format $fmt1 "PROBLEM" $IDctrlNode $IDctrlDOF [nodeDisp $IDctrlNod
e $IDctrlDOF]$LunitTXT]
} else {
puts [format $fmt1 "DONE" ..."
(file "Ex4.Portal2D.analyze.Static.Cycle.tcl" line 80)
invoked from within
"source Ex4.Portal2D.analyze.Static.Cycle.tcl"

[silvia]

you should read the message from the system carefully:

expected floating-point number but got " -0.00000000000034106051mm"
while executing
"format $fmt1 "DONE" $IDctrlNode $IDctrlDOF [nodeDisp $IDctrlNode $IDctrlDOF]$L
unitTXT"
invoked from within
"if {$ok != 0 } {
puts [format $fmt1 "PROBLEM" $IDctrlNode $IDctrlDOF [nodeDisp $IDctrlNod
e $IDctrlDOF]$LunitTXT]

look at the input and you'll see that, somehow, you lost a space between ] and $. i checked the file you are referring to, and the space seems to be there.

[lightegg]

Hi! zhangchuanchao.I'm from china too. My QQ number is 59228664. I think we can talk about how to build SRC section in opensees.

[eda000]

hi! I am a Chinese too!

My QQ number is 376596305. Can we dicuss some questions about static reversed cyclic analysis controled by deformation?

thank you very much!