public void SolveNewmark()
{
List <double> aConstants = CalculateIntegrationConstantsNewmark();
double[,] hatStiffnessMatrixNewmark = CalculateHatKMatrixNewmark(aConstants);
explicitSolution.Add(0, InitialValues.InitialDisplacementVector);
explicitVelocity.Add(0, InitialValues.InitialVelocityVector);
explicitAcceleration.Add(0, CalculateInitialAccelerationsNewmark());
TimeAtEachStep.Add(0, 0.0);
double[] nextSolution;
for (int i = 1; i < timeStepsNumber; i++)
{
double time = i * timeStep + InitialValues.InitialTime;
if (ActivateNonLinearSolution == false)
{
double[] hatRVectorNewmark = CalculateHatRVectorNewmark(i, aConstants);
nextSolution = LinearSolver.Solve(hatStiffnessMatrixNewmark, hatRVectorNewmark);
Console.WriteLine("Solution for Load Step {0} is:", i);
VectorOperations.PrintVector(nextSolution);
}
else
{
nextSolution = NewtonIterationsNewmark(ExternalForcesVector, i, aConstants);
Console.WriteLine("Solution for Load Step {0} is:", i);
VectorOperations.PrintVector(nextSolution);
}
explicitSolution.Add(i, nextSolution);
explicitAcceleration.Add(i, CalculateAccelerationNewmark(i, aConstants));
explicitVelocity.Add(i, CalculateVelocityNewmark(i, aConstants));
TimeAtEachStep.Add(i, time);
}
ExportToFile.ExportExplicitResults(explicitSolution, TimeAtEachStep, 1, 1);
//ShowToGUI.ShowResults(explicitSolution, TimeAtEachStep, 1, 1);
}
public void SolveExplicit()
{
double[,] hatMassMatrix = CalculateHatMMatrix();
explicitSolution.Add(-1, CalculatePreviousDisplacementVector());
explicitSolution.Add(0, InitialValues.InitialDisplacementVector);
explicitAcceleration.Add(0, CalculateInitialAccelerations());
TimeAtEachStep.Add(-1, -1 * timeStep + InitialValues.InitialTime);
TimeAtEachStep.Add(0, 0.0);
double[] nextSolution;
for (int i = 1; i < timeStepsNumber; i++)
{
double time = i * timeStep + InitialValues.InitialTime;
if (ActivateNonLinearSolution == false)
{
double[] hatRVector = CalculateHatRVector(i);
nextSolution = LinearSolver.Solve(hatMassMatrix, hatRVector);
Console.WriteLine("Solution for Load Step {0} is:", i);
VectorOperations.PrintVector(nextSolution);
}
else
{
double[] hatRVector = CalculateHatRVectorNL(i);
nextSolution = LinearSolver.Solve(hatMassMatrix, hatRVector);
//nextSolution = NewtonIterations(hatRVector);
Console.WriteLine("Solution for Load Step {0} is:", i);
VectorOperations.PrintVector(nextSolution);
}
explicitSolution.Add(i, nextSolution);
explicitAcceleration.Add(i, CalculateAccelerations());
TimeAtEachStep.Add(i, time);
}
ExportToFile.ExportExplicitResults(explicitSolution, TimeAtEachStep, 1, 1);
}