manolas
/
ReducedModelOptimization
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Testing the new axial and shear scenarios

This commit is contained in:
Iason 2021-03-02 20:32:43 +02:00
parent 4a03c435c2
commit ae636c3af5
3 changed files with 61 additions and 61 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
CMakeLists.txt
View File

@ -38,14 +38,14 @@ download_project(PROJ MYSOURCES
endif() endif()
##Polyscope ##Polyscope
#download_project(PROJ POLYSCOPE download_project(PROJ POLYSCOPE
# GIT_REPOSITORY https://github.com/nmwsharp/polyscope.git GIT_REPOSITORY https://github.com/nmwsharp/polyscope.git
# GIT_TAG master GIT_TAG master
# PREFIX ${CMAKE_CURRENT_SOURCE_DIR}/build/external/ PREFIX ${CMAKE_CURRENT_SOURCE_DIR}/build/external/
# ${UPDATE_DISCONNECTED_IF_AVAILABLE} ${UPDATE_DISCONNECTED_IF_AVAILABLE}
#) )
#add_subdirectory(${POLYSCOPE_SOURCE_DIR}) add_subdirectory(${POLYSCOPE_SOURCE_DIR})
#add_compile_definitions(POLYSCOPE_DEFINED) add_compile_definitions(POLYSCOPE_DEFINED)
#dlib #dlib
download_project(PROJ DLIB download_project(PROJ DLIB
@ -88,5 +88,5 @@ target_include_directories(${PROJECT_NAME}
) )
link_directories(${MYSOURCES_SOURCE_DIR}/boost_graph/libs) link_directories(${MYSOURCES_SOURCE_DIR}/boost_graph/libs)
#target_link_libraries(${PROJECT_NAME} polyscope Eigen3::Eigen matplot dlib::dlib) target_link_libraries(${PROJECT_NAME} polyscope Eigen3::Eigen matplot dlib::dlib)
target_link_libraries(${PROJECT_NAME} -static Eigen3::Eigen matplot dlib::dlib) #target_link_libraries(${PROJECT_NAME} -static Eigen3::Eigen matplot dlib::dlib)

4
src/main.cpp
View File

@ -52,8 +52,8 @@ int main(int argc, char *argv[]) {
settings_optimization.numberOfFunctionCalls = settings_optimization.numberOfFunctionCalls =
input_numberOfFunctionCallsDefined ? std::atoi(argv[3]) : 100; input_numberOfFunctionCallsDefined ? std::atoi(argv[3]) : 100;
settings_optimization.normalizationStrategy = settings_optimization.normalizationStrategy =
ReducedModelOptimizer::Settings::NormalizationStrategy::Epsilon; ReducedModelOptimizer::Settings::NormalizationStrategy::NonNormalized;
settings_optimization.normalizationParameter = 0.003; settings_optimization.normalizationParameter = 0.0003;
settings_optimization.solutionAccuracy = 0.01; settings_optimization.solutionAccuracy = 0.01;
// Optimize pair // Optimize pair

98
src/reducedmodeloptimizer.cpp
View File

@ -784,31 +784,31 @@ ReducedModelOptimizer::createScenarios(
//// Axial //// Axial
SimulationScenario scenarioName = SimulationScenario::Axial; SimulationScenario scenarioName = SimulationScenario::Axial;
// NewMethod // NewMethod
// for (auto viPairIt = m_fullPatternOppositeInterfaceViMap.begin(); for (auto viPairIt = m_fullPatternOppositeInterfaceViMap.begin();
// viPairIt != m_fullPatternOppositeInterfaceViMap.end(); viPairIt++) { viPairIt != m_fullPatternOppositeInterfaceViMap.end(); viPairIt++) {
// if (viPairIt != m_fullPatternOppositeInterfaceViMap.begin()) { if (viPairIt != m_fullPatternOppositeInterfaceViMap.begin()) {
// CoordType forceDirection(1, 0, 0); CoordType forceDirection(1, 0, 0);
// const auto viPair = *viPairIt; const auto viPair = *viPairIt;
// nodalForces[viPair.first] = nodalForces[viPair.first] =
// Vector6d({forceDirection[0], forceDirection[1], forceDirection[2], Vector6d({forceDirection[0], forceDirection[1], forceDirection[2], 0,
// 0, 0, 0}) *
// 0, 0}) * forceMagnitude * 10;
// forceMagnitude * 10; fixedVertices[viPair.second] =
// fixedVertices[viPair.second] = std::unordered_set<DoFType>{0, 1, 2, 3, 4, 5};
// std::unordered_set<DoFType>{0, 1, 2, 3, 4, 5}; }
// }
// }
// OldMethod
for (const auto &viPair : m_fullPatternOppositeInterfaceViMap) {
CoordType forceDirection =
(pMesh->vert[viPair.first].cP() - pMesh->vert[viPair.second].cP())
.Normalize();
nodalForces[viPair.first] = Vector6d({forceDirection[0], forceDirection[1],
forceDirection[2], 0, 0, 0}) *
forceMagnitude * 10;
fixedVertices[viPair.second] =
std::unordered_set<DoFType>{0, 1, 2, 3, 4, 5};
} }
// OldMethod
// for (const auto &viPair : m_fullPatternOppositeInterfaceViMap) {
// CoordType forceDirection =
// (pMesh->vert[viPair.first].cP() - pMesh->vert[viPair.second].cP())
// .Normalize();
// nodalForces[viPair.first] = Vector6d({forceDirection[0],
// forceDirection[1],
// forceDirection[2], 0, 0, 0}) *
// forceMagnitude * 10;
// fixedVertices[viPair.second] =
// std::unordered_set<DoFType>{0, 1, 2, 3, 4, 5};
// }
scenarios[scenarioName] = std::make_shared<SimulationJob>( scenarios[scenarioName] = std::make_shared<SimulationJob>(
SimulationJob(pMesh, simulationScenarioStrings[scenarioName], SimulationJob(pMesh, simulationScenarioStrings[scenarioName],
fixedVertices, nodalForces, {})); fixedVertices, nodalForces, {}));
@ -818,32 +818,32 @@ ReducedModelOptimizer::createScenarios(
fixedVertices.clear(); fixedVertices.clear();
nodalForces.clear(); nodalForces.clear();
// NewMethod // NewMethod
// for (auto viPairIt = m_fullPatternOppositeInterfaceViMap.begin(); for (auto viPairIt = m_fullPatternOppositeInterfaceViMap.begin();
// viPairIt != m_fullPatternOppositeInterfaceViMap.end(); viPairIt++) { viPairIt != m_fullPatternOppositeInterfaceViMap.end(); viPairIt++) {
// if (viPairIt != m_fullPatternOppositeInterfaceViMap.begin()) { if (viPairIt != m_fullPatternOppositeInterfaceViMap.begin()) {
// CoordType forceDirection(0, 1, 0); CoordType forceDirection(0, 1, 0);
// const auto viPair = *viPairIt; const auto viPair = *viPairIt;
// nodalForces[viPair.first] = nodalForces[viPair.first] =
// Vector6d({forceDirection[0], forceDirection[1], forceDirection[2], Vector6d({forceDirection[0], forceDirection[1], forceDirection[2], 0,
// 0, 0, 0}) *
// 0, 0}) * forceMagnitude * 10;
// forceMagnitude * 10; fixedVertices[viPair.second] =
// fixedVertices[viPair.second] = std::unordered_set<DoFType>{0, 1, 2, 3, 4, 5};
// std::unordered_set<DoFType>{0, 1, 2, 3, 4, 5}; }
// }
// }
// OldMethod
for (const auto &viPair : m_fullPatternOppositeInterfaceViMap) {
CoordType v =
(pMesh->vert[viPair.first].cP() - pMesh->vert[viPair.second].cP())
.Normalize();
CoordType forceDirection = (v ^ patternPlaneNormal).Normalize();
nodalForces[viPair.first] = Vector6d({forceDirection[0], forceDirection[1],
forceDirection[2], 0, 0, 0}) *
0.40 * forceMagnitude;
fixedVertices[viPair.second] =
std::unordered_set<DoFType>{0, 1, 2, 3, 4, 5};
} }
// OldMethod
// for (const auto &viPair : m_fullPatternOppositeInterfaceViMap) {
// CoordType v =
// (pMesh->vert[viPair.first].cP() - pMesh->vert[viPair.second].cP())
// .Normalize();
// CoordType forceDirection = (v ^ patternPlaneNormal).Normalize();
// nodalForces[viPair.first] = Vector6d({forceDirection[0],
// forceDirection[1],
// forceDirection[2], 0, 0, 0}) *
// 0.40 * forceMagnitude;
// fixedVertices[viPair.second] =
// std::unordered_set<DoFType>{0, 1, 2, 3, 4, 5};
// }
scenarios[scenarioName] = std::make_shared<SimulationJob>( scenarios[scenarioName] = std::make_shared<SimulationJob>(
SimulationJob(pMesh, simulationScenarioStrings[scenarioName], SimulationJob(pMesh, simulationScenarioStrings[scenarioName],
fixedVertices, nodalForces, {})); fixedVertices, nodalForces, {}));