Hi Michael,
Thank you for your response.
I agree that approach can weakly enforce Dirichlet BCs as required. There is an issue I foresee specifically for the kind of problems I'm working on.
1. I have a formulation that already uses Lagrange multipliers to impose certain conservation constraints.
2. The governing equations are posed in mixed form. So such constraints have to be imposed on both displacements and tractions (Neumann BCs related to stress at the boundary surfaces).
Using the penalty approach to apply constraints to a system with the above characteristics will likely affect the system of equations adversely. If possible, I'd prefer not to degrade the linear algebraic system with additional penalty constraints.
Is there a different approach to implement the constraint equations? (I am fine even if it's painful to implement
Thank you,
Anand