v0.14.0
SmallTransverselyIsotropic.hpp
Go to the documentation of this file.
1 /** \file SmallStrainTranverslyIsotropic.hpp
2  * \ingroup nonlinear_elastic_elem
3  * \brief Implementation of linear transverse isotropic material.
4  */
5 
6 
7 
8 #ifndef __SMALLSTRAINTRANVERSLYISOTROPIC_HPP__
9 #define __SMALLSTRAINTRANVERSLYISOTROPIC_HPP__
10 
11 #ifndef WITH_ADOL_C
12  #error "MoFEM need to be compiled with ADOL-C"
13 #endif
14 
15 /** \brief Hook equation
16  * \ingroup nonlinear_elastic_elem
17  */
18 template<typename TYPE>
19 struct SmallStrainTranverslyIsotropic: public NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI<TYPE> {
20 
21  SmallStrainTranverslyIsotropic(): NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI<TYPE>() {}
22 
23 
24 
25  ublas::matrix<TYPE> sTrain;
26  ublas::vector<TYPE> voightStrain;
27  TYPE tR;
28 
29  MoFEMErrorCode calculateStrain() {
30  MoFEMFunctionBeginHot;
31  sTrain.resize(3,3,false);
32  noalias(sTrain) = this->F;
33  for(int dd = 0;dd<3;dd++) {
34  sTrain(dd,dd) -= 1;
35  }
36  sTrain += trans(sTrain);
37  sTrain *= 0.5;
38  voightStrain.resize(6,false);
39  voightStrain[0] = sTrain(0,0);
40  voightStrain[1] = sTrain(1,1);
41  voightStrain[2] = sTrain(2,2);
42  voightStrain[3] = 2*sTrain(0,1);
43  voightStrain[4] = 2*sTrain(1,2);
44  voightStrain[5] = 2*sTrain(2,0);
45  MoFEMFunctionReturnHot(0);
46  }
47 
48  double nu_p, nu_pz, E_p, E_z, G_zp;
49  ublas::symmetric_matrix<TYPE,ublas::upper> localStiffnessMatrix;
50  MoFEMErrorCode calculateLocalStiffnesMatrix() {
51  MoFEMFunctionBeginHot;
52  double nu_zp=(nu_pz*E_z)/E_p;
53  double delta=((1+nu_p)*(1-nu_p-(2*nu_pz*nu_zp)))/(E_p*E_p*E_z);
54 
55  localStiffnessMatrix.resize(6);
56  localStiffnessMatrix.clear();
57  localStiffnessMatrix(0,0)=localStiffnessMatrix(1,1)=(1-nu_pz*nu_zp)/(E_p*E_z*delta);
58  localStiffnessMatrix(2,2)=(1-nu_p*nu_p)/(E_p*E_p*delta);
59 
60  localStiffnessMatrix(0,1)=localStiffnessMatrix(1,0)=(nu_p+nu_zp*nu_pz)/(E_p*E_z*delta);
61  localStiffnessMatrix(0,2)=localStiffnessMatrix(2,0)=localStiffnessMatrix(1,2)=localStiffnessMatrix(2,1)=(nu_zp+nu_p*nu_zp)/(E_p*E_z*delta);
62 
63  localStiffnessMatrix(3,3)=E_p/(2*(1+nu_p));
64  localStiffnessMatrix(4,4)=localStiffnessMatrix(5,5)=G_zp;
65 
66  MoFEMFunctionReturnHot(0);
67  }
68 
69  ublas::matrix<TYPE> aARotMat;
70  ublas::vector<TYPE> axVector;
71  TYPE axAngle;
72 
73  /**
74  * \brief Function to Calculate the Rotation Matrix at a given axis and angle of rotation
75 
76  * This function computes the rotational matrix for a given axis of rotation
77  * and angle of rotation about that angle <br>
78 
79  *\param axVector A vector representing the axis of rotation
80  *\param axAngle Angle of rotation along the axis (in radians)
81  */
82  MoFEMErrorCode calculateAxisAngleRotationalMatrix() {
83  MoFEMFunctionBeginHot;
84 
85  aARotMat.resize(3,3,false);
86  aARotMat.clear();
87 
88  TYPE norm = sqrt(pow(axVector[0],2) + pow(axVector[1],2) + pow(axVector[2],2));
89 
90  aARotMat(0,0) = 1-((1-cos(axAngle))*(pow(axVector[1],2)+pow(axVector[2],2))/pow(norm,2));
91  aARotMat(1,1) = 1-((1-cos(axAngle))*(pow(axVector[0],2)+pow(axVector[2],2))/pow(norm,2));
92  aARotMat(2,2) = 1-((1-cos(axAngle))*(pow(axVector[0],2)+pow(axVector[1],2))/pow(norm,2));
93 
94  aARotMat(0,1) = ((1-cos(axAngle))*axVector[0]*axVector[1]-norm*axVector[2]*sin(axAngle))/pow(norm,2);
95  aARotMat(1,0) = ((1-cos(axAngle))*axVector[0]*axVector[1]+norm*axVector[2]*sin(axAngle))/pow(norm,2);
96 
97  aARotMat(0,2) = ((1-cos(axAngle))*axVector[0]*axVector[2]+norm*axVector[1]*sin(axAngle))/pow(norm,2);
98  aARotMat(2,0) = ((1-cos(axAngle))*axVector[0]*axVector[2]-norm*axVector[1]*sin(axAngle))/pow(norm,2);
99 
100  aARotMat(1,2) = ((1-cos(axAngle))*axVector[1]*axVector[2]-norm*axVector[0]*sin(axAngle))/pow(norm,2);
101  aARotMat(2,1) = ((1-cos(axAngle))*axVector[1]*axVector[2]+norm*axVector[0]*sin(axAngle))/pow(norm,2);
102 
103  MoFEMFunctionReturnHot(0);
104  }
105 
106  ublas::matrix<TYPE> stressRotMat;
107 
108  /**
109  * \brief Function to Calculate Stress Transformation Matrix
110  * This function computes the stress transformation Matrix at a give axis and angle of rotation <br>
111  * One can also output the axis/angle rotational Matrix
112  */
113  MoFEMErrorCode stressTransformation() {
114  MoFEMFunctionBeginHot;
115 
116  stressRotMat.resize(6,6,false);
117  stressRotMat.clear();
118 
119  stressRotMat(0, 0) = aARotMat(0,0) * aARotMat(0,0);
120  stressRotMat(0, 1) = aARotMat(1,0) * aARotMat(1,0);
121  stressRotMat(0, 2) = aARotMat(2,0) * aARotMat(2,0);
122  stressRotMat(0, 3) = 2.0 * aARotMat(1,0) * aARotMat(0,0);
123  stressRotMat(0, 4) = 2.0 * aARotMat(2,0) * aARotMat(1,0);
124  stressRotMat(0, 5) = 2.0 * aARotMat(0,0) * aARotMat(2,0);
125 
126  stressRotMat(1, 0) = aARotMat(0,1) * aARotMat(0,1);
127  stressRotMat(1, 1) = aARotMat(1,1) * aARotMat(1,1);
128  stressRotMat(1, 2) = aARotMat(2,1) * aARotMat(2,1);
129  stressRotMat(1, 3) = 2.0 * aARotMat(1,1) * aARotMat(0,1);
130  stressRotMat(1, 4) = 2.0 * aARotMat(2,1) * aARotMat(1,1);
131  stressRotMat(1, 5) = 2.0 * aARotMat(0,1) * aARotMat(2,1);
132 
133  stressRotMat(2, 0) = aARotMat(0,2) * aARotMat(0,2);
134  stressRotMat(2, 1) = aARotMat(1,2) * aARotMat(1,2);
135  stressRotMat(2, 2) = aARotMat(2,2) * aARotMat(2,2);
136  stressRotMat(2, 3) = 2.0 * aARotMat(1,2) * aARotMat(0,2);
137  stressRotMat(2, 4) = 2.0 * aARotMat(2,2) * aARotMat(1,2);
138  stressRotMat(2, 5) = 2.0 * aARotMat(0,2) * aARotMat(2,2);
139 
140  stressRotMat(3, 0) = aARotMat(0,1) * aARotMat(0,0);
141  stressRotMat(3, 1) = aARotMat(1,1) * aARotMat(1,0);
142  stressRotMat(3, 2) = aARotMat(2,1) * aARotMat(2,0);
143  stressRotMat(3, 3) = ( aARotMat(1,1) * aARotMat(0,0) + aARotMat(0,1) * aARotMat(1,0) );
144  stressRotMat(3, 4) = ( aARotMat(2,1) * aARotMat(1,0) + aARotMat(1,1) * aARotMat(2,0) );
145  stressRotMat(3, 5) = ( aARotMat(0,1) * aARotMat(2,0) + aARotMat(2,1) * aARotMat(0,0) );
146 
147  stressRotMat(4, 0) = aARotMat(0,2) * aARotMat(0,1);
148  stressRotMat(4, 1) = aARotMat(1,2) * aARotMat(1,1);
149  stressRotMat(4, 2) = aARotMat(2,2) * aARotMat(2,1);
150  stressRotMat(4, 3) = ( aARotMat(1,2) * aARotMat(0,1) + aARotMat(0,2) * aARotMat(1,1) );
151  stressRotMat(4, 4) = ( aARotMat(2,2) * aARotMat(1,1) + aARotMat(1,2) * aARotMat(2,1) );
152  stressRotMat(4, 5) = ( aARotMat(0,2) * aARotMat(2,1) + aARotMat(2,2) * aARotMat(0,1) );
153 
154  stressRotMat(5, 0) = aARotMat(0,0) * aARotMat(0,2);
155  stressRotMat(5, 1) = aARotMat(1,0) * aARotMat(1,2);
156  stressRotMat(5, 2) = aARotMat(2,0) * aARotMat(2,2);
157  stressRotMat(5, 3) = ( aARotMat(1,0) * aARotMat(0,2) + aARotMat(0,0) * aARotMat(1,2) );
158  stressRotMat(5, 4) = ( aARotMat(2,0) * aARotMat(1,2) + aARotMat(1,0) * aARotMat(2,2) );
159  stressRotMat(5, 5) = ( aARotMat(0,0) * aARotMat(2,2) + aARotMat(2,0) * aARotMat(0,2) );
160 
161  MoFEMFunctionReturnHot(0);
162  }
163 
164  ublas::matrix<TYPE> strainRotMat;
165 
166  /**
167  * \brief Function to Calculate Strain Transformation Matrix<br>
168  * This function computes the strain transformation Matrix at a give axis and angle of rotation <br>
169  * One can also output the axis/angle rotational Matrix
170  */
171  MoFEMErrorCode strainTransformation() {
172  MoFEMFunctionBeginHot;
173 
174  strainRotMat.resize(6,6,false);
175  strainRotMat.clear();
176 
177  strainRotMat(0, 0) = aARotMat(0,0) * aARotMat(0,0);
178  strainRotMat(0, 1) = aARotMat(1,0) * aARotMat(1,0);
179  strainRotMat(0, 2) = aARotMat(2,0) * aARotMat(2,0);
180  strainRotMat(0, 3) = aARotMat(1,0) * aARotMat(0,0);
181  strainRotMat(0, 4) = aARotMat(2,0) * aARotMat(1,0);
182  strainRotMat(0, 5) = aARotMat(0,0) * aARotMat(2,0);
183 
184  strainRotMat(1, 0) = aARotMat(0,1) * aARotMat(0,1);
185  strainRotMat(1, 1) = aARotMat(1,1) * aARotMat(1,1);
186  strainRotMat(1, 2) = aARotMat(2,1) * aARotMat(2,1);
187  strainRotMat(1, 3) = aARotMat(1,1) * aARotMat(0,1);
188  strainRotMat(1, 4) = aARotMat(2,1) * aARotMat(1,1);
189  strainRotMat(1, 5) = aARotMat(0,1) * aARotMat(2,1);
190 
191  strainRotMat(2, 0) = aARotMat(0,2) * aARotMat(0,2);
192  strainRotMat(2, 1) = aARotMat(1,2) * aARotMat(1,2);
193  strainRotMat(2, 2) = aARotMat(2,2) * aARotMat(2,2);
194  strainRotMat(2, 3) = aARotMat(1,2) * aARotMat(0,2);
195  strainRotMat(2, 4) = aARotMat(2,2) * aARotMat(1,2);
196  strainRotMat(2, 5) = aARotMat(0,2) * aARotMat(2,2);
197 
198  strainRotMat(3, 0) = 2.0 * aARotMat(0,1) * aARotMat(0,0);
199  strainRotMat(3, 1) = 2.0 * aARotMat(1,1) * aARotMat(1,0);
200  strainRotMat(3, 2) = 2.0 * aARotMat(2,1) * aARotMat(2,0);
201  strainRotMat(3, 3) = ( aARotMat(1,1) * aARotMat(0,0) + aARotMat(0,1) * aARotMat(1,0) );
202  strainRotMat(3, 4) = ( aARotMat(2,1) * aARotMat(1,0) + aARotMat(1,1) * aARotMat(2,0) );
203  strainRotMat(3, 5) = ( aARotMat(0,1) * aARotMat(2,0) + aARotMat(2,1) * aARotMat(0,0) );
204 
205  strainRotMat(4, 0) = 2.0 * aARotMat(0,2) * aARotMat(0,1);
206  strainRotMat(4, 1) = 2.0 * aARotMat(1,2) * aARotMat(1,1);
207  strainRotMat(4, 2) = 2.0 * aARotMat(2,2) * aARotMat(2,1);
208  strainRotMat(4, 3) = ( aARotMat(1,2) * aARotMat(0,1) + aARotMat(0,2) * aARotMat(1,1) );
209  strainRotMat(4, 4) = ( aARotMat(2,2) * aARotMat(1,1) + aARotMat(1,2) * aARotMat(2,1) );
210  strainRotMat(4, 5) = ( aARotMat(0,2) * aARotMat(2,1) + aARotMat(2,2) * aARotMat(0,1) );
211 
212  strainRotMat(5, 0) = 2.0 * aARotMat(0,0) * aARotMat(0,2);
213  strainRotMat(5, 1) = 2.0 * aARotMat(1,0) * aARotMat(1,2);
214  strainRotMat(5, 2) = 2.0 * aARotMat(2,0) * aARotMat(2,2);
215  strainRotMat(5, 3) = ( aARotMat(1,0) * aARotMat(0,2) + aARotMat(0,0) * aARotMat(1,2) );
216  strainRotMat(5, 4) = ( aARotMat(2,0) * aARotMat(1,2) + aARotMat(1,0) * aARotMat(2,2) );
217  strainRotMat(5, 5) = ( aARotMat(0,0) * aARotMat(2,2) + aARotMat(2,0) * aARotMat(0,2) );
218 
219  MoFEMFunctionReturnHot(0);
220  }
221 
222  ublas::matrix<TYPE> dR;
223  ublas::matrix<TYPE> globalStiffnessMatrix;
224 
225  MoFEMErrorCode calculateGlobalStiffnesMatrix() {
226  MoFEMFunctionBeginHot;
227 
228  dR.resize(6,6,false);
229  noalias(dR) = prod(localStiffnessMatrix,strainRotMat);
230  globalStiffnessMatrix.resize(6,6,false);
231  noalias(globalStiffnessMatrix) = prod(stressRotMat,dR);
232 
233  MoFEMFunctionReturnHot(0);
234  }
235 
236  virtual MoFEMErrorCode calculateAngles() {
237  MoFEMFunctionBeginHot;
238  MoFEMFunctionReturnHot(0);
239  }
240 
241  ublas::vector<TYPE> voigtStress;
242 
243  /** \brief Calculate global stress
244 
245  This is small strain approach, i.e. Piola stress
246  is like a Cauchy stress, since configurations are notation
247  distinguished.
248 
249  */
250  virtual MoFEMErrorCode calculateP_PiolaKirchhoffI(
251  const NonlinearElasticElement::BlockData block_data,
252  boost::shared_ptr<const NumeredEntFiniteElement> fe_ptr
253  ) {
254  MoFEMFunctionBeginHot;
255  ierr = calculateAngles(); CHKERRG(ierr);
256  ierr = calculateStrain(); CHKERRG(ierr);
257  ierr = calculateLocalStiffnesMatrix(); CHKERRG(ierr);
258  ierr = calculateAxisAngleRotationalMatrix(); CHKERRG(ierr);
259  ierr = stressTransformation(); CHKERRG(ierr);
260  axAngle = -axAngle;
261  ierr = calculateAxisAngleRotationalMatrix(); CHKERRG(ierr);
262  ierr = strainTransformation(); CHKERRG(ierr);
263  ierr = calculateGlobalStiffnesMatrix(); CHKERRG(ierr);
264 
265  voigtStress.resize(6,false);
266  noalias(voigtStress) = prod(globalStiffnessMatrix,voightStrain);
267  this->P.resize(3,3,false);
268  this->P(0,0) = voigtStress[0];
269  this->P(1,1) = voigtStress[1];
270  this->P(2,2) = voigtStress[2];
271  this->P(0,1) = voigtStress[3];
272  this->P(1,2) = voigtStress[4];
273  this->P(0,2) = voigtStress[5];
274  this->P(1,0) = this->P(0,1);
275  this->P(2,1) = this->P(1,2);
276  this->P(2,0) = this->P(0,2);
277 
278  MoFEMFunctionReturnHot(0);
279  }
280 
281  /** \brief calculate density of strain energy
282  *
283  */
284  virtual MoFEMErrorCode calculateElasticEnergy(
285  const NonlinearElasticElement::BlockData block_data,
286  boost::shared_ptr<const NumeredEntFiniteElement> fe_ptr
287 ) {
288  MoFEMFunctionBeginHot;
289 
290  ierr = calculateAngles(); CHKERRG(ierr);
291  ierr = calculateStrain(); CHKERRG(ierr);
292  ierr = calculateLocalStiffnesMatrix(); CHKERRG(ierr);
293  ierr = calculateAxisAngleRotationalMatrix(); CHKERRG(ierr);
294  ierr = stressTransformation(); CHKERRG(ierr);
295  axAngle = -axAngle;
296  ierr = calculateAxisAngleRotationalMatrix(); CHKERRG(ierr);
297  ierr = strainTransformation(); CHKERRG(ierr);
298  ierr = calculateGlobalStiffnesMatrix(); CHKERRG(ierr);
299 
300  voigtStress.resize(6,false);
301  noalias(voigtStress) = prod(globalStiffnessMatrix,voightStrain);
302  this->eNergy = 0.5*inner_prod(voigtStress,voightStrain);
303  MoFEMFunctionReturnHot(0);
304  }
305 
306  VectorDouble normalizedPhi;
307  VectorDouble axVectorDouble;
308  double axAngleDouble;
309 
310  MoFEMErrorCode calculateFibreAngles() {
311  MoFEMFunctionBeginHot;
312 
313  try {
314 
315  int gg = this->gG; // number of integration point
316  MatrixDouble &phi = (this->commonDataPtr->gradAtGaussPts["POTENTIAL_FIELD"][gg]);
317  normalizedPhi.resize(3,false);
318  double nrm2_phi = sqrt(pow(phi(0,0),2)+pow(phi(0,1),2)+pow(phi(0,2),2));
319  for(int ii = 0;ii<3;ii++) {
320  normalizedPhi[ii] = -phi(0,ii)/nrm2_phi;
321  }
322 
323  axVectorDouble.resize(3,false);
324  const double zVec[3]={ 0.0,0.0,1.0 };
325  axVectorDouble[0] = normalizedPhi[1]*zVec[2]-normalizedPhi[2]*zVec[1];
326  axVectorDouble[1] = normalizedPhi[2]*zVec[0]-normalizedPhi[0]*zVec[2];
327  axVectorDouble[2] = normalizedPhi[0]*zVec[1]-normalizedPhi[1]*zVec[0];
328  double nrm2_ax_vector = norm_2(axVectorDouble);
329  const double eps = 1e-12;
330  if(nrm2_ax_vector<eps) {
331  axVectorDouble[0] = 0;
332  axVectorDouble[1] = 0;
333  axVectorDouble[2] = 1;
334  nrm2_ax_vector = 1;
335  }
336  axAngleDouble = asin(nrm2_ax_vector);
337 
338  } catch (const std::exception& ex) {
339  std::ostringstream ss;
340  ss << "throw in method: " << ex.what() << std::endl;
341  SETERRQ(PETSC_COMM_SELF,1,ss.str().c_str());
342  }
343 
344  MoFEMFunctionReturnHot(0);
345  }
346 
347 };
348 
349 
350 struct SmallStrainTranverslyIsotropicDouble: public SmallStrainTranverslyIsotropic<double> {
351 
352  SmallStrainTranverslyIsotropicDouble(): SmallStrainTranverslyIsotropic<double>() {}
353 
354  virtual MoFEMErrorCode calculateAngles() {
355  MoFEMFunctionBeginHot;
356 
357  try {
358 
359  ierr = calculateFibreAngles(); CHKERRG(ierr);
360  axVector.resize(3,false);
361  axVector[0] = axVectorDouble[0];
362  axVector[1] = axVectorDouble[1];
363  axVector[2] = axVectorDouble[2];
364  axAngle = axAngleDouble;
365 
366  } catch (const std::exception& ex) {
367  std::ostringstream ss;
368  ss << "throw in method: " << ex.what() << std::endl;
369  SETERRQ(PETSC_COMM_SELF,1,ss.str().c_str());
370  }
371  MoFEMFunctionReturnHot(0);
372  }
373 
374  virtual MoFEMErrorCode getDataOnPostProcessor(
375  std::map<std::string,std::vector<VectorDouble > > &field_map,
376  std::map<std::string,std::vector<MatrixDouble > > &grad_map
377  ) {
378  MoFEMFunctionBeginHot;
379  int nb_gauss_pts = grad_map["POTENTIAL_FIELD"].size();
380  this->commonDataPtr->gradAtGaussPts["POTENTIAL_FIELD"].resize(nb_gauss_pts);
381  for(int gg = 0;gg<nb_gauss_pts;gg++) {
382  this->commonDataPtr->gradAtGaussPts["POTENTIAL_FIELD"][gg].resize(1,3,false);
383  for(int ii = 0;ii<3;ii++) {
384  this->commonDataPtr->gradAtGaussPts["POTENTIAL_FIELD"][gg](0,ii) =
385  ((grad_map["POTENTIAL_FIELD"])[gg])(0,ii);
386  }
387  }
388  MoFEMFunctionReturnHot(0);
389  }
390 
391 };
392 
393 struct SmallStrainTranverslyIsotropicADouble: public SmallStrainTranverslyIsotropic<adouble> {
394 
395  SmallStrainTranverslyIsotropicADouble(): SmallStrainTranverslyIsotropic<adouble>() {}
396 
397  int nbActiveVariables0;
398 
399  virtual MoFEMErrorCode setUserActiveVariables(
400  int &nb_active_variables
401  ) {
402  MoFEMFunctionBeginHot;
403 
404  try {
405 
406  ierr = calculateFibreAngles(); CHKERRG(ierr);
407  axVector.resize(3,false);
408  axVector[0] <<= axVectorDouble[0];
409  axVector[1] <<= axVectorDouble[1];
410  axVector[2] <<= axVectorDouble[2];
411  axAngle <<= axAngleDouble;
412  nbActiveVariables0 = nb_active_variables;
413  nb_active_variables += 4;
414 
415  } catch (const std::exception& ex) {
416  std::ostringstream ss;
417  ss << "throw in method: " << ex.what() << std::endl;
418  SETERRQ(PETSC_COMM_SELF,1,ss.str().c_str());
419  }
420 
421  MoFEMFunctionReturnHot(0);
422  }
423 
424  virtual MoFEMErrorCode setUserActiveVariables(
425  VectorDouble &active_variables) {
426  MoFEMFunctionBeginHot;
427 
428  try {
429 
430  int shift = nbActiveVariables0; // is a number of elements in F
431  ierr = calculateFibreAngles(); CHKERRG(ierr);
432  active_variables[shift+0] = axVectorDouble[0];
433  active_variables[shift+1] = axVectorDouble[1];
434  active_variables[shift+2] = axVectorDouble[2];
435  active_variables[shift+3] = axAngleDouble;
436 
437  } catch (const std::exception& ex) {
438  std::ostringstream ss;
439  ss << "throw in method: " << ex.what() << std::endl;
440  SETERRQ(PETSC_COMM_SELF,1,ss.str().c_str());
441  }
442 
443  MoFEMFunctionReturnHot(0);
444  }
445 
446 };
447 
448 #endif //__SMALLSTRAINTRANVERSLYISOTROPIC_HPP__
SmallStrainTranverslyIsotropic::tR
TYPE tR
Definition: SmallTransverselyIsotropic.hpp:27
MoFEMFunctionReturnHot
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:460
SmallStrainTranverslyIsotropicADouble::SmallStrainTranverslyIsotropicADouble
SmallStrainTranverslyIsotropicADouble()
Definition: SmallTransverselyIsotropic.hpp:395
NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI::FunctionsToCalculatePiolaKirchhoffI
FunctionsToCalculatePiolaKirchhoffI()
Definition: NonLinearElasticElement.hpp:136
NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI::eNergy
TYPE eNergy
Definition: NonLinearElasticElement.hpp:152
SmallStrainTranverslyIsotropic::localStiffnessMatrix
ublas::symmetric_matrix< TYPE, ublas::upper > localStiffnessMatrix
Definition: SmallTransverselyIsotropic.hpp:49
SmallStrainTranverslyIsotropic::G_zp
double G_zp
Definition: SmallTransverselyIsotropic.hpp:48
NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI
Implementation of elastic (non-linear) St. Kirchhoff equation.
Definition: NonLinearElasticElement.hpp:79
SmallStrainTranverslyIsotropicADouble::nbActiveVariables0
int nbActiveVariables0
Definition: SmallTransverselyIsotropic.hpp:397
SmallStrainTranverslyIsotropic::normalizedPhi
VectorDouble normalizedPhi
Definition: SmallTransverselyIsotropic.hpp:306
NonlinearElasticElement
structure grouping operators and data used for calculation of nonlinear elastic element
Definition: HookeElement.hpp:27
SmallStrainTranverslyIsotropic::calculateLocalStiffnesMatrix
MoFEMErrorCode calculateLocalStiffnesMatrix()
Definition: SmallTransverselyIsotropic.hpp:50
SmallStrainTranverslyIsotropic::calculateP_PiolaKirchhoffI
virtual MoFEMErrorCode calculateP_PiolaKirchhoffI(const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
Calculate global stress.
Definition: SmallTransverselyIsotropic.hpp:250
phi
static double phi
Definition: poisson_2d_dis_galerkin.cpp:30
MoFEM::Exceptions::MoFEMErrorCode
PetscErrorCode MoFEMErrorCode
MoFEM/PETSc error code.
Definition: Exceptions.hpp:56
MoFEM::Types::MatrixDouble
UBlasMatrix< double > MatrixDouble
Definition: Types.hpp:77
SmallStrainTranverslyIsotropic::sTrain
ublas::matrix< TYPE > sTrain
Definition: SmallTransverselyIsotropic.hpp:25
SmallStrainTranverslyIsotropic::strainTransformation
MoFEMErrorCode strainTransformation()
Function to Calculate Strain Transformation Matrix This function computes the strain transformation ...
Definition: SmallTransverselyIsotropic.hpp:171
SmallStrainTranverslyIsotropicADouble::setUserActiveVariables
virtual MoFEMErrorCode setUserActiveVariables(int &nb_active_variables)
add additional active variables
Definition: SmallTransverselyIsotropic.hpp:399
SmallStrainTranverslyIsotropic::calculateAngles
virtual MoFEMErrorCode calculateAngles()
Definition: SmallTransverselyIsotropic.hpp:236
SmallStrainTranverslyIsotropic::nu_p
double nu_p
Definition: SmallTransverselyIsotropic.hpp:48
SmallStrainTranverslyIsotropic::axVector
ublas::vector< TYPE > axVector
Definition: SmallTransverselyIsotropic.hpp:70
SmallStrainTranverslyIsotropic::E_p
double E_p
Definition: SmallTransverselyIsotropic.hpp:48
SmallStrainTranverslyIsotropic::SmallStrainTranverslyIsotropic
SmallStrainTranverslyIsotropic()
Definition: SmallTransverselyIsotropic.hpp:21
SmallStrainTranverslyIsotropic
Hook equation.
Definition: SmallTransverselyIsotropic.hpp:19
SmallStrainTranverslyIsotropicDouble::SmallStrainTranverslyIsotropicDouble
SmallStrainTranverslyIsotropicDouble()
Definition: SmallTransverselyIsotropic.hpp:352
SmallStrainTranverslyIsotropic::dR
ublas::matrix< TYPE > dR
Definition: SmallTransverselyIsotropic.hpp:222
SmallStrainTranverslyIsotropic::calculateFibreAngles
MoFEMErrorCode calculateFibreAngles()
Definition: SmallTransverselyIsotropic.hpp:310
delta
static constexpr double delta
Definition: prism_elements_from_surface.cpp:18
SmallStrainTranverslyIsotropic::axAngle
TYPE axAngle
Definition: SmallTransverselyIsotropic.hpp:71
NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI::F
MatrixBoundedArray< TYPE, 9 > F
Definition: NonLinearElasticElement.hpp:147
double
SmallStrainTranverslyIsotropic::nu_pz
double nu_pz
Definition: SmallTransverselyIsotropic.hpp:48
SmallStrainTranverslyIsotropicADouble
Definition: SmallTransverselyIsotropic.hpp:393
SmallStrainTranverslyIsotropic::stressRotMat
ublas::matrix< TYPE > stressRotMat
Definition: SmallTransverselyIsotropic.hpp:106
SmallStrainTranverslyIsotropic::globalStiffnessMatrix
ublas::matrix< TYPE > globalStiffnessMatrix
Definition: SmallTransverselyIsotropic.hpp:223
NonlinearElasticElement::BlockData
data for calculation heat conductivity and heat capacity elements
Definition: HookeElement.hpp:32
SmallStrainTranverslyIsotropic::E_z
double E_z
Definition: SmallTransverselyIsotropic.hpp:48
SmallStrainTranverslyIsotropic::axAngleDouble
double axAngleDouble
Definition: SmallTransverselyIsotropic.hpp:308
SmallStrainTranverslyIsotropicDouble::getDataOnPostProcessor
virtual MoFEMErrorCode getDataOnPostProcessor(std::map< std::string, std::vector< VectorDouble > > &field_map, std::map< std::string, std::vector< MatrixDouble > > &grad_map)
Definition: SmallTransverselyIsotropic.hpp:374
SmallStrainTranverslyIsotropic::stressTransformation
MoFEMErrorCode stressTransformation()
Function to Calculate Stress Transformation Matrix This function computes the stress transformation M...
Definition: SmallTransverselyIsotropic.hpp:113
SmallStrainTranverslyIsotropic::aARotMat
ublas::matrix< TYPE > aARotMat
Definition: SmallTransverselyIsotropic.hpp:69
adouble
FTensor::dd
const Tensor2_symmetric_Expr< const ddTensor0< T, Dim, i, j >, typename promote< T, double >::V, Dim, i, j > dd(const Tensor0< T * > &a, const Index< i, Dim > index1, const Index< j, Dim > index2, const Tensor1< int, Dim > &d_ijk, const Tensor1< double, Dim > &d_xyz)
Definition: ddTensor0.hpp:33
SmallStrainTranverslyIsotropic::calculateStrain
MoFEMErrorCode calculateStrain()
Definition: SmallTransverselyIsotropic.hpp:29
eps
static const double eps
Definition: check_base_functions_derivatives_on_tet.cpp:11
MoFEM::Exceptions::ierr
static MoFEMErrorCodeGeneric< PetscErrorCode > ierr
Definition: Exceptions.hpp:76
SmallStrainTranverslyIsotropic::strainRotMat
ublas::matrix< TYPE > strainRotMat
Definition: SmallTransverselyIsotropic.hpp:164
MoFEM::Types::VectorDouble
UBlasVector< double > VectorDouble
Definition: Types.hpp:68
SmallStrainTranverslyIsotropic::voigtStress
ublas::vector< TYPE > voigtStress
Definition: SmallTransverselyIsotropic.hpp:241
SmallStrainTranverslyIsotropic::calculateElasticEnergy
virtual MoFEMErrorCode calculateElasticEnergy(const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
calculate density of strain energy
Definition: SmallTransverselyIsotropic.hpp:284
MoFEMFunctionBeginHot
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:453
SmallStrainTranverslyIsotropicDouble::calculateAngles
virtual MoFEMErrorCode calculateAngles()
Definition: SmallTransverselyIsotropic.hpp:354
SmallStrainTranverslyIsotropicDouble
Definition: SmallTransverselyIsotropic.hpp:350
SmallStrainTranverslyIsotropic::calculateGlobalStiffnesMatrix
MoFEMErrorCode calculateGlobalStiffnesMatrix()
Definition: SmallTransverselyIsotropic.hpp:225
SmallStrainTranverslyIsotropic::voightStrain
ublas::vector< TYPE > voightStrain
Definition: SmallTransverselyIsotropic.hpp:26
NonlinearElasticElement::CommonData::gradAtGaussPts
std::map< std::string, std::vector< MatrixDouble > > gradAtGaussPts
Definition: NonLinearElasticElement.hpp:108
CHKERRG
#define CHKERRG(n)
Check error code of MoFEM/MOAB/PETSc function.
Definition: definitions.h:496
SmallStrainTranverslyIsotropicADouble::setUserActiveVariables
virtual MoFEMErrorCode setUserActiveVariables(VectorDouble &active_variables)
Add additional independent variables More complex physical models depend on gradient of defamation an...
Definition: SmallTransverselyIsotropic.hpp:424
SmallStrainTranverslyIsotropic::axVectorDouble
VectorDouble axVectorDouble
Definition: SmallTransverselyIsotropic.hpp:307
NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI::gG
int gG
Gauss point number.
Definition: NonLinearElasticElement.hpp:154
NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI::P
MatrixBoundedArray< TYPE, 9 > P
Definition: NonLinearElasticElement.hpp:147
NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI::commonDataPtr
CommonData * commonDataPtr
Definition: NonLinearElasticElement.hpp:155
SmallStrainTranverslyIsotropic::calculateAxisAngleRotationalMatrix
MoFEMErrorCode calculateAxisAngleRotationalMatrix()
Function to Calculate the Rotation Matrix at a given axis and angle of rotation.
Definition: SmallTransverselyIsotropic.hpp:82
Generated by Doxygen 1.8.17 and hosted at University of Glasgow