def Stress_TR03(elem):
mat_id = elem["MaterialId"] - 1
material = globalvars.data["Materials"][mat_id]
Dmat = dmat_Solids2D(material)
Eps_vec = Strain_TR03(elem)
Sigma_vec = np.matmul(Dmat, Eps_vec)
return Sigma_vec
def stiffness_QU04(elem):
ProblemType = globalvars.data["ProblemType"]
mat_id = elem["MaterialId"] - 1
material = globalvars.data["Materials"][mat_id]
Dmat = dmat_Solids2D(material)
thickness = 1.0
if(material["Plane_Type"] == "Plane_Stress"):
thickness = material["Thickness"]
Dmat *= thickness
return stiffnessMatCalc(elem, "QU04", ProblemType, Dmat)
def Stress_Solid(elem, ElemType, Eps_vec):
mat_id = elem["MaterialId"] - 1
material = globalvars.data["Materials"][mat_id]
if (ElemType == "BAR02" or ElemType == "BAR03" or ElemType == "TRUSS02"):
Young = material["Young"]
Sigma_vec = Young * Eps_vec
else:
Dmat = dmat_Solids2D(material)
Sigma_vec = np.matmul(Dmat, Eps_vec)
return Sigma_vec
def stiffness_TR03(elem):
mat_id = elem["MaterialId"] - 1
material = globalvars.data["Materials"][mat_id]
Dmat = dmat_Solids2D(material)
Bmat, Area = Bmat_TR03(elem)
DmatxB = np.matmul(Dmat, Bmat)
thickness = 1.0
if(material["Plane_Type"] == "Plane_Stress"):
thickness =material["Thickness"]
rigimat = np.matmul(Bmat.transpose(), DmatxB) * Area * thickness
return rigimat