Recent content by famatto

I'm using an E=2000000 PSI which is common for wood. Obviously I do know that things have some level of margin for error here to rein in. There are maybe hierarchal methods to use, that I am working towards... Mostly the talk about the fudge factor is so far unconstructive, and its a very common...

One thing lately I have a quandry. There are hermitian shape functions: \begin{align*} H_1(\xi) &= \frac{1}{4}\xi^3 - \frac{3}{4}\xi + \frac{1}{2} \\ H_2(\xi) &= \frac{1}{4}\xi^3 - \frac{1}{4}\xi^2 - \frac{1}{4}\xi + \frac{1}{4} \\ H_3(\xi) &= -\frac{1}{4}\xi^3 + \frac{3}{4}\xi + \frac{1}{2} \\...

``` def strain(u): return 0.5 * (ufl.grad(u) + ufl.grad(u).T) epsilon = strain(u) def stress(epsilon, mu, lambda_): return 2 * mu * epsilon + lambda_ * ufl.tr(epsilon) * ufl.Identity(3) sigma = stress(epsilon, mu, lambda_) ``` The above is how stress and strain are being described. u...

So I had to put a scale factor so far of -8.66071e10, however it seems to be working... Anyone can tell me how to get things going without the scaling? It sounds like there is maybe a Jacobian to apply, maybe? So far the mesh can be exported with dolfinx 0.8.0... the rest of the code will run...