damage variables (Bilin Material,ModIMKPeakOriented Materi)

[kampix01]

Hi. I'm currently working with the Modified Ibarra-Medina-Krawinkler Deterioration Model with Bilinear Hysteretic Response (Bilin Material) and the Modified Ibarra-Medina-Krawinkler Deterioration Model with Peak-Oriented Hysteretic Response (ModIMKPeakOriented Material), respectively.
Is there any way to get a little closer look into the single damage parameter? In detail I want e.g. a timehistory of the damage parameter \beta_i, the reduced strength F_i or the reduced stiffness K_i for a better understanding. Or is it possible to record the dissipated energies of an excursion and then calculation the required variables 'manually'?

[fmk]

not w/o modifying the code and putting in the setResponse/getResponse methods. i will forward this to the author of the code.

[kampix01]

thx, would help me a lot!

[kostast88]

I have added setResponse/getResponse methods to this Material recording beta parameters and Energy.However I had trouble recording the functions local variables(beta only),so I duplicated them.No performance impact though.

[kampix01]

Would you be so kind to post a minimal example where you query these parameter. I want to to a time history analysis and check in each timestep the behaviour of one spring like:

# HYSTERESIS OF BASE SPRING
lappend basespring(M) [format %.2f [lindex [eleResponse 10177 "force"] 2]];
lappend basespring(phi) [format %.6f [eleResponse 10177 "deformation"]];
lappend basespring(dres) [eleResponse 10177 "material 1 dres"];

with the material '1':
# BaseSpring
uniaxialMaterial Bilin 1 197287690.432592 0.03 0.03 1463365.50 -1463365.5 0.5 0.5 0.5 0.5 1.0 1.0 1.0 1.0 7.410002 7.410002 2.289757 2.289757 0.1 0.1 970.717646 970.717646 1.0 1.0;

and the corresponding zerolength element:
element zeroLength 10177 1 2 -mat 1 -dir 6;

The results for the Moment and the rotation seem ok. But the result for the 'dres' is always empty, like:
basespring(M) = 103.30 488.54 1175.00 2068.03 3097.21 .....
basespring(dres) = {} {} {} {} {} .....
basespring(phi) = -0.000001 -0.000002 -0.000006 -0.000010 -0.000016 ....

.... and, instead of 'dres' I'm allowed to enter ANY keyword. E.g. if I query for 'anicetrain' the same thing happens (with version 2.4.3). BUT if I use version 2.4.1 i get an error message (Which seems intuitive).

can you please help me. Thanks a lot.
regards

[fmk]

kostast88 created his own material, a duplicate of the original with the setResponse/getResponse methods capable of returning the parameters he wanted.

http://opensees.berkeley.edu/WebSVN/fil ... iented.cpp

[kampix01]

Thanks for your answer,
If I get you right, this means that in the current OpenSees version 2.4.3. it is only possible in the ModIMKPeakOriented material to query for 'dres' but not in the 'Bilin' material? If yes, which parameter is represented by 'dres'. If not, though I'm not into C++ programming, would it be possible to add e.g. an output for beta_i parameter in the next OpenSees version?
Best regards

[fmk]

i will forward the request to the developer.