Skip to content

Unify computation of deformation gradient in new solid elements. #70

New issue
New issue
Open
#1565
Open
Unify computation of deformation gradient in new solid elements.#70
#1565
Labels
priority: lowIssues with low prioritystatus: unspecifiedIssues where the problem definition or possible solutions are unspecifiedtype: enhancementA new feature or enhancement to be implemented
@amgebauer

Description

@amgebauer
amgebauer
opened on Nov 29, 2024

Copied from https://gitlab.lrz.de/baci/baci/-/issues/1253

Description

There are mainly two ways to compute the deformation gradient:

  • A: using displacements (used for all except hex8)
$$\boldsymbol{F}=\boldsymbol{I} + \boldsymbol{\hat{u}}^T \frac{\partial \boldsymbol{N}}{\partial \boldsymbol{X}}^T$$
  • B: using current coordinates (used for hex8)
$$\boldsymbol{F}=(\boldsymbol{\hat{X}}+\boldsymbol{\hat{u}})^T \frac{\partial \boldsymbol{N}}{\partial \boldsymbol{X}}^T$$

with

  • $\boldsymbol{\hat{X}}$: reference nodal coordinates (num_nodes x num_dim)
  • $\boldsymbol{\hat{u}}$: nodal displacements (num_nodes x num_dim)
  • $\frac{\partial \boldsymbol{N}}{\partial \boldsymbol{X}}$: Derivative of the shape functions w.r.t. the reference coordinate $\boldsymbol{X}$ (num_dim x num_nodes)

Currently, using A does not work for some contact tests with hex8, B does not work for some ssi/ssti with tet4. Let's try to figure out which one is the better in floating point arithmetics and then try to fix either the ssi/ssti tests or the hex8 contact tests and consistently use the same way.

A seems to be the preferred solution, so we should try to migrate to A.

Possible Solution and Definition of Done

  • Decide which one is better from the perspective of floating point arithmetics
  • Consistently use the same way (i.e., either fix ssi/ssti tests or contact tests)

Additional Information

Note, in the failing tests, the difference between the deformation gradients is less than 1e-14

Error in SSI/SSTI tests 
2392: WARNING: PROC 1 has 92 nodes with zero/neg. concentration values for species 0
2392: WARNING: PROC 2 has 2 nodes with zero/neg. concentration values for species 0
2392: WARNING: PROC 0 has 103 nodes with zero/neg. concentration values for species 0
2392: WARNING: intercalation fraction X = c/c_max is violating prescribed bounds of ocp calculation model. Calculated values might therefore not be reasonable!
2392: X: 0.00109595 lower bound is: 0.3 upper bound is: 0.9825
2392: 
2392: 
2392: 
2392: =========================================================================
2392: 
2392: 
2392: =========================================================================
2392: 
2392: 
2392: =========================================================================
2392: PROC 2 ERROR in /builds/baci/baci/src/scatra_ele/4C_scatra_ele_boundary_calc_elch_electrode_sti_thermo.cpp, line 445:
2392: Temperature is non-positive!
2392: ------------------
2392:  0# 0x0000555F9A223541 in /builds/baci/build/4C
2392:  1# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalcElchElectrode<(FourC::Core::FE::CellType)5, 3>::evaluate_s2_i_coupling(FourC::Core::Elements::FaceElement const*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  2# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalc<(FourC::Core::FE::CellType)5, 3>::evaluate_action(FourC::Core::Elements::FaceElement*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::ScaTra::BoundaryAction, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  3# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalcElch<(FourC::Core::FE::CellType)5, 3>::evaluate_action(FourC::Core::Elements::FaceElement*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::ScaTra::BoundaryAction, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  4# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalcElchElectrodeSTIThermo<(FourC::Core::FE::CellType)5, 3>::evaluate_action(FourC::Core::Elements::FaceElement*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::ScaTra::BoundaryAction, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  5# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalc<(FourC::Core::FE::CellType)5, 3>::evaluate(FourC::Core::Elements::FaceElement*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  6# FourC::Discret::ELEMENTS::TransportBoundary::evaluate(Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  7# FourC::Core::FE::Discretization::evaluate_condition(Teuchos::ParameterList&, FourC::Core::FE::AssembleStrategy&, std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, int) in /builds/baci/build/4C
2392:  8# FourC::Core::FE::Discretization::evaluate_condition(Teuchos::ParameterList&, Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<Epetra_Vector>, Teuchos::RCP<Epetra_Vector>, Teuchos::RCP<Epetra_Vector>, std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, int) in /builds/baci/build/4C
2392:  9# FourC::ScaTra::MeshtyingStrategyS2I::EvaluateMeshtying() in /builds/baci/build/4C
2392: 10# FourC::ScaTra::MeshtyingStrategyS2IElch::EvaluateMeshtying() in /builds/baci/build/4C
2392: 11# FourC::ScaTra::ScaTraTimIntImpl::evaluate_solution_depending_conditions(Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<Epetra_Vector>) in /builds/baci/build/4C
2392: 12# FourC::ScaTra::ScaTraTimIntElch::evaluate_solution_depending_conditions(Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<Epetra_Vector>) in /builds/baci/build/4C
2392: 13# FourC::ScaTra::ScaTraTimIntImpl::assemble_mat_and_rhs() in /builds/baci/build/4C
2392: 14# virtual thunk to FourC::ScaTra::ScaTraTimIntElch::assemble_mat_and_rhs() in /builds/baci/build/4C
2392: 15# FourC::ScaTra::ScaTraTimIntImpl::PrepareLinearSolve() in /builds/baci/build/4C
2392: 16# FourC::SSTI::SSTIMono::evaluate_subproblems() in /builds/baci/build/4C
2392: 17# FourC::SSTI::SSTIMono::newton_loop() in /builds/baci/build/4C
2392: 18# FourC::SSTI::SSTIMono::Timeloop() in /builds/baci/build/4C
2392: 19# FourC::ssti_drt() in /builds/baci/build/4C
2392: 20# ntacal() in /builds/baci/build/4C
2392: 21# ntam(int, char**) in /builds/baci/build/4C
2392: 22# main in /builds/baci/build/4C
2392: 23# __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
2392: 24# _start in /builds/baci/build/4C
2392: 
2392: =========================================================================
2392: 
2392: 
2392: --------------------------------------------------------------------------
2392: MPI_ABORT was invoked on rank 2 in communicator MPI_COMM_WORLD
2392: with errorcode 1.
2392: 
2392: NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
2392: You may or may not see output from other processes, depending on
2392: exactly when Open MPI kills them.
2392: --------------------------------------------------------------------------
2392: PROC 0 ERROR in /builds/baci/baci/src/scatra_ele/4C_scatra_ele_boundary_calc_elch_electrode_sti_thermo.cpp, line 445:
2392: Temperature is non-positive!
2392: ------------------
2392:  0# 0x0000562C361A2541 in /builds/baci/build/4C
2392:  1# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalcElchElectrode<(FourC::Core::FE::CellType)5, 3>::evaluate_s2_i_coupling(FourC::Core::Elements::FaceElement const*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  2# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalc<(FourC::Core::FE::CellType)5, 3>::evaluate_action(FourC::Core::Elements::FaceElement*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::ScaTra::BoundaryAction, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  3# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalcElch<(FourC::Core::FE::CellType)5, 3>::evaluate_action(FourC::Core::Elements::FaceElement*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::ScaTra::BoundaryAction, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  4# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalcElchElectrodeSTIThermo<(FourC::Core::FE::CellType)5, 3>::evaluate_action(FourC::Core::Elements::FaceElement*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::ScaTra::BoundaryAction, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  5# FourC::Discret::ELEMENTS::ScaTraEleBoundaryCalc<(FourC::Core::FE::CellType)5, 3>::evaluate(FourC::Core::Elements::FaceElement*, Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  6# FourC::Discret::ELEMENTS::TransportBoundary::evaluate(Teuchos::ParameterList&, FourC::Core::FE::Discretization&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
2392:  7# FourC::Core::FE::Discretization::evaluate_condition(Teuchos::ParameterList&, FourC::Core::FE::AssembleStrategy&, std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, int) in /builds/baci/build/4C
2392:  8# FourC::Core::FE::Discretization::evaluate_condition(Teuchos::ParameterList&, Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<Epetra_Vector>, Teuchos::RCP<Epetra_Vector>, Teuchos::RCP<Epetra_Vector>, std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, int) in /builds/baci/build/4C
2392:  9# FourC::ScaTra::MeshtyingStrategyS2I::EvaluateMeshtying() in /builds/baci/build/4C
2392: 10# FourC::ScaTra::MeshtyingStrategyS2IElch::EvaluateMeshtying() in /builds/baci/build/4C
2392: 11# FourC::ScaTra::ScaTraTimIntImpl::evaluate_solution_depending_conditions(Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<Epetra_Vector>) in /builds/baci/build/4C
2392: 12# FourC::ScaTra::ScaTraTimIntElch::evaluate_solution_depending_conditions(Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<Epetra_Vector>) in /builds/baci/build/4C
2392: 13# FourC::ScaTra::ScaTraTimIntImpl::assemble_mat_and_rhs() in /builds/baci/build/4C
2392: 14# virtual thunk to FourC::ScaTra::ScaTraTimIntElch::assemble_mat_and_rhs() in /builds/baci/build/4C
2392: 15# FourC::ScaTra::ScaTraTimIntImpl::PrepareLinearSolve() in /builds/baci/build/4C
2392: 16# FourC::SSTI::SSTIMono::evaluate_subproblems() in /builds/baci/build/4C
2392: 17# FourC::SSTI::SSTIMono::newton_loop() in /builds/baci/build/4C
2392: 18# FourC::SSTI::SSTIMono::Timeloop() in /builds/baci/build/4C
2392: 19# FourC::ssti_drt() in /builds/baci/build/4C
2392: 20# ntacal() in /builds/baci/build/4C
2392: 21# ntam(int, char**) in /builds/baci/build/4C
2392: 22# main in /builds/baci/build/4C
2392: 23# __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
2392: 24# _start in /builds/baci/build/4C
2392: 
2392: =========================================================================
2392: 
2392: 
 235/2973 Test #2392: ssti_mono_3D_tet4_batt_with_anode_plate_and_current_collectors_elch_s2i_butlervolmer_BGS_AMG_11x11-p3 ..................................................................................***Failed    3.87
Error in contact tests 
440: Time for parallel redistribution...............6.810473e-03 secs
440: 
440: === Structural predictor: TangDis ===
440:  #It      ||F||_2         step    ||dx||_2
440: ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
440:    0  2.40129e+01  0.00000e+00  0.00000e+00
440:    1  7.73444e+08  1.00000e+00  9.13631e+00
440:    2  2.12544e+04  1.00000e+00  5.08254e+02
440:    3  7.85330e+07  1.00000e+00  5.19853e+02
440:    4  4.07291e+07  1.00000e+00  1.35340e+02
440:    5  6.45159e+10  1.00000e+00  7.09519e+03
440:    6  8.08302e+08  1.00000e+00  8.69232e+03
440:    7  8.76255e+06  1.00000e+00  6.24933e+03
440:    8  2.85326e+08  1.00000e+00  2.79859e+03
440:    9  1.14312e+08  1.00000e+00  8.59438e+02
440:   10  9.96095e+08  1.00000e+00  1.43515e+03
440:   11  4.01394e+08  1.00000e+00  1.55589e+03
440:   12  1.60703e+08  1.00000e+00  6.26128e+02
440:   13  4.22156e+09  1.00000e+00  1.69728e+02
440:   14  1.68442e+09  1.00000e+00  6.03183e+03
440:   15  6.72191e+08  1.00000e+00  5.98337e+03
440:   16  2.69234e+08  1.00000e+00  6.44689e+02
440:   17  1.07735e+08  1.00000e+00  3.30567e+02
440:   18  5.58101e+08  1.00000e+00  3.09539e+02
440:   19  1.36285e+10  1.00000e+00  5.00240e+02
440: 
440: 
440: =========================================================================
440: PROC 1 ERROR in /builds/baci/baci/apps/global_full/4C_global_full_main.cpp, line 195:
440: 4C produced a FE_OVERFLOW floating point exception.
440: ------------------
440:  0# (anonymous namespace)::sigfpe_handler(int) in /builds/baci/build/4C
440:  1# 0x00007F5E55BF2420 in /lib/x86_64-linux-gnu/libpthread.so.0
440:  2# FourC::Core::LinAlg::Voigt::VoigtUtils<(FourC::Core::LinAlg::Voigt::NotationType)1>::inverse_tensor(FourC::Core::LinAlg::Matrix<6u, 1u, double> const&, FourC::Core::LinAlg::Matrix<6u, 1u, double>&) in /builds/baci/build/4C
440:  3# FourC::Mat::ElastHyperEvaluate(FourC::Core::LinAlg::Matrix<3u, 3u, double> const&, FourC::Core::LinAlg::Matrix<6u, 1u, double> const&, Teuchos::ParameterList&, FourC::Core::LinAlg::Matrix<6u, 1u, double>&, FourC::Core::LinAlg::Matrix<6u, 6u, double>&, int, int, std::vector<Teuchos::RCP<FourC::Mat::Elastic::Summand>, std::allocator<Teuchos::RCP<FourC::Mat::Elastic::Summand> > > const&, FourC::Mat::SummandProperties const&, bool) in /builds/baci/build/4C
440:  4# FourC::Mat::ElastHyper::evaluate(FourC::Core::LinAlg::Matrix<3u, 3u, double> const*, FourC::Core::LinAlg::Matrix<6u, 1u, double> const*, Teuchos::ParameterList&, FourC::Core::LinAlg::Matrix<6u, 1u, double>*, FourC::Core::LinAlg::Matrix<6u, 6u, double>*, int, int) in /builds/baci/build/4C
440:  5# FourC::Discret::ELEMENTS::SolidEleCalc<(FourC::Core::FE::CellType)7, FourC::Discret::ELEMENTS::DisplacementBasedFormulation<(FourC::Core::FE::CellType)7> >::evaluate_nonlinear_force_stiffness_mass(FourC::Core::Elements::Element const&, FourC::Mat::So3Material&, FourC::Core::FE::Discretization const&, std::vector<int, std::allocator<int> > const&, Teuchos::ParameterList&, FourC::Core::LinAlg::SerialDenseVector*, FourC::Core::LinAlg::SerialDenseMatrix*, FourC::Core::LinAlg::SerialDenseMatrix*) in /builds/baci/build/4C
440:  6# _ZNSt8__detail9__variant17__gen_vtable_implILb1ENS0_12_Multi_arrayIPFvOZN5FourC7Discret8ELEMENTS5Solid8evaluateERN7Teuchos13ParameterListERNS3_4Core2FE14DiscretizationERSt6vectorIiSaIiEERNSA_6LinAlg17SerialDenseMatrixESK_RNSI_17SerialDenseVectorESM_SM_EUlRT_E2_RSt7variantIJNS5_7Details11VariantItemINS5_12SolidEleCalcILNSB_8CellTypeE7ENS5_28DisplacementBasedFormulationILSV_7EEEEEEENST_INSU_ILSV_9ENSW_ILSV_9EEEEEEENST_INSU_ILSV_10ENSW_ILSV_10EEEEEEENST_INSU_ILSV_11ENSW_ILSV_11EEEEEEENST_INSU_ILSV_28ENSW_ILSV_28EEEEEEENST_INSU_ILSV_12ENSW_ILSV_12EEEEEEENST_INSU_ILSV_13ENSW_ILSV_13EEEEEEENST_INSU_ILSV_14ENSW_ILSV_14EEEEEEENST_INSU_ILSV_16ENSW_ILSV_16EEEEEEENST_INSU_ILSV_7ENS5_44DisplacementBasedLinearKinematicsFormulationILSV_7EEEEEEENST_INSU_ILSV_9ENS1O_ILSV_9EEEEEEENST_INSU_ILSV_10ENS1O_ILSV_10EEEEEEENST_INSU_ILSV_11ENS1O_ILSV_11EEEEEEENST_INSU_ILSV_28ENS1O_ILSV_28EEEEEEENST_INSU_ILSV_12ENS1O_ILSV_12EEEEEEENST_INSU_ILSV_13ENS1O_ILSV_13EEEEEEENST_INSU_ILSV_14ENS1O_ILSV_14EEEEEEENST_INSU_ILSV_16ENS1O_ILSV_16EEEEEEENST_INSU_ILSV_7ENS5_15FBarFormulationILSV_7EEEEEEENST_INSU_ILSV_16ENS2G_ILSV_16EEEEEEENST_INS5_15SolidEleCalcEasILSV_7ELNS3_3STR8ELEMENTS7EasTypeE1ELNS3_5Inpar3STR9KinemTypeE2EEEEENST_INS2N_ILSV_7ELS2Q_2ELS2T_2EEEEENST_INS2N_ILSV_7ELS2Q_1ELS2T_1EEEEENST_INS2N_ILSV_7ELS2Q_2ELS2T_1EEEEENST_INSU_ILSV_7ENS5_15MulfFormulationILSV_7EEEEEEENST_INSU_ILSV_9ENS32_ILSV_9EEEEEEENST_INSU_ILSV_10ENS32_ILSV_10EEEEEEENST_INSU_ILSV_11ENS32_ILSV_11EEEEEEENST_INSU_ILSV_28ENS32_ILSV_28EEEEEEENST_INSU_ILSV_12ENS32_ILSV_12EEEEEEENST_INSU_ILSV_13ENS32_ILSV_13EEEEEEENST_INSU_ILSV_14ENS32_ILSV_14EEEEEEENST_INSU_ILSV_16ENS32_ILSV_16EEEEEEENST_INSU_ILSV_7ENS5_19MulfFBarFormulationILSV_7EEEEEEENST_INSU_ILSV_16ENS3U_ILSV_16EEEEEEEEEEJEEESt5tupleIJS42_EESt16integer_sequenceImJLm0EEEE14__visit_invokeESQ_S42_ in /builds/baci/build/4C
440:  7# FourC::Discret::ELEMENTS::Solid::evaluate(Teuchos::ParameterList&, FourC::Core::FE::Discretization&, std::vector<int, std::allocator<int> >&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
440:  8# std::_Function_handler<void (FourC::Core::Elements::Element&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&), FourC::Core::FE::Discretization::evaluate(Teuchos::ParameterList&, FourC::Core::FE::AssembleStrategy&)::{lambda(FourC::Core::Elements::Element&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&)#1}>::_M_invoke(std::_Any_data const&, FourC::Core::Elements::Element&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&) in /builds/baci/build/4C
440:  9# FourC::Core::FE::Discretization::evaluate(Teuchos::ParameterList&, FourC::Core::FE::AssembleStrategy&, std::function<void (FourC::Core::Elements::Element&, FourC::Core::Elements::Element::LocationArray&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseMatrix&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&, FourC::Core::LinAlg::SerialDenseVector&)> const&) in /builds/baci/build/4C
440: 10# FourC::Core::FE::Discretization::evaluate(Teuchos::ParameterList&, FourC::Core::FE::AssembleStrategy&) in /builds/baci/build/4C
440: 11# FourC::Core::FE::Discretization::evaluate(Teuchos::ParameterList&, Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>, Teuchos::RCP<Epetra_Vector>, Teuchos::RCP<Epetra_Vector>, Teuchos::RCP<Epetra_Vector>) in /builds/baci/build/4C
440: 12# FourC::STR::MODELEVALUATOR::Structure::evaluate_internal(Teuchos::ParameterList&, Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>*, Teuchos::RCP<Epetra_Vector>*) in /builds/baci/build/4C
440: 13# FourC::STR::MODELEVALUATOR::Structure::evaluate_internal(Teuchos::RCP<FourC::Core::LinAlg::SparseOperator>*, Teuchos::RCP<Epetra_Vector>*) in /builds/baci/build/4C
440: 14# FourC::STR::MODELEVALUATOR::Structure::apply_force_internal() in /builds/baci/build/4C
440: 15# FourC::STR::ModelEvaluator::apply_force(Epetra_Vector const&, Epetra_Vector&, double const&) const in /builds/baci/build/4C
440: 16# FourC::STR::IMPLICIT::Statics::apply_force(Epetra_Vector const&, Epetra_Vector&) in /builds/baci/build/4C
440: 17# FourC::STR::TimeInt::NoxInterface::computeF(Epetra_Vector const&, Epetra_Vector&, NOX::Epetra::Interface::Required::FillType) in /builds/baci/build/4C
440: 18# FourC::NOX::Nln::Group::computeF() in /builds/baci/build/4C
440: 19# NOX::Solver::LineSearchBased::step() in /builds/baci/build/4C
440: 20# NOX::Solver::LineSearchBased::solve() in /builds/baci/build/4C
440: 21# FourC::STR::Nln::SOLVER::Nox::Solve() in /builds/baci/build/4C
440: 22# FourC::STR::TimeInt::Implicit::Solve() in /builds/baci/build/4C
440: 23# FourC::Adapter::StructureTimeLoop::Integrate() in /builds/baci/build/4C
440: 24# FourC::dyn_nlnstructural_drt() in /builds/baci/build/4C
440: 25# FourC::caldyn_drt() in /builds/baci/build/4C
440: 26# ntacal() in /builds/baci/build/4C
440: 27# ntam(int, char**) in /builds/baci/build/4C
440: 28# main in /builds/baci/build/4C
440: 29# __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
440: 30# _start in /builds/baci/build/4C
440: 
440: =========================================================================
440: 
440: 
440: --------------------------------------------------------------------------
440: MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
440: with errorcode 1.
440: 
440: NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
440: You may or may not see output from other processes, depending on
440: exactly when Open MPI kills them.
440: --------------------------------------------------------------------------
 119/2973 Test  #440: contact3D_nitsche_hex8_sym_coulomb_new_struc-p2 ........................................................................................................................................***Failed    2.31 sec

Interested Parties

Metadata

Metadata

Assignees

No one assigned

    Labels

    priority: lowIssues with low prioritystatus: unspecifiedIssues where the problem definition or possible solutions are unspecifiedtype: enhancementA new feature or enhancement to be implemented

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions