void GetEIF____00_21_UV(CMLoad_21 Load, float fL, out float fA, out float fB, out float fMa, out float fMb)
{
// Temporary
fA = 0.0f;
fB = 0.0f;
fMa = 0.0f;
fMb = 0.0f;
}
// Uniform Load
void GetEIF_0__00_21_UV(CMLoad_21 Load, float fL, out float fA, out float fB, out float fMa, out float fMb)
{
// !!! Signs
fA = 3/8f * Load.Fq * fL;
fB = 5/8f*Load.Fq * fL;
fMa = 0f;
fMb = Load.Fq * fL * fL / 8f;
}
// Uniform Load
void GetEIF_0__0__21_UV(CMLoad_21 Load, float fL, out float fA, out float fB, out float fMa, out float fMb)
{
// Whole Member
fA = 0.5f * Load.Fq * fL;
fB = fA;
fMa = fMb = 0f;
// IF distribution formulas
float fx = 0; // Temp
float fx_M_max = fL / 2.0f;
float fM_max = Load.Fq * MathF.Pow2(fL);
float fM_x = Load.Fq * fx / 2.0f * (fL - fx);
// Deflection formulas
float fE = 0.0f, fI = 0.0f;
float fx_Delta_max = fL / 2.0f;
float fDelta_max = 5.0f * Load.Fq * MathF.Pow4(fL) / (384.0f * fE * fI);
}
// Uniform Load
void GetEIF_00_00_21_UV(CMLoad_21 Load, float fL, out float fA, out float fB, out float fMa, out float fMb)
{
fA = fB = Load.Fq * fL / 2f;
fMa = -Load.Fq * fL * fL / 12f;
fMb = -fMa; /*Load.Fq * fL * fL / 12f;*/
}
void GetMLoadPart_21(CMLoad_21 Load, CE_1D_BASE Member, FEM_CALC_BASE.Enums.EElemSuppType2D eMType, out float fA, out float fB, out float fMa, out float fMb)
{
switch (eMType)
{
case Enums.EElemSuppType2D.eEl_00_00: // Both Side restrained
{
GetEIF_00_00_21_UV(Load, Member.FLength, out fA, out fB, out fMa, out fMb);
return;
}
case Enums.EElemSuppType2D.eEl_00_0_:
{
GetEIF_00_0__21_UV(Load, Member.FLength, out fA, out fB, out fMa, out fMb);
return;
}
case Enums.EElemSuppType2D.eEl_0__00:
{
GetEIF_0__00_21_UV(Load, Member.FLength, out fA, out fB, out fMa, out fMb);
return;
}
case Enums.EElemSuppType2D.eEl_0__0_:
{
GetEIF_0__0__21_UV(Load, Member.FLength, out fA, out fB, out fMa, out fMb);
return;
}
case Enums.EElemSuppType2D.eEl_00___:
{
GetEIF_00____21_UV(Load, Member.FLength, out fA, out fB, out fMa, out fMb);
return;
}
case Enums.EElemSuppType2D.eEl____00:
{
GetEIF____00_21_UV(Load, Member.FLength, out fA, out fB, out fMa, out fMb);
return;
}
default:
{
// Exception
fA = fB = fMa = fMb = 0.0f;
return;
}
}
}
void GetMLoadPart_21(CMLoad_21 Load, CE_1D_BASE Member, FEM_CALC_BASE.Enums.EElemSuppType2D eMType, out float fA, out float fB)
{
// Temporary
fA = fB = 0.0f;
}
public CExample_2D_12()
{
m_eSLN = ESLN.e2DD_1D; // 1D members in 2D model
m_eNDOF = (int)ENDOF.e2DEnv; // DOF in 2D
m_eGCS = EGCS.eGCSLeftHanded; // Global coordinate system
m_arrNodes = new BaseClasses.CNode[5];
m_arrMembers = new CMember[4];
m_arrMat = new CMat_00[1];
m_arrCrSc = new CRSC.CCrSc[3];
m_arrNSupports = new BaseClasses.CNSupport[1];
m_arrMLoads = new BaseClasses.CMLoad[1];
m_arrLoadCases = new BaseClasses.CLoadCase[1];
m_arrLoadCombs = new BaseClasses.CLoadCombination[1];
// Auxiliary
// Use basic SI units
// Sobota, J. > Statika stavebnych konstrukcii 2
// Load
float fq = 20f; // Unit [N/m]
// Geometry
float fa = 3.0f,
fb = 4.0f,
fc = 1.0f,
fd = 2.0f,
fe = 2.5f;// Unit [m]
// Materials
// Materials List - Materials Array - Fill Data of Materials Array
m_arrMat[0] = new CMat_00();
// Auxiliary, dopocitat
m_arrMat[0].m_fE = 10000000f; // Unit [Pa]
m_arrMat[0].m_fNu = 0.3f; // Unit [-]
m_arrMat[0].m_fG = 4000000; // Unit [Pa]
// Cross-sections
// CrSc List - CrSc Array - Fill Data of Cross-sections Array
// Cross-section
m_arrCrSc[0] = new CCrSc_0_00(); //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
m_arrCrSc[0].FA_g = 1.18095238095E-02f; // Unit [m^2]
m_arrCrSc[0].FI_y = 9.52380952381E-05f; // Unit [m^4]
m_arrCrSc[0].m_Mat = m_arrMat[0]; // Set CrSc Material
m_arrCrSc[1] = new CCrSc_0_00(); //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
m_arrCrSc[1].FA_g = 1.48571428571E-02f; // Unit [m^2]
m_arrCrSc[1].FI_y = 1.90476190476E-04f; // Unit [m^4]
m_arrCrSc[1].m_Mat = m_arrMat[0]; // Set CrSc Material
m_arrCrSc[2] = new CCrSc_0_00(); //!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
m_arrCrSc[2].FA_g = 1.35238095238E-02f; // Unit [m^2]
m_arrCrSc[2].FI_y = 1.42857142857E-07f; // Unit [m^4]
m_arrCrSc[2].m_Mat = m_arrMat[0]; // Set CrSc Material
// Nodes
// Nodes List - Nodes Array
// Node 1
m_arrNodes[0] = new CNode();
m_arrNodes[0].ID = 1;
m_arrNodes[0].X = 0f;
m_arrNodes[0].Y = 0f;
m_arrNodes[0].Z = 0f;
// Node 2
m_arrNodes[1] = new CNode();
m_arrNodes[1].ID = 2;
m_arrNodes[1].X = 0f;
m_arrNodes[1].Y = -fc - fd;
m_arrNodes[1].Z = 0f;
// Node 3
m_arrNodes[2] = new CNode();
m_arrNodes[2].ID = 3;
m_arrNodes[2].X = fa;
m_arrNodes[2].Y = -fc - fd - fe;
m_arrNodes[2].Z = 0f;
// Node 4
m_arrNodes[3] = new CNode();
m_arrNodes[3].ID = 4;
m_arrNodes[3].X = fa + fb;
m_arrNodes[3].Y = -fc - fd - fe;
m_arrNodes[3].Z = 0f;
// Node 5
m_arrNodes[4] = new CNode();
m_arrNodes[4].ID = 5;
m_arrNodes[4].X = fa + fb;
m_arrNodes[4].Y = fc;
m_arrNodes[4].Z = 0f;
// Sort by ID
//Array.Sort(m_arrNodes, new BaseClasses.CCompare_NodeID());
// Members
// Members List - Members Array
// Member 1 - 1-2
m_arrMembers[0] = new CMember();
m_arrMembers[0].ID = 1;
m_arrMembers[0].NodeStart = m_arrNodes[0];
m_arrMembers[0].NodeEnd = m_arrNodes[1];
m_arrMembers[0].CrScStart = m_arrCrSc[0];
// Member 2 - 2-3
m_arrMembers[1] = new CMember();
m_arrMembers[1].ID = 2;
m_arrMembers[1].NodeStart = m_arrNodes[1];
m_arrMembers[1].NodeEnd = m_arrNodes[2];
m_arrMembers[1].CrScStart = m_arrCrSc[1];
// Member 3 - 3-4
m_arrMembers[2] = new CMember();
m_arrMembers[2].ID = 3;
m_arrMembers[2].NodeStart = m_arrNodes[2];
m_arrMembers[2].NodeEnd = m_arrNodes[3];
m_arrMembers[2].CrScStart = m_arrCrSc[1];
// Member 4 - 4-5
m_arrMembers[3] = new CMember();
m_arrMembers[3].ID = 4;
m_arrMembers[3].NodeStart = m_arrNodes[3];
m_arrMembers[3].NodeEnd = m_arrNodes[4];
m_arrMembers[3].CrScStart = m_arrCrSc[2];
//Sort by ID
//Array.Sort(m_arrMembers, new BaseClasses.CCompare_MemberID());
// Nodal Supports - fill values
// Support 1 - NodeIDs: 1,5
m_arrNSupports[0] = new CNSupport(m_eNDOF);
m_arrNSupports[0].ID = 1;
m_arrNSupports[0].m_bRestrain[0] = true; // true - 1 restraint (infinity) / false - 0 - free (zero rigidity)
m_arrNSupports[0].m_bRestrain[1] = true;
m_arrNSupports[0].m_bRestrain[2] = true;
m_arrNSupports[0].m_iNodeCollection = new int[2];
m_arrNSupports[0].m_iNodeCollection[0] = 1;
m_arrNSupports[0].m_iNodeCollection[1] = 5;
// Sort by ID
Array.Sort(m_arrNSupports, new BaseClasses.CCompare_NSupportID());
// Member loads
// Load 1 - MemberIDs: 3
CMLoad_21 MLoad_q = new CMLoad_21(fq);
MLoad_q.ID = 1;
MLoad_q.MLoadTypeDistr = EMLoadTypeDistr.eMLT_QUF_W_21;
MLoad_q.MLoadType = EMLoadType.eMLT_F;
MLoad_q.EDirPPC = EMLoadDirPCC1.eMLD_PCC_FYU_MZV;
MLoad_q.IMemberCollection = new int[1];
MLoad_q.IMemberCollection[0] = 3;
m_arrMLoads[0] = MLoad_q;
// Load Cases
// Load Case 1
CLoadCase LoadCase0 = new CLoadCase();
LoadCase0.ID = 1;
m_arrLoadCases[0] = LoadCase0;
// Load Combinations
// Load Combination 1
CLoadCombination LoadComb0 = new CLoadCombination();
LoadComb0.ID = 1;
m_arrLoadCombs[0] = LoadComb0;
}
public CExample_2D_11()
{
m_eSLN = ESLN.e2DD_1D; // 1D members in 2D model
m_eNDOF = (int)ENDOF.e2DEnv; // DOF in 2D
m_eGCS = EGCS.eGCSLeftHanded; // Global coordinate system
m_arrNodes = new BaseClasses.CNode[5];
m_arrMembers = new CMember[4];
m_arrMat = new CMat_00[1];
m_arrCrSc = new CRSC.CCrSc[1];
m_arrNSupports = new BaseClasses.CNSupport[3];
m_arrNReleases = new BaseClasses.CNRelease[1];
m_arrNLoads = new BaseClasses.CNLoad[1];
m_arrMLoads = new BaseClasses.CMLoad[4];
m_arrLoadCases = new BaseClasses.CLoadCase[1];
m_arrLoadCombs = new BaseClasses.CLoadCombination[1];
// Auxiliary
// Load
float fF1 = 45000f; // Unit [N]
float fF2 = 55000f; // Unit [N]
float fM = 60000000f; // Unit [Nm]
float fq = 22f; // Unit [N/m]
// Geometry
float fa = 2.8f,
fb = 5.6f,
fc = 4.2f; // Unit [m]
// Materials
// Materials List - Materials Array - Fill Data of Materials Array
m_arrMat[0] = new CMat_00();
m_arrMat[0].m_fE = 10000000f; // Unit [Pa]
m_arrMat[0].m_fNu = 0.3f; // Unit [-]
m_arrMat[0].m_fG = 4000000; // Unit [Pa]
// Cross-sections
// CrSc List - CrSc Array - Fill Data of Cross-sections Array
m_arrCrSc[0] = new CRSC.CCrSc_0_00();
m_arrCrSc[0].FA_g = 0.12f; // Unit [m^2]
m_arrCrSc[0].FI_y = 0.0016f; // Unit [m^4]
m_arrCrSc[0].FI_z = 0.0016f; // Unit [m^4]
m_arrCrSc[0].m_Mat = m_arrMat[0]; // Set CrSc Material
// Nodes
// Nodes List - Nodes Array
// Nodes
// Node 1
m_arrNodes[0] = new CNode();
m_arrNodes[0].ID = 1;
m_arrNodes[0].X = fa + fb;
m_arrNodes[0].Y = fc;
m_arrNodes[0].Z = 0f;
// Node 2
m_arrNodes[1] = new CNode();
m_arrNodes[1].ID = 2;
m_arrNodes[1].X = fb;
m_arrNodes[1].Y = 0f;
m_arrNodes[1].Z = 0f;
// Node 3
m_arrNodes[2] = new CNode();
m_arrNodes[2].ID = 3;
m_arrNodes[2].X = 0f;
m_arrNodes[2].Y = 0f;
m_arrNodes[2].Z = 0f;
// Node 4
m_arrNodes[3] = new CNode();
m_arrNodes[3].ID = 4;
m_arrNodes[3].X = 0f;
m_arrNodes[3].Y = -fa;
m_arrNodes[3].Z = 0f;
// Node 5
m_arrNodes[4] = new CNode();
m_arrNodes[4].ID = 5;
m_arrNodes[4].X = 0f;
m_arrNodes[4].Y = fc;
m_arrNodes[4].Z = 0f;
// Sort by ID
//Array.Sort(m_arrNodes, new BaseClasses.CCompare_NodeID());
// Members
// Members List - Members Array
// Member 1 - 1-2
m_arrMembers[0] = new CMember();
m_arrMembers[0].ID = 1;
m_arrMembers[0].NodeStart = m_arrNodes[0];
m_arrMembers[0].NodeEnd = m_arrNodes[1];
m_arrMembers[0].CrScStart = m_arrCrSc[0];
// Member 2 - 2-3
m_arrMembers[1] = new CMember();
m_arrMembers[1].ID = 2;
m_arrMembers[1].NodeStart = m_arrNodes[1];
m_arrMembers[1].NodeEnd = m_arrNodes[2];
m_arrMembers[1].CrScStart = m_arrCrSc[0];
// Member 3 - 3-4
m_arrMembers[2] = new CMember();
m_arrMembers[2].ID = 3;
m_arrMembers[2].NodeStart = m_arrNodes[2];
m_arrMembers[2].NodeEnd = m_arrNodes[3];
m_arrMembers[2].CrScStart = m_arrCrSc[0];
// Member 4 - 3-5
m_arrMembers[3] = new CMember();
m_arrMembers[3].ID = 4;
m_arrMembers[3].NodeStart = m_arrNodes[2];
m_arrMembers[3].NodeEnd = m_arrNodes[4];
m_arrMembers[3].CrScStart = m_arrCrSc[0];
//Sort by ID
//Array.Sort(m_arrMembers, new BaseClasses.CCompare_MemberID());
// Nodal Supports - fill values
// Support 1 - NodeIDs: 1,4
m_arrNSupports[0] = new CNSupport(m_eNDOF);
m_arrNSupports[0].ID = 1;
m_arrNSupports[0].m_bRestrain[0] = true; // true - 1 restraint (infinity) / false - 0 - free (zero rigidity)
m_arrNSupports[0].m_bRestrain[1] = true;
m_arrNSupports[0].m_bRestrain[2] = true;
m_arrNSupports[0].m_iNodeCollection = new int[2];
m_arrNSupports[0].m_iNodeCollection[0] = 1;
m_arrNSupports[0].m_iNodeCollection[1] = 4;
// Support 2 - NodeIDs: 3
m_arrNSupports[1] = new CNSupport(m_eNDOF);
m_arrNSupports[1].ID = 2;
m_arrNSupports[1].m_bRestrain[0] = true; // true - 1 restraint (infinity) / false - 0 - free (zero rigidity) - support in GCS X-axis
m_arrNSupports[1].m_bRestrain[1] = false;
m_arrNSupports[1].m_bRestrain[2] = false;
m_arrNSupports[1].m_iNodeCollection = new int[1];
m_arrNSupports[1].m_iNodeCollection[0] = 3;
// Support 3 - NodeIDs: 5
m_arrNSupports[2] = new CNSupport(m_eNDOF);
m_arrNSupports[2].ID = 3;
m_arrNSupports[2].m_bRestrain[0] = true; // true - 1 restraint (infinity) / false - 0 - free (zero rigidity) - support in GCS X and Y-axis
m_arrNSupports[2].m_bRestrain[1] = true;
m_arrNSupports[2].m_bRestrain[2] = false;
m_arrNSupports[2].m_iNodeCollection = new int[1];
m_arrNSupports[2].m_iNodeCollection[0] = 5;
// Sort by ID
Array.Sort(m_arrNSupports, new BaseClasses.CCompare_NSupportID());
// Member releases
bool?[] bRelTemp = new bool?[m_eNDOF];
bRelTemp[0] = true; // true - 1 restraint (infinity rigidity) / false - 0 - free (zero rigidity) - support in LCS X and Y-axis
bRelTemp[1] = true;
bRelTemp[2] = false;
CNRelease NRelease1 = new CNRelease(m_eNDOF, m_arrNodes[1], m_arrMembers[0], bRelTemp, 0);
NRelease1.m_iNodeCollection = new int[1];
NRelease1.m_iNodeCollection[0] = 2;
NRelease1.m_iMembCollection = new int[2];
NRelease1.m_iMembCollection[0] = 1;
NRelease1.m_iMembCollection[1] = 2;
m_arrNReleases[0] = NRelease1;
m_arrMembers[0].CnRelease2 = m_arrNReleases[0]; // Release at end node
m_arrMembers[1].CnRelease1 = m_arrNReleases[0]; // Release at start node
// Nodal loads
// Load 1 - NodeIDs: 2
CNLoadSingle NLoad0 = new CNLoadSingle();
NLoad0.INLoad_ID = 1;
NLoad0.NLoadType = ENLoadType.eNLT_Fy;
NLoad0.INodeCollection = new int[1];
NLoad0.INodeCollection[0] = 2;
NLoad0.Value = fF2; // Positive
m_arrNLoads[0] = NLoad0;
// Member loads
// Load 1 and 2 - MemberIDs: 1
float fAlpha_1 = (236.3099f / 360f) * 2 * MATH.MathF.fPI; // Radians
float fF1x = Math.Abs(fF1 * (float)Math.Cos(fAlpha_1)); // Force in local coordinate system of member + possitive orientation in x-axis
float fF1y = fF1 * (float)Math.Sin(fAlpha_1); // Force in local coordinate system of member + negative orientation in y-axis
CMLoad_11 MLoad_F1x = new CMLoad_11(fF1x, 0.5f * m_arrMembers[0].FLength);
MLoad_F1x.ID = 1;
MLoad_F1x.MLoadTypeDistr = EMLoadTypeDistr.eMLT_FS_H_12;
MLoad_F1x.MLoadType = EMLoadType.eMLT_F;
MLoad_F1x.EDirPPC = EMLoadDirPCC1.eMLD_PCC_FXX_MXX;
MLoad_F1x.IMemberCollection = new int[1];
MLoad_F1x.IMemberCollection[0] = 1;
m_arrMLoads[0] = MLoad_F1x;
CMLoad_11 MLoad_F1y = new CMLoad_11(fF1y, 0.5f * m_arrMembers[0].FLength);
MLoad_F1y.ID = 2;
MLoad_F1y.MLoadTypeDistr = EMLoadTypeDistr.eMLT_FS_H_12;
MLoad_F1y.MLoadType = EMLoadType.eMLT_F;
MLoad_F1y.EDirPPC = EMLoadDirPCC1.eMLD_PCC_FYU_MZV;
MLoad_F1y.IMemberCollection = new int[1];
MLoad_F1y.IMemberCollection[0] = 1;
m_arrMLoads[1] = MLoad_F1y;
// Load 3 - MemberIDs: 3
CMLoad_21 MLoad_q = new CMLoad_21(fq);
MLoad_q.ID = 3;
MLoad_q.MLoadTypeDistr = EMLoadTypeDistr.eMLT_QUF_W_21;
MLoad_q.MLoadType = EMLoadType.eMLT_F;
MLoad_q.EDirPPC = EMLoadDirPCC1.eMLD_PCC_FYU_MZV;
MLoad_q.IMemberCollection = new int[1];
MLoad_q.IMemberCollection[0] = 2;
m_arrMLoads[2] = MLoad_q;
// Load 4 - MemberIDs: 4
CMLoad_11 MLoad_M = new CMLoad_11(-fM, 0.5f * m_arrMembers[3].FLength);
MLoad_M.ID = 4;
MLoad_M.MLoadTypeDistr = EMLoadTypeDistr.eMLT_FS_H_12;
MLoad_M.MLoadType = EMLoadType.eMLT_M;
MLoad_M.EDirPPC = EMLoadDirPCC1.eMLD_PCC_FYU_MZV;
MLoad_M.IMemberCollection = new int[1];
MLoad_M.IMemberCollection[0] = 4;
m_arrMLoads[3] = MLoad_M;
// Load Cases
// Load Case 1
m_arrLoadCases[0] = new CLoadCase();
m_arrLoadCases[0].ID = 1;
// Load Combinations
// Load Combination 1
m_arrLoadCombs[0] = new CLoadCombination();
m_arrLoadCombs[0].ID = 1;
}