Method update_3D

tools/Assembler/Materials/ElementWiseOrthotropicMaterial.cpp:90–130

}

void ElementWiseOrthotropicMaterial::update_3D() {
    const size_t num_voxels = m_material_mesh->get_num_voxels();
    const VectorF& young_x = m_material_mesh->get_attribute(m_young_field_names[0]);
    const VectorF& young_y = m_material_mesh->get_attribute(m_young_field_names[1]);
    const VectorF& young_z = m_material_mesh->get_attribute(m_young_field_names[2]);
    const VectorF& poisson_yz = m_material_mesh->get_attribute(m_poisson_field_names[0]);
    const VectorF& poisson_zy = m_material_mesh->get_attribute(m_poisson_field_names[1]);
    const VectorF& poisson_zx = m_material_mesh->get_attribute(m_poisson_field_names[2]);
    const VectorF& poisson_xz = m_material_mesh->get_attribute(m_poisson_field_names[3]);
    const VectorF& poisson_xy = m_material_mesh->get_attribute(m_poisson_field_names[4]);
    const VectorF& poisson_yx = m_material_mesh->get_attribute(m_poisson_field_names[5]);
    const VectorF& shear_yz = m_material_mesh->get_attribute(m_shear_field_names[0]);
    const VectorF& shear_zx = m_material_mesh->get_attribute(m_shear_field_names[1]);
    const VectorF& shear_xy = m_material_mesh->get_attribute(m_shear_field_names[2]);

    assert(young_x.size() == num_voxels);
    assert(young_y.size() == num_voxels);
    assert(young_z.size() == num_voxels);
    assert(poisson_yz.size() == num_voxels);
    assert(poisson_zy.size() == num_voxels);
    assert(poisson_zx.size() == num_voxels);
    assert(poisson_xz.size() == num_voxels);
    assert(poisson_xy.size() == num_voxels);
    assert(poisson_yx.size() == num_voxels);
    assert(shear_yz.size() == num_voxels);
    assert(shear_zx.size() == num_voxels);
    assert(shear_xy.size() == num_voxels);

    m_materials.resize(num_voxels);
    for (size_t i=0; i<num_voxels; i++) {
        Vector3F elem_young(young_x[i], young_y[i], young_z[i]);
        VectorF elem_poisson(6);
        elem_poisson << poisson_yz[i], poisson_zy[i],
                        poisson_zx[i], poisson_xz[i],
                        poisson_xy[i], poisson_yx[i];
        Vector3F elem_shear(shear_yz[i], shear_zx[i], shear_xy[i]);

        m_materials[i] = Material::create_orthotropic(m_density,
                elem_young, elem_poisson, elem_shear);
    }
}