From ba6b338d521e106403267192d65032845556e598 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?=E7=A6=85=E5=85=83=E5=A4=A9=E9=81=93?= Date: Tue, 25 Jul 2023 14:36:10 +0800 Subject: [PATCH] =?UTF-8?q?=E7=A7=BB=E9=99=A4=E4=B8=8D=E9=9C=80=E8=A6=81?= =?UTF-8?q?=E7=9A=84=E4=BB=A3=E7=A0=81?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .../MSHdataExchange.cpp | 520 ------------------ src/PluginMeshDataExchange/MSHdataExchange.h | 23 +- 2 files changed, 1 insertion(+), 542 deletions(-) diff --git a/src/PluginMeshDataExchange/MSHdataExchange.cpp b/src/PluginMeshDataExchange/MSHdataExchange.cpp index 454d884..5886acc 100644 --- a/src/PluginMeshDataExchange/MSHdataExchange.cpp +++ b/src/PluginMeshDataExchange/MSHdataExchange.cpp @@ -42,51 +42,6 @@ namespace MeshData { _stream = nullptr; } - QString MSHdataExchange::readLine() - { - while(_threadRuning) { - if(_stream->atEnd()) { - //_threadRuning = false; - return QString(); - ; - } - QString line = _stream->readLine().simplified(); - if(line.isEmpty()) - continue; - return line; - } - return QString(); - } - - bool MSHdataExchange::readHeader() - { - QFileInfo tempinfo(_fileName); - if(!tempinfo.exists()) - return false; - QString name = tempinfo.fileName(); - QString path = tempinfo.filePath(); - QFile tempfile(_fileName); - if(!tempfile.open(QIODevice::ReadOnly)) - return false; - // QTextStream* tempStream=new QTextStream(&tempfile); - QTextStream tempStream(&tempfile); - while(!tempStream.atEnd()) { - if(!_threadRuning) - return false; - QString line = tempStream.readLine().simplified(); - if(line.contains("GAMBIT")) { - mMeshType = typeGambit; - tempfile.close(); - return true; - } else if(line.contains("Fluent")) { - mMeshType = typeFluent; - tempfile.close(); - return true; - } - } - // delete tempStream; - return false; - } bool MSHdataExchange::write() { @@ -114,12 +69,6 @@ namespace MeshData { bool MSHdataExchange::read() { - if(!readHeader()) { - // QMessageBox::warning(nullptr, tr("Prompt"), tr("The file format could not parse the - // amount!"), QMessageBox::Ok); - return false; - } - _file = new QFile(); QFileInfo info(_fileName); if(!info.exists()) @@ -236,475 +185,6 @@ namespace MeshData { // return true; } - void MSHdataExchange::readPoints10(vtkUnstructuredGrid* dataset, QString info) - { - vtkSmartPointer points = dataset->GetPoints(); - if(points == nullptr) { - points = vtkSmartPointer::New(); - dataset->SetPoints(points); - } - int type = 0; - int index = 0; - int setid = 0; - int startIndex = 0; - int endIndex = 0; - QStringList infos = info.split(" "); - if(infos.length() > 5) { - if(!_threadRuning) - return; - index = strToInt(infos.at(0), true); - setid = strToInt(infos.at(1), true); - startIndex = strToInt(infos.at(2), true); - endIndex = strToInt(infos.at(3), true); - _ND = strToInt(infos.at(5)); - MeshSet* set = new MeshSet; - set->setType(SetType::Node); - set->setName("Point"); - MeshSet* s = _idsetList.value(setid); - if(s != nullptr) { - delete s; - _idsetList.remove(setid); - } - // _idsetList[setid] = set; - - // vtkIdTypeArray* idlist = vtkIdTypeArray::New(); - double coordinate[3] = { 0.0, 0.0, 0.0 }; - // fluent文件会多出一行( - if(mMeshType == typeFluent) { - QString line = this->readLine(); - } - // - for(int i = startIndex; i <= endIndex; i++) { - if(!_threadRuning) - return; - // idlist->InsertNextValue(i); - set->appendTempMem(i); - - QString scoor = this->readLine(); - QStringList scoors = scoor.split(" "); - // qDebug() << i << endl; - coordinate[0] = scoors.at(0).toDouble(); - coordinate[1] = scoors.at(1).toDouble(); - if(_ND == 3) - coordinate[2] = scoors.at(2).toDouble(); - points->InsertNextPoint(coordinate); - } - // set->setIDList(idlist); - } - } - - void MSHdataExchange::readPoints130(vtkUnstructuredGrid* dataset, QString info) - { - vtkSmartPointer points = dataset->GetPoints(); - if(points == nullptr) { - points = vtkSmartPointer::New(); - dataset->SetPoints(points); - } - QStringList infos = info.split(" "); - const int setid = infos.at(1).toInt(); - const int startIndex = infos.at(2).toInt(); - const int endIndex = infos.at(3).toInt(); - _ND = infos.at(5).toInt(); - MeshSet* set = new MeshSet; - set->setType(SetType::Node); - set->setName("001"); - MeshSet* s = _idsetList.value(setid); - if(s != nullptr) { - delete s; - _idsetList.remove(setid); - } - _idsetList[setid] = set; - - // vtkIdTypeArray* idlist = vtkIdTypeArray::New(); - double coordinate[3] = { 0.0, 0.0, 0.0 }; - for(int i = startIndex; i <= endIndex; ++i) { - if(!_threadRuning) - return; - // idlist->InsertNextValue(i); - set->appendTempMem(i); - QString scoor = this->readLine(); - QStringList scoors = scoor.split(" "); - coordinate[0] = scoors.at(0).toDouble(); - coordinate[1] = scoors.at(1).toDouble(); - if(_ND == 3) - coordinate[2] = scoors.at(2).toDouble(); - points->InsertNextPoint(coordinate); - } - // set->setIDList(idlist); - } - - void MSHdataExchange::readCells120(vtkUnstructuredGrid* dataset, QString info) - { - QStringList infos = info.split(" "); - const int setid = infos.at(1).toInt(); - const int startIndex = infos.at(4).toInt(); - const int endIndex = infos.at(5).toInt(); - MeshSet* set = new MeshSet; - set->setType(SetType::Element); - set->setName("002"); - MeshSet* s = _idsetList.value(setid); - if(s != nullptr) { - delete s; - _idsetList.remove(setid); - } - _idsetList[setid] = set; - - // vtkIdTypeArray* idlist = vtkIdTypeArray::New(); - double coordinate[3] = { 0.0, 0.0, 0.0 }; - for(int i = startIndex; i <= endIndex; ++i) { - if(!_threadRuning) - return; - // idlist->InsertNextValue(i); - set->appendTempMem(i); - QString sele = this->readLine(); - QStringList selelist = sele.split(" "); - QList member; - vtkSmartPointer idlist = vtkSmartPointer::New(); - for(int i = 0; i < selelist.size(); ++i) { - if(!_threadRuning) - return; - int d = selelist.at(i).toInt(); - if(!member.contains(d)) { - member.append(d); - idlist->InsertNextId(d); - } - } - if(member.size() == 3) { - dataset->InsertNextCell(VTK_TRIANGLE, idlist); - } else if(member.size() == 4) { - dataset->InsertNextCell(VTK_QUAD, idlist); - } - } - // set->setIDList(idlist); - } - - void MSHdataExchange::readFluentCells12(vtkUnstructuredGrid* dataset, QString info) - { - QStringList infos = info.split(" "); - const int setid = strToInt(infos.at(1), true); - const int startIndex = strToInt(infos.at(2), true); - const int endIndex = strToInt(infos.at(3), true); - MeshSet* set = new MeshSet; - set->setType(SetType::Element); - set->setName("002"); - MeshSet* s = _idsetList.value(setid); - if(s != nullptr) { - delete s; - _idsetList.remove(setid); - } - _idsetList[setid] = set; - - // vtkIdTypeArray* idlist = vtkIdTypeArray::New(); - double coordinate[3] = { 0.0, 0.0, 0.0 }; - for(int i = startIndex; i <= endIndex; ++i) { - // idlist->InsertNextValue(i); - set->appendTempMem(i); - QString sele = this->readLine(); - QStringList selelist = sele.split(" "); - QList member; - vtkSmartPointer idlist = vtkSmartPointer::New(); - for(int i = 0; i < selelist.size(); ++i) { - int d = selelist.at(i).toInt(); - if(!member.contains(d)) { - member.append(d); - idlist->InsertNextId(d); - } - } - if(member.size() == 3) { - dataset->InsertNextCell(VTK_TRIANGLE, idlist); - } else if(member.size() == 4) { - dataset->InsertNextCell(VTK_QUAD, idlist); - } - } - // set->setIDList(idlist); - } - - void MSHdataExchange::readZone45(QString zone) - { - // qDebug() << zone; - QStringList zoneinfo = zone.split(" "); - const int setid = zoneinfo.at(1).toInt(); - MeshSet* set = _idsetList.value(setid); - if(set != nullptr) - set->setName(zoneinfo.at(3)); - } - - bool MSHdataExchange::readFace13(vtkUnstructuredGrid* dataset, QString info) - { - int index; - - int bctype; - int facenumber; - int facetype; - int elemtype; - int startIndex; - int endIndex; - int setid; - char* ch; - // qDebug() << info; - QStringList infos = info.split(" "); - - setid = strToInt(infos.at(1), true); - - startIndex = strToInt(infos.at(2), true); - endIndex = strToInt(infos.at(3), true); - bctype = strToInt(infos.at(4), true); - facetype = strToInt(infos.at(5)); - - MeshSet* set = new MeshSet; - set->setType(SetType::Element); - set->setName("002"); - MeshSet* s = _idsetList.value(setid); - if(s != nullptr) { - delete s; - _idsetList.remove(setid); - } - _idsetList[setid] = set; - // 获取面的数量 - facenumber = endIndex - startIndex + 1; - // vtkIdTypeArray* idlist = vtkIdTypeArray::New(); - double coordinate[3] = { 0.0, 0.0, 0.0 }; - for(int i = 0; i < facenumber; ++i) { - if(!_threadRuning) - return false; - QList member; - // idlist->InsertNextValue(_staticid); - set->appendTempMem(_staticid); - _staticid++; - QString sele = this->readLine(); - QStringList selelist = sele.split(" "); - int firstFaceType = 0; - if(mMeshType == typeGambit) { - // 取每一个面对应的类型，不是头中的类型 - if(selelist.count() > 1) { - firstFaceType = strToInt(selelist.at(0)); - } - } else if(mMeshType == typeFluent) // fluent类型文件，面中没有类型，这里取的是头中的类型 - { - firstFaceType = facetype; - } - - vtkSmartPointer pointIdList = vtkSmartPointer::New(); - // 第一个字符串是单元类型，不需要保存,最后两个字符串不需要保存，一个是cellid,另外一个为0 - // 4 1 77 32b0 29 0 1 - int j = 0; - if(mMeshType == typeFluent) { - j = 0; - } else if(mMeshType == typeGambit) { - j = 1; - } - for(j; j <= selelist.count() - 3; j++) { - if(!_threadRuning) - return false; - int pointIndex = strToInt(selelist.at(j), true); - if(pointIndex > 0) { - pointIndex = pointIndex - 1; - if(!member.contains(pointIndex)) { - member.append(pointIndex); - pointIdList->InsertNextId(pointIndex); - } - } - } - - if(firstFaceType == 3) { - dataset->InsertNextCell(VTK_TRIANGLE, pointIdList); - } else if(firstFaceType == 4) { - dataset->InsertNextCell(VTK_QUAD, pointIdList); - } else if(firstFaceType == 2) { - dataset->InsertNextCell(VTK_LINE, pointIdList); - } else if(firstFaceType == 5) { - dataset->InsertNextCell(VTK_POLYGON, pointIdList); - } - } - // set->setIDList(idlist); - return true; - } - - QString MSHdataExchange::getNextLineInfo(QString line) - { - QString nextLine; - nextLine = line.remove("(").remove(")"); - return nextLine; - } - - bool MSHdataExchange::isHex(QString line) - { - bool isHex = false; - for(int i = 0; i < line.length(); i++) { - if(line.at(i) > '9') { - isHex = true; - return isHex; - } - } - return isHex; - } - - int MSHdataExchange::strToInt(QString line, bool isHex) - { - if(isHex) { - return line.toInt(0, 16); - } else { - return line.toInt(); - } - } - - bool MSHdataExchange::readGambitFile(vtkUnstructuredGrid* dataset) - { - QString line; - do { - if(!_threadRuning) - return false; - line = this->readLine(); - char* ch; - QByteArray ba = line.toLatin1(); - ch = ba.data(); - if(line.endsWith(")(")) { - line = line.remove("(").remove(")"); - if(line.isEmpty()) - continue; - QStringList slist = line.split(" "); - - int index; - // 取出数字标识 - sscanf(ch, "(%d", &index); - if(slist.at(0) == "0") { - _describe = slist.at(1); - _describe.remove("\""); - } - - if(index == 130) { - readPoints130(dataset, line); - } else if(slist.at(0) == "120") { - readCells120(dataset, line); - } - // 10为GAMBIT to Fluent File - else if(index == 10) { - // 读取16进制 - readPoints10(dataset, line); - } - // 读取Face - if(index == 13) { - //(13(0 1 1a7c9 0)) - // 13如果为13标志，在13后插入空格，和其他格式对应上 - QString replaceLine; - replaceLine = line; - replaceLine.insert(2, " "); - readFace13(dataset, replaceLine); - } - } else if(line.endsWith("))") && line.startsWith("(")) { - line = line.remove("(").remove(")"); - if(line.isEmpty()) - continue; - QStringList slist = line.split(" "); - if(slist.at(0) == "45") - readZone45(line); - } - - } while(!_stream->atEnd()); - _file->close(); - QList setList = _idsetList.values(); - if(dataset != nullptr) { - MeshKernal* k = new MeshKernal; - k->setName(_baseFileName); - k->setPath(_filePath); - k->setMeshData((vtkDataSet*)dataset); - _meshData->appendMeshKernal(k); - int kid = k->getID(); - - for(int i = 0; i < setList.size(); ++i) { - if(!_threadRuning) - return false; - MeshSet* set = setList.at(i); - // set->setDataSet(dataset); - set->setKeneralID(kid); - _meshData->appendMeshSet(set); - } - return true; - } - - for(int i = 0; i < setList.size(); ++i) { - MeshSet* set = setList.at(i); - delete set; - } - _idsetList.clear(); - return false; - } - - bool MSHdataExchange::readFluentFile(vtkUnstructuredGrid* dataset) - { - QString line; - do { - if(!_threadRuning) - return false; - line = this->readLine(); - char* ch; - QByteArray ba = line.toLatin1(); - ch = ba.data(); - line = line.remove("(").remove(")"); - if(line.isEmpty()) - continue; - QStringList slist = line.split(" "); - - int startIndex = 0; - int index; - // 取出数字标识 - sscanf(ch, "(%d", &index); - if(slist.at(0) == "0") { - _describe = slist.at(1); - _describe.remove("\""); - } - // 10为GAMBIT to Fluent File - else if(index == 10) { - startIndex = strToInt(slist.at(1), true); - // 如果startindex<1应该是区的标识，不读取 - if(startIndex < 1) - continue; - // 读取16进制 - readPoints10(dataset, line); - } - // 读取Face - else if(index == 13) { - // face面总的信息不读取 - startIndex = strToInt(slist.at(1), true); - if(startIndex < 1) { - continue; - } - readFace13(dataset, line); - } - // 读取区域名称，包含边界 - if(index == 39) { - QStringList slist = line.split(" "); - readZone45(line); - } - } while(!_stream->atEnd()); - _file->close(); - QList setList = _idsetList.values(); - if(dataset != nullptr) { - MeshKernal* k = new MeshKernal; - k->setName(_baseFileName); - k->setPath(_filePath); - k->setMeshData((vtkDataSet*)dataset); - _meshData->appendMeshKernal(k); - - int kid = k->getID(); - - for(int i = 0; i < setList.size(); ++i) { - MeshSet* set = setList.at(i); - // set->setDataSet(dataset); - set->setKeneralID(kid); - _meshData->appendMeshSet(set); - } - return true; - } - - for(int i = 0; i < setList.size(); ++i) { - MeshSet* set = setList.at(i); - delete set; - } - _idsetList.clear(); - return false; - } - bool CFDMshFileParser::setFile(const std::string& file) { m_File = file; diff --git a/src/PluginMeshDataExchange/MSHdataExchange.h b/src/PluginMeshDataExchange/MSHdataExchange.h index 76955b1..0eb9370 100644 --- a/src/PluginMeshDataExchange/MSHdataExchange.h +++ b/src/PluginMeshDataExchange/MSHdataExchange.h @@ -49,33 +49,12 @@ namespace MeshData ~MSHdataExchange(); bool read(); - QString readLine(); - //读取文件头判断是GAMBIT to Flent文件还是Flent生成文件 - bool readHeader(); bool write(); void run() ; private: meshType mMeshType{ typeNone }; int _totalNumber; - //读取16进制点 - void readPoints10(vtkUnstructuredGrid* dataset, QString info); - void readPoints130(vtkUnstructuredGrid* dataset, QString info); - void readCells120(vtkUnstructuredGrid* dataset, QString info); - void readFluentCells12(vtkUnstructuredGrid* dataset, QString info); - void readZone45(QString zone); - //wzy添加用于读取面 - bool readFace13(vtkUnstructuredGrid* dataset, QString info); - //去掉头标志下一处中的括号等信息 - QString getNextLineInfo(QString line); - //判断头信息中是否有16进制字符串 - bool isHex(QString line); - //判断输入字符串，是16进制还是10进制，调用不同转换方法 - int strToInt(QString line,bool isHex=false); - //解析Gambit生成fluent文件 - bool readGambitFile(vtkUnstructuredGrid* dataset); - //解析Flent文件 - bool readFluentFile(vtkUnstructuredGrid* dataset); private: QFile *_file; @@ -90,7 +69,7 @@ namespace MeshData int _ND{ 3 }; int _staticid{ 0 }; int _modelId; - MeshOperation _operation; + MeshOperation _operation; }; /** * @brief msh文件解析类。解析后所得信息中索引均以1开始，使用vtk渲染时注意索引的转换。