Re: [Getfem-users] Shrink fit problem

[Konstantinos Poulios]

Hi Miikka,

some questions/comments/recommendations:

- Which version of getfem++ do you use, which operating system?

- Why don't you use "add_integral_contact_between_nonmatching_meshes_brick"? It is normally superior to nodal contact.

- Why don't you use a direct solver like superlu or even better mumps? How many degrees of freedom do you have?

- The easiest way of accounting for the interference fit is by modifying the meshes so that they incorporate the fit. You can modify the mesh in your mesher software or within getfem, e.g. by using the functions pts() and set_pts() of the mesh object in python.

Best regards

Kostas

On Tue, Aug 19, 2014 at 2:59 PM, Miikka Väntänen <address@hidden> wrote:

Hi all,

 

I’m relatively new to GetFem, but I’ll try to do my best in explaining what I’m struggling with…

 

I’m trying to solve a linear elastic problem that consists of two parts, a shaft and a bush. Between the two parts I use a contact “add_nodal_contact_between_nonmatching_meshes_brick” (in python api). I’m using a linear gmres solver with either ILU or ILUT preconditioning (probably ilut works better with saddle point problems?). There’s also some loading implemented to the shaft, but I reckon that’s irrelevant at this point.

 

The real problem I’m currently facing is with the shrink fit. The actual shaft-bush contact should have a shrink fit with 5 mm radial difference between the two parts. Is there any way in getfem to implement the shrink fit? Should I, for example, use temperature change to shrink/expand the parts in an individual step before calculating the actual static problem? Or should I implement a radial displacement to either of the faces in contact? Is there any way to do any of this?

Br.

Miikka V.

_______________________________________________
Getfem-users mailing list
address@hidden
https://mail.gna.org/listinfo/getfem-users