imcs-compsim / pinns_for_comp_mech

Implementation of Physics-Informed Neural Networks for Computational Mechanics based on the DeepXDE package.
MIT License
20 stars 3 forks source link

Beam contac #19

Open HumbertHumbert7 opened 5 months ago

HumbertHumbert7 commented 5 months ago

IMG_5340

Is this problem implemented? How can I solve it? Thank you

permas4edu commented 4 months ago

I guess you ask for a static contact analysis. This problem can be solved in PERMAS. Can you provide the geometrical (i.e. cross-section, length) and material properties (Young's modulus) of the elastic beam? Is it a distributed load? How large is the load?

tsarikahin commented 4 months ago

@HumbertHumbert7 Sorry for the late answer.

This problem is not implemented. But by applying Karush-Kuhn-Tucker inequalities, this can be solved. However, KKT constraints are formulated only for 2D domains. You need to modify them based on the 1D beam example.

Steps:

To enforce KKT constraints, you can use the Fisher-Burmeister function, but still, you have to calculate the gap and normal pressure (Normal pressure would be in this case the pressure in vertical direction, so the shear stress)

Thats basically all.

HumbertHumbert7 commented 4 months ago

Hi, thank you very much for your responses. it doesn't matter if the problem is not implemented, I have considered some other contact mechanics problem for which I have proceeded as indicated by tsarikahin and I have founded very useful the paper you have published.

tsarikahin commented 4 months ago

@HumbertHumbert7 Thank you for your kind words. Please feel free to make a PR to add your examples.

permas4edu commented 4 months ago

beam_in_contact_with_ground

permas4edu commented 4 months ago

The clamped-clamped beam is loaded by a line load. A node to ground contact with an initial gap is used here. Moreover, a hard contact is used in PERMAS.