71 lines
2.8 KiB
C++
71 lines
2.8 KiB
C++
/***************************************************************************
|
|
* Copyright (c) 2023 Ondsel, Inc. *
|
|
* *
|
|
* This file is part of OndselSolver. *
|
|
* *
|
|
* See LICENSE file for details about copyright. *
|
|
***************************************************************************/
|
|
|
|
/*********************************************************************
|
|
* @file MbDCode.cpp
|
|
*
|
|
* @brief Program to assemble a piston crank system.
|
|
*********************************************************************/
|
|
|
|
#include <filesystem>
|
|
#include "../OndselSolver/CADSystem.h"
|
|
#include "../OndselSolver/CREATE.h"
|
|
#include "../OndselSolver/GESpMatParPvPrecise.h"
|
|
#include "../OndselSolver/ASMTAssembly.h"
|
|
#include "../OndselSolver/MBDynSystem.h"
|
|
#include "../OndselSolver/MomentOfInertiaSolver.h"
|
|
|
|
using namespace MbD;
|
|
void runSpMat();
|
|
|
|
int main()
|
|
{
|
|
////ASMTAssembly::readWriteFile("piston.asmt");
|
|
//MBDynSystem::runFile("MBDynCaseJose.mbd");
|
|
//MBDynSystem::runFile("MBDynCase.mbd");
|
|
//MBDynSystem::runFile("crank_slider.mbd");
|
|
////ASMTAssembly::runSinglePendulumSuperSimplified(); //Mass is missing
|
|
////ASMTAssembly::runSinglePendulumSuperSimplified2(); //DOF has infinite acceleration due to zero mass and inertias
|
|
ASMTAssembly::runSinglePendulumSimplified();
|
|
ASMTAssembly::runSinglePendulum();
|
|
ASMTAssembly::runFile("../testapp/piston.asmt");
|
|
ASMTAssembly::runFile("../testapp/00backhoe.asmt");
|
|
//ASMTAssembly::runFile("circular.asmt"); //Needs checking
|
|
//ASMTAssembly::runFile("cirpendu.asmt"); //Under constrained. Testing ICKine.
|
|
//ASMTAssembly::runFile("engine1.asmt"); //Needs checking
|
|
ASMTAssembly::runFile("../testapp/fourbar.asmt");
|
|
//ASMTAssembly::runFile("fourbot.asmt"); //Very large but works
|
|
ASMTAssembly::runFile("../testapp/wobpump.asmt");
|
|
|
|
auto cadSystem = std::make_shared<CADSystem>();
|
|
cadSystem->runOndselSinglePendulum();
|
|
cadSystem->runOndselDoublePendulum();
|
|
//cadSystem->runOndselPiston(); //For debugging
|
|
cadSystem->runPiston();
|
|
runSpMat();
|
|
MomentOfInertiaSolver::example1();
|
|
}
|
|
|
|
void runSpMat() {
|
|
auto spMat = std::make_shared<SparseMatrix<double>>(3, 3);
|
|
spMat->atijput(0, 0, 1.0);
|
|
spMat->atijput(0, 1, 1.0);
|
|
spMat->atijput(1, 0, 1.0);
|
|
spMat->atijput(1, 1, 1.0);
|
|
spMat->atijput(1, 2, 1.0);
|
|
spMat->atijput(2, 1, 1.0);
|
|
spMat->atijput(2, 2, 1.0);
|
|
auto fullCol = std::make_shared<FullColumn<double>>(3);
|
|
fullCol->atiput(0, 1.0);
|
|
fullCol->atiput(1, 2.0);
|
|
fullCol->atiput(2, 3.0);
|
|
auto matSolver = CREATE<GESpMatParPvPrecise>::With();
|
|
auto answer = matSolver->solvewithsaveOriginal(spMat, fullCol, true);
|
|
auto aAx = spMat->timesFullColumn(answer);
|
|
}
|