OSG开发笔记（三十）：OSG加载动力学仿真K模型文件以及测试Demo

已于 2024-09-18 19:42:42 修改
阅读量902 收藏 11
点赞数 26
分类专栏： Qt开发 文章标签： k文件 osg打开k文件 .k文件 三维重建
于 2024-09-18 19:26:43 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq21497936/article/details/142340138
版权

​
若该文为原创文章，未经允许不得转载
本文章博客地址：https://blog.csdn.net/qq21497936/article/details/142340138
各位读者，知识无穷而人力有穷，要么改需求，要么找专业人士，要么自己研究

长沙红胖子Qt（长沙创微智科）博文大全：开发技术集合（包含Qt实用技术、树莓派、三维、OpenCV、OpenGL、ffmpeg、OSG、单片机、软硬结合等等）持续更新中…

OSG开发专栏（点击传送门）

上一篇：《OSG开发笔记（二十九）：OSG加载模型文件、加载3DMax三维型文件Demo》
下一篇： 持续补充中…


前言

  Osg需要打开模型文件，但是遇到显示动力学仿真的K模型文件，.k文件是一种描述材料属性的文件，比如密度、弹性模量等，该模型文件不是常规中间开放格式，无法直接支持，需要自定义解析并且重建三维模型。


Demo

  实际非常流程，因为视频转gif导致部分看起来不行：
  请添加图片描述

  请添加图片描述

  请添加图片描述


交互流畅性测试

  实际研发需要用不同的策略进行在不影响视觉观感的前提下，进行模型展示优化处理，本次测试直接copy多份模型加载：

  • 172万多个四边形，丝滑
      在这里插入图片描述

  • 344万多个四边形，丝滑
      在这里插入图片描述

  • 518万多个四边形，开始稍微有点卡顿了
      在这里插入图片描述

  • 691万多个四边形，开始更不利索了
      此时翻倍就在加载场景崩溃，这个崩溃确认代码没问题，是运行内存较大了，加大运行时内存配置：

QMAKE_LFLAGS_WINDOWS += /LARGEADDRESSAWARE

  运行成功，交互更卡一点。
  在这里插入图片描述

  更大又不行了：
  在这里插入图片描述


注意

  分析的是本文件解析，可能包含数据格式关键字不全。
  在这里插入图片描述

  实际遇到的为：*TITLE，*PART，*ELEMENT_SHELL，*NODE， *END。
  （PS：以后项目上遇到了再另起篇章补充。）


K文件

概览

  .k文件是一种描述材料属性的文件，比如密度、弹性模量等。
  打开其的软件LS-DYNA和LSPREPOST

  • LS-DYNA
      LS-DYNA程序（最新版本17.2版）是功能齐全的几何非线性（大位移、大转动和大应变）、材料非线性（140多种材料动态模型）和接触非线性（50多种）程序。它以Lagrange算法为主，兼有ALE和Euler算法；以显式求解为主，兼有隐式求解功能；以结构分析为主，兼有热分析、流体-结构耦合功能；以非线性动力分析为主，兼有静力分析功能（如动力分析前的预应力计算和薄板冲压成型后的回弹计算）；军用和民用相结合的通用结构分析非线性有限元程序，是显式动力学程序的鼻祖和先驱。
  • LSPREPOST
      lsprepost是一套专为LS-DYNA开发的高级有限元前后处理软件，上线时间是1976年1月1日。
      文件内容概览
      文件整体概览如下：
      在这里插入图片描述

  以上所知，分为六个关键字，分割为五个区域，归属关系如下：
  在这里插入图片描述

关键字：*KEYWORK

关键字：*TITLE

关键字：*PART

  PART关键字来定义部件及其相关特性
  在这里插入图片描述

  第一行为title属性，为标题，没有实际意义。
  在这里插入图片描述

关键字：*ELEMENT_SHELL

  SECTION系列关键字指定所采用的单元算法、积分分规则及各种几何参数(壳单元的厚度、梁单元的截面信息参数等)，对不同的单元类型，或者需要相应的SECTION关键字段来定义单元的算法和参数。
  *SECTION_SEHLL定义单壳体单元的特性。这些属性对于模拟如汽车碰撞、航空航天结构、以及任何涉及薄壳结构动态响应的模拟都至关重要。然而，具体到ELEMENT_SHELL属性在K文件中的详细内容和格式，可能因LS-DYNA的版本和具体应用而有所不同。
  在这里插入图片描述

  在这里插入图片描述

关键字：*NODE

  *NODE关键字来定义结构模型的每一节点及其在总体直角坐标系中的坐标及受约束情况，该关键字段包含如下的变量信息：
  在这里插入图片描述
  在这里插入图片描述

  然后部件与node的对应关系，通过索引对应：
  在这里插入图片描述

  解析文件即可。

关键字：*END

  表示文件结束。


测试加载K文件

  在这里插入图片描述

  在这里插入图片描述


Demo关键源码

模型结构体

// 模型结构体
    struct Element_Shell    // *ELEMENT_SHELL
    {
        Element_Shell() {
        }
        qint64 eid;         // 单元id
        qint64 pid;         // 材料id
        qint64 n1;          // 节点1，定义几何形状
        qint64 n2;          // 节点2，定义几何形状
        qint64 n3;          // 节点3，定义几何形状
        qint64 n4;          // 节点4，定义几何形状
        qint64 n5;          // 厚度，额外的节点在标准的LS-DYNA四边形壳单元定义中是没有意义的。
        qint64 n6;          // 积分点数，额外的节点在标准的LS-DYNA四边形壳单元定义中是没有意义的。
        qint64 n7;          // 额外的节点在标准的LS-DYNA四边形壳单元定义中是没有意义的。
        qint64 n8;          // 额外的节点在标准的LS-DYNA四边形壳单元定义中是没有意义的。
    };
    struct Part             // *PART
    {
        Part() {
        }
        qint64 pid;         // 部件的id号，唯一
        qint64 secid;       // 有*section关键字定义的section的id号
        QList<Element_Shell> listElementShell;  // 部件片元
        qint64 mid;         // 部件的材料号
        qint64 eosid;       // 部件所属材料涉及的状态方程号，由*EOS关键字定义
        qint64 hgid;        // 沙漏或体积粘性参数编号，由*HOURGLASS关键字定义，取0表示将采用默认的数值：
        qint64 grav;        // 仅对实体单元有效，取0表示对所有PART进行重力初始化，取1表示仅对当前材料初始化
        qint64 adpopt;      // 标识该部件是否采用自适应网格划分，取0表示不采用
        qint64 tmid;        // 标识该部件是否采用自适应网格划分，取0表示不采用
    };
    struct Node {
        Node() {
        }
        qint64 nid;         // 结点号，唯一
        double x;           // 三维x坐标（全局）
        double y;           // 三维y坐标（全局）
        double z;           // 三维z坐标（全局）
        int tc;             // 平动自由度受约束状态，枚举值：0-无平动约束,1-X方向平动约束,2-Y方向平动约束
        int rc;             // 转动自由度收约束状态，枚举值：0-无转动约束,1-X方向转动约束,2-Y方向转动约束
    };
    struct K_Mode
    {
        K_Mode() {}
        QList<Part> listPart;
        QList<Node> listNode;
        QHash<int, Node> hashNid2Node;
    };

绘制部分

osg::ref_ptr<osg::Group> pGroup = new osg::Group;
// 绘图
{
    for(int partIndex = 0; partIndex < kMode.listPart.size(); partIndex++)
    {
        // 创建一个用户保存几何信息的对象
        osg::ref_ptr<osg::Geometry> pGeometry = new osg::Geometry;
        // 创建四个顶点的数组
        osg::ref_ptr<osg::Vec3Array> pVec3Array = new osg::Vec3Array;
        // 添加四个顶点
        pGeometry->setVertexArray(pVec3Array.get());

        // 创建四种颜色的数据
        osg::ref_ptr<osg::Vec4Array> pVec4Array = new osg::Vec4Array;
        // 添加四种颜色
        pGeometry->setColorArray(pVec4Array.get());
        // 绑定颜色
        pGeometry->setColorBinding(osg::Geometry::BIND_PER_VERTEX);

        double r, g, b;
        r = qrand() % 100 * 1.0f / 100;
        g = qrand() % 100 * 1.0f / 100;
        b = qrand() % 100 * 1.0f / 100;
        for(int elementShellIndex = 0; elementShellIndex < kMode.listPart.at(partIndex).listElementShell.size(); elementShellIndex++)
        {
            //                               x     y     z
#if 0
            pVec3Array->push_back(osg::Vec3( 1.0, 0.0, 0.0));
            pVec3Array->push_back(osg::Vec3( 3.0, 0.0, 0.0));
            pVec3Array->push_back(osg::Vec3( 3.0, 0.0, 1.0));
            pVec3Array->push_back(osg::Vec3( 1.0, 0.0, 1.0));
#endif
            pVec3Array->push_back(osg::Vec3(kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n1).x,
                                            kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n1).y,
                                            kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n1).z));
            pVec3Array->push_back(osg::Vec3(kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n2).x,
                                            kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n2).y,
                                            kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n2).z));
            pVec3Array->push_back(osg::Vec3(kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n3).x,
                                            kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n3).y,
                                            kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n3).z));
            pVec3Array->push_back(osg::Vec3(kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n4).x,
                                            kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n4).y,
                                            kMode.hashNid2Node.value(kMode.listPart.at(partIndex).listElementShell.at(elementShellIndex).n4).z));


            //                               r    g    b    a(a设置无效，估计需要其他属性配合)
            pVec4Array->push_back(osg::Vec4(r, g, b, 1.0));
            pVec4Array->push_back(osg::Vec4(r, g, b, 1.0));
            pVec4Array->push_back(osg::Vec4(r, g, b, 1.0));
            pVec4Array->push_back(osg::Vec4(r, g, b, 1.0));

        }
        // 注意：此处若不绑定画笔，则表示使用之前绑定的画笔

        // 为唯一的法线创建一个数组    法线: normal
        osg::ref_ptr<osg::Vec3Array> pVec3ArrayNormal = new osg::Vec3Array;
        pGeometry->setNormalArray(pVec3ArrayNormal.get());
        pGeometry->setNormalBinding(osg::Geometry::BIND_OVERALL);
        pVec3ArrayNormal->push_back(osg::Vec3(0.0, -1.0, 0.0));

        // 由保存的数据绘制四个顶点的多边形
        pGeometry->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::QUADS, 0, kMode.listPart.at(partIndex).listElementShell.size() * 4));
//            pGeometry->addPrimitiveSet(new osg::DrawArrays(osg::PrimitiveSet::QUADS, 0, 4));

        // 向Geode类添加几何体(Drawable)
        osg::ref_ptr<osg::Geode> pGeode = new osg::Geode;
        pGeode->addDrawable(pGeometry.get());

        {
            osg::ref_ptr<osg::StateSet> pStateSet = pGeometry->getOrCreateStateSet();
            osg::ref_ptr<osg::PolygonMode> pPolygonMode = new osg::PolygonMode(osg::PolygonMode::FRONT_AND_BACK, osg::PolygonMode::LINE);
            pStateSet->setAttribute(pPolygonMode);
        }

        pGroup->addChild(pGeode);
    }
}

工程模板：对应版本号v1.32.0

  在这里插入图片描述


入坑

入坑一：绑定颜色不出来

问题

  去掉纹理后颜色不出来，显示一直灰色。

原理

  开启光照，就颜色不出来，需要关闭

解决

// 始终是灰色，这里需要设置关闭光照：OFF
osg::StateSet *pStateSet = pGroup->getOrCreateStateSet();
pStateSet->setMode(GL_LIGHTING, osg::StateAttribute::OFF);
//pStateSet->setMode(GL_LIGHTING, osg::StateAttribute::ON);

  在这里插入图片描述

  在这里插入图片描述

  
在这里插入图片描述

入坑二：运行起来就崩溃

问题

  数据都加载绘制时，运行就崩溃

解决过程

  • 步骤一：先减少面数量，发现小于一档数量可以没问题，初步怀疑绘制太多，但是osg本身做深度测试相关很不错，除非时内存爆了，所以应然需要继续。
  • 步骤二：使用绘制四边形的，绘制10000个也会崩，加重了内存怀疑。
  • 步骤三：由于之前对于顶点颜色没有进行全部点绑定，其他的看纹理的方式都是重复的，所以这块就直接只设置了4个点，但本次去掉颜色之后，发现是可以的。
  • 步骤四：将颜色给他赋值上，一个部件一个纯色，把对应绑定点都补上，发现没问题

原因

  就是因为绑定颜色角点，只绑定了四个，按道理来说就算其越界也只是读取的乱值，颜色是乱的，但是点多了就不行了，所以这里还是需要一一对应，不偷懒。

解决

  颜色单独，每一个部件一个颜色，然后循环内一一对应：
  在这里插入图片描述


上一篇：《OSG开发笔记（二十九）：OSG加载模型文件、加载3DMax三维型文件Demo》
下一篇： 持续补充中…


本文章博客地址：https://blog.csdn.net/qq21497936/article/details/142340138

OSG开发笔记（二十八）：OSG模型固定路径动画
08-20 3万+
若该文为原创文章，未经允许不得转载 原博主博客地址：https://blog.csdn.net/qq21497936 本文章博客地址：https://blog.csdn.net/qq21497936/article/details/99816214 目录 前言 Demo效果 路径动画 动画关键点 动画路径类AnimationPath 飞机模型 飞机模型动画 模型路径动画构架图 ...
OSG开发笔记（八）：OSG模型文件存储与读取
热门推荐
07-10 4万+
若该文为原创文章，未经允许不得转载 原博主博客地址：https://blog.csdn.net/qq21497936 本文章博客地址：https://blog.csdn.net/qq21497936/article/details/95303864 目录 前言 读取文件功能 存入文件功能 演示 步骤一：代码存入1.osg 步骤二：存储文件 步骤三：使用打开文件加载 步骤四(补充)...
OSG开发笔记（二十九）：OSG加载模型文件、加载3DMax三维型文件Demo
10-27 6847
Osg深入之后需要打开模型文件，这些模型文件是已有的模型文件，加载入osg之后可以在常见中展示模型文件，该节点可以操作，多个逼真的模型的节点就实现了基本的场景构建。自定义文件插件主要是自定义一个插件读写类，继承osgDB::ReaderWriter类，然后根据需求重写如readNode等函数方法即可。
OSG —— 笔记2 - 加载模型（附源码）
wx微信公众号“码农总动员”，获取更多编程知识
03-23 2403
效果     关注 笔者 - jxd 微信公众号搜索 “码农总动员” 或 微信扫描下方二维码，了解更多你不知道的XX，O(∩_∩)O
OSG开发笔记（二十三）：Qt使用QOpenGLWidget渲染OSG和地球仪
08-07 4万+
若该文为原创文章，未经允许不得转载 原博主博客地址：https://blog.csdn.net/qq21497936 本文章博客地址：https://blog.csdn.net/qq21497936/article/details/98730328 目录 前言 地球控件 运行效果 运行Demo下载地址 QGLWidget与QOpenGLWidget详解 QOpenGLWidget放置...
OSG开发笔记（十四）：OSG交互
07-18 3万+
若该文为原创文章，未经允许不得转载 原博主博客地址：https://blog.csdn.net/qq21497936 本文章博客地址：https://blog.csdn.net/qq21497936/article/details/96422466 目录 前言 交互概述 交互实现方式与步骤 交互动作 模型多种类型动作的理解 代码 地球代码 事件代码 操作过程 工程模板：对应版...
OSG开发笔记（十）：OSG模型的变换之平移、旋转和缩放
07-13 3万+
若该文为原创文章，未经允许不得转载 原博主博客地址：https://blog.csdn.net/qq21497936 本文章博客地址：https://blog.csdn.net/qq21497936/article/details/95725133 目录 前言 变换模型的添加与删除 添加变换模型 删除变换模型 平移 缩放 旋转 扩展：OSG中正确的矩阵变换级联顺序 操作过程 ...
OSG开发笔记（十五）：OSG光照
07-19 3万+
该文为原创文章，未经允许不得转载 原博主博客地址：https://blog.csdn.net/qq21497936 本文章博客地址：https://blog.csdn.net/qq21497936/article/details/96473289 目录 前言 OSG光照 地球光照效果 设置光照后的效果 白光 红光 光照实现方式与步骤 光照代码 入坑 入坑一：光照后的效果是黑色...
OSG开发笔记（十三）：OSG三维纹理映射(体渲染)
07-17 3万+
若该文为原创文章，未经允许不得转载 原博主博客地址：https://blog.csdn.net/qq21497936 本文章博客地址：https://blog.csdn.net/qq21497936/article/details/96285960 目录 前言 代码 总结步骤 步骤一：绘制几何类型 步骤二：加载图片添加纹理 步骤三：渲染设置，结点完成 入坑 入坑一：2D纹理贴图未...
OSG数据加载性能优化：纹理压缩
欧特GO
07-06 3105
1、理论知识 我们都知道，纹理的存储是RGBA，这四个数字假如不压缩的话其占用的存储空间是相当大的，每个通道按0~255 一个字节来计算，一个像素RGBA是4个字节（byte)，一个1024x1024这么大的图片会占用1024x1024x4=4M。如果不经压缩，一个GPU显存在4GB的显卡，也只能存储1024张1024x1024的图片。这显然不合常理。有时候我们在处理tiff图像时，指定不使用压缩，则随便一个大小出来就是几十个GB，那就是什么显卡也存不了了。 所以无论是CPU还是GPU都支持对纹理进行压缩
OSG载入地形和模型文件
09-20
在OpenSceneGraph (OSG) 中，加载地形和模型文件是一项关键任务，它涉及图形学、数据处理和用户交互等多个方面。OpenSceneGraph是一个高性能的3D图形库，广泛用于实时三维可视化应用，如游戏开发、虚拟现实和科学...
osg模型文件精简工具V1.00_osg简化_osg格式文件_模型简化_osg_图像化界面_
10-03
**osg模型文件精简工具V1.00** 在3D建模和渲染领域，OpenSceneGraph（OSG）是一种高性能的开源图形库，用于创建复杂的实时3D应用程序。OSG支持多种3D模型格式，其中包括它自己的专有格式。然而，大型3D模型在处理和...
report ETL .ffff
09-17
report ETL .ffff
基于ssm的共享充电宝管理系统设计与实现.docx
09-17
基于ssm的共享充电宝管理系统设计与实现.docx
操作系统大作业_Linux_命令行_OSWork.zip
最新发布
09-17
操作系统大作业_Linux_命令行_OSWork
Java项目-基于SSM+Jsp的社区管理与服务系统的设计与实现（源码+数据库脚本+部署视频+代码讲解视频+全套软件）
09-17
【基于SSM+Jsp的社区管理与服务系统的设计与实现】高分通过项目，已获导师指导。 本项目是一套基于SSM+Jsp的社区管理与服务系统，主要针对计算机相关专业的正在做毕设的学生和需要项目实战练习的Java学习者。也可作为课程设计、期末大作业 包含：项目源码、数据库脚本、开发说明文档、部署视频、代码讲解视频、全套软件等，该项目可以直接作为毕设使用。 项目都经过严格调试，确保可以运行！
吉他谱_Show Me How To Live - Audioslave.pdf
09-17
初级入门吉他谱 guitar tab
Linux系统初始化脚本，包括设置时区并同步时间、禁用selinux、关闭防火墙
09-17
Linux系统初始化脚本，包括设置时区并同步时间、禁用selinux、关闭防火墙、历史命令显示操作时_xiaopenglinux_shell02
基于ssm的即动运动网站设计与实现设计.docx
09-17
基于ssm的即动运动网站设计与实现设计.docx
ffmpeg windows
07-27
FFmpeg 是一个开源的跨平台音视频处理工具，可以在 Windows 上使用。你可以按照以下步骤在 Windows 上安装和使用 FFmpeg： 1. 访问 FFmpeg 官方网站：https://www.ffmpeg.org/ ，点击 "Download" 进入下载页面。 2. 在下载页面中，找到 "Windows" 标题下的链接，选择 "Static Builds"。 3. 在 "Static Builds" 页面中，你可以选择根据你的系统架构下载对应的版本。如果你的 Windows 是 64 位的，点击 "64-bit" 下载；如果是 32 位的，点击 "32-bit" 下载。 4. 下载完成后，解压缩下载的压缩包到一个你喜欢的目录。 5. 打开解压缩后的文件夹，在 "bin" 文件夹下可以找到可执行文件 "ffmpeg.exe"。 6. 可以将该文件夹添加到系统的环境变量中，这样就可以在任意目录下使用 FFmpeg。 7. 现在你可以在命令行或者其他应用程序中直接使用 FFmpeg 命令进行音视频处理了。 请注意，FFmpeg 的命令行非常强大和灵活，你可以参考官方文档或者其他教程来学习如何使用 FFmpeg 进行各种音视频处理操作。
评论 3
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

长沙红胖子Qt（长沙创微智科）

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值