alias.h

Finite element program for mechanical problem. It can solve various problem in solid problem

#ifndef ALIAS_H
#define ALIAS_H
#include "kwdset.h"

// aliases for type of solved problem
enum problemtype {linear_statics=1, eigen_dynamics=2, forced_dynamics=3,
                  linear_stability=5,mat_nonlinear_statics=10,geom_nonlinear_statics=11,
		  mech_timedependent_prob=15,growing_mech_structure=17,
		  earth_pressure=20,var_stiff_method=25,layered_linear_statics=30,
		  lin_floating_subdomain=35,nonlin_floating_subdomain=36,
                  nonlinear_dynamics=40};
const enumstr problemtypestr[14] = {{"linear_statics",1}, {"eigen_dynamics",2}, {"forced_dynamics",3},
                                    {"linear_stability",5}, {"mat_nonlinear_statics",10}, {"geom_nonlinear_statics",11},
                                    {"mech_timedependent_prob",15}, {"growing_mech_structure",17}, {"earth_pressure",20},
                                    {"var_stiff_method",25}, {"layered_linear_statics",30}, {"lin_floating_subdomain",35},
                                    {"nonlin_floating_subdomain",36}, {"nonlinear_dynamics",40}};
const kwdset problemtype_kwdset(sizeof(problemtypestr)/sizeof(*problemtypestr), problemtypestr);



// aliases for forced dynamics solvers
enum forcedsolver {newmark=1,findiff=2,modal_analysis=51};
const enumstr forcedsolverstr[3] = {{"newmark",1}, {"findiff",2}, {"modal_analysis",51}};
const kwdset forcedsolver_kwdset(sizeof(forcedsolverstr)/sizeof(*forcedsolverstr), forcedsolverstr);



//aliases for nonlinear statics solvers
enum nonlinsolvertype {arcl=1,newton=2,arclrv1=3,newtonrv1=4,newtonrestart=5,displctrl=6,displctrlrv=7,arclrv=8,adaptram=10};
const enumstr nonlinsolvertypestr[9] = {{"arcl",1}, {"newton",2}, {"arclrv1",3}, {"newtonrv1",4}, {"newtonrestart",5}, 
                                        {"displctrl",6}, {"displctrlrv",7}, {"arclrv",8}, {"adaptram",10}};
const kwdset nonlinsolvertype_kwdset(sizeof(nonlinsolvertypestr)/sizeof(*nonlinsolvertypestr), nonlinsolvertypestr);

//aliases for nonlinear time-dependent solvers
enum nonlintimesolvertype {fullnewton=1,modnewton=2};


// aliases for earth pressure solvers
enum epsolvertype {gep_sol=1, gepvarsup_sol=2, epplast_sol=3};
const enumstr epsolvertypestr[3] = {{"gep_sol",1}, {"gepvarsup_sol",2}, {"epplast_sol",3}};
const kwdset epsolvertypes_kwdset(sizeof(epsolvertypestr)/sizeof(*epsolvertypestr), epsolvertypestr);



// aliases for type of displacement norm computation
enum displacementnorm {alldofs=1,seldofs=2,seldofscoord=3,selecnodes=6,nodesdistincr=8};
const enumstr displacementnormstr[5] = {{"alldofs",1}, {"seldofs",2}, {"seldofscoord",3}, {"selecnodes",6}, {"nodesdistincr",8}};
const kwdset displacementnorm_kwdset(sizeof(displacementnormstr)/sizeof(*displacementnormstr), displacementnormstr);



// aliases for types fo stress return algorithms
enum stressretalgtype {nostressretalg=0,cp=1,gsra=10};
const enumstr stressretalgtypestr[3] = {{"nostressretalg",0}, {"cp",1}, {"gsra",10}};
const kwdset stressretalgtype_kwdset(sizeof(stressretalgtypestr)/sizeof(*stressretalgtypestr), stressretalgtypestr);


// aliases for types of stiffness matrix computation
enum stiffmatrix {initial_stiff=1,tangent_stiff=2,secant_stiff=3};
const enumstr stiffmatrixstr[3] = {{"initial_stiff",1}, {"tangent_stiff",2}, {"secant_stiff",3}};
const kwdset stiffmatrix_kwdset(sizeof(stiffmatrixstr)/sizeof(*stiffmatrixstr), stiffmatrixstr);



// aliases for damping types
enum damping {nodamp=0,massdamp=1,stiffdamp=2,rayleighdamp=3};
const enumstr dampingstr[4] = {{"nodamp",0}, {"massdamp",1}, {"stiffdamp",2}, {"rayleighdamp",3}};
const kwdset damping_kwdset(sizeof(dampingstr)/sizeof(*dampingstr), dampingstr);



// aliases used for properties in preprocessor
enum proptype {ndofs=1, boundarycond=2, loadnodes=3, crosssec=4, localsystem=5,
               eltype=6, matel=7, loadelems=8, dloadnodes=9, dloadelems=10,
               initcond=11, comcodnum=12, loadedge=13, sscomp=14, nodetemp=15, 
               timefun=16,surfaceload=17, elsurfaceload=18, nodtimefunc=19, 
               eltimefunc=20};



// aliases for element types
enum elemtype {bar2d=1,beam2d=2,bar3d=3,beam3d=4,beamg3d=5,barq2d=6,barq3d=7,subsoilbeam=8,beam2dsp=9,
               spring_1 = 10, spring_2 = 11, spring_3 = 12,
               spring_4 = 13, spring_5 = 14, spring_6 = 15,
               planeelementlt=20,planeelementqt=21,planeelementrotlt=22,
               planeelementlq=23,planeelementqq=24,planeelementrotlq=25,
               planeelementsubqt=30,planequadcontact=35,
               cctel=41,dktel=42,dstel=43,q4plateel=45,subsoilplatetr=50,subsoilplateq=51,
               axisymmlt=60,axisymmlq=63,axisymmqq=64,
               shelltrelem=80,shellqelem=81,
               lineartet=100,linearhex=102,quadrhex=103,linearwed=105,quadrwed=106,
               particleelem=200};
const enumstr elemtypestr[40] = {{"bar2d",1}, {"beam2d",2}, {"bar3d",3}, {"beam3d",4}, {"beamg3d",5}, 
                                 {"barq2d",6},{"barq3d",7}, {"subsoilbeam",8}, {"beam2dsp",9}, {"spring_1",10}, 
                                 {"spring_2",11}, {"spring_3",12}, {"spring_4",13}, {"spring_5",14}, {"spring_6",15},
                                 {"planeelementlt",20}, {"planeelementqt",21}, {"planeelementrotlt",22},
                                 {"planeelementlq",23}, {"planeelementqq",24}, {"planeelementrotlq",25},
                                 {"planeelementsubqt",30}, {"planequadcontact",35}, {"cctel",41},
                                 {"dktel",42}, {"dstel",43}, {"q4plateel",45},
                                 {"subsoilplatetr",50}, {"subsoilplateq",51}, {"axisymmlt",60},
                                 {"axisymmlq",63}, {"axisymmqq",64}, {"shelltrelem",80}, 
                                 {"shellqelem",81}, {"lineartet",100}, {"linearhex",102},
                                 {"quadrhex",103}, {"linearwed",105}, {"quadrwed",106},
                                 {"particleelem",200}};
const kwdset elemtype_kwdset(sizeof(elemtypestr)/sizeof(*elemtypestr), elemtypestr);



// aliases for strain/stress state
enum strastrestate {bar=1,plbeam=2,spacebeam=5,
                    planestress=10,planestrain=11,plate=15,
                    axisymm=20,shell=25,
		    spacestress=30};
const enumstr strastrestatestr[9] = {{"bar",1}, {"plbeam",2}, {"spacebeam",5},
                                     {"planestress",10}, {"planestrain",11}, {"plate",15},
                                     {"axisymm",20}, {"shell",25}, {"spacestress",30}};
const kwdset strastrestate_kwdset(sizeof(strastrestatestr)/sizeof(*strastrestatestr), strastrestatestr);



// aliases for  material model classes
enum materialmodel {local=0,nonlocal=1};
const enumstr materialmodelstr[2] = {{"local",0}, {"nonlocal",1}};
const kwdset materialmodel_kwdset(sizeof(materialmodelstr)/sizeof(*materialmodelstr), materialmodelstr);



// aliases for material models
enum mattype {elisomat=1,elgmat3d=2,
              simplas1d=10,jflow=11,jflow2=12,mohrcoulparab=22,mohrcoul=24,
              boermaterial=25,druckerprager=26,modcamclaymat=30,
	      shefpl=40,chenplast=42,
              microplaneM4=50,microsimp=51,microfibro=52,
              simvisplas=70,simvisc=71,lemaitr=75,
              scaldamage=100,scaldamagecc=101,glasgowdamage=102,scaldamage_time=103, anisodamage=104,
	      contmat=150,
              aci=200,cebfip=202,
              nonlocplastmat=310,nonlocdamgmat=320,nonlocalmicroM4=340,
              graphm=400,geoelast=401,graphmtime=402,
              creepbaz=500,creepb3=502,creepdpl=503,creeprs=504,winklerpasternak=550,consolidation=600,
	      glasgowmechmat=700,
	      rcmatmodnorm=800,
              therisodilat=900, therisodilattime=910,
              damage_plasticity=1000, viscoplasticity=1010, elasttimemat=1020, 
              creep_damage=1030, time_switchmat=1040,
              lenjonespot=2000};
const enumstr mattypestr[48] = {{"elisomat",1}, {"elgmat3d",2}, {"simplas1d",10}, 
                                {"jflow",11}, {"jflow2",12}, {"mohrcoulparab",22}, 
                                {"mohrcoul",24}, {"boermaterial",25}, {"druckerprager",26}, 
                                {"modcamclaymat",30}, {"shefpl",40}, {"chenplast",42}, 
                                {"microplaneM4",50}, {"microsimp",51}, {"microfibro",52}, 
                                {"simvisplas",70}, {"simvisc",71}, {"lemaitr",75}, 
                                {"scaldamage",100}, {"scaldamagecc",101}, {"glasgowdamage",102}, 
                                {"scaldamage_time",103}, {"anisodamage",104}, {"contmat",150}, 
                                {"aci",200}, {"cebfip",202}, {"nonlocplastmat",310}, 
                                {"nonlocdamgmat",320}, {"nonlocalmicroM4",340}, {"graphm",400}, 
                                {"geoelast",401}, {"graphmtime",402}, {"creepbaz",500}, 
                                {"creepb3",502}, {"creepdpl",503}, {"creeprs",504}, 
                                {"winklerpasternak",550}, {"consolidation",600}, {"glasgowmechmat",700}, 
                                {"rcmatmodnorm",800}, {"therisodilat",900}, {"therisodilattime",910}, 
                                {"damage_plasticity",1000}, {"viscoplasticity",1010}, {"elasttimemat",1020},  
                                {"creep_damage",1030}, {"time_switchmat",1040}, {"lenjonespot",2000}};