Pinching Limit State Material: Difference between revisions

From OpenSeesWiki
Jump to navigation Jump to search
No edit summary
mNo edit summary
 
(5 intermediate revisions by the same user not shown)
Line 30: Line 30:
|-
|-
|  '''$driftAxis''' || integer to indicate the drift axis in which lateral-strength degradation will occur. This axis should be orthogonal to the axis of measured rotation (see $rotAxis in Rotation Shear Curve definition)
|  '''$driftAxis''' || integer to indicate the drift axis in which lateral-strength degradation will occur. This axis should be orthogonal to the axis of measured rotation (see $rotAxis in Rotation Shear Curve definition)
driftAxis = 1 – Drift along the x-axis
|-
|  || driftAxis = 1 – Drift along the x-axis
driftAxis = 2 – Drift along the y-axis
driftAxis = 2 – Drift along the y-axis
driftAxis = 3 – Drift along the z-axis
driftAxis = 3 – Drift along the z-axis
Line 37: Line 38:
|-
|-
|  '''$crvTyp''' || integer flag to indicate the type of limit curve associated with this material.
|  '''$crvTyp''' || integer flag to indicate the type of limit curve associated with this material.
crvTyp = 0 – No limit curve
|-
|  || crvTyp = 0 – No limit curve
crvTyp = 1 – axial limit curve
crvTyp = 1 – axial limit curve
crvTyp = 2 – [[RotationShearCurve]]
crvTyp = 2 – [[RotationShearCurve]]
Line 44: Line 46:
|-
|-
|  '''$YpinchUPN''' || floating point unloading force pinching factor for loading in the negative direction
|  '''$YpinchUPN''' || floating point unloading force pinching factor for loading in the negative direction
Note: This value must be between zero and unity
|-
|  || Note: This value must be between zero and unity
|-
|-
|  '''$YpinchRPN''' || floating point reloading force pinching factor for loading in the negative direction
|  '''$YpinchRPN''' || floating point reloading force pinching factor for loading in the negative direction
Note: This value must be between negative one and unity
|-
|  || Note: This value must be between negative one and unity
|-
|-
|  '''$XpinchRPN''' || floating point reloading displacement pinching factor for loading in the negative direction
|  '''$XpinchRPN''' || floating point reloading displacement pinching factor for loading in the negative direction
Note: This value must be between negative one and unity
|-
|  || Note: This value must be between negative one and unity
|-
|-
|  '''$YpinchUNP''' || floating point unloading force pinching factor for loading in the positive direction  
|  '''$YpinchUNP''' || floating point unloading force pinching factor for loading in the positive direction  
Note: This value must be between zero and unity  
|-
|  || Note: This value must be between zero and unity  
|-
|-
|  '''$YpinchRNP''' || floating point reloading force pinching factor for loading in the positive direction  
|  '''$YpinchRNP''' || floating point reloading force pinching factor for loading in the positive direction  
Note: This value must be between negative one and unity
|-
|  || Note: This value must be between negative one and unity
|-
|-
|  '''$XpinchRNP''' || floating point reloading displacement pinching factor for loading in the positive direction
|  '''$XpinchRNP''' || floating point reloading displacement pinching factor for loading in the positive direction
Note: This value must be between negative one and unity
|-
|-
|  '''$dmgStrsLimE''' || floating point force limit for elastic stiffness damage (typically defined as the lowest of shear strength or shear at flexrual yielding). This value is used to compute the maximum deformation at flexural yield (δmax Eq. 1) and using the initial elastic stiffness (Kelastic) the monotonic energy (Emono Eq. 1) to yield. Input 1 if this type of damage is not required and set $dmgE1, $dmgE2, $dmgE3, $dmgE4, and $dmgELim to zero
|  || Note: This value must be between negative one and unity
|-
|  '''$dmgStrsLimE''' || floating point force limit for elastic stiffness damage (typically defined as the lowest of shear strength or shear at flexrual yielding).  
|-
|  || This value is used to compute the maximum deformation at flexural yield (δmax Eq. 1) and using the initial elastic stiffness (Kelastic) the monotonic energy (Emono Eq. 1) to yield. Input 1 if this type of damage is not required and set $dmgE1, $dmgE2, $dmgE3, $dmgE4, and $dmgELim to zero
|-
|  '''$dmgDispMax''' || floating point for ultimate drift at failure (δmax Eq. 1) and is used for strength and stiffness damage.
|-
|-
'''$dmgDispMax''' || floating point for ultimate drift at failure (δmax Eq. 1) and is used for strength and stiffness damage. This value is used to compute the monotonic energy at axial failure (Emono Eq. 2) by computing the area under the backbone in the positive loading direction up to δmax. Input 1 if this type of damage is not required and set $dmgR1, $dmgR2, $dmgR3, $dmgR4, and $dmgRLim to zero for reloading stiffness damage. Similarly set $dmgS1, $dmgS2, $dmgS3, $dmgS4, and $dmgSLim to zero if reloading strength damage is not required
|  || This value is used to compute the monotonic energy at axial failure (Emono Eq. 2) by computing the area under the backbone in the positive loading direction up to δmax. Input 1 if this type of damage is not required and set $dmgR1, $dmgR2, $dmgR3, $dmgR4, and $dmgRLim to zero for reloading stiffness damage. Similarly set $dmgS1, $dmgS2, $dmgS3, $dmgS4, and $dmgSLim to zero if reloading strength damage is not required
|-
|-
|  '''$dmgE1 $dmgE2 $dmgE3 $dmgE4''' || floating point elastic stiffness damage factors ''α1,α2,α3,α4'' shown in Eq. 1
|  '''$dmgE1 $dmgE2 $dmgE3 $dmgE4''' || floating point elastic stiffness damage factors ''α1,α2,α3,α4'' shown in Eq. 1
Line 84: Line 96:




== '''MODE 2: Calibrated for Shear Critical Concrete Columns''' ==
== '''MODE 2: Calibrated Model for Shear-Critical Concrete Columns''' ==




{|  
{|  
| style="background:yellow; color:black; width:800px" | '''$matTag $nodeT $nodeB $driftAxis $Kelas $crvTyp $crvTag $eleTag  
| style="background:yellow; color:black; width:800px" | '''uniaxialMaterial PinchingLimitStateMaterial $matTag $nodeT $nodeB $driftAxis $Kelas $crvTyp $crvTag $eleTag $b $d $h $a $st $As $Acc $ld $db $rhot $f'c $fy $fyt'''
$b $d $h $a $st $As $Acc $ld $db $rhot $f'c $fy $fyt'''
|}
|}


Line 102: Line 113:
|-
|-
|  '''$driftAxis''' || integer to indicate the drift axis in which lateral-strength degradation will occur. This axis should be orthogonal to the axis of measured rotation (see $rotAxis in Rotation Shear Curve definition)
|  '''$driftAxis''' || integer to indicate the drift axis in which lateral-strength degradation will occur. This axis should be orthogonal to the axis of measured rotation (see $rotAxis in Rotation Shear Curve definition)
driftAxis = 1 – Drift along the x-axis
|-
|  || driftAxis = 1 – Drift along the x-axis
driftAxis = 2 – Drift along the y-axis
driftAxis = 2 – Drift along the y-axis
driftAxis = 3 – Drift along the z-axis
driftAxis = 3 – Drift along the z-axis
|-
|-
|  '''$Kelas''' || floating point value to define the shear stiffness (Kelastic) of the shear spring prior to shear failure
|  '''$Kelas''' || floating point value to define the shear stiffness (Kelastic) of the shear spring prior to shear failure
|-
|  ||
Kelas = -4 – Shear stiffness calculated assuming double curvature and shear springs at both column element ends
Kelas = -4 – Shear stiffness calculated assuming double curvature and shear springs at both column element ends
Kelas = -3 – Shear  stiffness calculated assuming double curvature and a shear spring at one column element end
Kelas = -3 – Shear  stiffness calculated assuming double curvature and a shear spring at one column element end
Kelas = -2 – Shear stiffness calculated assuming single curvature and shear springs at both column element ends
Kelas = -2 – Shear stiffness calculated assuming single curvature and shear springs at both column element ends
Kelas = -1 – Shear stiffness calculated assuming single curvature and a shear spring at one column element end
Kelas = -1 – Shear stiffness calculated assuming single curvature and a shear spring at one column element end
Kelas > 0 – Shear stiffness is the input value
Kelas > 0 – Shear stiffness is the input value
Note: integer inputs allow the model to know whether column height equals the shear span (cantelever) or twice the shear span (double curvature). For columns in frames, input the value for the case that best approximates column end conditions or manually input shear stiffness (typically double curvature better estimates framed column behavior)
Note: integer inputs allow the model to know whether column height equals the shear span (cantelever) or twice the shear span (double curvature). For columns in frames, input the value for the case that best approximates column end conditions or manually input shear stiffness (typically double curvature better estimates framed column behavior)
|-
|-
Line 130: Line 149:
|  '''$As''' || floating point total area (inches squared) of longitudinal steel bars in section
|  '''$As''' || floating point total area (inches squared) of longitudinal steel bars in section
|-
|-
|  '''$Acc''' || floating point gross confined concrete area (inches squared) bounded by the transverse reinforcement in column section.
|  '''$Acc''' || floating point gross confined concrete area (inches squared) bounded by the transverse reinforcement in column section
|-
|-
|  '''$ld''' || floating point development length (inches) of longitudinal bars using ACI 318-11 Eq. 12-1 and Eq. 12-2.
|  '''$ld''' || floating point development length (inches) of longitudinal bars using ACI 318-11 Eq. 12-1 and Eq. 12-2
|-
|-
|  '''$db''' || floating point diameter (inches) of longitudinal bars in column section
|  '''$db''' || floating point diameter (inches) of longitudinal bars in column section
|-
|-
|  '''$rhot''' || floating point transverse reinforcement ratio (Ast/st•db)
|  '''$rhot''' || floating point transverse reinforcement ratio (Ast/st.db)
|-
|-
|  '''$f'c''' || floating point concrete compressive strength (ksi)
|  '''$f'c''' || floating point concrete compressive strength (ksi)
Line 143: Line 162:
|-
|-
|  '''$fyt''' || floating point transverse steel yield strength (ksi)
|  '''$fyt''' || floating point transverse steel yield strength (ksi)
|}


----
----


NOTES:
DESCRIPTION:


[[Image:PinchingLimitStateMaterial2.png|550px]]
[[Image:PinchingLimitStateMaterial2.png|550px]]
Line 173: Line 193:


[[PinchingLimitStateMaterial Example]]  
[[PinchingLimitStateMaterial Example]]  
----
DESCRIPTION:




Line 188: Line 201:
1. LeBorgne M. R., 2012, "Modeling the Post Shear Failure Behavior of Reinforced Concrete Columns." Austin, Texas: University of Texas at Austin, PhD, 301.
1. LeBorgne M. R., 2012, "Modeling the Post Shear Failure Behavior of Reinforced Concrete Columns." Austin, Texas: University of Texas at Austin, PhD, 301.


2. LeBorgne M. R. , Ghannoum W. M., 2013, "Analytical Element for Simulating Lateral-Strength Degradation in Reinforced Concrete Columns and Other Frame Members," Journal of Structural Engineering, V. doi: 10.1061/(ASCE)ST.1943-541X.0000925
2. LeBorgne, M.R., Ghannoum, W.M., 2014, "Analytical Element for Simulating Lateral-Strength Degradation in Reinforced Concrete Columns and Other Frame Members," Journal of Structural Engineering, V. 140, No. 7, pp. 04014038 1-12.
 
3.      LeBorgne, M.R., Ghannoum, W.M., 2014, "Calibrated Analytical Element for Lateral-Strength Degradation of Reinforced Concrete Columns," Engineering Structures, V. 81, pp. 35-48.


3. Ghannoum W. M., Moehle J. P., 2012, "Rotation-Based Shear Failure Model for Lightly Confined Reinforced Concrete Columns," Journal of Structural Engineering, V. 138, No. 10, 1267-78.
3. Ghannoum W. M., Moehle J. P., 2012, "Rotation-Based Shear Failure Model for Lightly Confined Reinforced Concrete Columns," Journal of Structural Engineering, V. 138, No. 10, 1267-78.

Latest revision as of 14:49, 10 July 2015




This command is used to construct a uniaxial material that simulates a pinched load-deformation response and exhibits degradation under cyclic loading. This material works with the RotationShearCurve limit surface that can monitor a key deformation and/or a key force in an associated frame element and trigger a degrading behavior in this material when a limiting value of the deformation and/or force are reached. The material can be used in two modes: 1) direct input mode, where pinching and damage parameters are directly input; and 2) calibrated mode for shear-critical concrete columns, where only key column properties are input for model to fully define pinching and damage parameters.


MODE 1: Direct Input

uniaxialMaterial PinchingLimitStateMaterial $matTag $nodeT $nodeB $driftAxis $Kelas $crvTyp $crvTag

$YpinchUPN $YpinchRPN $XpinchRPN $YpinchUNP $YpinchRNP $XpinchRNP $dmgStrsLimE $dmgDispMax $dmgE1 $dmgE2 $dmgE3 $dmgE4 $dmgELim $dmgR1 $dmgR2 $dmgR3 $dmgR4 $dmgRLim $dmgRCyc $dmgS1 $dmgS2 $dmgS3 $dmgS4 $dmgSLim $dmgSCyc


$matTag unique material object integer tag
$nodeT integer node tag to define the first node at the extreme end of the associated flexural frame member (L3 or D5 in Figure)
$nodeB integer node tag to define the last node at the extreme end of the associated flexural frame member (L2 or D2 in Figure)
$driftAxis integer to indicate the drift axis in which lateral-strength degradation will occur. This axis should be orthogonal to the axis of measured rotation (see $rotAxis in Rotation Shear Curve definition)
driftAxis = 1 – Drift along the x-axis

driftAxis = 2 – Drift along the y-axis driftAxis = 3 – Drift along the z-axis

$Kelas floating point value to define the initial material elastic stiffness (Kelastic); Kelas > 0
$crvTyp integer flag to indicate the type of limit curve associated with this material.
crvTyp = 0 – No limit curve

crvTyp = 1 – axial limit curve crvTyp = 2 – RotationShearCurve

$crvTag integer tag for the unique limit curve object associated with this material
$YpinchUPN floating point unloading force pinching factor for loading in the negative direction
Note: This value must be between zero and unity
$YpinchRPN floating point reloading force pinching factor for loading in the negative direction
Note: This value must be between negative one and unity
$XpinchRPN floating point reloading displacement pinching factor for loading in the negative direction
Note: This value must be between negative one and unity
$YpinchUNP floating point unloading force pinching factor for loading in the positive direction
Note: This value must be between zero and unity
$YpinchRNP floating point reloading force pinching factor for loading in the positive direction
Note: This value must be between negative one and unity
$XpinchRNP floating point reloading displacement pinching factor for loading in the positive direction
Note: This value must be between negative one and unity
$dmgStrsLimE floating point force limit for elastic stiffness damage (typically defined as the lowest of shear strength or shear at flexrual yielding).
This value is used to compute the maximum deformation at flexural yield (δmax Eq. 1) and using the initial elastic stiffness (Kelastic) the monotonic energy (Emono Eq. 1) to yield. Input 1 if this type of damage is not required and set $dmgE1, $dmgE2, $dmgE3, $dmgE4, and $dmgELim to zero
$dmgDispMax floating point for ultimate drift at failure (δmax Eq. 1) and is used for strength and stiffness damage.
This value is used to compute the monotonic energy at axial failure (Emono Eq. 2) by computing the area under the backbone in the positive loading direction up to δmax. Input 1 if this type of damage is not required and set $dmgR1, $dmgR2, $dmgR3, $dmgR4, and $dmgRLim to zero for reloading stiffness damage. Similarly set $dmgS1, $dmgS2, $dmgS3, $dmgS4, and $dmgSLim to zero if reloading strength damage is not required
$dmgE1 $dmgE2 $dmgE3 $dmgE4 floating point elastic stiffness damage factors α1,α2,α3,α4 shown in Eq. 1
$dmgELim floating point elastic stiffness damage limit Dlim shown in Eq. 1; Note: This value must be between zero and unity
$dmgR1 $dmgR2 $dmgR3 $dmgR4 floating point reloading stiffness damage factors α1,α2,α3,α4 shown in Eq. 1
$dmgRLim floating point reloading stiffness damage limit Dlim shown in Eq. 1; Note: This value must be between zero and unity
$dmgRCyc floating point cyclic reloading stiffness damage index; Note: This value must be between zero and unity
$dmgS1 $dmgS2 $dmgS3 $dmgS4 floating point backbone strength damage factors α1,α2,α3,α4 shown in Eq. 1
$dmgSLim floating point backbone strength damage limit Dlim shown in Eq. 1; Note: This value must be between zero and unity
$dmgSCyc floating point cyclic backbone strength damage index; Note: This value must be between zero and unity


MODE 2: Calibrated Model for Shear-Critical Concrete Columns

uniaxialMaterial PinchingLimitStateMaterial $matTag $nodeT $nodeB $driftAxis $Kelas $crvTyp $crvTag $eleTag $b $d $h $a $st $As $Acc $ld $db $rhot $f'c $fy $fyt

$matTag unique material object integer tag
$nodeT integer node tag to define the first node at the extreme end of the associated flexural frame member (L3 or D5 in Figure)
$nodeB integer node tag to define the last node at the extreme end of the associated flexural frame member (L2 or D2 in Figure)
$driftAxis integer to indicate the drift axis in which lateral-strength degradation will occur. This axis should be orthogonal to the axis of measured rotation (see $rotAxis in Rotation Shear Curve definition)
driftAxis = 1 – Drift along the x-axis

driftAxis = 2 – Drift along the y-axis driftAxis = 3 – Drift along the z-axis

$Kelas floating point value to define the shear stiffness (Kelastic) of the shear spring prior to shear failure

Kelas = -4 – Shear stiffness calculated assuming double curvature and shear springs at both column element ends

Kelas = -3 – Shear stiffness calculated assuming double curvature and a shear spring at one column element end

Kelas = -2 – Shear stiffness calculated assuming single curvature and shear springs at both column element ends

Kelas = -1 – Shear stiffness calculated assuming single curvature and a shear spring at one column element end

Kelas > 0 – Shear stiffness is the input value

Note: integer inputs allow the model to know whether column height equals the shear span (cantelever) or twice the shear span (double curvature). For columns in frames, input the value for the case that best approximates column end conditions or manually input shear stiffness (typically double curvature better estimates framed column behavior)

$crvTag integer tag for the unique limit curve object associated with this material
$eleTag integer element tag to define the associated beam-column element used to extract axial load
$b floating point column width (inches)
$d floating point column depth (inches)
$h floating point column height (inches)
$a floating point shear span length (inches)
$st floating point transverse reinforcement spacing (inches) along column height
$As floating point total area (inches squared) of longitudinal steel bars in section
$Acc floating point gross confined concrete area (inches squared) bounded by the transverse reinforcement in column section
$ld floating point development length (inches) of longitudinal bars using ACI 318-11 Eq. 12-1 and Eq. 12-2
$db floating point diameter (inches) of longitudinal bars in column section
$rhot floating point transverse reinforcement ratio (Ast/st.db)
$f'c floating point concrete compressive strength (ksi)
$fy floating point longitudinal steel yield strength (ksi)
$fyt floating point transverse steel yield strength (ksi)

DESCRIPTION:


The material model coupled with the RotationShearCurve limit surface: 1) has the ability to continually monitor forces and deformations in the flexural elements for conditions that trigger lateral-strength degradation, 2) has a built-in function that compensates for flexural deformation offsets that arise from the degrading behavior of the material in shear springs, and 3) is able to trigger lateral-strength degradation through either a limiting lateral force or element deformations (whichever is reached first). The material introduces several functionalities that give users a high degree of control over the triggering of strength degradation and the ensuing cyclic degrading behavior. Damage algorithms are implemented to control the degrading behavior through elastic stiffness, reloading stiffness, and backbone strength degradation (Fig. 2). The rate of damage accumulation can be controlled by energy-, displacement-, and cycle-based damage computation algorithms.

During the degrading behavior, the model automatically adjusts reloading stiffness to achieve a symmetric global-element lateral load-vs lateral displacement behavior. The model does so by automatically adjusting the reloading stiffness and backbone curve of the material model to compensate for dissymmetry introduced by the unloading of the flexural elements in series with shear springs governed by the model.


DAMAGE:

Damage accumulations effects based on numbers of cycles can be introduced to reloading stiffness and backbone strength through the simple parameters $dmgRCyc and $dmgSCyc with values ranging from 0 to 1.

Elastic stiffness, reloading stiffness, and strength can be adjusted using the following energy and displacement damage model (from Mitra and Lowes (2007)):



EXAMPLE:

PinchingLimitStateMaterial Example



REFERENCES:

1. LeBorgne M. R., 2012, "Modeling the Post Shear Failure Behavior of Reinforced Concrete Columns." Austin, Texas: University of Texas at Austin, PhD, 301.

2. LeBorgne, M.R., Ghannoum, W.M., 2014, "Analytical Element for Simulating Lateral-Strength Degradation in Reinforced Concrete Columns and Other Frame Members," Journal of Structural Engineering, V. 140, No. 7, pp. 04014038 1-12.

3. LeBorgne, M.R., Ghannoum, W.M., 2014, "Calibrated Analytical Element for Lateral-Strength Degradation of Reinforced Concrete Columns," Engineering Structures, V. 81, pp. 35-48.

3. Ghannoum W. M., Moehle J. P., 2012, "Rotation-Based Shear Failure Model for Lightly Confined Reinforced Concrete Columns," Journal of Structural Engineering, V. 138, No. 10, 1267-78.

4. Mitra Nilanjan, Lowes Laura N., 2007, "Evaluation, Calibration, and Verification of a Reinforced Concrete Beam--Column Joint Model," Journal of Structural Engineering, V. 133, No. 1, 105-20.




Code Developed by: Matthew Leborgne and Wassim M. Ghannoum, University of Texas at Austin