Elmer Discussion Forum

Hello, I am writing a Solver, I have a 2D body with different materials and I need to define the boundary conditions between the materials. For that purpose I need the normal vector. In "ElmerSolver Manual" is written that the normal vector "is pointing towards the body with lower den...

Hi Peter,

Here is attached the small example. It presents a cantilever beam excited with harmonic force on the free edge. I compared the results with Timoshenko beam model and they are in agreement.

Thank you for your help!

Best regards,
Stan

Hello, I am using Elmer version 6.1. on Ubuntu 64 bits. I have a problem when I use "Calculate Velocity = Logical True" in the Solution section. The velocity is written at the .ep file but the following error appears at the end of the program: ElmerSolver: *** Elmer Solver: ALL DONE *** **...

Hello, I think that for second order problems, like the Linear Elasticity problem, the displacement and the velocity are necessary to continue from existing solution. Thus the vector of velocity has to be written in the Restart File, i.e. one should use "Calculate Velocity = Logical True" ...

Hello Elmer users, I think if someone needs to have the mass and the stiffness matrices, they can be accessed through the pointer Matrix_t: TYPE(Matrix_t), POINTER :: A A => Solver % Matrix A % BulkValues ! pointer to the values of the stiffness matrix A % MassValues ! pointer to the values of the m...

Hi Peter, There is a small problem. I am using the vector of solutions, i.e. the values of the nodal points. I take it from the solver, after calling DefaultSolve(), by vector => Solver % Variable % Values When I generate the mesh using Elmer the values are the same as in .ep file. Bu if I generate ...

Thank you again Peter!

It was not difficult to implement such a function by using the values of the shape functions at that point and the corresponding values of the nodal coordinates.

Regards,
Stan

Hello, I need to obtain the solution of a problem at an arbitrary point inside the domain. I did not find an appropriate subroutine for this purpose. Is there a subroutine which gives you this value of a point, if you give as input the coordinates of the point? Something like: GetResult(x, y, z, res...

Thank you Peter,

both functions, GetElementNOFNodes and GetElementNOFDOFs, work well. I will use them instead of Nodes % NumberOfNodes.

Regards,
Stan

Hello, I am using the subroutine GetElementNodes( Nodes ), where Nodes is of type Nodes_t and after calling this subroutine it gets the nodes of the active element. I suppose that there is a problem with the type Nodes_t or with the subroutine GetElementNodes because the attribute NumberOfNodes is z...