|
From: | Konstantinos Poulios |
Subject: | Re: rigid contact with arbitrary surface |
Date: | Tue, 19 Apr 2022 13:16:17 +0200 |
Dear all
I have tried the following rigid contact option for an arbitrary 2D surface obstacle and a 3D mesh in python scripting
It looks quite successful
However it takes a lot of time for me to finish a viscoelastic computation (28h)
How can I speed it up? By reducing the number of Newton iteration for example?
I can read in the documentation that there are three other options
How to they compare in terms of speed and accuracy?
How do they compare to a nonlinear stiffness penalty method between the rigid contact surface considered as master and the 3D solid mesh considered as slave nodes sets?
Thank you
Regards
Anne-Cecile
ldof = mflambda.dof_on_region(BRAIN_BOUND)
mflambda_partial = gf.MeshFem('partial', mflambda, ldof)
md.add_fem_variable('lambda_n', mflambda_partial)
md.add_initialized_data('r', [r])
md.add_initialized_fem_data('obstacle', mfpb, levelset)
md.add_integral_contact_with_rigid_obstacle_brick(mim_friction, 'ub', 'lambda_n',
'obstacle', 'r', BRAIN_BOUND, 1);
The information contained in this e-mail message may be privileged, confidential, and/or protected from disclosure. This e-mail message may contain protected health information (PHI); dissemination of PHI should comply with applicable federal and state laws. If you are not the intended recipient, or an authorized representative of the intended recipient, any further review, disclosure, use, dissemination, distribution, or copying of this message or any attachment (or the information contained therein) is strictly prohibited. If you think that you have received this e-mail message in error, please notify the sender by return e-mail and delete all references to it and its contents from your systems.
[Prev in Thread] | Current Thread | [Next in Thread] |