Hello Ana Regina,

I expect the reason your model has very little stress in the artery is that there is no initial interference between the rod and the artery. The geometry was constructed with the same diameter for the rod and artery ID, so I am not surprised that there is no stress.

When a stent is expanded in an artery by a balloon, it stretches the artery. When the balloon deflates, the stretched artery applies a pressure to the stent, squeezing it. This is the state you need the model to be in before you start extracting the rod. There are several ways to put the model in that state.

The first way is to assign an initial state of stress on the artery. There is a command called INISTATE. This method is quite technical to implement, so let’s leave that aside for now.

A second way is to construct the geometry so that the stent (rod) diameter is larger than the ID of the artery. This would be a 2-step solution. In step 1, frictional contact between the artery and stent would be resolved. The contact algorithm would detect the interference and ramp up the contact force until the interference is driven to practically zero, stretching the artery in the process. In step 2, you pull the stent out of the artery. This can work for small values of interference.

The third way is to use a 2-step solution, but start with the geometry you have. In step 1, apply a thermal condition to raise the temperature on the rod, causing it to grow in diameter. The rod material will be given a Secant Coefficient of Thermal Expansion (CTE) material property using a number you will calculate to give the desired interference between the stent and the artery for the temperature load you apply. You can use orthotropic properties instead of isotropic properties so only the diameter gets larger, but the length doesn’t change because the CTE in the axial direction is 0. In step 1, the rod temperature increases by 1 degree C, the rod grows in diameter, stretching the artery. In step 2, the rod temperature stays the same and the stent is pulled out of the artery by the remote displacement.

In all these cases, delete the Fixed Supports from all the faces of the artery. You need a Remote Displacement, Behavior = Flexible, applied to the flat end face of the artery on the end opposite the direction the rod is being pulled. Set all six DOF to zero.

You also need a Remote Displacement, Behavior = Flexible on the flat end face of the rod on the end that is in the direction the rod will be pulled.  Set five DOF to zero, leaving the axial direction Free. Promote that Remote Displacement to a Remote Point. Apply another Remote Displacement, but scope it to the Remote Point you just created. This Remote Displacement has a Tabular input for the 2-step solution. The axial displacement is 0 in step 1, which holds the rod in place while the artery expands in step 1. The pull distance is applied in step 2.

For this initial simple model, I suggest you make the rod longer than the artery. Edit the geometry of the rod and put a generous blend radius on the sharp edge at the end face. The reason is that as the rod is removed, the artery, which is stretched over the rod, will want to contract to its original diameter. It will be numerically very difficult for the nodes of the artery to slide off a sharp edge of the rod.  Make that end of the rod have a dome shape, or at least a blend that leaves a very small flat.

I hope that gets you started. I will check back for any follow-up questions.

Kind regards,