|
From: | Miikka Väntänen |
Subject: | Re: [Getfem-users] Shrink fit problem |
Date: | Wed, 20 Aug 2014 09:44:47 +0300 |
KostasBest regards- 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.- Why don't you use a direct solver like superlu or even better mumps? How many degrees of freedom do you have?- Why don't you use "add_integral_contact_between_nonmatching_meshes_brick"? It is normally superior to nodal contact.- Which version of getfem++ do you use, which operating system?Hi Miikka,some questions/comments/recommendations:
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
[Prev in Thread] | Current Thread | [Next in Thread] |