1 | /*! \file Penta.h
|
---|
2 | * \brief: header file for penta object
|
---|
3 | */
|
---|
4 |
|
---|
5 | #ifndef _PENTA_H
|
---|
6 | #define _PENTA_H
|
---|
7 |
|
---|
8 | class Object;
|
---|
9 | class Node;
|
---|
10 | class Hook;
|
---|
11 | class DataSet;
|
---|
12 | struct IoModel;
|
---|
13 |
|
---|
14 |
|
---|
15 | #include "./Object.h"
|
---|
16 | #include "./Element.h"
|
---|
17 | #include "./Matpar.h"
|
---|
18 | #include "./Numpar.h"
|
---|
19 | #include "./Matice.h"
|
---|
20 | #include "./Tria.h"
|
---|
21 | #include "./Hook.h"
|
---|
22 | #include "../ModelProcessorx/IoModel.h"
|
---|
23 | #include "./Node.h"
|
---|
24 |
|
---|
25 | class Penta: public Element{
|
---|
26 |
|
---|
27 | private:
|
---|
28 |
|
---|
29 | int id;
|
---|
30 | Hook hnodes; //hook to 6 nodes
|
---|
31 | Hook hmatice; //hook to 1 matice
|
---|
32 | Hook hmatpar; //hook to 1 matpar
|
---|
33 | Hook hnumpar; //hook to 1 numpar
|
---|
34 |
|
---|
35 | Inputs* inputs;
|
---|
36 |
|
---|
37 | public:
|
---|
38 |
|
---|
39 | /*FUNCTION constructors, destructors {{{1*/
|
---|
40 | Penta();
|
---|
41 | Penta(int penta_id,int* penta_node_ids, int penta_matice_id, int penta_matpar_id, int penta_numpar_id);
|
---|
42 | Penta(int penta_id,Hook* penta_hnodes, Hook* penta_hmatice, Hook* penta_hmatpar, Hook* penta_hnumpar, Inputs* inputs);
|
---|
43 | Penta(int i, IoModel* iomodel);
|
---|
44 | ~Penta();
|
---|
45 | /*}}}*/
|
---|
46 | /*FUNCTION object management {{{1*/
|
---|
47 | void Configure(void* loads,void* nodes,void* materials,void* parameters);
|
---|
48 | Object* copy();
|
---|
49 | void DeepEcho();
|
---|
50 | void Demarshall(char** pmarshalled_dataset);
|
---|
51 | void Echo();
|
---|
52 | int Enum();
|
---|
53 | int GetId();
|
---|
54 | char* GetName();
|
---|
55 | void Marshall(char** pmarshalled_dataset);
|
---|
56 | int MarshallSize();
|
---|
57 | int MyRank();
|
---|
58 | void* SpawnTria(int g0, int g1, int g2);
|
---|
59 | void UpdateFromDakota(void* inputs);
|
---|
60 | void UpdateInputs(double* solution, int analysis_type, int sub_analysis_type);
|
---|
61 | void SetClone(int* minranks);
|
---|
62 |
|
---|
63 | /*}}}*/
|
---|
64 | /*FUNCTION element numerical routines {{{1*/
|
---|
65 | void CreateKMatrix(Mat Kgg,int analysis_type,int sub_analysis_type);
|
---|
66 | void CreateKMatrixDiagnosticHoriz( Mat Kgg, int analysis_type,int sub_analysis_type);
|
---|
67 | void CreateKMatrixDiagnosticVert( Mat Kgg, int analysis_type,int sub_analysis_type);
|
---|
68 | void CreatePVector(Vec pg, int analysis_type,int sub_analysis_type);
|
---|
69 | void GetDofList(int* doflist,int* pnumberofdofs);
|
---|
70 | void GetDofList1(int* doflist);
|
---|
71 | void* GetMatPar();
|
---|
72 | bool GetShelf();
|
---|
73 | void GetNodes(void** nodes);
|
---|
74 | bool GetOnBed();
|
---|
75 | void Du(Vec du_g,int analysis_type,int sub_analysis_type);
|
---|
76 | void Gradj(Vec grad_g,int analysis_type,int sub_analysis_type,char* control_type);
|
---|
77 | void GradjDrag(Vec grad_g,int analysis_type,int sub_analysis_type);
|
---|
78 | void GradjB(Vec grad_g,int analysis_type,int sub_analysis_type);
|
---|
79 | double Misfit(int analysis_type,int sub_analysis_type);
|
---|
80 | double SurfaceArea(int analysis_type,int sub_analysis_type);
|
---|
81 | double CostFunction(int analysis_type,int sub_analysis_type);
|
---|
82 |
|
---|
83 | void GetThicknessList(double* thickness_list);
|
---|
84 | void GetBedList(double* bed_list);
|
---|
85 |
|
---|
86 | void GetStrainRate(double* epsilon, double* velocity, double* xyz_list, double* gauss_coord);
|
---|
87 | void GetB(double* pB, double* xyz_list, double* gauss_coord);
|
---|
88 | void GetBPrime(double* B, double* xyz_list, double* gauss_coord);
|
---|
89 | void GetB_vert(double* B, double* xyz_list, double* gauss_coord);
|
---|
90 | void GetBPrime_vert(double* B, double* xyz_list, double* gauss_coord);
|
---|
91 | void GetJacobianDeterminant(double* Jdet, double* xyz_list,double* gauss_coord);
|
---|
92 | void GetNodalFunctionsDerivatives(double* dh1dh6,double* xyz_list, double* gauss_coord);
|
---|
93 | void GetJacobian(double* J, double* xyz_list,double* gauss_coord);
|
---|
94 | void GetNodalFunctionsDerivativesReference(double* dl1dl6,double* gauss_coord);
|
---|
95 | void GetJacobianInvert(double* Jinv, double* xyz_list,double* gauss_coord);
|
---|
96 | void CreatePVectorDiagnosticHoriz( Vec pg, int analysis_type,int sub_analysis_type);
|
---|
97 | void CreatePVectorDiagnosticVert( Vec pg, int analysis_type,int sub_analysis_type);
|
---|
98 | void GetParameterValue(double* pvalue, double* v_list,double* gauss_coord);
|
---|
99 | void GetParameterDerivativeValue(double* p, double* p_list,double* xyz_list, double* gauss_coord);
|
---|
100 | void GetNodalFunctions(double* l1l6, double* gauss_coord);
|
---|
101 | void FieldExtrude(Vec field,double* field_serial,char* field_name, int iscollapsed);
|
---|
102 | void ComputeBasalStress(Vec sigma_b,int analysis_type,int sub_analysis_type);
|
---|
103 | void ComputePressure(Vec p_g,int analysis_type,int sub_analysis_type);
|
---|
104 | void ComputeStrainRate(Vec eps,int analysis_type,int sub_analysis_type);
|
---|
105 | void CreateKMatrixSlopeCompute(Mat Kgg,int analysis_type,int sub_analysis_type);
|
---|
106 | void CreatePVectorSlopeCompute( Vec pg, int analysis_type,int sub_analysis_type);
|
---|
107 | void CreateKMatrixPrognostic(Mat Kgg,int analysis_type,int sub_analysis_type);
|
---|
108 | void CreatePVectorPrognostic( Vec pg, int analysis_type,int sub_analysis_type);
|
---|
109 |
|
---|
110 | void CreateKMatrixDiagnosticStokes( Mat Kgg, int analysis_type,int sub_analysis_type);
|
---|
111 | void CreatePVectorDiagnosticStokes( Vec pg, int analysis_type,int sub_analysis_type);
|
---|
112 | void ReduceMatrixStokes(double* Ke_reduced, double* Ke_temp);
|
---|
113 | void GetMatrixInvert(double* Ke_invert, double* Ke);
|
---|
114 | void SurfaceNormal(double* surface_normal, double xyz_list[3][3]);
|
---|
115 | void GetStrainRateStokes(double* epsilon, double* velocity, double* xyz_list, double* gauss_coord);
|
---|
116 | void GetBStokes(double* B, double* xyz_list, double* gauss_coord);
|
---|
117 | void GetBprimeStokes(double* B_prime, double* xyz_list, double* gauss_coord);
|
---|
118 | void GetLStokes(double* LStokes, double* gauss_coord_tria);
|
---|
119 | void GetLprimeStokes(double* LprimeStokes, double* xyz_list, double* gauss_coord_tria, double* gauss_coord);
|
---|
120 | void GetNodalFunctionsDerivativesStokes(double* dh1dh7,double* xyz_list, double* gauss_coord);
|
---|
121 | void GetNodalFunctionsDerivativesReferenceStokes(double* dl1dl7,double* gauss_coord);
|
---|
122 | void ReduceVectorStokes(double* Pe_reduced, double* Ke_temp, double* Pe_temp);
|
---|
123 | void GetNodalFunctionsStokes(double* l1l7, double* gauss_coord);
|
---|
124 | void CreateKMatrixThermal(Mat Kgg,int analysis_type,int sub_analysis_type);
|
---|
125 | void GetB_conduct(double* B_conduct, double* xyz_list, double* gauss_coord);
|
---|
126 | void GetB_advec(double* B_advec, double* xyz_list, double* gauss_coord);
|
---|
127 | void GetBprime_advec(double* Bprime_advec, double* xyz_list, double* gauss_coord);
|
---|
128 | void GetB_artdiff(double* B_artdiff, double* xyz_list, double* gauss_coord);
|
---|
129 | void CreateKMatrixMelting(Mat Kgg,int analysis_type,int sub_analysis_type);
|
---|
130 | void CreatePVectorThermal( Vec pg, int analysis_type,int sub_analysis_type);
|
---|
131 | void CreatePVectorMelting( Vec pg, int analysis_type,int sub_analysis_type);
|
---|
132 | void GetPhi(double* phi, double* epsilon, double viscosity);
|
---|
133 | double MassFlux(double* segment,double* ug);
|
---|
134 | /*}}}*/
|
---|
135 |
|
---|
136 |
|
---|
137 | };
|
---|
138 | #endif /* _PENTA_H */
|
---|