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

Corrected error in G

This commit is contained in:
Iason 2021-01-14 16:46:26 +02:00
parent a71efffe14
commit ab1206be9b
4 changed files with 12 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
beam.hpp
View File

@ -31,13 +31,13 @@ struct CylindricalBeamDimensions {
}; };
struct ElementMaterial { struct ElementMaterial {
float G; // poisson's ratio float poissonsRatio; // poisson's ratio
float E; // ym in pascal float E; // ym in pascal
ElementMaterial(const float &poissonsRatio, const float &youngsModulus) ElementMaterial(const float &poissonsRatio, const float &youngsModulus)
: G(poissonsRatio), E(youngsModulus) { : poissonsRatio(poissonsRatio), E(youngsModulus) {
assert(poissonsRatio <= 0.5 && poissonsRatio >= -1); assert(poissonsRatio <= 0.5 && poissonsRatio >= -1);
} }
ElementMaterial() : G(0.3), E(200 * 1e9) {} ElementMaterial() : poissonsRatio(0.3), E(200 * 1e9) {}
}; };
#endif // BEAM_HPP #endif // BEAM_HPP

5
beamformfinder.cpp
View File

@ -756,9 +756,8 @@ double FormFinder::computeTotalPotentialEnergy() const {
element.properties.A / (2 * element.initialLength); element.properties.A / (2 * element.initialLength);
const double theta1Diff = element.rotationalDisplacements_jplus1.theta1 - const double theta1Diff = element.rotationalDisplacements_jplus1.theta1 -
element.rotationalDisplacements_j.theta1; element.rotationalDisplacements_j.theta1;
const double torsionalPE = element.properties.material.G * const double torsionalPE = element.properties.G * element.properties.J *
element.properties.J * pow(theta1Diff, 2) / pow(theta1Diff, 2) / (2 * element.initialLength);
(2 * element.initialLength);
const double &theta2_j = element.rotationalDisplacements_j.theta2; const double &theta2_j = element.rotationalDisplacements_j.theta2;
const double &theta2_jplus1 = element.rotationalDisplacements_jplus1.theta2; const double &theta2_jplus1 = element.rotationalDisplacements_jplus1.theta2;
const double firstBendingPE_inBracketsTerm = 4 * pow(theta2_j, 2) + const double firstBendingPE_inBracketsTerm = 4 * pow(theta2_j, 2) +

10
elementalmesh.cpp
View File

@ -419,6 +419,7 @@ double computeAngle(const VectorType &vector0, const VectorType &vector1,
void Element::Properties::computeMaterialProperties( void Element::Properties::computeMaterialProperties(
const ElementMaterial &material) { const ElementMaterial &material) {
this->material = material; this->material = material;
this->G = material.E / (2 * (1 + material.poissonsRatio));
} }
void Element::Properties::computeDimensionsProperties( void Element::Properties::computeDimensionsProperties(
@ -459,7 +460,7 @@ Element::Properties::Properties(const CrossSectionType &dimensions,
Element::Properties::Properties(const std::array<double, 6> &arr) { Element::Properties::Properties(const std::array<double, 6> &arr) {
assert(arr.size() == 6); assert(arr.size() == 6);
material.E = arr[0]; material.E = arr[0];
material.G = arr[1]; G = arr[1];
A = arr[2]; A = arr[2];
I2 = arr[3]; I2 = arr[3];
I3 = arr[4]; I3 = arr[4];
@ -467,14 +468,13 @@ Element::Properties::Properties(const std::array<double, 6> &arr) {
} }
std::array<double, 6> Element::Properties::toArray() const { std::array<double, 6> Element::Properties::toArray() const {
return std::array<double, 6>({material.E, material.G, A, return std::array<double, 6>(
static_cast<double>(I2), {material.E, G, A, static_cast<double>(I2), static_cast<double>(I3), J});
static_cast<double>(I3), J});
} }
void Element::updateRigidity() { void Element::updateRigidity() {
rigidity.axial = properties.material.E * properties.A / initialLength; rigidity.axial = properties.material.E * properties.A / initialLength;
rigidity.torsional = properties.material.G * properties.J / initialLength; rigidity.torsional = properties.G * properties.J / initialLength;
rigidity.firstBending = rigidity.firstBending =
2 * properties.material.E * properties.I2 / initialLength; 2 * properties.material.E * properties.I2 / initialLength;
rigidity.secondBending = rigidity.secondBending =

1
elementalmesh.hpp
View File

@ -56,6 +56,7 @@ struct Element {
struct Properties { struct Properties {
CrossSectionType dimensions; CrossSectionType dimensions;
ElementMaterial material; ElementMaterial material;
double G{0};
double A{0}; // cross sectional area double A{0}; // cross sectional area
double I2{0}; // second moment of inertia double I2{0}; // second moment of inertia
double I3{0}; // third moment of inertia double I3{0}; // third moment of inertia