public void GenerateMesh1D(CModel TopoModel)
{
// Generate FEM nodes
// Nodes
for (int i = 0; i < TopoModel.m_arrNodes.Length; i++)
{
// Create FEM node
CFemNode TempNode = new CFemNode(TopoModel.m_arrNodes[i]);
m_arrFemNodes[i] = TempNode;
// Set FEM node DOF
// Get nodal support
// Search if node is in list of supported nodes for each nodal support
for (int i2 = 0; i2 < TopoModel.m_arrNSupports.Length; i2++) // Check all nodal supports
{
for (int j = 0; j < TopoModel.m_arrNSupports[i2].m_iNodeCollection.Length; j++) // Check list of nodes (Nodes IDs collection)
{
if (m_arrFemNodes[i].ID == TopoModel.m_arrNSupports[i2].m_iNodeCollection[j])
{
// DOF of nodes are free - zero rigidity of restraints false as default
m_arrFemNodes[i].m_ArrNodeDOF[(int)e3D_DOF.eUX] = TopoModel.m_arrNSupports[i2].m_bRestrain[(int)e3D_DOF.eUX]; // !!! 3D Environment enum
m_arrFemNodes[i].m_ArrNodeDOF[(int)e3D_DOF.eUY] = TopoModel.m_arrNSupports[i2].m_bRestrain[(int)e3D_DOF.eUY]; // !!! 3D Environment enum
m_arrFemNodes[i].m_ArrNodeDOF[(int)e3D_DOF.eUZ] = TopoModel.m_arrNSupports[i2].m_bRestrain[(int)e3D_DOF.eUZ]; // !!! 3D Environment enum
m_arrFemNodes[i].m_ArrNodeDOF[(int)e3D_DOF.eRX] = TopoModel.m_arrNSupports[i2].m_bRestrain[(int)e3D_DOF.eRX]; // !!! 3D Environment enum
m_arrFemNodes[i].m_ArrNodeDOF[(int)e3D_DOF.eRY] = TopoModel.m_arrNSupports[i2].m_bRestrain[(int)e3D_DOF.eRY]; // !!! 3D Environment enum
m_arrFemNodes[i].m_ArrNodeDOF[(int)e3D_DOF.eRZ] = TopoModel.m_arrNSupports[i2].m_bRestrain[(int)e3D_DOF.eRZ]; // !!! 3D Environment enum
}
}
}
}
// Generate FEM members
// Members
for (int i = 0; i < TopoModel.m_arrMembers.Length; i++)
{
CE_1D TempMember = new CE_1D(TopoModel.m_arrMembers[i], m_arrFemNodes);
m_arrFemMembers[i] = TempMember;
// Set FEM Member DOF
if (TopoModel.m_arrNReleases != null) // Some releases exist
{
for (int j = 0; j < TopoModel.m_arrNReleases.Length; j++) // Check each release
{
for (int k = 0; k < TopoModel.m_arrNReleases[j].m_iMembCollection.Length; k++) // Check each member in collection
{
if (TopoModel.m_arrNReleases[j].m_iMembCollection[k] == m_arrFemMembers[i].ID) // if release exists on member
{
if (TopoModel.m_arrMembers[i].CnRelease1 != null)
m_arrFemMembers[i].CnRelease1 = TopoModel.m_arrMembers[i].CnRelease1;
if (TopoModel.m_arrMembers[i].CnRelease2 != null)
m_arrFemMembers[i].CnRelease2 = TopoModel.m_arrMembers[i].CnRelease2;
}
}
}
}
}
////////////////////////////////////////////////////////////////////////////////////////////////
// Additional data of nodes depending on generated members
for (int i = 0; i < m_arrFemNodes.Length; i++)
{
m_arrFemNodes[i].m_iMemberCollection = new System.Collections.ArrayList(); // Allocate collection memory
for (int j = 0; j < m_arrFemMembers.Length; j++)
{
if (m_arrFemNodes[i].ID == m_arrFemMembers[j].NodeStart.ID || m_arrFemNodes[i].ID == m_arrFemMembers[j].NodeEnd.ID) // Is node ID same as member start or end node ID
{
m_arrFemNodes[i].m_iMemberCollection.Add(m_arrFemMembers[j].ID); // Add FEMmember ID to the FEM node list
}
}
}
// Releases !!!! - copied from 2D solution, check implementation in 3D !!!
// If only two members are connected in one node and if release exists at that node, copy this release from one member to the another
for (int i = 0; i < m_arrFemNodes.Length; i++)
{
if (m_arrFemNodes[i].m_iMemberCollection != null && m_arrFemNodes[i].m_iMemberCollection.Count == 2) // Node is connected to two members
{
// We know member ID, so we can get index of members in list
int iMember1_index = -1;
int iMember2_index = -1;
for (int j = 0; j < m_arrFemMembers.Length; j++)
{
// 1st member index
if (iMember1_index < 0 && m_arrFemNodes[i].m_iMemberCollection.Contains(m_arrFemMembers[j].ID)) // if Member ID is in the list
{
iMember1_index = j; // Set 1st
}
if (iMember1_index > -1) // Index was defined, we can break cycle
break;
}
// 2nd member index
for (int k = iMember1_index + 1; k < m_arrFemMembers.Length; k++) // Search for second only in interval between first founded member and last member
{
if (iMember2_index < 0 && m_arrFemNodes[i].m_iMemberCollection.Contains(m_arrFemMembers[k].ID)) // if Member ID is in the list interval
{
iMember2_index = k;
}
if (iMember2_index > -1) // Index was defined, we can break cycle
break;
}
// If relases exist, they are neccesary to define DOF of both members, therefore copy release of one member to the another one
if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember1_index].NodeStart.ID && m_arrFemMembers[iMember1_index].CnRelease1 != null)
{
if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember2_index].NodeStart.ID && m_arrFemMembers[iMember2_index].CnRelease1 == null)
m_arrFemMembers[iMember2_index].CnRelease1.m_bRestrain = m_arrFemMembers[iMember1_index].CnRelease1.m_bRestrain;
else if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember2_index].NodeEnd.ID && m_arrFemMembers[iMember2_index].CnRelease2 == null)
m_arrFemMembers[iMember2_index].CnRelease2.m_bRestrain = m_arrFemMembers[iMember1_index].CnRelease1.m_bRestrain;
}
else if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember1_index].NodeEnd.ID && m_arrFemMembers[iMember1_index].CnRelease2 != null)
{
if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember2_index].NodeStart.ID && m_arrFemMembers[iMember2_index].CnRelease1 == null)
m_arrFemMembers[iMember2_index].CnRelease1.m_bRestrain = m_arrFemMembers[iMember1_index].CnRelease2.m_bRestrain;
else if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember2_index].NodeEnd.ID && m_arrFemMembers[iMember2_index].CnRelease2 == null)
m_arrFemMembers[iMember2_index].CnRelease2.m_bRestrain = m_arrFemMembers[iMember1_index].CnRelease2.m_bRestrain;
}
else if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember2_index].NodeStart.ID && m_arrFemMembers[iMember2_index].CnRelease1 != null)
{
if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember1_index].NodeStart.ID && m_arrFemMembers[iMember1_index].CnRelease1 == null)
m_arrFemMembers[iMember1_index].CnRelease1.m_bRestrain = m_arrFemMembers[iMember2_index].CnRelease1.m_bRestrain;
else if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember1_index].NodeEnd.ID && m_arrFemMembers[iMember1_index].CnRelease2 == null)
m_arrFemMembers[iMember1_index].CnRelease2.m_bRestrain = m_arrFemMembers[iMember2_index].CnRelease1.m_bRestrain;
}
else if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember2_index].NodeEnd.ID && m_arrFemMembers[iMember2_index].CnRelease2 != null)
{
if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember1_index].NodeStart.ID && m_arrFemMembers[iMember1_index].CnRelease1 == null)
m_arrFemMembers[iMember1_index].CnRelease1.m_bRestrain = m_arrFemMembers[iMember2_index].CnRelease2.m_bRestrain;
else if (m_arrFemNodes[i].ID == m_arrFemMembers[iMember1_index].NodeEnd.ID && m_arrFemMembers[iMember1_index].CnRelease2 == null)
m_arrFemMembers[iMember1_index].CnRelease2.m_bRestrain = m_arrFemMembers[iMember2_index].CnRelease2.m_bRestrain;
}
}
}
// Additional data of members
// Fill Members stifeness matrices
for (int i = 0; i < m_arrFemMembers.Length; i++)
{
m_arrFemMembers[i].FillBasic3_StiffMatrices();
}
//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// External Loads
// Fill vectors
// Nodal loads
// Fill vector of external load for each node in list
// Node could be included in all lists only once, more than one nodal load in one node is not allowed
if (!object.ReferenceEquals(null, TopoModel.m_arrNLoads))
{
for (int i = 0; i < TopoModel.m_arrNLoads.Length; i++) // Each load
{
if (TopoModel.m_arrNLoads[i].INodeCollection != null) // Check if some nodes are in list
{
for (int j = 0; j < TopoModel.m_arrNLoads[i].INodeCollection.Length; j++) // Each node in collection
{
// Set load vector values
for (int k = 0; k < m_arrFemNodes.Length; k++) // Neefektivne prechadzat zbytocne vsetky uzly kedze pozname konkretne ID zatazenych
{
if (TopoModel.m_arrNLoads[i].INodeCollection.Contains(TopoModel.m_arrNodes[k].ID)) // If node ID is same as collection item
{
// Check object class
if (TopoModel.m_arrNLoads[i] is BaseClasses.CNLoadAll)
{
CNLoadAll oTemp = new CNLoadAll(); // Create new object of necessary type
oTemp = (CNLoadAll)(TopoModel.m_arrNLoads[i]); // Change object type
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eFX] = oTemp.Value_FX;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eFY] = oTemp.Value_FY;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eFZ] = oTemp.Value_FZ;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eMX] = oTemp.Value_MX;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eMY] = oTemp.Value_MY;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eMZ] = oTemp.Value_MZ;
/*
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eFX] = TopoModel.m_arrNLoads[i].Value_FX;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eFY] = TopoModel.m_arrNLoads[i].Value_FY;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eFZ] = TopoModel.m_arrNLoads[i].Value_FZ;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eMX] = TopoModel.m_arrNLoads[i].Value_MX;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eMY] = TopoModel.m_arrNLoads[i].Value_MY;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eMZ] = TopoModel.m_arrNLoads[i].Value_MZ;
*/
}
else
{
CNLoadSingle oTemp = new CNLoadSingle(); // Create new object of necessary type
oTemp = (CNLoadSingle)(TopoModel.m_arrNLoads[i]); // Change object type
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eFX] = oTemp.NLoadType == ENLoadType.eNLT_Fx ? oTemp.Value : 0.0f;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eFY] = oTemp.NLoadType == ENLoadType.eNLT_Fy ? oTemp.Value : 0.0f;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eFZ] = oTemp.NLoadType == ENLoadType.eNLT_Fz ? oTemp.Value : 0.0f;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eMX] = oTemp.NLoadType == ENLoadType.eNLT_Mx ? oTemp.Value : 0.0f;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eMY] = oTemp.NLoadType == ENLoadType.eNLT_My ? oTemp.Value : 0.0f;
m_arrFemNodes[k].m_VDirNodeLoad.FVectorItems[(int)e3D_E_F.eMZ] = oTemp.NLoadType == ENLoadType.eNLT_Mz ? oTemp.Value : 0.0f;
}
}
}
}
}
}
}
// Member loads
// Set primary end forces only due to member loads in local coordinate system LCS
// Summation of loads applied on one member
// There can by more loads on one member, member could be in various loads lists (but only once in one list)
if (!object.ReferenceEquals(null, TopoModel.m_arrMLoads))
{
for (int i = 0; i < TopoModel.m_arrMLoads.Length; i++) // Each load
{
if (TopoModel.m_arrMLoads[i].IMemberCollection != null) // Check if some members are in list
{
for (int j = 0; j < TopoModel.m_arrMLoads[i].IMemberCollection.Length; j++) // Each node in collection
{
// Set end forces due to member load
for (int k = 0; k < m_arrFemMembers.Length; k++) // Neefektivne prechadzat zbytocne vsetky pruty kedze pozname konkretne ID zatazenych
{
// Temporary value of end forces due to particular external force
float fTemp_A_UXRX = 0f, fTemp_B_UXRX = 0f, fTemp_Ma_UXRX = 0f, fTemp_Mb_UXRX = 0f;
float fTemp_A_UYRZ = 0f, fTemp_B_UYRZ = 0f, fTemp_Ma_UYRZ = 0f, fTemp_Mb_UYRZ = 0f;
float fTemp_A_UZRY = 0f, fTemp_B_UZRY = 0f, fTemp_Ma_UZRY = 0f, fTemp_Mb_UZRY = 0f;
// Fill end forces due to external forces vectors for member particular load index i and member index k
// Depends on load dirrection and member support type
// Use member PRINCIPAL coordinate system
CMLoadPart objAux = new CMLoadPart(TopoModel, m_arrFemMembers, i, k,
out fTemp_A_UXRX, out fTemp_Ma_UXRX, out fTemp_A_UYRZ, out fTemp_Ma_UYRZ, out fTemp_A_UZRY, out fTemp_Ma_UZRY,
out fTemp_B_UXRX, out fTemp_Mb_UXRX, out fTemp_B_UYRZ, out fTemp_Mb_UYRZ, out fTemp_B_UZRY, out fTemp_Mb_UZRY
);
// Add values of temporary end forces due to particular load to the end forces items of vector
// Primary end forces due member loading in local coordinate system LCS
// Start Node
m_arrFemMembers[k].m_VElemPEF_LCS_StNode.FVectorItems[(int)e3D_E_F.eFX] += fTemp_A_UXRX;
m_arrFemMembers[k].m_VElemPEF_LCS_StNode.FVectorItems[(int)e3D_E_F.eMX] += fTemp_Ma_UXRX;
m_arrFemMembers[k].m_VElemPEF_LCS_StNode.FVectorItems[(int)e3D_E_F.eFY] += fTemp_A_UYRZ; // !!! Signs - nutne skontrolovat znamienka podla smeru lokalnych osi a orientacie zatazenia
m_arrFemMembers[k].m_VElemPEF_LCS_StNode.FVectorItems[(int)e3D_E_F.eMZ] += fTemp_Ma_UYRZ;
m_arrFemMembers[k].m_VElemPEF_LCS_StNode.FVectorItems[(int)e3D_E_F.eFZ] += fTemp_A_UZRY; // !!! Signs - nutne skontrolovat znamienka podla smeru lokalnych osi a orientacie zatazenia
m_arrFemMembers[k].m_VElemPEF_LCS_StNode.FVectorItems[(int)e3D_E_F.eMY] += fTemp_Ma_UZRY;
// End Node
m_arrFemMembers[k].m_VElemPEF_LCS_EnNode.FVectorItems[(int)e3D_E_F.eFX] += fTemp_B_UXRX;
m_arrFemMembers[k].m_VElemPEF_LCS_EnNode.FVectorItems[(int)e3D_E_F.eMX] += fTemp_Mb_UXRX;
m_arrFemMembers[k].m_VElemPEF_LCS_EnNode.FVectorItems[(int)e3D_E_F.eFY] += fTemp_B_UYRZ; //-fTemp_B_UYRZ; // Zmena znamienka pre silu Vb na konci pruta, znamienko je opacne nez u reakcie, toto su vsak reakcie
m_arrFemMembers[k].m_VElemPEF_LCS_EnNode.FVectorItems[(int)e3D_E_F.eMZ] += fTemp_Mb_UYRZ;
m_arrFemMembers[k].m_VElemPEF_LCS_EnNode.FVectorItems[(int)e3D_E_F.eFZ] += fTemp_B_UZRY; //-fTemp_B_UZRY; // Zmena znamienka pre silu Vb na konci pruta, znamienko je opacne nez u reakcie, toto su vsak reakcie
m_arrFemMembers[k].m_VElemPEF_LCS_EnNode.FVectorItems[(int)e3D_E_F.eMY] += fTemp_Mb_UZRY;
}
}
}
}
}
// Set primary end forces only due to member loads in global coordinate system GCS
foreach (CE_1D Elem in m_arrFemMembers)
{
Elem.SetGetPEF_GCS();
}
}
// public CCrSc_0_20 objCrScHollow = new CCrSc_0_20(200, 10);
// public CCrSc_0_22 objCrScHollow = new CCrSc_0_22(200, 10);
// public CCrSc_0_23 objCrScHollow = new CCrSc_0_23(200, 100, 5);
// public CCrSc_0_24 objCrScHollow = new CCrSc_0_24(100, 10);
// public CCrSc_0_25 objCrScHollow = new CCrSc_HL(200, 100, 10, 5);
// public CCrSc_0_26 objCrScHollow = new CCrSc_0_26(200,10);
// public CCrSc_0_27 objCrScHollow = new CCrSc_0_27(200, 10);
// public CCrSc_0_28 objCrScHollow = new CCrSc_0_28(200, 10);
// public CCrSc_0_50 objCrScSolid = new CCrSc_0_50(200,100,10,5);
// public CCrSc_0_50 objCrScSolid = new CCrSc_0_50(200, 100, 150, 10, 15, 5, 70);
// public CCrSc_0_52 objCrScSolid = new CCrSc_0_52(200, 100, 10, 5, 30);
// public CCrSc_0_52 objCrScSolid = new CCrSc_0_52(200, 100, 150, 10, 15, 5, 30, 70);
// public CCrSc_0_54 objCrScSolid = new CCrSc_0_54(200, 100, 10, 5,30,50);
// public CCrSc_0_56 objCrScSolid = new CCrSc_0_56(200, 100, 10, 5,130);
// public CCrSc_0_56 objCrScSolid = new CCrSc_0_56(200, 100, 30, 10, 5, 20, 130);
// public CCrSc_0_58 objCrScSolid = new CCrSc_0_58(200, 100, 10, 5);
// public CCrSc_0_58 objCrScSolid = new CCrSc_0_58(200, 100, 150, 10, 15, 5, 30, 70);
// public CCrSc_0_60 objCrScSolid = new CCrSc_0_60(200,100,10);
// public CCrSc_0_61 objCrScSolid = new CCrSc_0_61(100,10);
// public CCrSc_3_00 objCrScSolid = new CCrSc_3_00(0,200,90,11.3f,7.5f,7.5f,4.5f,159.1f);
// public CCrSc_3_00 objCrScSolid = new CCrSc_3_00(1, 200, 90, 11.3f, 7.5f, 4.5f, 159.1f);
// public CCrSc_3_00 objCrScSolid = new CCrSc_3_00(2, 200, 90, 11.3f, 7.5f, 7.5f, 159.1f);
// public CCrSc_3_01 objCrScSolid = new CCrSc_3_01(0, 200, 90, 120, 11.3f, 11.3f, 7.5f, 7.5f, 4.5f, 159.1f);
// public CCrSc_3_01 objCrScSolid = new CCrSc_3_01(1, 200, 90, 120, 11.3f, 11.3f, 7.5f, 4.5f, 80.0f);
// public CCrSc_3_01 objCrScSolid = new CCrSc_3_01(2, 200, 80, 100, 10, 8, 12, 10, 5, 5, 80);
// public CCrSc_3_02 objCrScSolid = new CCrSc_3_02(0,300, 100, 16, 10, 10, 8, 232, 27);
// public CCrSc_3_02 objCrScSolid = new CCrSc_3_02(2,300, 100, 12.5f, 10, 10, 250, 30);
// public CCrSc_3_03 objCrScSolid = new CCrSc_3_03(150, 10, 15, 10, 50);
// public CCrSc_3_04 objCrScSolid = new CCrSc_3_04(250, 150, 20, 15, 20, 50, 75);
// public CCrSc_3_05 objCrScHollow = new CCrSc_3_05(200,10);
// public CCrSc_3_06 objCrScHollow = new CCrSc_3_06(200, 100, 20);
// public CCrSc_3_07 objCrScHollow = new CCrSc_3_07(0, 200, 100, 10, 30); // Both radii, coincident centres
// public CCrSc_3_07 objCrScHollow = new CCrSc_3_07(1, 200, 100, 10, 5, 30); // Both radii, incoincident centres
// public CCrSc_3_07 objCrScHollow = new CCrSc_3_07(2, 500, 300, 20, 50); // Outside radius = 0
// public CCrSc_3_07 objCrScHollow = new CCrSc_3_07(3, 500, 300, 30); // Outside radius = 0, coincident centres
// public CCrSc_3_07 objCrScHollow = new CCrSc_3_07(0, 400, 150, 20); // Both radii, coincident centres
// public CCrSc_3_07 objCrScHollow = new CCrSc_3_07(2, 500, 300, 20); // Outside radius = 0
// public CCrSc_3_07 objCrScHollow = new CCrSc_3_07(3, 500, 300, 30); // Inside radius = 0, coincident centres
// public CCrSc_3_07 objCrScHollow = new CCrSc_3_07(5, 400, 150, 20); // No radii, Outside radius = 0, Inside radius = 0
// public CCrSc_3_08 objCrScSolid = new CCrSc_3_08(0, 120, 100, 12.5f, 15, 100, 10, 10, 5, 90);
// public CCrSc_3_08 objCrScSolid = new CCrSc_3_08(1, 120, 100, 12.5f, 15, 100, 10, 10, 5, 95);
// public CCrSc_3_08 objCrScSolid = new CCrSc_3_08(2, 120, 100, 12.5f, 15, 100, 10, 5, 105);
// public CCrSc_3_08 objCrScSolid = new CCrSc_3_08(3, 120, 100, 12.5f, 15, 100, 10, 5, 95);
public CTest1()
{
// !!!!!! Vytvarat len ODKAZY na objekty "ref", aby sa zbytocne nevytvarali lokalne kopie
// Materials
// Materials List - Materials Array - Fill Data of Materials Array
arrMat[0] = new CMat_00(); // Vytvarat priamo pre konkretne typy materialov (uzivatelsky, ocel, beton, drevo, hlinik)
arrMat[1] = new CMat_02_00();
arrMat[2] = new CMat_03_00();
arrMat[3] = new CMat_05_00();
arrMat[4] = new CMat_09_00();
// CrSc
// CrSc List - CrSc Array - Fill Data of Cross-sections Array
arrCrSc[0] = new CRSC.CCrSc_3_07(3, 500, 300, 30);
arrCrSc[1] = new CRSC.CCrSc_0_05(100, 50);
arrCrSc[2] = new CRSC.CCrSc_3_04(250, 150, 20, 15, 20, 50, 75);
arrCrSc[0].m_Mat = arrMat[0];
arrCrSc[1].m_Mat = arrMat[1];
arrCrSc[2].m_Mat = arrMat[4];
// Nodes
// Nodes List - Nodes Array
arrNodes[0] = new BaseClasses.CNode(1, 500, 0, 2500, 0);
arrNodes[1] = new BaseClasses.CNode(2, 2500, 0, 2500, 0);
arrNodes[2] = new BaseClasses.CNode(3, 5500, 0, 2500, 0);
arrNodes[3] = new BaseClasses.CNode(4, 500, 0, 500, 0);
arrNodes[4] = new BaseClasses.CNode(5, 2500, 0, 500, 0);
arrNodes[5] = new BaseClasses.CNode(6, 5500, 0, 500, 0);
// Sort by ID
//Array.Sort(arrNodes, new BaseClasses.CCompare_NodeID());
// Members
// Members List - Members Array
arrMembers[0] = new BaseClasses.CMember(1, arrNodes[0], arrNodes[1], arrCrSc[0], 0);
arrMembers[1] = new BaseClasses.CMember(2, arrNodes[1], arrNodes[2], arrCrSc[0], 0);
arrMembers[2] = new BaseClasses.CMember(3, arrNodes[0], arrNodes[3], arrCrSc[1], 0);
arrMembers[3] = new BaseClasses.CMember(4, arrNodes[1], arrNodes[4], arrCrSc[0], 0);
arrMembers[4] = new BaseClasses.CMember(5, arrNodes[2], arrNodes[5], arrCrSc[0], 0);
arrMembers[5] = new BaseClasses.CMember(6, arrNodes[3], arrNodes[4], arrCrSc[2], 0);
arrMembers[6] = new BaseClasses.CMember(7, arrNodes[4], arrNodes[5], arrCrSc[0], 0);
arrMembers[7] = new BaseClasses.CMember(8, arrNodes[1], arrNodes[3], arrCrSc[0], 0);
arrMembers[8] = new BaseClasses.CMember(9, arrNodes[1], arrNodes[5], arrCrSc[1], 0);
//Sort by ID
//Array.Sort(arrMembers, new BaseClasses.CCompare_MemberID());
// Nodal Supports - fill values
// Set values
bool[] bSupport1 = { true, false, true, false, true, false };
bool[] bSupport2 = { false, false, true, false, true, false };
bool[] bSupport3 = { true, false, false, false, false, false };
// Create Support Objects
// Pozn. Jednym z parametrov by malo byt pole ID uzlov v ktorych je zadefinovana tato podpora
// objekt podpory bude len jeden a dotknute uzly budu vediet ze na ich podpora existuje a ake je konkretne ID jej nastaveni
arrNSupports[0] = new BaseClasses.CNSupport(6,1, arrNodes[0], bSupport1, 0);
arrNSupports[1] = new BaseClasses.CNSupport(6,2, arrNodes[2], bSupport2, 0);
arrNSupports[2] = new BaseClasses.CNSupport(6,3, arrNodes[5], bSupport3, 0);
// Sort by ID
Array.Sort(arrNSupports, new BaseClasses.CCompare_NSupportID());
// Member Releases / hinges - fill values
// Pozn. Jednym z parametrov by malo byt pole ID prutov v ktorych je zadefinovany klb
// objekt klbu bude len jeden a dotknute pruty budu vediet ze na ich klb existuje, v ktorom uzle (start 1 / End 2) a ake je konkretne ID jeho nastaveni
// Set values
//bool[] bMembRelase1 = { true, false, false, false, false, false };
//bool[] bMembRelase2 = { false, false, true, false, false, false };
//bool[] bMembRelase3 = { false, false, false, false, true, false };
bool? [] bMembRelase4 = { true, false, true, false, true, false };
// Create Release / Hinge Objects
arrMembers[7].CnRelease1 = new BaseClasses.CNRelease(6,bMembRelase4, 0);
arrMembers[8].CnRelease2 = new BaseClasses.CNRelease(6,bMembRelase4, 0);
// Nodal Forces - fill values
// Pozn. Jednym z parametrov by malo byt pole ID uzlov v ktorych je zadefinovane uzlove zatazenie
// objekt zatazenia bude len jeden a dotknute uzly budu vediet ze na ich dane zatazenie existuje a ake je konkretne ID jeho vlastnosti
arrNLoads[0] = new BaseClasses.CNLoadSingle(arrNodes[1], ENLoadType.eNLT_Fx, 40.0f, true, 0);
arrNLoads[1] = new BaseClasses.CNLoadSingle(arrNodes[4], ENLoadType.eNLT_Fx, -60.0f, true, 0);
arrNLoads[2] = new BaseClasses.CNLoadSingle(arrNodes[5], ENLoadType.eNLT_Fz, 80.0f, true, 0);
}
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;
}
