beam_bending_thermal_stress_evaluation.cpp

Go to the documentation of this file.

/*
 * MAST: Multidisciplinary-design Adaptation and Sensitivity Toolkit
 * Copyright (C) 2013-2020  Manav Bhatia and MAST authors
 *
 * This library is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * This library is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with this library; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */


// BOOST includes
#include <boost/test/unit_test.hpp>


// MAST includes
#include "examples/structural/beam_bending_thermal_stress_with_offset/beam_bending_thermal_stress.h"
#include "tests/base/check_sensitivity.h"



BOOST_FIXTURE_TEST_SUITE  (Structural1DBeamBending,
                           MAST::BeamBendingThermalStress)

//BOOST_AUTO_TEST_CASE   (BeamBendingWithThermalStressSolution) {
//    
//    this->solve();
//    
//    // check the solution
//    // iterate over each node, and compare the nodal solution with the
//    // expected anlaytical value
//    unsigned int
//    dof_num = 0;
//    libMesh::MeshBase::const_node_iterator
//    it     =  _mesh->local_nodes_begin(),
//    end    =  _mesh->local_nodes_end();
//    
//    Real
//    temp        = (*_temp)(),
//    th_y        = (*_thy)(),
//    x           = 0.,
//    xi          = 0.,
//    eta         = 0.,
//    Eval        = (*_E)(),
//    alpha       = (*_alpha)(),
//    analytical  = 0.,
//    numerical   = 0.;
//    
//    // analytical solution to the simply supported problem is
//    // w(x)    = -alpha*T*A*th_y/2/Izz*(x^2/2 - L*x/2);
//    // dwdx(x) = -alpha*T*A*th_y/2/Izz*(x - L/2);
//    
//    BOOST_TEST_MESSAGE("  ** v-displacement and theta-z rotation **");
//    for ( ; it!=end; it++) {
//        const libMesh::Node* node = *it;
//        x            = (*node)(0);
//        
//        // v-displacement
//        analytical   = -alpha*temp/th_y*3./2.*(pow(x,2)/2. - _length*x/2.);
//        
//        dof_num      = node->dof_number(_sys->number(),
//                                        _structural_sys->vars()[1], // v-displ.
//                                        0);
//        numerical    =   _sys->solution->el(dof_num);
//        BOOST_CHECK(MAST::compare_value(analytical, numerical, tol));
//        
//        // theta-z rotation
//        analytical   = -alpha*temp/th_y*3./2.*(x - _length/2.);
//        
//        dof_num      = node->dof_number(_sys->number(),
//                                        _structural_sys->vars()[5], // tz-rotation
//                                        0);
//        numerical    =   _sys->solution->el(dof_num);
//        BOOST_CHECK(MAST::compare_value(analytical, numerical, tol));
//        
//    }
//    
//    
//    // make sure that each stress object has a single stored value
//    for (unsigned int i=0; i<_outputs.size(); i++) {
//        BOOST_CHECK(_outputs[i]->n_elem_in_storage() == 1);
//    }
//    
//    // now check the stress value in each element, which should be the same as
//    // the pressure value specified for the problem
//    BOOST_TEST_MESSAGE("  ** Stress **");
//    for (unsigned int i=0; i<_outputs.size(); i++) {
//        
//        // get the element and the nodes to evaluate the stress
//        const libMesh::Elem& e  = **(_outputs[i]->get_elem_subset().begin());
//        
//        const std::vector<MAST::StressStrainOutputBase::Data*>&
//        data = _outputs[i]->get_stress_strain_data_for_elem(&e);
//        
//        // find the location of quadrature point
//        for (unsigned int j=0; j<data.size(); j++) {
//            
//            // logitudinal strain for this location
//            numerical = data[j]->stress()(0);
//            
//            xi   = data[j]->point_location_in_element_coordinate()(0);
//            eta  = data[j]->point_location_in_element_coordinate()(1);
//            
//            // assuming linear Lagrange interpolation for elements
//            x =  e.point(0)(0) * (1.-xi)/2. +  e.point(1)(0) * (1.+xi)/2.;
//            // stress is a combination of the bending and compressive stress.
//            analytical   = Eval*alpha*temp*(0.75*(eta+1.) - 1.);
//            
//            
//            BOOST_CHECK(MAST::compare_value(analytical, numerical, tol));
//        }
//    }
//}


BOOST_AUTO_TEST_CASE   (BeamBendingWithThermalStressSensitivity) {
    
    MAST::check_sensitivity(*this);
    
}


BOOST_AUTO_TEST_SUITE_END()