MAST
Multidisciplinary-design Adaptation and Sensitivity Toolkit (MAST)
structural_element_2d.h
Go to the documentation of this file.
1 /*
2  * MAST: Multidisciplinary-design Adaptation and Sensitivity Toolkit
3  * Copyright (C) 2013-2020 Manav Bhatia and MAST authors
4  *
5  * This library is free software; you can redistribute it and/or
6  * modify it under the terms of the GNU Lesser General Public
7  * License as published by the Free Software Foundation; either
8  * version 2.1 of the License, or (at your option) any later version.
9  *
10  * This library is distributed in the hope that it will be useful,
11  * but WITHOUT ANY WARRANTY; without even the implied warranty of
12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13  * Lesser General Public License for more details.
14  *
15  * You should have received a copy of the GNU Lesser General Public
16  * License along with this library; if not, write to the Free Software
17  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
18  */
19 
20 #ifndef __mast__structural_element_2d__
21 #define __mast__structural_element_2d__
22 
23 // C++ includes
24 #include <memory>
25 
26 
27 // MAST includes
28 #include "elasticity/bending_structural_element.h"
29 
30 
31 
32 
33 namespace MAST {
34 
35  // Forward declerations
36  class BendingOperator2D;
37  class BoundaryCondition;
38  class FEMOperatorMatrix;
39 
40 
41 
42  class StructuralElement2D:
43  public MAST::BendingStructuralElem {
44 
45  public:
46  StructuralElement2D(MAST::SystemInitialization& sys,
47  const MAST::GeomElem& elem,
48  const MAST::ElementPropertyCardBase& p);
49 
50  virtual ~StructuralElement2D();
51 
56  virtual unsigned int n_direct_strain_components() {
57  return 3;
58  }
59 
63  virtual unsigned int n_von_karman_strain_components() {
64  return 2;
65  }
66 
71  virtual bool internal_residual(bool request_jacobian,
72  RealVectorX& f,
73  RealMatrixX& jac);
74 
79  virtual bool internal_residual_sensitivity(const MAST::FunctionBase& p,
80  bool request_jacobian,
81  RealVectorX& f,
82  RealMatrixX& jac);
86  virtual bool
87  internal_residual_jac_dot_state_sensitivity (RealMatrixX& jac);
88 
94  virtual void
95  internal_residual_boundary_velocity(const MAST::FunctionBase& p,
96  const unsigned int s,
97  const MAST::FieldFunction<RealVectorX>& vel_f,
98  bool request_jacobian,
99  RealVectorX& f,
100  RealMatrixX& jac);
101 
106  virtual bool prestress_residual (bool request_jacobian,
107  RealVectorX& f,
108  RealMatrixX& jac);
109 
114  virtual bool prestress_residual_sensitivity (const MAST::FunctionBase& p,
115  bool request_jacobian,
116  RealVectorX& f,
117  RealMatrixX& jac);
118 
123  virtual bool if_incompatible_modes() const {
124  return false;
125  }
126 
130  virtual bool calculate_stress(bool request_derivative,
131  const MAST::FunctionBase* p,
132  MAST::StressStrainOutputBase& output);
133 
138  virtual void
139  calculate_stress_boundary_velocity(const MAST::FunctionBase& p,
140  MAST::StressStrainOutputBase& output,
141  const unsigned int s,
142  const MAST::FieldFunction<RealVectorX>& vel_f);
143 
144 
145  virtual void
146  calculate_stress_temperature_derivative(MAST::FEBase& fe_thermal,
147  MAST::StressStrainOutputBase& output);
148 
149  virtual void
150  thermal_residual_temperature_derivative (const MAST::FEBase& fe_thermal,
151  RealMatrixX& m);
152 
153  protected:
154 
155 
159  virtual bool
160  surface_traction_residual(bool request_jacobian,
161  RealVectorX& f,
162  RealMatrixX& jac,
163  const unsigned int side,
164  MAST::BoundaryConditionBase& bc);
165 
166 
171  virtual bool
172  surface_traction_residual_sensitivity(const MAST::FunctionBase& p,
173  bool request_jacobian,
174  RealVectorX& f,
175  RealMatrixX& jac,
176  const unsigned int side,
177  MAST::BoundaryConditionBase& bc);
178 
179 
184  virtual bool
185  surface_traction_residual_shifted_boundary(bool request_jacobian,
186  RealVectorX& f,
187  RealMatrixX& jac,
188  const unsigned int side,
189  MAST::BoundaryConditionBase& bc);
190 
195  virtual bool
196  surface_traction_residual_shifted_boundary_sensitivity(const MAST::FunctionBase& p,
197  bool request_jacobian,
198  RealVectorX& f,
199  RealMatrixX& jac,
200  const unsigned int side,
201  MAST::BoundaryConditionBase& bc);
202 
203 
207  bool surface_pressure_residual(bool request_jacobian,
208  RealVectorX &f,
209  RealMatrixX &jac,
210  const unsigned int side,
211  MAST::BoundaryConditionBase& bc);
212 
216  virtual bool
217  surface_pressure_residual_sensitivity(const MAST::FunctionBase& p,
218  bool request_jacobian,
219  RealVectorX& f,
220  RealMatrixX& jac,
221  const unsigned int side,
222  MAST::BoundaryConditionBase& bc);
223 
224 
228  virtual bool thermal_residual(bool request_jacobian,
229  RealVectorX& f,
230  RealMatrixX& jac,
231  MAST::BoundaryConditionBase& bc);
232 
237  virtual bool thermal_residual_sensitivity(const MAST::FunctionBase& p,
238  bool request_jacobian,
239  RealVectorX& f,
240  RealMatrixX& jac,
241  MAST::BoundaryConditionBase& bc);
242 
248  virtual void
249  thermal_residual_boundary_velocity(const MAST::FunctionBase& p,
250  const unsigned int s,
251  const MAST::FieldFunction<RealVectorX>& vel_f,
252  MAST::BoundaryConditionBase& bc,
253  bool request_jacobian,
254  RealVectorX& f,
255  RealMatrixX& jac);
256 
264  virtual bool piston_theory_residual(bool request_jacobian,
265  RealVectorX &f,
266  RealMatrixX& jac_xdot,
267  RealMatrixX& jac,
268  MAST::BoundaryConditionBase& bc);
269 
277  virtual bool
278  piston_theory_residual_sensitivity(const MAST::FunctionBase& p,
279  bool request_jacobian,
280  RealVectorX &f,
281  RealMatrixX& jac_xdot,
282  RealMatrixX& jac,
283  MAST::BoundaryConditionBase& bc);
284 
291  virtual bool piston_theory_residual(bool request_jacobian,
292  RealVectorX &f,
293  RealMatrixX& jac_xdot,
294  RealMatrixX& jac,
295  const unsigned int side,
296  MAST::BoundaryConditionBase& bc);
297 
298 
305  virtual bool
306  piston_theory_residual_sensitivity(const MAST::FunctionBase& p,
307  bool request_jacobian,
308  RealVectorX &f,
309  RealMatrixX& jac_xdot,
310  RealMatrixX& jac,
311  const unsigned int side,
312  MAST::BoundaryConditionBase& bc);
313 
314 
319  virtual bool
320  linearized_frequency_domain_surface_pressure_residual
321  (bool request_jacobian,
322  ComplexVectorX& f,
323  ComplexMatrixX& jac,
324  const unsigned int side,
325  MAST::BoundaryConditionBase& bc) {
326 
327  libmesh_error(); // to be implemented
328  }
329 
330 
331 
336  virtual bool
337  linearized_frequency_domain_surface_pressure_residual_sensitivity
338  (const MAST::FunctionBase& p,
339  bool request_jacobian,
340  ComplexVectorX& f,
341  ComplexMatrixX& jac,
342  const unsigned int side,
343  MAST::BoundaryConditionBase& bc) {
344 
345  libmesh_error(); // to be implemented
346  }
347 
348 
355  virtual void
356  initialize_von_karman_strain_operator(const unsigned int qp,
357  const MAST::FEBase& fe,
358  RealVectorX& vk_strain,
359  RealMatrixX& vk_dwdxi_mat,
360  MAST::FEMOperatorMatrix& Bmat_vk);
361 
367  virtual void
368  initialize_von_karman_strain_operator_sensitivity(const unsigned int qp,
369  const MAST::FEBase& fe,
370  RealMatrixX& vk_dwdxi_mat_sens);
371 
372 
373  virtual void
374  initialize_green_lagrange_strain_operator(const unsigned int qp,
375  const MAST::FEBase& fe,
376  const RealVectorX& local_disp,
377  RealVectorX& epsilon,
378  RealMatrixX& mat_x,
379  RealMatrixX& mat_y,
380  MAST::FEMOperatorMatrix& Bmat_lin,
381  MAST::FEMOperatorMatrix& Bmat_nl_x,
382  MAST::FEMOperatorMatrix& Bmat_nl_y,
383  MAST::FEMOperatorMatrix& Bmat_nl_u,
384  MAST::FEMOperatorMatrix& Bmat_nl_v);
385 
386 
387  virtual void
388  initialize_strain_operator_gradient(const unsigned int qp,
389  const MAST::FEBase& fe,
390  const RealVectorX& local_disp,
391  RealMatrixX& epsilon_grad,
392  std::vector<MAST::FEMOperatorMatrix>& dBmat_lin);
393 
399  void _compute_stress(MAST::FEBase& fe, unsigned int qp,
400  RealVectorX& stress,
401  RealMatrixX* dstress_dX);
402 
403 
404  void _compute_stress_sensitivity(const MAST::FunctionBase& f,
405  MAST::FEBase& fe, unsigned int qp,
406  RealVectorX& stress);
407 
414  void _compute_stress_gradient(MAST::FEBase& fe,
415  unsigned int qp,
416  RealMatrixX& stress_grad,
417  std::vector<RealMatrixX>* dstress_grad_dX);
418 
419  void _compute_stress_gradient_sensitivity(const MAST::FunctionBase& f,
420  MAST::FEBase& fe,
421  unsigned int qp,
422  RealMatrixX& stress_grad);
423 
428  void _surface_normal_voigt_notation(const RealVectorX& normal,
429  RealMatrixX& normal_mat);
430 
431 
437  virtual void
438  _internal_residual_operation(bool if_vk,
439  const unsigned int n2,
440  const unsigned int qp,
441  const MAST::FEBase& fe,
442  const std::vector<Real>& JxW,
443  bool request_jacobian,
444  RealVectorX& local_f,
445  RealMatrixX& local_jac,
446  RealVectorX& local_disp,
447  RealVectorX& strain_mem,
448  MAST::BendingOperator2D* bend,
449  FEMOperatorMatrix& Bmat_lin,
450  FEMOperatorMatrix& Bmat_nl_x,
451  FEMOperatorMatrix& Bmat_nl_y,
452  FEMOperatorMatrix& Bmat_nl_u,
453  FEMOperatorMatrix& Bmat_nl_v,
454  MAST::FEMOperatorMatrix& Bmat_bend,
455  MAST::FEMOperatorMatrix& Bmat_vk,
456  RealMatrixX& mat_x,
457  RealMatrixX& mat_y,
458  RealMatrixX& stress,
459  RealMatrixX& vk_dwdxi_mat,
460  RealMatrixX& material_A_mat,
461  RealMatrixX& material_B_mat,
462  RealMatrixX& material_D_mat,
463  RealVectorX& vec1_n1,
464  RealVectorX& vec2_n1,
465  RealVectorX& vec3_n2,
466  RealVectorX& vec4_2,
467  RealVectorX& vec5_2,
468  RealVectorX& vec6_n2,
469  RealMatrixX& mat1_n1n2,
470  RealMatrixX& mat2_n2n2,
471  RealMatrixX& mat3,
472  RealMatrixX& mat4_2n2,
473  RealMatrixX& mat5_3n2);
474 
475 
479  void _convert_prestress_A_mat_to_vector(const RealMatrixX& mat,
480  RealVectorX& vec) const;
481 
485  void _convert_prestress_B_mat_to_vector(const RealMatrixX& mat,
486  RealVectorX& vec) const;
487 
488  };
489 }
490 
491 
492 
493 #endif // __mast__structural_element_2d__
MAST::StructuralElement2D::calculate_stress
virtual bool calculate_stress(bool request_derivative, const MAST::FunctionBase *p, MAST::StressStrainOutputBase &output)
Calculates the stress tensor.
Definition: structural_element_2d.cpp:307
MAST::StructuralElement2D::surface_traction_residual_shifted_boundary_sensitivity
virtual bool surface_traction_residual_shifted_boundary_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)
Calculates the sensitivity of force vector and Jacobian due to surface traction and sensitiity due to...
Definition: structural_element_2d.cpp:2663
MAST::StructuralElement2D::_convert_prestress_B_mat_to_vector
void _convert_prestress_B_mat_to_vector(const RealMatrixX &mat, RealVectorX &vec) const
converts the prestress stress tensor to a vector representation
Definition: structural_element_2d.cpp:2039
MAST::StructuralElement2D::_surface_normal_voigt_notation
void _surface_normal_voigt_notation(const RealVectorX &normal, RealMatrixX &normal_mat)
creates a matrix that can be multiplied with the Voigt notation of stress to compute the surface norm...
Definition: structural_element_2d.cpp:5053
MAST::StressStrainOutputBase
Data structure provides the mechanism to store stress and strain output from a structural analysis...
Definition: stress_output_base.h:61
MAST::StructuralElement2D::surface_traction_residual_sensitivity
virtual bool surface_traction_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)
Calculates the sensitivity of element vector and matrix quantities for surface traction boundary cond...
Definition: structural_element_2d.cpp:2431
MAST::StructuralElement2D::n_direct_strain_components
virtual unsigned int n_direct_strain_components()
row dimension of the direct strain matrix, also used for the bending operator row dimension ...
Definition: structural_element_2d.h:56
MAST::StructuralElement2D::initialize_strain_operator_gradient
virtual void initialize_strain_operator_gradient(const unsigned int qp, const MAST::FEBase &fe, const RealVectorX &local_disp, RealMatrixX &epsilon_grad, std::vector< MAST::FEMOperatorMatrix > &dBmat_lin)
Definition: structural_element_2d.cpp:255
MAST::StructuralElement2D::_internal_residual_operation
virtual void _internal_residual_operation(bool if_vk, const unsigned int n2, const unsigned int qp, const MAST::FEBase &fe, const std::vector< Real > &JxW, bool request_jacobian, RealVectorX &local_f, RealMatrixX &local_jac, RealVectorX &local_disp, RealVectorX &strain_mem, MAST::BendingOperator2D *bend, FEMOperatorMatrix &Bmat_lin, FEMOperatorMatrix &Bmat_nl_x, FEMOperatorMatrix &Bmat_nl_y, FEMOperatorMatrix &Bmat_nl_u, FEMOperatorMatrix &Bmat_nl_v, MAST::FEMOperatorMatrix &Bmat_bend, MAST::FEMOperatorMatrix &Bmat_vk, RealMatrixX &mat_x, RealMatrixX &mat_y, RealMatrixX &stress, RealMatrixX &vk_dwdxi_mat, RealMatrixX &material_A_mat, RealMatrixX &material_B_mat, RealMatrixX &material_D_mat, RealVectorX &vec1_n1, RealVectorX &vec2_n1, RealVectorX &vec3_n2, RealVectorX &vec4_2, RealVectorX &vec5_2, RealVectorX &vec6_n2, RealMatrixX &mat1_n1n2, RealMatrixX &mat2_n2n2, RealMatrixX &mat3, RealMatrixX &mat4_2n2, RealMatrixX &mat5_3n2)
performs integration at the quadrature point for the provided matrices.
Definition: structural_element_2d.cpp:1767
MAST::StructuralElement2D::prestress_residual_sensitivity
virtual bool prestress_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac)
Calculates the internal residual vector and Jacobian due to strain energy coming from a prestress...
Definition: structural_element_2d.cpp:2205
ComplexVectorX
Matrix< Complex, Dynamic, 1 > ComplexVectorX
Definition: mast_data_types.h:37
MAST::StructuralElement2D::internal_residual_jac_dot_state_sensitivity
virtual bool internal_residual_jac_dot_state_sensitivity(RealMatrixX &jac)
calculates d[J]/d{x} .
Definition: structural_element_2d.cpp:1572
bending_structural_element.h
MAST::StructuralElement2D::_compute_stress
void _compute_stress(MAST::FEBase &fe, unsigned int qp, RealVectorX &stress, RealMatrixX *dstress_dX)
Computes the stress at the specified quadrature point of the FE data structure.
Definition: structural_element_2d.cpp:4178
MAST::StructuralElement2D::surface_traction_residual_shifted_boundary
virtual bool surface_traction_residual_shifted_boundary(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)
Calculates the sensitivity of force vector and Jacobian due to surface traction and sensitiity due to...
Definition: structural_element_2d.cpp:2510
MAST::StructuralElement2D::_compute_stress_sensitivity
void _compute_stress_sensitivity(const MAST::FunctionBase &f, MAST::FEBase &fe, unsigned int qp, RealVectorX &stress)
Definition: structural_element_2d.cpp:4477
MAST::StructuralElement2D::linearized_frequency_domain_surface_pressure_residual
virtual bool linearized_frequency_domain_surface_pressure_residual(bool request_jacobian, ComplexVectorX &f, ComplexMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)
Calculates the force vector and Jacobian due to small perturbation surface pressure.
Definition: structural_element_2d.h:321
MAST::StructuralElement2D::thermal_residual_temperature_derivative
virtual void thermal_residual_temperature_derivative(const MAST::FEBase &fe_thermal, RealMatrixX &m)
Definition: structural_element_2d.cpp:3625
MAST::StructuralElement2D::internal_residual
virtual bool internal_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac)
Calculates the internal residual vector and Jacobian due to strain energy.
Definition: structural_element_2d.cpp:1076
MAST::StructuralElement2D::calculate_stress_temperature_derivative
virtual void calculate_stress_temperature_derivative(MAST::FEBase &fe_thermal, MAST::StressStrainOutputBase &output)
Definition: structural_element_2d.cpp:958
MAST::StructuralElement2D::piston_theory_residual
virtual bool piston_theory_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac, MAST::BoundaryConditionBase &bc)
Calculates the force vector and Jacobian due to piston-theory based surface pressure on the entire el...
Definition: structural_element_2d.cpp:3818
MAST::StructuralElement2D::thermal_residual_boundary_velocity
virtual void thermal_residual_boundary_velocity(const MAST::FunctionBase &p, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f, MAST::BoundaryConditionBase &bc, bool request_jacobian, RealVectorX &f, RealMatrixX &jac)
calculates the term on side s: .
Definition: structural_element_2d.cpp:3421
RealMatrixX
Matrix< Real, Dynamic, Dynamic > RealMatrixX
Definition: mast_data_types.h:44
ComplexMatrixX
Matrix< Complex, Dynamic, Dynamic > ComplexMatrixX
Definition: mast_data_types.h:46
MAST::StructuralElement2D::thermal_residual_sensitivity
virtual bool thermal_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac, MAST::BoundaryConditionBase &bc)
Calculates the sensitivity of residual vector and the Jacobian due to thermal stresses.
Definition: structural_element_2d.cpp:3208
MAST::BendingOperator2D
Bending strain operator for 1D element.
Definition: bending_operator.h:157
MAST::StructuralElement2D::if_incompatible_modes
virtual bool if_incompatible_modes() const
Definition: structural_element_2d.h:123
MAST::StructuralElement2D::prestress_residual
virtual bool prestress_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac)
Calculates the internal residual vector and Jacobian due to strain energy coming from a prestress...
Definition: structural_element_2d.cpp:2054
RealVectorX
Matrix< Real, Dynamic, 1 > RealVectorX
Definition: mast_data_types.h:35
MAST::StructuralElement2D::StructuralElement2D
StructuralElement2D(MAST::SystemInitialization &sys, const MAST::GeomElem &elem, const MAST::ElementPropertyCardBase &p)
Definition: structural_element_2d.cpp:39
MAST::StructuralElement2D::initialize_von_karman_strain_operator
virtual void initialize_von_karman_strain_operator(const unsigned int qp, const MAST::FEBase &fe, RealVectorX &vk_strain, RealMatrixX &vk_dwdxi_mat, MAST::FEMOperatorMatrix &Bmat_vk)
initialze the von Karman strain in vK_strain, the operator matrices needed for Jacobian calculation...
Definition: structural_element_2d.cpp:57
MAST::StructuralElement2D::initialize_von_karman_strain_operator_sensitivity
virtual void initialize_von_karman_strain_operator_sensitivity(const unsigned int qp, const MAST::FEBase &fe, RealMatrixX &vk_dwdxi_mat_sens)
initialze the sensitivity of von Karman operator matrices needed for Jacobian calculation.
Definition: structural_element_2d.cpp:107
MAST::GeomElem
This class acts as a wrapper around libMesh::Elem for the purpose of providing a uniform interface fo...
Definition: geom_elem.h:59
MAST::StructuralElement2D::surface_pressure_residual
bool surface_pressure_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)
Calculates the force vector and Jacobian due to surface pressure.
Definition: structural_element_2d.cpp:2857
MAST::StructuralElement2D::internal_residual_sensitivity
virtual bool internal_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac)
Calculates the sensitivity internal residual vector and Jacobian due to strain energy.
Definition: structural_element_2d.cpp:1234
MAST::ElementBase::elem
const MAST::GeomElem & elem() const
Definition: elem_base.h:109
MAST::StructuralElement2D::surface_traction_residual
virtual bool surface_traction_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)
Calculates the force vector and Jacobian due to surface traction.
Definition: structural_element_2d.cpp:2358
MAST::StructuralElement2D::piston_theory_residual_sensitivity
virtual bool piston_theory_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac, MAST::BoundaryConditionBase &bc)
Calculates the sensitivity of force vector and Jacobian due to piston-theory based surface pressure o...
Definition: structural_element_2d.cpp:3996
MAST::StructuralElement2D::_convert_prestress_A_mat_to_vector
void _convert_prestress_A_mat_to_vector(const RealMatrixX &mat, RealVectorX &vec) const
converts the prestress stress tensor to a vector representation
Definition: structural_element_2d.cpp:2025
MAST::StructuralElement2D::surface_pressure_residual_sensitivity
virtual bool surface_pressure_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)
Calculates the force vector and Jacobian due to surface pressure.
Definition: structural_element_2d.cpp:2932
MAST::StructuralElement2D::_compute_stress_gradient
void _compute_stress_gradient(MAST::FEBase &fe, unsigned int qp, RealMatrixX &stress_grad, std::vector< RealMatrixX > *dstress_grad_dX)
computes the gradient of stress in Voigt notation in stress_grad where the three columns are the deri...
Definition: structural_element_2d.cpp:4657
MAST::StructuralElement2D::linearized_frequency_domain_surface_pressure_residual_sensitivity
virtual bool linearized_frequency_domain_surface_pressure_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, ComplexVectorX &f, ComplexMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)
Calculates the sensitivity of force vector and Jacobian due to small is applicable for perturbation s...
Definition: structural_element_2d.h:338
MAST::StructuralElement2D::calculate_stress_boundary_velocity
virtual void calculate_stress_boundary_velocity(const MAST::FunctionBase &p, MAST::StressStrainOutputBase &output, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f)
Calculates the boundary velocity term contributions to the sensitivity of stress at the specified bou...
Definition: structural_element_2d.cpp:716
MAST::StructuralElement2D::initialize_green_lagrange_strain_operator
virtual void initialize_green_lagrange_strain_operator(const unsigned int qp, const MAST::FEBase &fe, const RealVectorX &local_disp, RealVectorX &epsilon, RealMatrixX &mat_x, RealMatrixX &mat_y, MAST::FEMOperatorMatrix &Bmat_lin, MAST::FEMOperatorMatrix &Bmat_nl_x, MAST::FEMOperatorMatrix &Bmat_nl_y, MAST::FEMOperatorMatrix &Bmat_nl_u, MAST::FEMOperatorMatrix &Bmat_nl_v)
Definition: structural_element_2d.cpp:143
MAST::StructuralElement2D::internal_residual_boundary_velocity
virtual void internal_residual_boundary_velocity(const MAST::FunctionBase &p, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f, bool request_jacobian, RealVectorX &f, RealMatrixX &jac)
calculates the term on side s: .
Definition: structural_element_2d.cpp:1403
MAST::StructuralElement2D::_compute_stress_gradient_sensitivity
void _compute_stress_gradient_sensitivity(const MAST::FunctionBase &f, MAST::FEBase &fe, unsigned int qp, RealMatrixX &stress_grad)
Definition: structural_element_2d.cpp:4873
MAST::StructuralElement2D::thermal_residual
virtual bool thermal_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to thermal stresses.
Definition: structural_element_2d.cpp:3012
MAST::StructuralElement2D::n_von_karman_strain_components
virtual unsigned int n_von_karman_strain_components()
row dimension of the von Karman strain matrix
Definition: structural_element_2d.h:63