Elmer Discussion Forum

Dear All, I am looking in the development of the Finite Elasticity model and the applied linearization (for the Newton's method), and I wish to ask if someone can confirm or correct my understanding of the code. In the definition of the LocalMatrix at file ElasticSolve.src, the force vestor is defin...

Hi,

on my experience, you may use umfpack for 3D problems if the DOF are about (or less) 400 000.
If the system is begger, use MUMPS with parallel processors. It works on my case for 2-3 millions DOF.

Regards,
Stan

Hi Peter, Thank you for the suggestion of using block preconditioner. Unfotionately I couldn't obtain convergence in the sequential version, using your example. What I cannot understand is why this example works for 2D probelems (the solver example is similar to BlockLinElast3), but it does not work...

Hello Elmer team, I am modeling elastic structures with geometrical nonlinearity and I use ElasticSolver. I obtain a solution when I use different iterative methods on single processor, but when I go to parallel computations I cannot obtain a solution, i.e. the method does not converge. Here is an e...

Yes, I am thinking to implement the modifications and write to the developers, so if they agree they might upload the source.

Thank you Peter, yes, in the case of solver, I may compile and use it through the sif file. But if I use the StressSolver and I need to change, for example, the SUBROUTINE Add2ndOrderTime, which is called from StreesSolver through Default2ndOrderTime, then I would need to compile SolverUtilis. But s...

Hello Elmer team, In some of my old posts I described about a problem that in 2nd order time-depenent systems, the solution obtained by restarting from previous solution differs from solution obtained without restarting. I noticed this when I used it for elasticity problems, such as StressSolve, Ela...

Hello, If I need to modify some of the source files from /fem/src, how to update the changes in the libraries? One way is to run the configure file and reinstall the complete Elmer. But there should be a way to compile only the modified file and add it to the library /lib/libelmersolver.so. Could so...

Hello again, I am back to the problem with continuation from previous solution. Still I cannot obtain the same solution if I restart from some time step (the difference is small, but there should be no difference). For example, if I find the solutions for time steps from 0 to 10, then if I restart f...

Hi Juha, Thank you for the answer! Yes, one option is to use "Normal Target Body = Integer n", but I wish to do general program which works for complex domains and to avoid specifying the direction of the normal vector manually. Thus, I would prefer to use the default normal, i.e. as it po...