ARCHIVE/pdd_forming_arch/include/pdh_plast_flow_supg.h

Go to the documentation of this file.

/************************************************************************
File pdh_plast_flow.h - problem module's types and functions
 BY ASSUMPTION NOT CALLED BY OTHER MODULES' FUNCTIONS
 (functions called by other modules are in pdh_plast_flow_weakform.h)

Contains problem module defines (see below)


Contains declarations of routines:
 BY ASSUMPTION NOT CALLED BY OTHER MODULES' FUNCTIONS
 (functions called by other modules are in pdh_plast_flow_problem.h)
  pdr_plast_flow_time_integration - time integration driver (in time_integration)
  pdr_plast_flow_ZZ_error - to compute estimated norm of error based on 
                 recovered first derivatives - the notorious ZZ error estimate
  pdr_plast_flow_err_indi_ZZ - to return error indicator for an element,
        based on ZZ first derivative recovery 
  pdr_plast_flow_adapt - to enforce adaptation strategy for ns problem (in adapt)
  pdr_plast_flow_dump_data - dump data to files (in input_output)
  pdr_plast_flow_write_paraview - to write graphics data to file (in input_output)

utilities:
  pdr_plast_flow_post_process
  pdr_plast_flow_write_profile
  pdr_plast_flow_initial_condition

------------------------------
History:
    initial version - Krzysztof Banas (pobanas@cyf-kr.edu.pl)
*************************************************************************/

#ifndef PDH_PLAST_FLOW
#define PDH_PLAST_FLOW

#include <stdio.h>

/* problem dependent interface with the PDEs  */
#include "pdh_intf.h"

/* types and functions related to problem structures */
#include "pdh_plast_flow_supg_problem.h" 
// bc and material header files are included in problem header files
#include "aph_intf.h"


#ifdef __cplusplus
extern "C" 
{
#endif

/**************************************/
/* DEFINES                            */
/**************************************/
/* Rules:
/* - always uppercase */
/* - name starts with PDC_ */

#define PDC_PLAST_FLOW_ID  1

/**************************************/
/* TYPES                              */
/**************************************/
/* Rules:
/* - type name starts always witn pdt_ */


/**************************************/
/* GLOBAL VARIABLES                   */
/**************************************/
/* Rules:
/* - name always begins with pdv_ */
/* - constants always uppercase and start with PDC_ */


// ID of the current problem
extern int pdv_plast_flow_current_problem_id;   /* ID of the current problem */
// problem structure for plast_flow module
extern pdt_plast_flow_problem pdv_plast_flow_problem;

#ifdef PARALLEL
  extern int pdv_exchange_table_index; // for a single problem
#endif

/**************************************/
/* INTERNAL PROCEDURES                */
/**************************************/
/* Rules:
/* - name always begins with pdr_ */
/* - argument names start uppercase */


extern void pdr_plast_flow_time_integration(
  char* Work_dir,
  FILE *Interactive_input, 
  FILE *Interactive_output);


extern int pdr_plast_flow_adapt(
  char* Work_dir,
  FILE *Interactive_input, 
  FILE *Interactive_output);

double pdr_plast_flow_ZZ_error(     
       /* returns  - Zienkiewicz-Zhu error for the whole mesh */
  char* Work_dir,
  FILE *Interactive_input, 
  FILE *Interactive_output
);

extern double pdr_plast_flow_err_indi_ZZ(       
                        /* returns error indicator for an element */
  int Problem_id,   /* in: data structure to be used  */
  int El    /* in: element number */
    );

double pdr_plast_flow_sol_mesh(/* returns total volume of the mesh */
        int Problem_id,                 /* in: data structure to be used  */
        double *sol_previous,           /* in: averaged previous solution (21) */
        double *sol_averaged            /* out: averaged current solution (21) */
);

/*-------------------------------------------------------------
pdr_plast_flow_sol_el - to return volume of element and
    averaged: strain intensity (1), flow stress (1), strain rate
    intensity (1), strain tensor (9), strain rate tensor (9)
    for an element (=21)
-------------------------------------------------------------*/
extern double pdr_plast_flow_sol_el(/* returns volume of element */
        int Problem_id,                 /* in: data structure to be used  */
        int El,                         /* in: element number */
        double *sol_current             /* out: averaged solution (21) */
);

extern int pdr_plast_flow_dump_data(
  char* Work_dir,
  FILE *Interactive_input,
  FILE *Interactive_output
);

/*---------------------------------------------------------
pdr_write_strain_tensor - to dump tensor, derived field
---------------------------------------------------------*/
int pdr_write_strain_tensor(
  int Field_id,     // derived field ID
  int* Dofs_write,  // dofs to write: dofs_write[0] - number of dofs
                    //                dofs_write[i] - IDs of dofs to write
  double dt
);

/*---------------------------------------------------------
pdr_write_stress_tensor - to dump tensor, derived field
---------------------------------------------------------*/
int pdr_write_stress_tensor(
  int Field_id,     // derived field ID
  int* Dofs_write,  // dofs to write: dofs_write[0] - number of dofs
                    //                dofs_write[i] - IDs of dofs to write
  double dt
);

extern int pdr_plast_flow_write_paraview(
  char* Work_dir,
  FILE *Interactive_input, 
  FILE *Interactive_output
  );

/*---------------------------------------------------------
---------------------------------------------------------*/
void pdr_write_paraview_mesh_nodes_displacement(
  int Mesh_id, int Mesh_id_initial, FILE  *File);


  // IN: MAIN/PDS_PLAST_FLOW_UTIL.C

double pdr_plast_flow_post_process(
  char* Work_dir,
  FILE *Interactive_input,
  FILE *Interactive_output
);

int pdr_plast_flow_write_profile(
  char* Work_dir,
  FILE *Interactive_input,
  FILE *Interactive_output
);

double pdr_plast_flow_initial_condition(
  int Field_id, // field_id - each problem should know its field id
  double *Coor,   // point coordinates
  int Sol_comp_id // solution component
);

double pdr_strain_initial_condition(
  int Field_id, // field_id - each problem should know its field id
  double *Coor,   // point coordinates
  int Sol_comp_id // solution component
);

#ifdef __cplusplus
}

#endif

#endif

---