Modeling balloon expansion inside coronary artery
Modeling balloon expansion inside coronary artery
(OP)
As suggested by the title, I am attempting to model a balloon expansion inside a coronary artery. The balloon is a uniform cylinder, whereas the artery is a non-uniform cylinder (See attached picture). I want to apply a given pressure inside of the balloon, causing it to expand until it touches the vessel. Then, through contact, it should cause the vessel itself to expand with it.
I have specified a contact interaction between the outer surface of the balloon and the inner surface of the vessel. However, when I apply a Pressure Load on the inner surface of the balloon, it ends up ignoring any contact with the vessel and simply expanding as if it were alone.
When I try to apply a Body Force Load instead of Pressure Load, the balloon moves in the specified direction and interacts with the vessel as it should, causing the vessel to also move in the same direction through contact. However, this is not ideal since I cannot make a Body Force act in all directions (as the Pressure Load does).
Does anybody have any tips? Either in how to make the pressure load work or how to apply another type of force that would act in all directions?
I would be happy to provide more details if needed. Thanks!
I have specified a contact interaction between the outer surface of the balloon and the inner surface of the vessel. However, when I apply a Pressure Load on the inner surface of the balloon, it ends up ignoring any contact with the vessel and simply expanding as if it were alone.
When I try to apply a Body Force Load instead of Pressure Load, the balloon moves in the specified direction and interacts with the vessel as it should, causing the vessel to also move in the same direction through contact. However, this is not ideal since I cannot make a Body Force act in all directions (as the Pressure Load does).
Does anybody have any tips? Either in how to make the pressure load work or how to apply another type of force that would act in all directions?
I would be happy to provide more details if needed. Thanks!





RE: Modeling balloon expansion inside coronary artery
and how have you defined the interaction between the balloon/artery?
also, can you describe how you are modelling the balloon - they are usually very thin and folded about themselves?
RE: Modeling balloon expansion inside coronary artery
The interaction properties are hard contact for normal behavior and penalty with specified friction coefficient for tangential behavior.
As of now I am modelling the balloon as a uniform, thin, isotropic elastic cylinder with specified parameters.
RE: Modeling balloon expansion inside coronary artery
RE: Modeling balloon expansion inside coronary artery
RE: Modeling balloon expansion inside coronary artery
- Check if you have positive results for CPRESS. If yes, then the contact works and your penalty stiffness might be too low.
RE: Modeling balloon expansion inside coronary artery
By the way, it may be that the perspective is making me think the geometry has sharp transitions so take this with a grain of salt. You might want to consider smoothening the geometry a bit. It will help the contact solver. If I were doing this analysis, I'd start with a much smaller and simpler geometry to have the workflow down first.
By the way, no one seems to take into account the moving elastic foundation of the heart in to account in these analyses. Perhaps, BCs are challenging to account for, I do not know. However, it may be an interesting avenue to explore given the rapid advances tagged MRI seems to have made in the last few years.
*********************************************************
Are you new to this forum? If so, please read these FAQs:
http://www.eng-tips.com/faqs.cfm?fid=376
http://www.eng-tips.com/faqs.cfm?fid=1083
RE: Modeling balloon expansion inside coronary artery
Like Mustaine suggested, I would check the contact output and diagnostics to figure out your issue. You would also see severe discontinuity iterations in the .sta file in any increments where contact constraints change state (open/closed).
@IceBreakerSours
this paper incorporates the movement of the cardiac tissue. Only one I'm aware of:
https://doi.org/10.1016/j.jbiomech.2014.01.007
RE: Modeling balloon expansion inside coronary artery
It turns out that only the first section that comes into contact has a positive CPRESS value. Everywhere else, the balloon goes right through the vessel and the respective CPRESS value is 0. My guess is the vessel geometry needs to be smoothed as @IceBreakerSours mentioned.
RE: Modeling balloon expansion inside coronary artery
@Andre1212
Contact solver likes the opposing surfaces to be "similar". The closer in "similarity" they are, the better it is for the solver. By "similar", I mean, in stiffness (material and geometric contributions) and discretization (mesh). The more "dissimilar" the opposing surfaces are, the trickier the contact problem starts to become. Metallic stents coming to contact with a soft vessel falls in the latter category. When you introduce geometric features such as sharp edges at the interface, you are just making life a hell for the contact solver and, consequently, for yourself :)
By the way, I am not entirely convinced if finite sliding is required while you are still getting your model to work. Friction coefficients seem to appear out of thin air in these analyses so I'd rather get the small sliding model running first to have something meaningful to report and then consider whether finite sliding is a more appropriate interaction or not.
*********************************************************
Are you new to this forum? If so, please read these FAQs:
http://www.eng-tips.com/faqs.cfm?fid=376
http://www.eng-tips.com/faqs.cfm?fid=1083
RE: Modeling balloon expansion inside coronary artery
Is NLGEOM active?
What happens with General Contact?
And no, I don't think that Small Sliding is applicable here - independent of friction. The displacements and deformations will be too large and also be relative to each other.
I also guess, that the overall goal of the simulation might be easier to achieve with quasi-static in A/Explicit. With the complex contact situation, especially when both expand while being in contact, it will be hard to get it running smooth in an implicit solver.
What happens in the z-directions? Maybe it is possible to run it in 2D...
RE: Modeling balloon expansion inside coronary artery
RE: Modeling balloon expansion inside coronary artery
Assuming surfaces and meshes are smooth enough, /Standard can get the job done. While I will admit it is not easy, it is do-able and I mention this because of the benefits of running these jobs in implicit. That said, if you can not get over the hump, I would follow the advice above and switch over to /Explicit.
*********************************************************
Are you new to this forum? If so, please read these FAQs:
http://www.eng-tips.com/faqs.cfm?fid=376
http://www.eng-tips.com/faqs.cfm?fid=1083