Search found 10 matches
- 10 Aug 2022, 13:18
- Forum: ElmerSolver
- Topic: Project material properties over mortar boundary
- Replies: 0
- Views: 1013
Project material properties over mortar boundary
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...
- 06 Aug 2021, 19:53
- Forum: ElmerSolver
- Topic: Sliding contact along a curved boundary
- Replies: 6
- Views: 2720
Re: Sliding contact along a curved boundary
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...
- 06 Aug 2021, 18:20
- Forum: ElmerSolver
- Topic: Create a 2D functional
- Replies: 4
- Views: 1373
Re: Create a 2D functional
You already have B defined on line 22
Code: Select all
$ B= 0.5 ! tesla
- 06 Aug 2021, 17:55
- Forum: ElmerSolver
- Topic: Create a 2D functional
- Replies: 4
- Views: 1373
Re: Create a 2D functional
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) ) " .......
- 05 Aug 2021, 13:24
- Forum: ElmerSolver
- Topic: Contact and Periodic BCs
- Replies: 9
- Views: 3248
Re: Contact and Periodic BCs
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...
- 30 Jul 2021, 18:48
- Forum: ElmerSolver
- Topic: Sliding contact along a curved boundary
- Replies: 6
- Views: 2720
Re: Sliding contact along a curved boundary
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...
- 29 Jul 2021, 19:50
- Forum: ElmerSolver
- Topic: Sliding contact along a curved boundary
- Replies: 6
- Views: 2720
Re: Sliding contact along a curved boundary
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...
- 28 Jul 2021, 18:58
- Forum: ElmerSolver
- Topic: Sliding contact along a curved boundary
- Replies: 6
- Views: 2720
Re: Sliding contact along a curved boundary
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...
- 20 Jul 2021, 19:09
- Forum: ElmerSolver
- Topic: Sliding contact along a curved boundary
- Replies: 6
- Views: 2720
Re: Sliding contact along a curved boundary
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". ...
- 19 Jul 2021, 17:51
- Forum: ElmerSolver
- Topic: Sliding contact along a curved boundary
- Replies: 6
- Views: 2720
Sliding contact along a curved boundary
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...