uniaxialMaterial ModIMKPinching

[mariateresa]

Hello,
I need to use the Modified Ibarra-Medina-Krawinkler Deterioration Model with Pinched Hysteretic Response
(uniaxialMaterial ModIMKPInching) in order to include pinching effect in my model.
I've downloaded the latest OpenSees version a few hours ago, but the uniaxialMaterial ModIMKPInching doesn't work.
A warning message appears: "WARNING could not create uniaxialMaterial ModIMKPinching".
What could be the problem?
This is the script that I used:


model BasicBuilder -ndm 3 -ndf 6

# define joint nodes in (0,0)
node 1 0 0 0
node 2 0 0 0

# define restrains
fix 1 1 1 1 1 1 1
fix 2 1 1 1 1 0 1

################
set matTagBilinPinch 1

################
set theta_cr 4.7360e-004;
set theta_y 0.0092
set theta_u 0.0910
set M_cr 19.7877;
set M_y 98.4101;
set M_u 98.4111;
################
set ni 0.0556
set rosh 0.0033
set s_d 0.3704
set fc 25
set a_sl 1;

################ set parameters for Bilin "Modified Ibarra-Krawinkler Model with Bilinear Hysteretic Response"
set K0 [expr $M_y/$theta_y]; #elastic stiffness
set as_Plus [expr ((1.25*pow(0.89,$ni)*pow(0.91,(0.01*$fc))*$M_y - $M_y)/(0.13*(1+0.55*$a_sl)*pow(0.13,$ni)*pow((0.02+40*$rosh),0.65)*pow(0.57,(0.01*$fc))))/$K0]; #strain hardening ratio for positive loading direction 1.25*pow(0.89,$ni)*pow(0.91,(0.01*$fc))
set as_Neg [expr ((1.25*pow(0.89,$ni)*pow(0.91,(0.01*$fc))*$M_y - $M_y)/(0.13*(1+0.55*$a_sl)*pow(0.13,$ni)*pow((0.02+40*$rosh),0.65)*pow(0.57,(0.01*$fc))))/$K0]; #strain hardening ratio for negative loading direction
set My_Plus $M_y; #effective yield strength for positive loading direction
set My_Neg -$M_y; #effective yield strength for negative loading direction (negative value)
set Lamda_S [expr 170.7*pow(0.27,$ni)*pow(0.1,$s_d)*(0.13*(1+0.55*$a_sl)*pow(0.13,$ni)*pow((0.02+40*$rosh),0.65)*pow(0.57,(0.01*$fc)))]; #Cyclic deterioration parameter for strength deterioration [expr 170.7*pow(0.27,$ni)*pow(0.1,$s_d)]
set Lamda_C [expr 170.7*pow(0.27,$ni)*pow(0.1,$s_d)*(0.13*(1+0.55*$a_sl)*pow(0.13,$ni)*pow((0.02+40*$rosh),0.65)*pow(0.57,(0.01*$fc)))]; #Cyclic deterioration parameter for post-capping strength deterioration
set Lamda_A [expr 170.7*pow(0.27,$ni)*pow(0.1,$s_d)*(0.13*(1+0.55*$a_sl)*pow(0.13,$ni)*pow((0.02+40*$rosh),0.65)*pow(0.57,(0.01*$fc)))]; #Cyclic deterioration parameter for acceleration reloading stiffness deterioration (is not a deterioration mode for a component with Bilinear hysteretic response).
set Lamda_K [expr 170.7*pow(0.27,$ni)*pow(0.1,$s_d)*(0.13*(1+0.55*$a_sl)*pow(0.13,$ni)*pow((0.02+40*$rosh),0.65)*pow(0.57,(0.01*$fc)))]; #Cyclic deterioration parameter for unloading stiffness deterioration
set c_S 1; #rate of strength deterioration. The default value is 1.0.
set c_C 1; #rate of post-capping strength deterioration. The default value is 1.0.
set c_A 1; #rate of accelerated reloading deterioration. The default value is 1.0.
set c_K 1; #rate of unloading stiffness deterioration. The default value is 1.0.
set theta_p_Plus [expr 0.13*(1+0.55*$a_sl)*pow(0.13,$ni)*pow((0.02+40*$rosh),0.65)*pow(0.57,(0.01*$fc))]; #pre-capping rotation for positive loading direction (often noted as plastic rotation capacity)
set theta_p_Neg [expr 0.13*(1+0.55*$a_sl)*pow(0.13,$ni)*pow((0.02+40*$rosh),0.65)*pow(0.57,(0.01*$fc))]; #pre-capping rotation for negative loading direction (often noted as plastic rotation capacity) (positive value)
set theta_pc_Plus [expr 0.76*pow(0.031,$ni)*pow((0.02+40*$rosh),1.02)]; #post-capping rotation for positive loading direction
set theta_pc_Neg [expr 0.76*pow(0.031,$ni)*pow((0.02+40*$rosh),1.02)]; #post-capping rotation for negative loading direction (positive value)
set Res_Pos 0; #residual strength ratio for positive loading direction
set Res_Neg 0; #residual strength ratio for negative loading direction (positive value)
set theta_u_Plus $theta_u; #ultimate rotation capacity for positive loading direction
set theta_u_Neg $theta_u; #ultimate rotation capacity for negative loading direction (positive value)
set D_Plus 1; #rate of cyclic deterioration in the positive loading direction
set D_Neg 1; #rate of cyclic deterioration in the negative loading direction

################ set new parameters for Bilin "Modified Ibarra-Krawinkler Model with Pinched Hysteretic Response"
set FprPos 1;
set FprNeg 1;
set A_pinch 1;

##### define uniaxialMaterial
uniaxialMaterial ModIMKPinching $matTagBilinPinch $K0 $as_Plus $as_Neg $My_Plus $My_Neg $FprPos $FprNeg $A_pinch $Lamda_S $Lamda_K $Lamda_A $Lamda_C $c_S $c_K $c_A $c_C $theta_p_Plus $theta_p_Neg $theta_pc_Plus $theta_pc_Neg $Res_Pos $Res_Neg $theta_u_Plus $theta_u_Neg $D_Plus $D_Neg

#define elements
#element zeroLength $eleTag $iNode $jNode -mat $matTag1 $matTag2 ... -dir $dir1 $dir2 ... <-orient $x1 $x2 $x3 $yp1 $yp2 $yp3>
element zeroLength $secTag 1 2 -mat $matTagBilinPinch -dir 5 -orient 1 0 0 0 1 0

#define recorders
recorder Node -file rotazione.out -time -node 2 -dof 5 disp
recorder Node -file momento.out -time -node 1 -dof 5 reaction

#define analyses
set P 125000;# -80000.0
pattern Plain 1 "Constant" {
load 2 $P 0.0 0.0 0.0 0.0 0.0
}
integrator LoadControl 0 1

system BandGeneral; #SparseGeneral -piv
test NormUnbalance 1.0e-3 10
numberer Plain
constraints Plain
algorithm NewtonLineSearch -type Secant -tol 0.1 -maxIter 20;
analysis Static
analyze 1
pattern Plain 2 "Linear" {
load 2 0.0 0.0 0.0 0.0 1 0.0
}

#integrator DisplacementControl $node $dof $incr <$numIter $?Umin $?Umax>
set dK [expr $theta_cr/10];
integrator DisplacementControl 2 5 $dK 1
analyze 10
integrator DisplacementControl 2 5 -$dK 1
analyze 20
integrator DisplacementControl 2 5 $dK 1
analyze 10



Thank you in advance!

[vesna]

This material will be in version 2.4.0 of OpenSees which is not yet released. It is expected to be released by the end of next week.

[vesna]

Would you send me an e-mail with this message?
VESNA AT BERKELEY DOTT EDU

[saijbz]

Hi， mariateresa
I also need to use the Modified Ibarra-Medina-Krawinkler Deterioration Model with Pinched Hysteretic Response in order to include pinching effect in my model. can you share your experimece when using uniaxialMaterial ModIMKPInching（I mean how to define the parameters of the model) and if possible the lastest model . my email is saijbz@126.com. Looking forward to your reply.