#mshFileName = 'quadCurvedGeoStrut.msh' mshFileName = 'quadRealGeoStrut.msh' #mshFileName = 'quadCurvedGeo.msh' #mshFileName = 'quad1.msh' #mshFileName = 'quad2.msh' #mshFileName = 'quad3.msh' #mshFileName = 'quad4.msh' #mshFileName = 'quad5.msh' pathMSHFile = searchMSH.Find(mshFileName) if pathMSHFile == 'File not found': sys.exit() mesh = importMSH.Quad2D(pathMSHFile, mshFileName) numNodes = mesh.numNodes numElements = mesh.numElements x = mesh.x y = mesh.y IEN = mesh.IEN boundaryEdges = mesh.boundaryEdges boundaryNodes = mesh.boundaryNodes neighborsNodes = mesh.neighborsNodes neighborsNodesALE = mesh.neighborsNodesALE neighborsElements = mesh.neighborsElements minLengthMesh = mesh.minLengthMesh FreedomDegree = mesh.FreedomDegree numPhysical = mesh.numPhysical
sys.exit()
mesh = importMSH.Mini2D(pathMSHFile, mshFileName, equation_number)
mesh.coord()
mesh.ien()
# Quad Element
elif polynomial_option == 3:
mshFileName = 'quadHalfPoiseuille_quad.msh'
equation_number = 3
pathMSHFile = searchMSH.Find(mshFileName)
if pathMSHFile == 'File not found':
sys.exit()
mesh = importMSH.Quad2D(pathMSHFile, mshFileName, equation_number)
mesh.coord()
mesh.ien()
# Cubic Element
elif polynomial_option == 4:
mshFileName = 'cubicHalfPoiseuille_cubic.msh'
equation_number = 3
pathMSHFile = searchMSH.Find(mshFileName)
if pathMSHFile == 'File not found':
sys.exit()
mesh = importMSH.Cubic2D(pathMSHFile, mshFileName, equation_number)
mesh.coord()
mesh.ien()