Elmer Discussion Forum

Hi, I'm currently trying to implement a partition function that splits the share of heat flux between two meshes in contact via the mortar contact BC condition. The function should split the flux based on the density, conductivity and specific heat properties of the two bodies, and it is possible fo...

Hi, By dropping the number of nodal points on meshing face curve from 50 to 30, it now solves pretty nicely. In both cases using second order elements from gmsh. I'll run with this mesh configuration for now. But I am still interested in any help to make it more robust to meshing as my plan is to pr...

You already have B defined on line 22

In MATC? You can refer to other defined MATC functions in your MATC. $ function B(y) { _B = 3*y+1 } $ function h(y) { _h = 4*y+6 } $ function V(x,y) { _V = cosh(B(y)*h(y)) + sinh(x) } .... Boundary Condition .... var = Variable Coodinate 1, Coordinate 2 Real MATC " V( tx(0), tx(1) ) " .......

Hi Paul, Take this with a pinch of salt, as I don't really know what I am doing, I've just been digging around to try and solve my own contact BC problems, but: It looks like in the SolverUtils.F90 file, a function called DetermineContact is called if 'apply contact bcs = true', and this function wi...

Hi, So running with the idea above, I wrote a short solver, copied below, to be called after each timestep that checks if a BC has a Plane Projector Normal, and if so removes it. The idea being that then the subroutine PlaneInterfaceMeshes in MeshUtils.F90 will rerun the part that generates the plan...

Hi, So having played with the files some more, in both steady state analysis and in smaller subsets of angles the gears rotate through. If I change the number of time steps, the model always starts to go wrong around the same percentage though the simulation. For the mesh attached previously, the an...

Hi, I have updated the boundaries on the pinion and removed RigidMeshMapper. However, the solution still becomes nonsensical from about 80% of the way though the timesteps, please see attached video and image of peak stress compared with calculated Hertzian pressure. The gentler transition in stress...

Hi Pete, Sorry, I removed the results files to get the file size small enough to post. An image of the setup is shown below. case.jpg The keyword "Linear System Residual Mode = True" seems to have made a big difference, it almost works perfectly even with "displace mesh = true". ...

Hi, I've been attempting to model the contact between a set of gear teeth. So far, using the mortar method in slide contact mode I've been able to produce static cases that are in good agreement with the Hertzian analytical solution. However, when I add the "rigidmeshmapper" solver to prod...