自研有限元软件与ANSYS精度对比-Bar3D2Node三维杆单元模型-央视大裤衩实例

目录

1、“央视大裤衩”自研有限元软件求解

1.1、选择单元类型

1.2、导入“央视大裤衩”工程

1.3、节点坐标定义

1.4、单元连接关系、材料定义

1.5、约束定义

1.6、外载定义

1.7、矩阵求解

1.8、变形云图展示

1.9、节点位移

1.10、单元应力

1.11、节点支反力

2、“央视大裤衩”Ansys软件求解

2.1、变形云图

2.2、节点位移

2.3、单元应力

2.4、节点支反力

3、结果对比

3.1、节点位移对比

3.2、单元应力对比

3.3、节点支反力对比


本系列文章致力于实现“手搓有限元,干翻Ansys的目标”,基本框架为前端显示使用QT实现交互,后端计算采用Visual Studio C++。

基于前文的准备工作,已具备初步的有限元分析的全流程计算,为了验证计算的准确性,将自研的有限元软件与商业软件ANSYS的计算结果进行对比。

手搓有限元 干翻Ansys-央视大裤衩有限元分析

(也可以点击此处进行视频跳转)

1、“央视大裤衩”自研有限元软件求解

使用自研有限元软件对“央视大裤衩”进行求解。模型中所有单元的杨氏模量E = 7e10Pa,横截面积A = 2e-4m^2。模型状态如下图所示(这里借用了Ansys的模型展示状态):

在模型底端全部使用铰支座进行X、Y、Z三个方向固定,除与地面接触的节点外,所有节点受到X正方向大小为1000N的力,模拟为风荷载。

1.1、选择单元类型

这里选择Bar_3D2Node单元(点击这里跳转到Bar3D2Node介绍)。

1.2、导入“央视大裤衩”工程

这里导入提前准备好的“央视大裤衩”模型。

1.3、节点坐标定义

这里显示的是模型中节点的ID、X坐标、Y坐标、Z坐标。

1.4、单元连接关系、材料定义

这里显示的是单元的ID、起始点ID、结束点ID、杨氏模量E、横截面积A。

1.5、约束定义

这里显示的是模型节点X方向是否约束(0为自由,1为约束)、X方向初始位移、Y方向是否约束(0为自由,1为约束)、Y方向初始位移、Z方向是否约束(0为自由,1为约束)、Z方向初始位移。所有方向的初始位移均设置为0.

1.6、外载定义

这里显示的是模型节点X方向初始荷载、Y方向初始荷载、Z方向初始荷载,依照荷载条件,与地面接触的节点无风荷载,除此之外所有的节点均受到X方向大小为1000N的初始荷载。

1.7、矩阵求解

控制台窗口打印的是每个单元的全局刚度矩阵。

1.8、变形云图展示

这里展示的是模型的变形云图,可以调整云图的色阶以及比例缩放。

1.9、节点位移

这里展示的是不同节点的位移。

1.10、单元应力

这里展示的是不同单元的应力。

1.11、节点支反力

这里展示的是不同节点的支反力。

2、“央视大裤衩”Ansys软件求解

2.1、变形云图

2.2、节点位移

这里展示的是局部节点位移。

2.3、单元应力

这里展示的是局部单元应力,由于与地面相连节点没有加风荷载,所以到27号单元应力才不为0.

2.4、节点支反力

在前12号节点均有荷载。

3、结果对比

在自研的有限元软件与Ansys软件的计算结果进行对比。

3.1、节点位移对比

节点位移对比
节点ID(从1开始计数)自研软件节点位移(m)Ansys节点位移(m)
1

X:0.000000e+00

Y:0.000000e+00

Z:0.000000e+00

X:0.000000e+00

Y:0.000000e+00

Z:0.000000e+00

13

X:1.145608e-03

Y:-5.583721e-05

Z:-1.500025e-04

X:0.11456E-02

Y:-0.55837E-04

Z:-0.15000E-03

44

X:3.593812e-02

Y:-1.132481e-02

Z:9.641387e-03

X:0.35938E-01

Y:-0.11325E-01

Z:0.96414E-02

56

X:4.634266e-02

Y:-1.434761e-02

Z:1.284828e-02

X:0.46343E-01

Y:-0.14348E-01

Z:0.12848E-01

80

X:7.366890e-02

Y:-3.727720e-02

Z:2.816048e-02

X:0.73669E-01

Y:-0.37277E-01

Z:0.28160E-01

这里只选取了部分节点进行数据对比。可以看出计算结果几乎一致,除了舍入误差保留精度之外。

3.2、单元应力对比

单元应力对比
单元ID(从1开始计数)自研软件单元应力(Pa)Ansys单元应力(Pa)
10.000000e+00

0.0000

273.000622e+060.30006E+07
88-2.147518e+07-0.21475E+08
1254.148358e+070.41484E+08
2665.679809e+060.56798E+07

这里只选取了部分单元进行数据对比。可以看出计算结果几乎一致,除了舍入误差保留精度之外。

3.3、节点支反力对比

节点支反力对比
节点ID(从1开始计数)自研软件节点支反力(N)Ansys节点支反力(N)
1

X:-1.598386e+04

Y:1.061466e+03

Z:7.351554e+03

X:-15984.

Y:1061.5

Z:7351.6

2

X:-7.255753e+03

Y:-2.371771e+04

Z:-1.842327e+04

X:-7255.8        

Y:-23718.

Z:-18423.

3

X:-4.658990e+03

Y:-1.039396e+04

Z:-4.412918e+03

X:-4659.0        

Y:-10394.

Z:-4412.9

4

X:-1.697886e+04

Y:-1.189124e+04

Z:-6.063854e+02

X:-16979.        

Y:-11891.

Z:-606.39

5

X:-1.135751e+04

Y:-5.277917e+04

Z:-1.239206e+04

X:-11358.        

Y:-52779.

Z:-12392.

6

X:4.755060e+02

Y:3.057443e+04

Z:1.210375e+04

X:475.51          

Y: 30574.

Z:12104.

7

X:-7.595668e+03

Y:-4.415977e+04

Z:9.528992e+03

X:-7595.7          

Y:-44160.

Z:9529.0

8

X:-2.386487e+03

Y:5.215669e+04

Z:-4.339480e+03

X:-2386.5            

Y:52157.

Z:-4339.5

9

X:-3.410433e+03

Y:-2.099796e+04

Z:-4.992923e+03

X:-3410.4        

Y:-20998.

Z:-4992.9

10

X:-2.388713e+03

Y:-9.652879e+03

Z:5.519742e+03

X:-2388.7          

Y:-9652.9

Z:5519.7

11

X:3.610844e+03

Y:3.707225e+04

Z:1.442215e+04

X:3610.8            

Y:37072.

Z:14422.

12

X:-7.007064e+01

Y:5.272786e+04

Z:-3.759157e+03

X:-70.071          

Y:52728.

Z: -3759.2

这里选取了全部固定节点进行数据对比。可以看出计算结果几乎一致,除了舍入误差保留精度之外。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/964367.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

家政预约小程序12服务详情

目录 1 修改数据源2 创建页面3 搭建轮播图4 搭建基本信息5 显示服务规格6 搭建服务描述7 设置过滤条件总结 我们已经在首页、分类页面显示了服务的列表信息,当点击服务的内容时候需要显示服务的详情信息,本篇介绍一下详情页功能的搭建。 1 修改数据源 在…

Multi-Scale Heterogeneous Text-Attributed Graph Datasets From Diverse Domains

Multi-Scale Heterogeneous Text-Attributed Graph Datasets From Diverse Domains WWW25 推荐指数:#paper/⭐⭐⭐#​ 代码地址:https://github.com/Cloudy1225/HTAG 作者主页:Yunhui Lius Homepage 一句话总结:提出了涵盖多…

day31-综合架构开篇(中)

轻松理解网站架构 运维内网维护架构

使用 SpringBoot+Thymeleaf 模板引擎进行 Web 开发

目录 一、什么是 Thymeleaf 模板引擎 二、Thymeleaf 模板引擎的 Maven 坐标 三、配置 Thymeleaf 四、访问页面 五、访问静态资源 六、Thymeleaf 使用示例 七、Thymeleaf 常用属性 前言 在现代 Web 开发中,模板引擎被广泛用于将动态内容渲染到静态页面中。Thy…

Verilog基础(一):基础元素

verilog基础 我先说,看了肯定会忘,但是重要的是这个过程,我们知道了概念,知道了以后在哪里查询。语法都是术,通用的概念是术。所以如果你有相关的软件编程经验,那么其实开启这个学习之旅,你会感…

保姆级教程Docker部署Kafka官方镜像

目录 一、安装Docker及可视化工具 二、单节点部署 1、创建挂载目录 2、运行Kafka容器 3、Compose运行Kafka容器 4、查看Kafka运行状态 三、集群部署 在Kafka2.8版本之前,Kafka是强依赖于Zookeeper中间件的,这本身就很不合理,中间件依赖…

yes镜像站群/PHP驱动的镜像站群架构实践

▍当前站群运维的三大技术困局 在近期与多个IDC服务商的交流中发现,传统站群系统普遍面临: 同步效率瓶颈:跨服务器内容同步耗时超过行业标准的42%SEO权重稀释:镜像站点重复率导致70%的站点无法进入百度前3页运维成本失控&#x…

知识库管理系统提升知识利用效率与企业创新能力的全面方案分析

内容概要 知识库管理系统作为现代企业信息化建设的重要组成部分,承担着整合、存储与分发知识资源的重任。其核心功能包括信息的分类与索引、知识的共享与协作,以及数据分析与挖掘等。这些功能不仅提高了企业内部信息流动的速度和效果,还为员…

STM32F103ZET6完整技术点(持续更新~)

①STM32②F③103④Z⑤E⑥T⑦6简介: ①基于ARM核心的32位微控制器,②通用类型,③增强型,④引脚数目144个 ⑤闪存存储器容量:512K字节,⑥封装:LQFP,⑦温度范围:工业级温度范围&#xf…

html基本结构和常见元素

html5文档基本结构 <!DOCTYPE html> <html lang"zh-CN"> <head><meta charset"UTF-8"><title>文档标题</title> </head> <body>文档正文部分 </body> </html> html文档可分为文档头和文档体…

Markdown转换器中间件

目录 需求 文本编码检测 Markdown→HTML 注意 实现 需求 Markdown是一种文本格式&#xff1b;不被浏览器支持&#xff1b;编写一个在服务器端把Markdown转换为HTML的中间件。我们开发的中间件是构建在ASP.NET Core内置的StaticFiles中间件之上&#xff0c;并且在它之前运…

数据降维技术研究:Karhunen-Loève展开与快速傅里叶变换的理论基础及应用

在现代科学计算和数据分析领域&#xff0c;数据降维与压缩技术对于处理高维数据具有重要意义。本文主要探讨两种基础而重要的数学工具&#xff1a;Karhunen-Love展开&#xff08;KLE&#xff09;和快速傅里叶变换&#xff08;FFT&#xff09;。通过分析这两种方法的理论基础和应…

使用LightGlue进行图像配准并提取图像重叠区域

发表日期&#xff1a;2023年6月23日 项目地址&#xff1a;https://github.com/cvg/LightGlue https://github.com/cvg/glue-factory/ LightGlue是一个在精度上媲美Superglue&#xff0c;但在速度上比Superglue快一倍的模型。通过博主实测&#xff0c;LightGlue的配准效果比Su…

小书包:让阅读更美的二次开发之作

小书包是在一款知名阅读软件的基础上进行二次开发的产品。在保留原有软件的基本功能和用户体验的同时&#xff0c;对其界面和视觉效果进行了精心美化&#xff0c;让阅读体验更加舒适和愉悦。 内置了171条书源&#xff0c;虽然数量不算多&#xff0c;但都是作者精挑细选出来的&a…

期末数据库课程设计基于Java+MySQL+JDBC+JavaSwing实现的图书进销管理系统源代码+数据库

期末数据库课程设计&#xff0c; 图书进销信息管理系统 直接将数据库文件导入就可以快速建表 效果图 系统登录弹窗 书库管理页 信息查询页 图书销售页 系统设置页 编码&#xff1a; GBK 开发环境 jdk12MySQL8.0 推荐用IDEA运行项目 补充 UI美化&#xff08;引用当前系统…

DeepSeek最新图像模型Janus-Pro论文阅读

目录 论文总结 摘要 1. 引言 2. 方法 2.1 架构 2.2 优化的训练策略 2.4 模型扩展 3. 实验 3.1 实施细节 3.2 评估设置 3.3 与最新技术的比较 3.4 定性结果 4. 结论 论文总结 Janus-Pro是DeepSeek最新开源的图像理解生成模型&#xff0c;Janus-Pro在多模态理解和文…

深入解析“legit”的地道用法——从俚语到正式表达:Sam Altman用来形容DeepSeek: legit invigorating(真的令人振奋)

深入解析“legit”的地道用法——从俚语到正式表达 一、引言 在社交媒体、科技圈甚至日常对话中&#xff0c;我们经常会看到或听到“legit”这个词。比如最近 Sam Altman 在 X&#xff08;原 Twitter&#xff09;上发的一条帖子中写道&#xff1a; we will obviously deliver …

物联网领域的MQTT协议,优势和应用场景

MQTT&#xff08;Message Queuing Telemetry Transport&#xff09;作为轻量级发布/订阅协议&#xff0c;凭借其低带宽消耗、低功耗与高扩展性&#xff0c;已成为物联网通信的事实标准。其核心优势包括&#xff1a;基于TCP/IP的异步通信机制、支持QoS&#xff08;服务质量&…

【Hadoop】Hadoop的HDFS

这里写目录标题 HDFS概述HDFS产出背景及定义HDFS产生背景HDFS定义 HDFS优缺点HDFS优点HDFS缺点 HDFS组成架构HDFS文件块大小 HDFS的Shell操作常用命令实操准备工作上传下载HDFS直接操作 HDFS的API操作客户端环境准备HDFS的API案例实操HDFS文件上传HDFS文件下载HDFS文件更名和移…

二、面向对象

一、结构体类型 结构体类型是一种自定义类型&#xff0c;用于创建我们游戏或者实际业务中的自定义类型. 代码中变量有通用的&#xff0c;可以使用结构体&#xff0c;包裹起来。 1、成员变量 /// <summary> /// 英雄结构体 /// </summary> struct Hero {//成员p…