I am trying to model an energy dissipater element that can undergo plastic deformation when contacted by a body, e.g. elastic perfectly plastic relationship
I have tried using a sensor as described in Article Number
000002688:
What I have so far:
-----------------------
Sensor event definition
.MODEL_1.Deformation - SENVAL(DeformationSensor)
Sensor event evaluation:
.MODEL_1.Deformation
SForce:
IF( .MODEL_1.Deformation: forceValue, 0, 0 )
Questions:
-------------
I am finding it difficult to figure out how to incorporate the plastic deformation. According to me plasticDeformation = SENVAL(DeformationSensor) therefore the force becomes (where force value: f= min(kx, plasticForceValue):
IF( .MODEL_1.Deformation + plasticDeformation: forceValue, 0, 0 ) but this is obviously incomplete and incorrect.
Is the sensor approach correct for this application, or do I need to look at doing a subroutine? Is there another approach that will work?
OK, here is an example for you. Rename the attachment to ElastoPlast.cmd and import into View.
Simulate for 10 s, 500 steps. Plot request .ElastoPlast.Last_Run.Ground_Data.Ground_Force vs .ElastoPlast.Last_Run.Ground_Data.Deformation
This model has three parameters:
PlastDef0: The deformation that defines the "knee" of the elasto-plastic curve. Set to zero for pure plastic behavior.
K_elast: Stiffness of the elastic part (as well as off-loading behavior). Set to very high for pure plastic behavior.
K_plast: Stiffness of the plastic deformation. Set to a low value for true plastic.
You can probably change the model to be true constant plastic behavior by setting PlastDef0=0, K_elast very high (>1e6) and changing the solver variable Plast_Force to a constant value.
OK, here is an example for you. Rename the attachment to ElastoPlast.cmd and import into View.
Simulate for 10 s, 500 steps. Plot request .ElastoPlast.Last_Run.Ground_Data.Ground_Force vs .ElastoPlast.Last_Run.Ground_Data.Deformation
This model has three parameters:
PlastDef0: The deformation that defines the "knee" of the elasto-plastic curve. Set to zero for pure plastic behavior.
K_elast: Stiffness of the elastic part (as well as off-loading behavior). Set to very high for pure plastic behavior.
K_plast: Stiffness of the plastic deformation. Set to a low value for true plastic.
You can probably change the model to be true constant plastic behavior by setting PlastDef0=0, K_elast very high (>1e6) and changing the solver variable Plast_Force to a constant value.