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

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
    评论
好的,让我来回答你的问题。 首先,我们需要知道悬臂梁的几何参数和材料参数,例如长度、宽度、高度、截面面积、杨氏模量等。假设悬臂梁的长度为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、付费专栏及课程。

余额充值