Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Arc length method

Please login with a confirmed email address before reporting spam

Hi! I am modelling a snap-back behavior in a structural mechanics problem. I am wondering how I can implement the arc length method in COMSOL. Is there any built-in functionality in COMSOL that I can use or is there any way to make COMSOL compute interactively with Matlab? Thanks!

Regards, Tianyu


9 Replies Last Post 01.05.2021, 08:54 GMT-4
Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 4 years ago 20.04.2021, 10:35 GMT-4
Updated: 4 years ago 20.04.2021, 10:38 GMT-4

Here are two examples:

https://www.comsol.com/model/spherical-cap-with-central-point-load-67161

https://www.comsol.com/model/postbuckling-analysis-of-a-hinged-cylindrical-shell-10257

-------------------
Henrik Sönnerlind
COMSOL
Here are two examples:

Please login with a confirmed email address before reporting spam

Posted: 4 years ago 21.04.2021, 23:16 GMT-4
Updated: 4 years ago 21.04.2021, 23:17 GMT-4

Thank you sir! I have reviewed these 2 models and I found the second one similar with my problem. However, the second model is time-independent. Do we have arc length method implemented for time-dependent problems? Thanks.

Thank you sir! I have reviewed these 2 models and I found the second one similar with my problem. However, the second model is time-independent. Do we have arc length method implemented for time-dependent problems? Thanks.

Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 4 years ago 23.04.2021, 03:38 GMT-4

In a time dependent analysis, the inertia will stabilize the structure. You can just use a load that increases with time.

That is what happens during a snap-through event in real life: When the limit point is reached, the applied load is balanced by inertial forces. See the attached plot. The vibrations is what causes the characteristic sound emission.

-------------------
Henrik Sönnerlind
COMSOL
In a time dependent analysis, the inertia will stabilize the structure. You can just use a load that increases with time. That is what happens during a snap-through event in real life: When the limit point is reached, the applied load is balanced by inertial forces. See the attached plot. The vibrations is what causes the characteristic sound emission.


Please login with a confirmed email address before reporting spam

Posted: 4 years ago 25.04.2021, 08:36 GMT-4

In a time dependent analysis, the inertia will stabilize the structure. You can just use a load that increases with time.

That is what happens during a snap-through event in real life: When the limit point is reached, the applied load is balanced by inertial forces. See the attached plot. The vibrations is what causes the characteristic sound emission.

Thank you very much!

>In a time dependent analysis, the inertia will stabilize the structure. You can just use a load that increases with time. > >That is what happens during a snap-through event in real life: When the limit point is reached, the applied load is balanced by inertial forces. See the attached plot. The vibrations is what causes the characteristic sound emission. Thank you very much!

Please login with a confirmed email address before reporting spam

Posted: 4 years ago 26.04.2021, 03:23 GMT-4
Updated: 4 years ago 26.04.2021, 10:35 GMT-4

Hi! I have tried this method and also the parametric sweep and they work under the assumption that the loading type is force. However, I still have one question, what if the loading type is prescribed displacement? For example, in the first pitcure is an single edge notched specimen under tension and the loading condition is a prescribed displacement on the upper bound. In my time-dependent simulation, I used displacement control and the computation diverges as the second pitcure indicates and I think it is because equilibrium path has a snap-back behaviour. So how can I use arc length method in this example when the loading condition requires that the upper bound of the specimen moves identically in the y direction? Thank you again for your patience.

Hi! I have tried this method and also the parametric sweep and they work under the assumption that the loading type is force. However, I still have one question, what if the loading type is prescribed displacement? For example, in the first pitcure is an single edge notched specimen under tension and the loading condition is a prescribed displacement on the upper bound. In my time-dependent simulation, I used displacement control and the computation diverges as the second pitcure indicates and I think it is because equilibrium path has a snap-back behaviour. So how can I use arc length method in this example when the loading condition requires that the upper bound of the specimen moves identically in the y direction? Thank you again for your patience.


Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 4 years ago 28.04.2021, 09:44 GMT-4

I guess that you are working with a deteriorating material model, like a cohesive zone model. That is a much nastier type of 'negative stiffness' than the ordinary snap-through problem. The negative stiffness will occur locally in several elements, and it can be extremely difficult to find a solution.

If you do it in time domain, stress waves will start traveling from the crack tip, out to the boundaries, and then be reflected back.

In this case, you may need some regularization method, like the ones available in the Damage material model, or in Decohesion under Contact.

-------------------
Henrik Sönnerlind
COMSOL
I guess that you are working with a deteriorating material model, like a cohesive zone model. That is a much nastier type of 'negative stiffness' than the ordinary snap-through problem. The negative stiffness will occur locally in several elements, and it can be extremely difficult to find a solution. If you do it in time domain, stress waves will start traveling from the crack tip, out to the boundaries, and then be reflected back. In this case, you may need some regularization method, like the ones available in the Damage material model, or in Decohesion under Contact.

Please login with a confirmed email address before reporting spam

Posted: 4 years ago 28.04.2021, 10:23 GMT-4

I guess that you are working with a deteriorating material model, like a cohesive zone model. That is a much nastier type of 'negative stiffness' than the ordinary snap-through problem. The negative stiffness will occur locally in several elements, and it can be extremely difficult to find a solution.

If you do it in time domain, stress waves will start traveling from the crack tip, out to the boundaries, and then be reflected back.

In this case, you may need some regularization method, like the ones available in the Damage material model, or in Decohesion under Contact.

Thanks for your reply. I am doing quasi-static analysis of a cohesive zone model implemented with phase field method. The post-peak behaviour is very unstable and the author of the paper mentioned that CMOD control should be used to track the post-peak behaviour. I tried parametric sweep with force type loading condition and it works well but I don't know how to use this method with displacement type loading condition.

>I guess that you are working with a deteriorating material model, like a cohesive zone model. That is a much nastier type of 'negative stiffness' than the ordinary snap-through problem. The negative stiffness will occur locally in several elements, and it can be extremely difficult to find a solution. > >If you do it in time domain, stress waves will start traveling from the crack tip, out to the boundaries, and then be reflected back. > >In this case, you may need some regularization method, like the ones available in the Damage material model, or in Decohesion under Contact. Thanks for your reply. I am doing quasi-static analysis of a cohesive zone model implemented with phase field method. The post-peak behaviour is very unstable and the author of the paper mentioned that CMOD control should be used to track the post-peak behaviour. I tried parametric sweep with force type loading condition and it works well but I don't know how to use this method with displacement type loading condition.

Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 4 years ago 29.04.2021, 07:08 GMT-4

I have never seen an example where a prescribed displacement is controlled in this way. It is possible that there can be some subtle reason for why it would be problematic. If so, you can try replacing the prescribed displacement with a rigid connector, and control the force on the rigid connector from the CMOD. But that is just a guess from my side.

-------------------
Henrik Sönnerlind
COMSOL
I have never seen an example where a prescribed displacement is controlled in this way. It is possible that there can be some subtle reason for why it would be problematic. If so, you can try replacing the prescribed displacement with a rigid connector, and control the force on the rigid connector from the CMOD. But that is just a guess from my side.

Please login with a confirmed email address before reporting spam

Posted: 4 years ago 01.05.2021, 08:54 GMT-4

I have never seen an example where a prescribed displacement is controlled in this way. It is possible that there can be some subtle reason for why it would be problematic. If so, you can try replacing the prescribed displacement with a rigid connector, and control the force on the rigid connector from the CMOD. But that is just a guess from my side.

Got it. Thank you very much!

>I have never seen an example where a prescribed displacement is controlled in this way. It is possible that there can be some subtle reason for why it would be problematic. If so, you can try replacing the prescribed displacement with a rigid connector, and control the force on the rigid connector from the CMOD. But that is just a guess from my side. Got it. Thank you very much!

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.