File numerics/src/FrictionContact/gfc3d_Solvers.h

Go to the source code of this file

Subroutines for the resolution of contact problems with friction (3-dimensional case).

Typedefs

typedef void (*ComputeErrorGlobalPtr)(GlobalFrictionContactProblem *, double *, double *, double *, double, SolverOptions *, double, double, double *)

pointer to function used to update velocity and compute error

typedef void (*FreeSolverGlobalPtr)(GlobalFrictionContactProblem *)
typedef void (*PostSolverGlobalPtr)(int, double *)
typedef void (*SolverGlobalPtr)(int, int, double *, int *, double *)

Functions

void gfc3d_ACLMFixedPoint(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Solver based on the fixed-point iteration proposed by Cadoux for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • globalVelocity: global vector (m), in-out parameters

  • info: return 0 if the solution is found

  • options: the solver options : iparam[0] : Maximum iteration number iparam[4] ; local strategy dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_aclmfp_set_default(SolverOptions *options)
void gfc3d_ADMM(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
void gfc3d_ADMM_free(GlobalFrictionContactProblem *problem, SolverOptions *options)
void gfc3d_ADMM_init(GlobalFrictionContactProblem *problem, SolverOptions *options)
void gfc3d_admm_set_default(SolverOptions *options)
void gfc3d_admm_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
void gfc3d_AVI_gams_path(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, int *info, SolverOptions *options)

solver using PATH (via GAMS) for friction-contact 3D problem based on an AVI reformulation

Parameters
  • problem: the friction-contact 3D problem to solve

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • info: return 0 if the solution is found

  • options: the solver options

void gfc3d_AVI_gams_pathvi(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, int *info, SolverOptions *options)

solver using PATHVI (via GAMS) for friction-contact 3D problem based on an AVI reformulation

Parameters
  • problem: the friction-contact 3D problem to solve

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • info: return 0 if the solution is found

  • options: the solver options

int gfc3d_checkTrivialCaseGlobal(int dim, double *q, double *velocity, double *reaction, double *globalVelocity, SolverOptions *options)

Check for trivial solution in the friction-contact 3D problem.

Return

int =0 if a trivial solution has been found, else = -1

Parameters
  • dim: of the problem

  • q: global vector (n)

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • globalVelocity: the velocity in global coordinates

  • options: the pointer to the array of options to set

void gfc3d_DeSaxceFixedPoint_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Fixed Point iteration on De Saxe formulation solver with reformulation for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • globalVelocity: global vector (m), in-out parameters

  • info: return 0 if the solution is found

  • options: the solver options : iparam[0] : Maximum iteration number dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_nonsmooth_Newton_AlartCurnier(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
void gfc3d_nonsmooth_Newton_AlartCurnier_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
void gfc3d_nsgs(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Non-Smooth Gauss Seidel solver for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • globalVelocity: global vector (m), in-out parameters

  • info: return 0 if the solution is found

  • options: the solver options : iparam[0] : Maximum iteration number iparam[4] ; local strategy dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_nsgs_velocity_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Non-Smooth Gauss Seidel solver for friction-contact 3D problem with iteration on velocities.

Parameters
  • problem: the friction-contact 3D problem to solve

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • globalVelocity: global vector (m), in-out parameters

  • info: return 0 if the solution is found

  • options: the solver options : iparam[0] : Maximum iteration number dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_nsgs_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Non-Smooth Gauss Seidel solver with reformulation for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • globalVelocity: global vector (m), in-out parameters

  • info: return 0 if the solution is found

  • options: the solver options : iparam[0] : Maximum iteration number iparam[4] : localsolver choice 0: projection on Cone, 1: Newton/AlartCurnier, 2: projection on Cone with local iteration, 2: projection on Disk with diagonalization, dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_nsn_ac_set_default(SolverOptions *options)
void gfc3d_proximal_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Proximal point solver with reformulation for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • globalVelocity: global vector (m), in-out parameters

  • info: return 0 if the solution is found

  • options: the solver options : iparam[0] : Maximum iteration number iparam[4] : localsolver choice 0: projection on Cone, 1: Newton/AlartCurnier, 2: projection on Cone with local iteration, 2: projection on Disk with diagonalization, dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_set_internalsolver_tolerance(GlobalFrictionContactProblem *problem, SolverOptions *options, SolverOptions *internalsolver_options, double error)
void gfc3d_TrescaFixedPoint_wr(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)

Fied Point iteration on Tresca Friction Cylinder with reformulation for friction-contact 3D problem.

Parameters
  • problem: the friction-contact 3D problem to solve

  • velocity: global vector (n), in-out parameter

  • reaction: global vector (n), in-out parameters

  • globalVelocity: global vector (m), in-out parameters

  • info: return 0 if the solution is found

  • options: the solver options : iparam[0] : Maximum iteration number dparam[0] : tolerance dparam[2] : localtolerance dparam[1] : (out) error

void gfc3d_VI_ExtraGradient(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)
void gfc3d_VI_FixedPointProjection(GlobalFrictionContactProblem *problem, double *reaction, double *velocity, double *globalVelocity, int *info, SolverOptions *options)