有限元——ANSYS求解悬臂梁均布载荷问题

最新推荐文章于 2024-01-26 19:43:39 发布
最新推荐文章于 2024-01-26 19:43:39 发布
阅读量1.2w 收藏 58
点赞数 16
分类专栏： # ANSYS 文章标签： ANSYS
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mahoon411/article/details/112299694
版权

1. 摘要

对于矩形悬臂梁均布载荷问题，分别采用不同单元模型进行分析，对比考察分析结果的差异。

2. 问题描述

在这里插入图片描述

如图所示矩形截面悬臂梁，左端固定，材料为钢。材料参数为：弹性模量E=2E11Pa，泊松比u=0.28，密度ρ=7800kg/m^3。几何参数如图中所示。梁上作用有均布载荷。
分别采用三维实体单元、平面应力单元及梁单元进行分析；对比考察梁右端的位移和左端的应力分布；讨论理论计算值与有限元计算结果差异。

3. 技术路线

此问题属于结构分析范畴，借助ANSYS Mechanical APDL模块，通过软件界面操作方式实现。单位制：m、kg、s、N、Pa。
分别采用三维实体单元Brick 8 node 185、平面应力单元Quad 4 node 182及梁单元2 node 188进行分析。

4. 分析过程

4.1 三维实体单元分析

(一) 定义工作文件名
依次单击Utility Menu>File>Change Jobname，在出现的对话框中输入：beam，在“New Log and error files”处选中“yes”，单击“OK”。
(二) 定义单元类型

  1. 依次单击 MainMenu>Preprocessor>Element type>Add/Edit/Delete，出现对话框，单击“Add”，出现一个 “Library of Element Type”对话框。在“Library of Element Type”左面的列表栏中选择“Solid”，在右面的列表栏中选择“Brick 8 node 185”，单击“OK”。
    在这里插入图片描述

(三) 设置材料属性
依次单击Main Menu>Preprocessor>Material Props>Material Models，出现“Define Material Model Behavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural>Linear>Elastic>Isotropic”,出现对话框，输入弹性模量EX＝2e11，泊松比PRXY＝0.28，单击“OK”，单击“Material>Exit”。
在这里插入图片描述
(四) 建立几何模型

  1. 创建长方体
    拾取菜单 Main Menu>Preprocessor>Modeling>Create>Volumes>Block> By Dimensions。弹出如下图所示的对话框，在“X1,X2” 文本框中输入 0,1.6，在“Y1,Y2”文本框中输入 0,0.05， 在“Z1,Z2” 文本框中输入 0,0.06，单击 OK。
    在这里插入图片描述
    在这里插入图片描述

(五) 划分单元

  1. 拾取菜单Main Menu>Mesh>Mesh Tool。弹出对话框，在 Mesh 选项中选择Volums，在 Shape 选项中选择 Hex 和 Mapped，然后点下方的 Mesh，即可划分好。
    在这里插入图片描述

(六) 施加载荷和边界条件

  1. MainMenu>Solution>DefineLoads>Apply>Structural>Displacement>On Areas,拾取左端面，然后点击OK。并在弹出的对话框里选择 All DOF，点击 ok。
    在这里插入图片描述
    在这里插入图片描述
  2. 选择 MainMenu>Solution>DefineLoads>Apply>Structural>pressure>On Areas,拾取上端面，点击 OK。在弹出的对话框中输入200000。
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

(七) 求解模型
依次点击Main Menu>Solution>Solve>Current LS进行求解，当提示Solution is Done说明求解完成。
在这里插入图片描述
(八) 结果查看
依次单击Main Menu>GeneralPostproc>PlotResults>Contour Plot>Nodal Solu，弹出对话框，在“Item to be Contoured”栏中选择“Stress”，在其右面的栏中，选择“von Mises stress”后，单击“OK”，这时在图形窗口处显示出平面的应力云图，如图所示。
在这里插入图片描述
在这里插入图片描述
依次单击 Main Menu>GeneralPostproc>Plot Results>Contour Plot>Nodal Solu，弹出对话框，在“Item to be Contoured”栏中选择“DOF Solution”，在其右面的栏中，选择“Displacement vector sum”后，单击“OK”，这时在图形窗口处显示出三维应变图，如图所示。
在这里插入图片描述
在这里插入图片描述

4.2 平面应力单元分析

(一) 定义工作文件名
依次单击Utility Menu>File>Change Jobname，在出现的对话框中输入：beam，在“New Log and error files”处选中“yes”，单击“OK”。
(二) 定义单元类型

  1. 依次单击 MainMenu>Preprocessor>Element type>Add/Edit/Delete，出现对话框，单击“Add”，出现一个 “Library of Element Type”对话框。在“Library of Element Type”左面的列表栏中选择“Solid”，在右面的列表栏中选择“Quad 4 node 182”，单击“OK”。
    在这里插入图片描述

(三) 设置材料属性
依次单击Main Menu>Preprocessor>Material Props>Material Models，出现“Define Material Model Behavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural>Linear>Elastic>Isotropic”,出现对话框，输入弹性模量EX＝2e11，泊松比PRXY＝0.28，单击“OK”，单击“Material>Exit”。
在这里插入图片描述
(四) 建立几何模型

  1. 创建长方形
    拾取菜单MainMenu>Preprocessor>Modeling>Create>Areas>Rectangle>By 2 Corners。弹出如下图所示对话框，Width填入1.6，Height填入0.06。点击OK。
    在这里插入图片描述
    在这里插入图片描述

(五) 划分单元

  1. 拾取菜单 Main Menu>Mesh>Mesh Tool。在 Mesh 选项中选择 Lines旁边的 set，分别划分两条长边为 15 份，两条短边为 3 份；然后点击最下方的mesh，选择矩形，点击 OK。
    在这里插入图片描述

(六) 施加载荷和边界条件

  1. MainMenu>Solution>DefineLoads>Apply>Structural>Displacement>On Lines,拾取左端线，然后点击OK。并在弹出的对话框里选择 All DOF，点击 ok。
    在这里插入图片描述
    在这里插入图片描述
  2. 选择 MainMenu>Solution>DefineLoads>Apply>Structural>
    pressure>On Lines,拾取上端线，点击 OK。在弹出的对话框中输入200000。
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

(七) 求解模型
依次点击Main Menu>Solution>Solve>Current LS进行求解，当提示Solution is Done说明求解完成。
在这里插入图片描述
(八) 结果查看
依次单击Main Menu>GeneralPostproc>PlotResults>Contour Plot>Nodal Solu，弹出对话框，在“Item to be Contoured”栏中选择“Stress”，在其右面的栏中，选择“von Mises stress”后，单击“OK”，这时在图形窗口处显示出平面的应力云图，如图所示。
在这里插入图片描述
在这里插入图片描述
依次单击 Main Menu>GeneralPostproc>Plot Results>Contour Plot>Nodal Solu，弹出对话框，在“Item to be Contoured”栏中选择“DOF Solution”，在其右面的栏中，选择“Displacement vector sum”后，单击“OK”，这时在图形窗口处显示出三维应变图，如图所示。
在这里插入图片描述
在这里插入图片描述

4.3 梁单元分析

(一) 定义工作文件名
依次单击Utility Menu>File>Change Jobname，在出现的对话框中输入：beam，在“New Log and error files”处选中“yes”，单击“OK”。
(二) 定义单元类型

  1. 依次单击 MainMenu>Preprocessor>Element type>Add/Edit/Delete，出现对话框，单击“Add”，出现一个 “Library of Element Type”对话框。在“Library of Element Type”左面的列表栏中选择“Beam”，在右面的列表栏中选择“2 node 188”，单击“OK”。
    在这里插入图片描述

(三) 设置材料属性
依次单击Main Menu>Preprocessor>Material Props>Material Models，出现“Define Material Model Behavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural>Linear>Elastic>Isotropic”,出现对话框，输入弹性模量EX＝2e11，泊松比PRXY＝0.28，单击“OK”，单击“Material>Exit”。
在这里插入图片描述

(四) 建立几何模型

  1. 创建直线
    创建关键点1（0，0，0）和关键点2（1.6，0，0），再创建关键点 10（0，1， 0），这是截面的方向点。然后连接点 1 和 2 形成直线。
    在这里插入图片描述
  2. 定义截面属性
    在这里插入图片描述

(五) 划分单元

  1. 打开 Mesh tool，在上面的 set 中选 line;选择关键点 10 为截面方向点，选好之后点击中间 line 后的 set，划分成15份，然后点击最底下的 mesh，选择直线，点击 OK。
    在这里插入图片描述
    在这里插入图片描述

(六) 施加载荷和边界条件

  1. 对左端点施加全约束载荷。
    在这里插入图片描述
    在这里插入图片描述
  2. 选择 MainMenu>Solution>DefineLoads>Apply>Structural>pressure>On Beams，框选整根梁。
    在这里插入图片描述
    由于载荷作用在y轴负方向，因此在弹出窗口的第一个空中输入1，由于需要将题目中的均布在面上的载荷转换为均布在线上的载荷，因此第二个以及第三个空中输入10000(200000*0.05)。
    在这里插入图片描述
    在这里插入图片描述

(七) 求解模型
依次点击Main Menu>Solution>Solve>Current LS进行求解，当提示Solution is Done说明求解完成。
在这里插入图片描述
(八) 结果查看
依次单击Main Menu>GeneralPostproc>PlotResults>Contour Plot>Nodal Solu，弹出对话框，在“Item to be Contoured”栏中选择“Stress”，在其右面的栏中，选择“von Mises stress”后，单击“OK”，这时在图形窗口处显示出平面的应力云图，如图所示。
在这里插入图片描述
在这里插入图片描述
依次单击 Main Menu>GeneralPostproc>Plot Results>Contour Plot>Nodal Solu，弹出对话框，在“Item to be Contoured”栏中选择“DOF Solution”，在其右面的栏中，选择“Displacement vector sum”后，单击“OK”，这时在图形窗口处显示出三维应变图，如图所示。
在这里插入图片描述
在这里插入图片描述

5. 结果和讨论

有限元计算结果与理论计算值的对比如下表所示。
在这里插入图片描述
由表可以看出，采用三维实体单元与平面应力单元进行分析，得到的结果与理论计算结果相差较大，采用梁单元分析得到的结果与理论计算结果相差较小。产生的误差可能是由于模型误差导致的。
因此可得出结论，采用梁单元进行分析，与理论计算结果的差异更小。
附：理论计算过程
在这里插入图片描述
在这里插入图片描述

吮指原味张
关注
  • 16
    点赞
  • 踩
  • 58
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
ANSYS经典界面中梁单元实例全解析
python_AI_fans的博客
05-24 2万+
转载链接：http://www.jishulink.com/content/post/361389 文中红色部分为自己添加！！！ 转载的文章中提到的梁单元截面是常用的工字形，那么要是不规则的截面形状的话需要自己进行截面形状绘制，具体可以参考这篇文章：https://blog.csdn.net/qq_45769063/article/details/106316667 目录 导读： 一、问题描述 二、GUI步骤 1．定义单元属性 (1)定义单元类型： (2)设置材料属性： (3)定义截面.
ansys经典界面分析工字钢_经典仿真案例教程 | 01 分布荷载的应用
weixin_39799646的博客
12-29 1114
点击蓝字，关注我们斯姆勒数值仿真技术研究院经典仿真案例雅典娜仿真技术共享云平台斯姆勒数值仿真技术研究院协同雅典娜仿真技术共享云平台，共同推出经典仿真案例系列文章，定期为大家提供有限元仿真案例教学内容，以供学习交流。学仿真，上雅典娜，体系化学习！01分布荷载的应用本教程的目的是解释如何应用分布式荷载和使用图元表提取数据。分布荷载为1000N/m(1N/mm)将施加在矩形截面实心钢梁上，如下...
有限元法求解悬臂梁弯曲问题的matlab实现.zip
02-13
1.版本：matlab2014/2019a/2021a，内含运行结果，不会运行可私信 2.附赠案例数据可直接运行matlab程序。 3.代码特点：参数化编程、参数可方便更改、代码编程思路清晰、注释明细。 4.适用对象：计算机，电子信息工程、数学等专业的大学生课程设计、期末大作业和毕业设计。
【ansys workbench】18.实际问题：均布载荷的简化
糊涂懿的博客
02-01 9913
本篇博客是根据阅读公众号“机械人读书笔记”而来的学习笔记~
ANSYS 静力 分析实例 悬臂梁
热门推荐
Hulunbuir
06-30 1万+
本文以力学中最简单的悬臂梁受载为例，介绍了采用 ANSYS 进行有限元分析的基本操作流程，并结合基础的有限元知识，探讨了有限元近似解的精度水平。本文还讨论了单元类型、单元阶次、积分方案、网格尺寸等对计算结果的影响。
《有限元法》学习笔记2 求表面力产生的节点载荷
lijil168的专栏
10-03 694
有限元法》求表面力产生的节点载荷，公式很简单，但要掌握其物理意义并不轻松，理解了其中含义，就可以举一反三，处理其他形式的载荷。
悬臂梁分析,悬臂梁分析误差,ansys
09-10
ANSYS是一款强大的有限元分析软件，其APDL（ANSYS Parametric Design Language）命令流功能使得用户可以通过编程方式执行复杂分析任务。 在这个案例中，我们将通过ANSYS APDL来对实体悬臂梁进行分析，以理解其工作...
ansys 实验二 悬臂梁质量块的实体建模与仿真
12-28
在ANSYS软件的应用中，悬臂梁的质量块实体建模与仿真是一项基础且重要的任务，它涉及到工程领域中的许多实际问题，如结构力学分析、振动研究等。本实验以一个100mm长的悬臂梁为例，详细介绍如何利用ANSYS进行模型...
工字悬臂梁的有限元分析
01-11
在本实验中，我们使用 ANSYS 软件对工字悬臂梁进行有限元分析，并将结果与手算结果进行比较。实验步骤包括建立几何模型、离散化处理、选择位移模式、确定单元应变与位移、应力与应变的关系、施加约束和载荷、求解...
beam_梁固有频率_悬臂梁振型_平面曲线悬臂梁模态_
10-04
悬臂梁的振型是描述其在特定频率下振动形态的物理量，这对于预测梁在动态载荷下的行为至关重要。振型通常通过数学公式表示，可以是正弦函数或复数形式，它们描述了梁在各个位置上的位移分布。 平面曲线悬臂梁模态则...
梁单元均布荷载
03-16
3．/solu 加边界条件 设置求解选项 定义载荷步 求解载荷步 4./post1（通用后处理） 5./post26（时间历程后处理） 6.PLOTCONTROL菜单命令 7.参数化设计语言
悬臂梁分析,悬臂梁分析误差,ansys源码.zip
10-10
悬臂梁分析是工程力学中的一个经典问题，主要涉及到结构力学和有限元分析等领域。在实际工程中，悬臂梁通常被用作支撑结构的一部分，例如桥梁、建筑或机械部件，它们的一端固定，另一端自由，因此承受的载荷会引发...
Python与Ansys apdl有限元系列三：单层单跨梁单元受竖向均布力，水平集中力
瓜地考拉
04-17 1219
单层单跨梁单元受竖向均布力，水平集中力
02初尝有限元分析——悬臂梁案例
陶小桃Blog
09-02 6129
一、简述有限元 基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，单元上所作用的力等效到节点上，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，就是用叉值函数来近似代替 ，借助于变分原理或加权余量法，将微分方程离散求解。 有限元分析基本过程：几何建模->网格划分->模型求解->结果判定与输出 拓展阅读： 有限元方法的核心思想是什么？ 有限元法基本思想和分类 什么是有限元？方法、基本思想是.
ID CS6怎么加框_苹果手机id密码忘了怎么办，超级实用。
weixin_39618806的博客
11-04 159
在使用苹果手机的时候，必须要有iPhone id和密码才能正常的使用苹果手机，并在APP Store里面下载软件，如果在使用的时候，ID密码忘记了，就会影响到很多操作，给我们带来很大的麻烦，如果苹果id密码忘了怎么办呢?今天就让大家利用手机快速设置id密码。1.在手机中，点击【设置】，并在设置的页面点击如图位置，进入到apple ID的页面。2.在【apple ID】页面点击下方的【itunes ...
ansys经典界面分析工字钢_ansys工字钢实例分析.doc
weixin_30632267的博客
01-13 901
ansys工字钢实例分析.doc ANSYS工字钢分析学院机械工程学院专业材料成型及控制工程班级2013级1班姓名邓祥丰学号201310112101一，定义单元类型1选择MAINMENUPREPROCESSORELEMENTTYPEADD/EDIT/DELETE,在弹出的对话框中选择ADD，点击SOLID中的BRICK8NODE45选项。二，定义材料参数1选择MAINMENUPREPROCESSO...
Ansys APDL施加任意形式载荷（线性/非线性，通过命令流/GUI）
ThreeS_tones的博客
01-26 5729
Ansys APDL施加任意形式载荷（线性/非线性，通过命令流/GUI）
ansys经典界面分析工字钢_【干货】基于ANSYS经典界面的接触分析例子
weixin_34684705的博客
12-23 802
原标题：【干货】基于ANSYS经典界面的接触分析例子1．问题描述一个钢销插在一个钢块中的光滑销孔中。已知钢销的半径是0.5 units, 长是2.5 units，而钢块的宽是4 Units, 长4 Units,高为1 Units,方块中的销孔半径为0.49 units,是一个通孔。钢块与钢销的弹性模量均为36e6，泊松比为0.3.由于钢销的直径比销孔的直径要大，所以它们之间是过盈配合。现在要对该问...
LIBSVM参数实例详解.rar
最新发布
07-26
神经网络的matlab案例，本案例介绍如下： 技术深度：案例详细介绍了如何使用MATLAB进行深度学习模型的构建、训练和测试。 实际应用：通过具体的图像识别任务，展示模型的实际应用效果，让你直观感受神经网络的强大能力。 代码解析：提供完整的MATLAB代码，并对关键部分进行详细注释，帮助你理解每一步的工作原理。 优化策略：探讨不同的训练策略和参数调整方法，优化模型性能。
matlab四节点四边形单元悬臂梁均布载荷
06-10
好的，让我来回答你的问题。 首先，我们需要知道悬臂梁的几何参数和材料参数，例如长度、宽度、高度、截面面积、杨氏模量等。假设悬臂梁的长度为L，宽度为b，高度为h，截面面积为A，杨氏模量为E。 其次，我们需要定义四节点四边形单元的节点坐标和节点编号，以及单元的材料属性和截面属性。假设单元的节点坐标为(x1,y1),(x2,y2),(x3,y3),(x4,y4)，节点编号为1,2,3,4，材料属性为E，截面属性为A。 最后，我们可以使用有限元分析软件（如MATLAB）来求解悬臂梁的位移和应力。其中，均布载荷可以表示为单元节点上的等效力，即每个节点上的力大小相等，方向相反。 具体的步骤如下： 1. 定义节点坐标和节点编号； 2. 定义单元的材料属性和截面属性； 3. 定义均布载荷的等效力； 4. 求解单元的刚度矩阵和等效力向量； 5. 组装整个系统的刚度矩阵和等效力向量； 6. 求解位移向量和应力。 下面给出MATLAB的示例代码： ```matlab % 定义悬臂梁的几何参数和材料参数 L = 1; % 长度 b = 0.1; % 宽度 h = 0.2; % 高度 A = b*h; % 截面面积 E = 2e11; % 杨氏模量 % 定义四节点四边形单元的节点坐标和节点编号 x1 = 0; y1 = 0; x2 = L; y2 = 0; x3 = L; y3 = h; x4 = 0; y4 = h; node = [1, x1, y1; 2, x2, y2; 3, x3, y3; 4, x4, y4]; elem = [1, 1, 2, 3, 4]; % 定义单元的材料属性和截面属性 D = E*A/L; Ke = [ 12, 6*L, -12, 6*L; 6*L, 4*L^2, -6*L, 2*L^2; -12, -6*L, 12, -6*L; 6*L, 2*L^2, -6*L, 4*L^2]; fe = [0; -0.5*b*h*L; 0; -0.5*b*h*L]; % 定义均布载荷的等效力 f = [0; -0.5*b*h*L; 0; -0.5*b*h*L]; % 求解位移和应力 K = zeros(8,8); f_e = zeros(8,1); for i = 1:size(elem,1) n1 = elem(i,2); n2 = elem(i,3); n3 = elem(i,4); n4 = elem(i,5); x1 = node(n1,2); y1 = node(n1,3); x2 = node(n2,2); y2 = node(n2,3); x3 = node(n3,2); y3 = node(n3,3); x4 = node(n4,2); y4 = node(n4,3); Le = sqrt((x2-x1)^2 + (y2-y1)^2); cosx = (x2-x1)/Le; sinx = (y2-y1)/Le; T = [ cosx, sinx, 0, 0; -sinx, cosx, 0, 0; 0, 0, cosx, sinx; 0, 0, -sinx, cosx]; Ke_local = D*Ke/Le; fe_local = T*fe; K([n1,n2,n3,n4],[n1,n2,n3,n4]) = K([n1,n2,n3,n4],[n1,n2,n3,n4]) + Ke_local; f_e([n1,n2,n3,n4],1) = f_e([n1,n2,n3,n4],1) + fe_local; end f_eq = f_e - f; U = K\f_eq; stress = D*Ke*[U(2)-U(1); (U(2)-U(1))/L; U(4)-U(3); (U(4)-U(3))/L]; ``` 注意，在这个示例代码中，我们假设均布载荷是沿y方向的，大小为0.5bhL，即每个节点上的等效力为0.25bhL。 希望这个回答能够解决你的问题。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

吮指原味张

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

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

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

打赏作者

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

抵扣说明：

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

余额充值