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 <[email protected]> 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 > [email protected] > https://mail.gna.org/listinfo/getfem-users > >
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
