MAST
Multidisciplinary-design Adaptation and Sensitivity Toolkit (MAST)
structural_element_base.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_base__
21 #define __mast__structural_element_base__
22 
23 // C++ includes
24 #include <memory>
25 #include <map>
26 
27 // MAST includes
28 #include "base/elem_base.h"
29 
30 
31 namespace MAST {
32 
33  // Forward declerations
34  class GeomElem;
35  class ElementPropertyCardBase;
36  class BoundaryConditionBase;
37  class FEMOperatorMatrix;
38  class StressStrainOutputBase;
39  template <typename ValType> class FieldFunction;
40 
41 
42  class StructuralElementBase:
43  public MAST::ElementBase
44  {
45  public:
49  StructuralElementBase(MAST::SystemInitialization& sys,
50  const MAST::GeomElem& elem,
51  const MAST::ElementPropertyCardBase& p);
52 
53  virtual ~StructuralElementBase();
54 
55 
60  virtual void set_solution(const RealVectorX& vec,
61  bool if_sens = false);
62 
63 
68  virtual void set_perturbed_solution(const RealVectorX& vec,
69  bool if_sens = false);
70 
71 
76  virtual void set_velocity(const RealVectorX& vec,
77  bool if_sens = false);
78 
83  virtual void set_perturbed_velocity(const RealVectorX& vec,
84  bool if_sens = false);
85 
86 
91  virtual void set_acceleration(const RealVectorX& vec,
92  bool if_sens = false);
93 
94 
99  virtual void set_perturbed_acceleration(const RealVectorX& vec,
100  bool if_sens = false);
101 
102 
108  const RealVectorX& local_solution(bool if_sens = false) const {
109 
110  if (!if_sens)
111  return _local_sol;
112  else
113  return _local_sol_sens;
114  }
115 
119  const MAST::ElementPropertyCardBase& elem_property() {
120  return _property;
121  }
122 
123 
127  virtual bool internal_residual (bool request_jacobian,
128  RealVectorX& f,
129  RealMatrixX& jac) = 0;
130 
135  virtual bool
136  linearized_internal_residual (bool request_jacobian,
137  RealVectorX& f,
138  RealMatrixX& jac);
139 
145  virtual void
146  internal_residual_boundary_velocity(const MAST::FunctionBase& p,
147  const unsigned int s,
148  const MAST::FieldFunction<RealVectorX>& vel_f,
149  bool request_jacobian,
150  RealVectorX& f,
151  RealMatrixX& jac) = 0;
152 
156  virtual bool
157  internal_residual_jac_dot_state_sensitivity (RealMatrixX& jac) = 0;
158 
159 
160 
164  virtual bool damping_residual (bool request_jacobian,
165  RealVectorX& f,
166  RealMatrixX& jac_xdot,
167  RealMatrixX& jac) {
168 
169  // to be implemented nothing is done yet
170  return false;
171  }
172 
176  virtual bool inertial_residual (bool request_jacobian,
177  RealVectorX& f,
178  RealMatrixX& jac_xddot,
179  RealMatrixX& jac_xdot,
180  RealMatrixX& jac);
181 
186  virtual bool linearized_inertial_residual (bool request_jacobian,
187  RealVectorX& f,
188  RealMatrixX& jac_xddot,
189  RealMatrixX& jac_xdot,
190  RealMatrixX& jac);
191 
192 
199  bool side_external_residual (bool request_jacobian,
200  RealVectorX& f,
201  RealMatrixX& jac_xdot,
202  RealMatrixX& jac,
203  std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
204 
211  bool
212  linearized_side_external_residual
213  (bool request_jacobian,
214  RealVectorX& f,
215  RealMatrixX& jac_xdot,
216  RealMatrixX& jac,
217  std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
218 
219 
220 
227  bool
228  linearized_frequency_domain_side_external_residual
229  (bool request_jacobian,
230  ComplexVectorX& f,
231  ComplexMatrixX& jac,
232  std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
233 
234 
238  virtual bool prestress_residual (bool request_jacobian,
239  RealVectorX& f,
240  RealMatrixX& jac) = 0;
241 
248  bool volume_external_residual (bool request_jacobian,
249  RealVectorX& f,
250  RealMatrixX& jac_xdot,
251  RealMatrixX& jac,
252  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
253 
254 
258  void volume_external_residual_boundary_velocity
259  (const MAST::FunctionBase& p,
260  const unsigned int s,
261  const MAST::FieldFunction<RealVectorX>& vel_f,
262  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc,
263  bool request_jacobian,
264  RealVectorX& f,
265  RealMatrixX& jac);
266 
267 
274  bool
275  linearized_volume_external_residual
276  (bool request_jacobian,
277  RealVectorX& f,
278  RealMatrixX& jac_xdot,
279  RealMatrixX& jac,
280  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
281 
282 
288  bool
289  linearized_frequency_domain_volume_external_residual
290  (bool request_jacobian,
291  ComplexVectorX& f,
292  ComplexMatrixX& jac,
293  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
294 
295 
296 
300  virtual bool internal_residual_sensitivity (const MAST::FunctionBase& p,
301  bool request_jacobian,
302  RealVectorX& f,
303  RealMatrixX& jac) = 0;
304 
308  virtual bool damping_residual_sensitivity (const MAST::FunctionBase& p,
309  bool request_jacobian,
310  RealVectorX& f,
311  RealMatrixX& jac_xdot,
312  RealMatrixX& jac) {
313  return false;
314  }
315 
319  virtual bool inertial_residual_sensitivity (const MAST::FunctionBase& p,
320  bool request_jacobian,
321  RealVectorX& f,
322  RealMatrixX& jac_xddot,
323  RealMatrixX& jac_xdot,
324  RealMatrixX& jac);
325 
329  virtual void
330  inertial_residual_boundary_velocity (const MAST::FunctionBase& p,
331  const unsigned int s,
332  const MAST::FieldFunction<RealVectorX>& vel_f,
333  bool request_jacobian,
334  RealVectorX& f,
335  RealMatrixX& jac_xddot,
336  RealMatrixX& jac_xdot,
337  RealMatrixX& jac);
338 
342  bool side_external_residual_sensitivity (const MAST::FunctionBase& p,
343  bool request_jacobian,
344  RealVectorX& f,
345  RealMatrixX& jac_xdot,
346  RealMatrixX& jac,
347  std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
348 
352  virtual bool prestress_residual_sensitivity (const MAST::FunctionBase& p,
353  bool request_jacobian,
354  RealVectorX& f,
355  RealMatrixX& jac) = 0;
356 
360  bool volume_external_residual_sensitivity (const MAST::FunctionBase& p,
361  bool request_jacobian,
362  RealVectorX& f,
363  RealMatrixX& jac_xdot,
364  RealMatrixX& jac,
365  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
366 
367 
368 
373  virtual bool if_incompatible_modes() const = 0;
374 
375 
379  virtual unsigned int incompatible_mode_size() const {
380  // needs to be implemented only in inherited elements
381  // that use incompatible modes
382  libmesh_error();
383  return 0;
384  }
385 
386 
392  void set_incompatible_mode_solution(RealVectorX& vec) {
393  _incompatible_sol = &vec;
394  }
395 
396 
402  virtual void update_incompatible_mode_solution(const RealVectorX& dsol) {
403  // should be implemented in derived classes, if relevant
404  libmesh_error();
405  }
406 
407 
413  virtual bool calculate_stress(bool request_derivative,
414  const MAST::FunctionBase* f,
415  MAST::StressStrainOutputBase& output) = 0;
416 
417 
422  virtual void
423  calculate_stress_boundary_velocity(const MAST::FunctionBase& p,
424  MAST::StressStrainOutputBase& output,
425  const unsigned int s,
426  const MAST::FieldFunction<RealVectorX>& vel_f) = 0;
427 
428  virtual void
429  calculate_stress_temperature_derivative(MAST::FEBase& fe_thermal,
430  MAST::StressStrainOutputBase& output) = 0;
431 
432  virtual void
433  thermal_residual_temperature_derivative (const MAST::FEBase& fe_thermal,
434  RealMatrixX& m) = 0;
435 
439  bool follower_forces;
440 
441 
442  template <typename ValType>
443  void transform_matrix_to_global_system(const ValType& local_mat,
444  ValType& global_mat) const;
445 
446  template <typename ValType>
447  void transform_vector_to_local_system(const ValType& global_vec,
448  ValType& local_vec) const;
449 
450  template <typename ValType>
451  void transform_vector_to_global_system(const ValType& local_vec,
452  ValType& global_vec) const;
453 
454  protected:
455 
459  virtual bool
460  surface_traction_residual(bool request_jacobian,
461  RealVectorX& f,
462  RealMatrixX& jac,
463  const unsigned int side,
464  MAST::BoundaryConditionBase& bc) = 0;
465 
470  virtual bool
471  surface_traction_residual_sensitivity(const MAST::FunctionBase& p,
472  bool request_jacobian,
473  RealVectorX& f,
474  RealMatrixX& jac,
475  const unsigned int side,
476  MAST::BoundaryConditionBase& bc) = 0;
477 
478 
483  virtual bool
484  surface_traction_residual_shifted_boundary(bool request_jacobian,
485  RealVectorX& f,
486  RealMatrixX& jac,
487  const unsigned int side,
488  MAST::BoundaryConditionBase& bc) = 0;
489 
494  virtual bool
495  surface_traction_residual_shifted_boundary_sensitivity(const MAST::FunctionBase& p,
496  bool request_jacobian,
497  RealVectorX& f,
498  RealMatrixX& jac,
499  const unsigned int side,
500  MAST::BoundaryConditionBase& bc) = 0;
501 
505  virtual bool
506  surface_pressure_residual(bool request_jacobian,
507  RealVectorX& f,
508  RealMatrixX& jac,
509  const unsigned int side,
510  MAST::BoundaryConditionBase& bc) = 0;
511 
512 
518  virtual bool
519  surface_pressure_residual(bool request_jacobian,
520  RealVectorX& f,
521  RealMatrixX& jac,
522  MAST::BoundaryConditionBase& bc);
523 
529  virtual void
530  surface_pressure_boundary_velocity(const MAST::FunctionBase& p,
531  const unsigned int s,
532  const MAST::FieldFunction<RealVectorX>& vel_f,
533  MAST::BoundaryConditionBase& bc,
534  bool request_jacobian,
535  RealVectorX& f,
536  RealMatrixX& jac);
537 
543  virtual bool
544  linearized_surface_pressure_residual(bool request_jacobian,
545  RealVectorX& f,
546  RealMatrixX& jac,
547  MAST::BoundaryConditionBase& bc);
548 
549 
554  Real piston_theory_cp(const unsigned int order,
555  const Real vel_U,
556  const Real gamma,
557  const Real mach);
558 
559 
560 
565  Real piston_theory_dcp_dv(const unsigned int order,
566  const Real vel_U,
567  const Real gamma,
568  const Real mach);
569 
570 
578  virtual bool piston_theory_residual(bool request_jacobian,
579  RealVectorX &f,
580  RealMatrixX& jac_xdot,
581  RealMatrixX& jac,
582  MAST::BoundaryConditionBase& bc) = 0;
583 
590  virtual bool piston_theory_residual(bool request_jacobian,
591  RealVectorX &f,
592  RealMatrixX& jac_xdot,
593  RealMatrixX& jac,
594  const unsigned int side,
595  MAST::BoundaryConditionBase& bc) = 0;
596 
597 
598 
606  virtual bool piston_theory_residual_sensitivity(const MAST::FunctionBase& p,
607  bool request_jacobian,
608  RealVectorX &f,
609  RealMatrixX& jac_xdot,
610  RealMatrixX& jac,
611  MAST::BoundaryConditionBase& bc) = 0;
612 
619  virtual bool piston_theory_residual_sensitivity(const MAST::FunctionBase& p,
620  bool request_jacobian,
621  RealVectorX &f,
622  RealMatrixX& jac_xdot,
623  RealMatrixX& jac,
624  const unsigned int side,
625  MAST::BoundaryConditionBase& bc) = 0;
626 
627 
632  virtual bool
633  linearized_frequency_domain_surface_pressure_residual
634  (bool request_jacobian,
635  ComplexVectorX& f,
636  ComplexMatrixX& jac,
637  const unsigned int side,
638  MAST::BoundaryConditionBase& bc) = 0;
639 
640 
649  virtual bool
650  linearized_frequency_domain_surface_pressure_residual
651  (bool request_jacobian,
652  ComplexVectorX& f,
653  ComplexMatrixX& jac,
654  MAST::BoundaryConditionBase& bc);
655 
656 
661  virtual bool
662  linearized_frequency_domain_surface_pressure_residual_sensitivity
663  (const MAST::FunctionBase& p,
664  bool request_jacobian,
665  ComplexVectorX& f,
666  ComplexMatrixX& jac,
667  const unsigned int side,
668  MAST::BoundaryConditionBase& bc) = 0;
669 
670 
676  virtual bool
677  linearized_frequency_domain_surface_pressure_residual_sensitivity
678  (const MAST::FunctionBase& p,
679  bool request_jacobian,
680  ComplexVectorX& f,
681  ComplexMatrixX& jac,
682  MAST::BoundaryConditionBase& bc) {
683 
684  libmesh_error(); // to be implemented
685  }
686 
687 
691  virtual bool
692  surface_pressure_residual_sensitivity(const MAST::FunctionBase& p,
693  bool request_jacobian,
694  RealVectorX& f,
695  RealMatrixX& jac,
696  const unsigned int side,
697  MAST::BoundaryConditionBase& bc) = 0;
698 
703  virtual bool
704  surface_pressure_residual_sensitivity(const MAST::FunctionBase& p,
705  bool request_jacobian,
706  RealVectorX& f,
707  RealMatrixX& jac,
708  MAST::BoundaryConditionBase& bc);
709 
710 
711 
716  virtual bool thermal_residual(bool request_jacobian,
717  RealVectorX& f,
718  RealMatrixX& jac,
719  MAST::BoundaryConditionBase& bc) = 0;
720 
725  virtual bool thermal_residual_sensitivity(const MAST::FunctionBase& p,
726  bool request_jacobian,
727  RealVectorX& f,
728  RealMatrixX& jac,
729  MAST::BoundaryConditionBase& bc) = 0;
730 
735  virtual void thermal_residual_boundary_velocity(const MAST::FunctionBase& p,
736  const unsigned int s,
737  const MAST::FieldFunction<RealVectorX>& vel_f,
738  MAST::BoundaryConditionBase& bc,
739  bool request_jacobian,
740  RealVectorX& f,
741  RealMatrixX& jac) = 0;
742 
746  const MAST::ElementPropertyCardBase& _property;
747 
748 
752  RealVectorX _local_sol;
753 
754 
758  RealVectorX _local_delta_sol;
759 
760 
764  RealVectorX _local_sol_sens;
765 
766 
770  RealVectorX _local_delta_sol_sens;
771 
772 
776  RealVectorX _local_vel;
777 
778 
782  RealVectorX _local_delta_vel;
783 
784 
788  RealVectorX _local_vel_sens;
789 
790 
794  RealVectorX _local_delta_vel_sens;
795 
796 
800  RealVectorX _local_accel;
801 
802 
806  RealVectorX _local_delta_accel;
807 
808 
812  RealVectorX _local_accel_sens;
813 
814 
818  RealVectorX _local_delta_accel_sens;
819 
820 
824  RealVectorX* _incompatible_sol;
825 
826  };
827 
828 
832  std::unique_ptr<MAST::StructuralElementBase>
833  build_structural_element(MAST::SystemInitialization& sys,
834  const MAST::GeomElem& elem,
835  const MAST::ElementPropertyCardBase& p);
836 
837 }
838 
839 
840 #endif // __mast__structural_element_base__
MAST::StructuralElementBase::_local_accel_sens
RealVectorX _local_accel_sens
local acceleration sensitivity
Definition: structural_element_base.h:812
MAST::StructuralElementBase::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)=0
calculates the term on side s: .
MAST::StructuralElementBase::linearized_frequency_domain_volume_external_residual
bool linearized_frequency_domain_volume_external_residual(bool request_jacobian, ComplexVectorX &f, ComplexMatrixX &jac, std::multimap< libMesh::subdomain_id_type, MAST::BoundaryConditionBase *> &bc)
Calculates the frequency domain volume external force contribution to system residual.
Definition: structural_element_base.cpp:923
MAST::StructuralElementBase::surface_pressure_boundary_velocity
virtual void surface_pressure_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 force vector and Jacobian due to surface pressure applied on the entire element domain...
Definition: structural_element_base.cpp:1318
MAST::StructuralElementBase::thermal_residual
virtual bool thermal_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, MAST::BoundaryConditionBase &bc)=0
Calculates the force vector and Jacobian due to thermal stresses.
MAST::StructuralElementBase::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)=0
Calculates the force vector and Jacobian due to small perturbation surface pressure.
MAST::StructuralElementBase::thermal_residual_sensitivity
virtual bool thermal_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac, MAST::BoundaryConditionBase &bc)=0
Calculates the sensitivity of force vector and Jacobian due to thermal stresses.
MAST::StructuralElementBase::update_incompatible_mode_solution
virtual void update_incompatible_mode_solution(const RealVectorX &dsol)
updates the incompatible solution for this element.
Definition: structural_element_base.h:402
MAST::StressStrainOutputBase
Data structure provides the mechanism to store stress and strain output from a structural analysis...
Definition: stress_output_base.h:61
MAST::StructuralElementBase::_local_delta_vel_sens
RealVectorX _local_delta_vel_sens
local perturbed velocity sensitivity
Definition: structural_element_base.h:794
MAST::StructuralElementBase::inertial_residual_sensitivity
virtual bool inertial_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xddot, RealMatrixX &jac_xdot, RealMatrixX &jac)
sensitivity of the inertial force contribution to system residual
Definition: structural_element_base.cpp:403
MAST::StructuralElementBase::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)=0
Calculates the force vector and Jacobian due to piston-theory based surface pressure on the entire el...
MAST::StructuralElementBase::linearized_inertial_residual
virtual bool linearized_inertial_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xddot, RealMatrixX &jac_xdot, RealMatrixX &jac)
inertial force contribution to system residual of linerized problem
Definition: structural_element_base.cpp:360
MAST::StructuralElementBase::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)=0
Calculates the sensitivity of force vector and Jacobian due to thermal stresses.
MAST::StructuralElementBase::set_perturbed_acceleration
virtual void set_perturbed_acceleration(const RealVectorX &vec, bool if_sens=false)
stores vec as perturbed acceleration for element level calculations, or its sensitivity if if_sens is...
Definition: structural_element_base.cpp:197
MAST::StructuralElementBase::set_acceleration
virtual void set_acceleration(const RealVectorX &vec, bool if_sens=false)
stores vec as acceleration for element level calculations, or its sensitivity if if_sens is true...
Definition: structural_element_base.cpp:170
ComplexVectorX
Matrix< Complex, Dynamic, 1 > ComplexVectorX
Definition: mast_data_types.h:37
MAST::StructuralElementBase::_local_delta_sol_sens
RealVectorX _local_delta_sol_sens
local perturbed solution sensitivity
Definition: structural_element_base.h:770
MAST::StructuralElementBase::incompatible_mode_size
virtual unsigned int incompatible_mode_size() const
Definition: structural_element_base.h:379
MAST::StructuralElementBase::inertial_residual
virtual bool inertial_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xddot, RealMatrixX &jac_xdot, RealMatrixX &jac)
inertial force contribution to system residual
Definition: structural_element_base.cpp:257
MAST::StructuralElementBase::_local_delta_accel_sens
RealVectorX _local_delta_accel_sens
local perturbed acceleration sensitivity
Definition: structural_element_base.h:818
MAST::StructuralElementBase::damping_residual
virtual bool damping_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac)
damping force contribution to system residual
Definition: structural_element_base.h:164
MAST::StructuralElementBase::_local_accel
RealVectorX _local_accel
local acceleration
Definition: structural_element_base.h:800
MAST::StructuralElementBase::linearized_volume_external_residual
bool linearized_volume_external_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac, std::multimap< libMesh::subdomain_id_type, MAST::BoundaryConditionBase *> &bc)
volume external force contribution to system residual.
Definition: structural_element_base.cpp:862
elem_base.h
MAST::StructuralElementBase::internal_residual_sensitivity
virtual bool internal_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac)=0
sensitivity of the internal force contribution to system residual
Real
libMesh::Real Real
Definition: mast_data_types.h:32
MAST::StructuralElementBase::_local_sol_sens
RealVectorX _local_sol_sens
local solution sensitivity
Definition: structural_element_base.h:764
MAST::StructuralElementBase::if_incompatible_modes
virtual bool if_incompatible_modes() const =0
MAST::StructuralElementBase::transform_vector_to_local_system
void transform_vector_to_local_system(const ValType &global_vec, ValType &local_vec) const
Definition: structural_element_base.cpp:1608
MAST::StructuralElementBase::thermal_residual_temperature_derivative
virtual void thermal_residual_temperature_derivative(const MAST::FEBase &fe_thermal, RealMatrixX &m)=0
MAST::StructuralElementBase::piston_theory_cp
Real piston_theory_cp(const unsigned int order, const Real vel_U, const Real gamma, const Real mach)
Definition: structural_element_base.cpp:1711
MAST::StructuralElementBase::_local_delta_vel
RealVectorX _local_delta_vel
local perturbed velocity
Definition: structural_element_base.h:782
MAST::StructuralElementBase::set_velocity
virtual void set_velocity(const RealVectorX &vec, bool if_sens=false)
stores vec as velocity for element level calculations, or its sensitivity if if_sens is true...
Definition: structural_element_base.cpp:116
MAST::StructuralElementBase::_local_delta_sol
RealVectorX _local_delta_sol
local perturbed solution
Definition: structural_element_base.h:758
MAST::StructuralElementBase::follower_forces
bool follower_forces
flag for follower forces
Definition: structural_element_base.h:439
MAST::StructuralElementBase::elem_property
const MAST::ElementPropertyCardBase & elem_property()
returns a constant reference to the finite element object
Definition: structural_element_base.h:119
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::StructuralElementBase::prestress_residual_sensitivity
virtual bool prestress_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac)=0
sensitivity of the prestress force contribution to system residual
MAST::StructuralElementBase::piston_theory_dcp_dv
Real piston_theory_dcp_dv(const unsigned int order, const Real vel_U, const Real gamma, const Real mach)
Definition: structural_element_base.cpp:1742
MAST::StructuralElementBase::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)=0
Calculates the sensitivity of element vector and matrix quantities for surface traction boundary cond...
MAST::StructuralElementBase::side_external_residual
bool side_external_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac, std::multimap< libMesh::boundary_id_type, MAST::BoundaryConditionBase *> &bc)
side external force contribution to system residual.
Definition: structural_element_base.cpp:630
MAST::StructuralElementBase::set_solution
virtual void set_solution(const RealVectorX &vec, bool if_sens=false)
stores vec as solution for element level calculations, or its sensitivity if if_sens is true...
Definition: structural_element_base.cpp:57
MAST::StructuralElementBase::calculate_stress
virtual bool calculate_stress(bool request_derivative, const MAST::FunctionBase *f, MAST::StressStrainOutputBase &output)=0
Calculates the stress tensor.
RealVectorX
Matrix< Real, Dynamic, 1 > RealVectorX
Definition: mast_data_types.h:35
MAST::StructuralElementBase::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)=0
Calculates the sensitivity of force vector and Jacobian due to small is applicable for perturbation s...
MAST::StructuralElementBase::volume_external_residual
bool volume_external_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac, std::multimap< libMesh::subdomain_id_type, MAST::BoundaryConditionBase *> &bc)
volume external force contribution to system residual.
Definition: structural_element_base.cpp:805
MAST::build_structural_element
std::unique_ptr< MAST::StructuralElementBase > build_structural_element(MAST::SystemInitialization &sys, const MAST::GeomElem &elem, const MAST::ElementPropertyCardBase &p)
builds the structural element for the specified element type
Definition: structural_element_base.cpp:1681
MAST::StructuralElementBase::piston_theory_residual
virtual bool piston_theory_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac, MAST::BoundaryConditionBase &bc)=0
Calculates the force vector and Jacobian due to piston-theory based surface pressure on the entire el...
MAST::StructuralElementBase::_property
const MAST::ElementPropertyCardBase & _property
element property
Definition: structural_element_base.h:746
MAST::StructuralElementBase::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)=0
Calculates the sensitivity of force vector and Jacobian due to surface traction and sensitiity due to...
MAST::StructuralElementBase::StructuralElementBase
StructuralElementBase(MAST::SystemInitialization &sys, const MAST::GeomElem &elem, const MAST::ElementPropertyCardBase &p)
Constructor.
Definition: structural_element_base.cpp:38
MAST::StructuralElementBase::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)=0
Calculates the boundary velocity term contributions to the sensitivity of stress at the specified bou...
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::StructuralElementBase::set_perturbed_velocity
virtual void set_perturbed_velocity(const RealVectorX &vec, bool if_sens=false)
stores vec as perturbed velocity for element level calculations, or its sensitivity if if_sens is tru...
Definition: structural_element_base.cpp:143
MAST::StructuralElementBase::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)=0
Calculates the sensitivity of force vector and Jacobian due to surface traction and sensitiity due to...
MAST::StructuralElementBase::local_solution
const RealVectorX & local_solution(bool if_sens=false) const
Definition: structural_element_base.h:108
MAST::StructuralElementBase::_local_vel_sens
RealVectorX _local_vel_sens
local velocity sensitivity
Definition: structural_element_base.h:788
MAST::StructuralElementBase::linearized_surface_pressure_residual
virtual bool linearized_surface_pressure_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, MAST::BoundaryConditionBase &bc)
Calculates the force vector and Jacobian due to surface pressure applied on the entire element domain...
Definition: structural_element_base.cpp:1411
MAST::ElementBase::elem
const MAST::GeomElem & elem() const
Definition: elem_base.h:109
MAST::StructuralElementBase::calculate_stress_temperature_derivative
virtual void calculate_stress_temperature_derivative(MAST::FEBase &fe_thermal, MAST::StressStrainOutputBase &output)=0
MAST::StructuralElementBase::set_perturbed_solution
virtual void set_perturbed_solution(const RealVectorX &vec, bool if_sens=false)
stores vec as perturbed solution for element level calculations, or its sensitivity if if_sens is tru...
Definition: structural_element_base.cpp:86
MAST::StructuralElementBase::damping_residual_sensitivity
virtual bool damping_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac)
sensitivity of the damping force contribution to system residual
Definition: structural_element_base.h:308
MAST::StructuralElementBase::side_external_residual_sensitivity
bool side_external_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac, std::multimap< libMesh::boundary_id_type, MAST::BoundaryConditionBase *> &bc)
sensitivity of the side external force contribution to system residual
Definition: structural_element_base.cpp:972
MAST::StructuralElementBase::prestress_residual
virtual bool prestress_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac)=0
prestress force contribution to system residual
MAST::StructuralElementBase::linearized_side_external_residual
bool linearized_side_external_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac, std::multimap< libMesh::boundary_id_type, MAST::BoundaryConditionBase *> &bc)
side external force contribution to system residual.
Definition: structural_element_base.cpp:706
MAST::StructuralElementBase::internal_residual_jac_dot_state_sensitivity
virtual bool internal_residual_jac_dot_state_sensitivity(RealMatrixX &jac)=0
calculates d[J]/d{x} .
MAST::StructuralElementBase::transform_matrix_to_global_system
void transform_matrix_to_global_system(const ValType &local_mat, ValType &global_mat) const
Definition: structural_element_base.cpp:1568
MAST::StructuralElementBase::linearized_frequency_domain_side_external_residual
bool linearized_frequency_domain_side_external_residual(bool request_jacobian, ComplexVectorX &f, ComplexMatrixX &jac, std::multimap< libMesh::boundary_id_type, MAST::BoundaryConditionBase *> &bc)
Calculates the external force due to frequency domain side external force contribution to system resi...
Definition: structural_element_base.cpp:752
MAST::StructuralElementBase::linearized_internal_residual
virtual bool linearized_internal_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac)
internal force contribution to system residual of the linearized problem
Definition: structural_element_base.cpp:224
MAST::StructuralElementBase::_incompatible_sol
RealVectorX * _incompatible_sol
incompatible mode solution vector
Definition: structural_element_base.h:824
MAST::StructuralElementBase::_local_delta_accel
RealVectorX _local_delta_accel
local perturbed acceleration
Definition: structural_element_base.h:806
MAST::StructuralElementBase::_local_sol
RealVectorX _local_sol
local solution
Definition: structural_element_base.h:752
MAST::StructuralElementBase::surface_pressure_residual
virtual bool surface_pressure_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)=0
Calculates the force vector and Jacobian due to surface pressure.
MAST::StructuralElementBase::volume_external_residual_sensitivity
bool volume_external_residual_sensitivity(const MAST::FunctionBase &p, bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac, std::multimap< libMesh::subdomain_id_type, MAST::BoundaryConditionBase *> &bc)
sensitivity of the volume external force contribution to system residual
Definition: structural_element_base.cpp:1054
MAST::StructuralElementBase::_local_vel
RealVectorX _local_vel
local velocity
Definition: structural_element_base.h:776
MAST::StructuralElementBase::set_incompatible_mode_solution
void set_incompatible_mode_solution(RealVectorX &vec)
sets the pointer to the incompatible mode solution vector.
Definition: structural_element_base.h:392
MAST::StructuralElementBase::surface_traction_residual
virtual bool surface_traction_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int side, MAST::BoundaryConditionBase &bc)=0
Calculates the force vector and Jacobian due to surface traction.
MAST::StructuralElementBase::internal_residual
virtual bool internal_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac)=0
internal force contribution to system residual
MAST::StructuralElementBase::inertial_residual_boundary_velocity
virtual void inertial_residual_boundary_velocity(const MAST::FunctionBase &p, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f, bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xddot, RealMatrixX &jac_xdot, RealMatrixX &jac)
sensitivity of the inertial force contribution to system residual
Definition: structural_element_base.cpp:508
MAST::StructuralElementBase::volume_external_residual_boundary_velocity
void volume_external_residual_boundary_velocity(const MAST::FunctionBase &p, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f, std::multimap< libMesh::subdomain_id_type, MAST::BoundaryConditionBase *> &bc, bool request_jacobian, RealVectorX &f, RealMatrixX &jac)
boundary velocity contribution of volume external force.
Definition: structural_element_base.cpp:1113
MAST::StructuralElementBase::transform_vector_to_global_system
void transform_vector_to_global_system(const ValType &local_vec, ValType &global_vec) const
Definition: structural_element_base.cpp:1644
MAST::StructuralElementBase::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)=0
Calculates the force vector and Jacobian due to surface pressure.
MAST::ElementBase
This is the base class for elements that implement calculation of finite element quantities over the ...
Definition: elem_base.h:72