有限元悬臂梁仿真

有限元悬臂梁仿真
有限元悬臂梁仿真

有限元方法大作业

课程设计题目:

若干个质量不等的仪器要安装在均匀悬臂梁(或板)不同位置上,仪器间要有预留安全距离,试确定一种安装方法,使梁(或板)的变形最小或第一阶固有频率最高。

题目分析:

1 题目中没有给定梁的材料和形状、仪器的数量和质量,以及仪器的安全距离。在这里不妨假定,梁的材料为结构钢,其密度为、杨氏模量为Pa、泊松比为0.3,梁的形状为。仪器的数量为3个,均匀的安装在梁上,其质量及其组合如表1所示。

表1 仪器的质量、及其组合

2 本次采用solidworks建立梁的实体模型,并导入ansys workbench软件中进行计算。梁模型左端固定,仪器安装顺序依次从左到右。

3 在题目中,需要找到一种安装组合使得梁的变形最小或第一阶固有频率最高,这分别是静力学分析问题和模态分析问题。在静力学分析中,如图2-1所示,在梁上安装仪器的位置上,加上一个加力面(半径为20mm的圆)。在加力面上可以施加均布载荷,这里将仪器的质量换算成相应的均布载荷,施加到相应的加力面上,如图2-2所示。

图2-1 ansys workbench实体梁的静力学分析

图2-2 加力面和加力面上的均布载荷

采用solid187单元对模型进行网格划分,solid187单元是一个高阶3维10节点固体结构单元,如图2-3所示,单元通过10个节点来定义,每个节点有3个沿着xyz方向平移的自由度。并对加力面附近进行加密,如图2-4所示。进而进行静力学分析,得到梁的总体变形量(total-Deformation)。

图2-3 solid187单元

图2-4 梁模型网格划分和加力面加密

4 梁的固有频率可由无阻尼自由振动方程求解:

令:

得到:

当:

从而求的梁的自振频率。在ansys workbench中,将仪器的质量用质量点代替,并安置在相应的位置上,如图2-5所示。采用solid186单元对模型进行网格划分,其结果如图2-6所示,solid186是一个高阶3维20节点固体结构单元,如图2-7所示,单元通过20个节点来定义,每个节点有3个沿着xyz方向平移的自由度。然后,求解梁模型的前6阶的固有频率。

图2-5 梁模型上的质量点

图2-6 梁模型模态分析网格划分

图2-7 solid186单元

分析结果:

1 静力学分析

1)组合1:m1-m2-m3

图3-1 组合1:m1-m2-m3变形图

2)组合2:m1-m3-m2

图3-2 组合2:m1-m3-m2变形图

3)组合3:m2-m1-m3

图3-3 组合3:m2-m1-m3变形图

4)组合4:m2-m3-m1

图3-4 组合3:m2-m3-m1变形图

5)组合5:m3-m1-m2

图3-5 组合5:m3-m1-m2变形图

6)组合6:m3-m2-m1

图3-6 组合6:m3-m2-m1变形图2 模态分析

组合1 组合2 组合3

组合4 组合5 组合6

图3-7 各组合前6阶固有频率

结论

梁模型的静力学分析和模态分析结果如下表2所示。

梁模型左端固定,仪器安装顺序依次从左到右。有上表可知,组合6从左到右,仪器的质量依次减小,其变形量最小,一阶固有频率最大。并且,当质量越大的仪器越靠近固定端的时候,整个梁的变形量最小,其一阶固有频率越大。

悬臂梁分析报告

悬臂梁受力分析报告 高一博 2016.11.13 西安理工大学 机械与精密仪器工程学院

摘要 利用ANSYS对悬臂梁进行有限元静力学分析,得到悬臂梁的最大应力和挠度位移。从而校验结构强度和尺寸定义,从而对结构进行最优化设计修正。 关键词:悬臂梁,变形分析,应力分析

目录 一.问题描述: (4) 二.分析的目的和内容: (4) 三.分析方案和有限元建模方法: (4) 四.几何模型 (4) 五.有限元模型 (4) 六.计算结果: (5) 七.结果合理性的讨论、分析 (8) 八.结论 (8) 参考文献 (8)

一.问题描述: 现有一悬臂梁,长500MM,一端固定,另外一端施加一个竖直向下的集中力200N。 其截面20MMX20MM的矩形,现在要分析该梁的在集中力作用下产生的位移,应力和局部应力。 二.分析的目的和内容: 1.观察悬臂梁的变形情况; 2.观察分析悬臂梁的应力变化; 3.找出其最大变形和最大应力点,分析形成原因; 三.分析方案和有限元建模方法: 1.使用ANSYS-modeling-create-volumes-block建模, 2.对梁进行材料定义,网格划分。 3.一端固定,另外一端施加一个向下的200N的力。 4.后处理中查看梁的应力和变形情况。 四.几何模型 500X20X20的梁在在ANSYS中进行绘制.由于结构简单规则,无需简化。 五.有限元模型 单元类型:solid brick8node45 材料参数:弹性模量2e+11pa,泊松比0.3 边界条件:一端固定,一端施加载荷 载荷:F=200N 划分网格后的悬臂梁模型

ANSYS悬臂梁的自由端受力的有限元计算[1]

悬臂梁自由端受力的有限元计算 任柳杰10110290005 一、计算目的 1、掌握ANSYS软件的基本几何形体构造、网格划分、边界条件施加等方法。 2、熟悉有限元建模、求解及结果分析步骤和方法。 3、利用ANSYS软件对梁结构进行有限元计算。 4、梁的变形、挠曲线等情况的分析。 5、一维梁单元,二维壳单元,三维实体单元对计算结果的影响。 6、载荷施加在不同的节点上对结果的影响。 二、计算设备 PC,ANSYS软件(版本为11.0) 三、计算内容 悬臂梁受力模型 如上图所示,一段长100[mm]的梁,一端固定,另一段受到平行于梁截面的集中力F的作用,F=100[N]。梁的截面为正方形,边长为10[mm]。梁所用的材料:弹性模量E=2.0 105[MPa],泊松比0.3。 四、计算步骤(以梁单元为例) 1、分析问题。 分析该物理模型可知,截面边长/梁长度=0.1是一个较小的值,我们可以用梁单元来分析这样的模型。当然,建立合适的壳单元模型和实体单元模型也是可以的。故拟采用这三种不同的 方式建立模型。以下主要阐述采用梁单元的模型的计算步骤。 2、建立有限元模型。 a)创建工作文件夹并添加标题; 在个人的工作目录下创建一个文件夹,命名为beam,用于保存分析过程中生成的各种文件。 启动ANSYS后,使用菜单“File”——“Change Directory…”将工作目录指向beam 文件夹;使用/FILNAME,BEAM命令将文件名改为BEAM,这样分析过程中生成的文件均 以BEAM为前缀。 偏好设定为结构分析,操作如下: GUI: Main Menu > Preferences > Structural b)选择单元; 进入单元类型库,操作如下: GUI: Main Menu > Preprocessor > Element Type > Add/Edit/Delete > Add… 对话框左侧选择Beam选项,在右侧列表中选择2D elastic 3选项,然后单击OK按钮。

有限元悬臂梁仿真

有限元方法大作业 课程设计题目: 若干个质量不等的仪器要安装在均匀悬臂梁(或板)不同位置上,仪器间要有预留安全距离,试确定一种安装方法,使梁(或板)的变形最小或第一阶固有频率最高。 题目分析: 1 题目中没有给定梁的材料和形状、仪器的数量和质量,以及仪器的安全距离。在这里不妨假定,梁的材料为结构钢,其密度为、杨氏模量为Pa、泊松比为0.3,梁的形状为。仪器的数量为3个,均匀的安装在梁上,其质量及其组合如表1所示。 表1 仪器的质量、及其组合 2 本次采用solidworks建立梁的实体模型,并导入ansys workbench软件中进行计算。梁模型左端固定,仪器安装顺序依次从左到右。 3 在题目中,需要找到一种安装组合使得梁的变形最小或第一阶固有频率最高,这分别是静力学分析问题和模态分析问题。在静力学分析中,如图2-1所示,在梁上安装仪器的位置上,加上一个加力面(半径为20mm的圆)。在加力面上可以施加均布载荷,这里将仪器的质量换算成相应的均布载荷,施加到相应的加力面上,如图2-2所示。 图2-1 ansys workbench实体梁的静力学分析

图2-2 加力面和加力面上的均布载荷 采用solid187单元对模型进行网格划分,solid187单元是一个高阶3维10节点固体结构单元,如图2-3所示,单元通过10个节点来定义,每个节点有3个沿着xyz方向平移的自由度。并对加力面附近进行加密,如图2-4所示。进而进行静力学分析,得到梁的总体变形量(total-Deformation)。 图2-3 solid187单元 图2-4 梁模型网格划分和加力面加密 4 梁的固有频率可由无阻尼自由振动方程求解: 令: 得到: 当: 从而求的梁的自振频率。在ansys workbench中,将仪器的质量用质量点代替,并安置在相应的位置上,如图2-5所示。采用solid186单元对模型进行网格划分,其结果如图2-6所示,solid186是一个高阶3维20节点固体结构单元,如图2-7所示,单元通过20个节点来定义,每个节点有3个沿着xyz方向平移的自由度。然后,求解梁模型的前6阶的固有频率。

有限元分析及应用报告-利用ANSYS软件分析带孔悬臂梁

有限元分析及应用报告 题目:利用ANSY软件分析带孔悬臂梁 姓名:xxx 学号:xxx 班级:机械xxx 学院: 机械学院 指导老师:xxx 二零一五年一月

问题概述 图示为一隧道断面,其内受均布水压力q,外受土壤均 布压力p;试采用不同单元计算断面内的位移及应力,并分别分析q=0或p=0时的位移和应力分布情况。(材料为钢,隧道几何尺寸和压力大小自行确定) 本例假定内圆半径为1m,外圆半径为2m,外受均布 压力p=10000pa ,内受均布压力为q=20000pa 。 问题分析 由题目可知,隧道的的长度尺寸远远大于截面尺寸,并且压力在长度方向上均匀分布,因此本问题可以看作为平面应变问题。由于在一个截面内,压力沿截面四周均匀分布,且截面是对称的圆环,所以可以只取截面1/4进行有限元建模分析,这样不仅简化了建模分析过程,也能保证得到精确的结果。由以上分析,可以选取单元类型plane42进行有限元分析,在option中选择K3 为plane strain。

三.有限元建模 1.设置计算类型 由问题分析可知本问题属于平面静应力问题,所以选择preferences 为structure 。 2.单元类型选定 选取平面四节点常应变单元plane42,来计算分析隧道截面的位移和应力。由于此问题为平面应变问题,在设置element type的K3时将其设置为plane strain。 3.材料参数 隧道的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比(T =0.3 4.几何建模 按照题目所给尺寸利用ansys的modeling依次建立keypoint : 1(0,0),2(1,0),3(2,0),4(0,2),5(0,1) , create LINES 依次连接keypoint 2、3和4、5即可创建两条直线,使用create article 的By cent & radius 创建两条圆弧。create AREAS依次选择四条线即建立了所需的1/4截面。 5.网格划分

悬臂梁ansys有限元分析求最大挠度

(一) 悬臂梁ansys 有限元分析求最大挠度 问题:悬臂梁长1000mm ,宽50mm ,高10mm ,左端固定,求其在自重作用下的最大挠度? 解:弯矩方程: 221) ()(x l q x M --= 微分方程: 22 1'')(x l q y EI z -= 积分求解:D Cx qx qlx x ql y EI C qx qlx x ql y EI z z +++-=++-=4322322'24 1 6125.06 1 5.05.0 由边界条件:0; 0, 0' ' ====A A A y y x θ 得:C=0, D=0 I=1/12*h^3*b,h 为梁截面的高,b 为梁截面的宽。 q=ρ*g*a*h*l 材料力学公式求:Y=EI 85 gahl^ρ=5.733mm L

ANSYS 模拟求:Y=5.5392mm,详细见下步骤 ANSYS 软件设置及其具体过程如下: 步骤1:建立一个模型,在model下creat一个长1,宽0.05,高0.01的长方体实体。(单位默认为m) 步骤2:材料属性设置。密度:7800,杨氏模量:2E11,泊松比0.3。

步骤3:划分网格。设置网格单元为structure solid brick 8node 185,mesh tool中设置网格大小为0.002,HEX下点击mesh。

步骤4:施加载荷;在preprocessor中inertia中设置重力加速度Y方向为9.8。在左面施加固定约束(三个方向固定)

步骤5::求解。在solve下solve current LS。 步骤6:后处理查看。在result中plot result,查看nodes displacement。List查看文本,观察nodes的最大位移点。

ansys-二维悬臂梁有限元分析

1 研究目的与问题阐述 1.1 基本研究目的 (1) 掌握ANSYS软件的基本几何形体构造、网格划分、边界条件施加等方法。 (2) 熟悉有限元建模、求解及结果分析步骤和方法。 (3) 利用ANSYS软件对梁结构进行有限元计算。 (4) 研究不同泊松比对同一位置应力的影响。 1.2 基本问题提出 图1.1 模型示意图 如图1.1所示,当EX=3.01e6,F=5000N,悬臂梁杆一端固定,另一端为自由端。当悬臂梁的泊松比u为:0.2、0.25、0.3、0.35、0.4时,确定同一位置的应力分布,得出分布云图。 采用二维模型,3*0.09m。

2 软件知识学习 2.1 软件的使用与介绍 软件介绍: ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件。由世界上最大的有限元分析软件公司之一的美国ANSYS开发,它能与多数CAD软件接口,实现数据的共享和交换,如Pro/Engineer, NASTRAN, Alogor, I-DEAS, AutoCAD等,是现代产品设计中的高级CAE工具之一。 ANSYS有限元软件包是一个多用途的有限元法计算机设计程序,可以用来求解结构、流体、电力、电磁场及碰撞等问题。因此它可应用于以下工业领域:航空航天、汽车工业、生物医学、桥梁、建筑、电子产品、重型机械、微机电系统、运动器械等。 软件主要包括三个部分:前处理模块,分析计算模块和后处理模块。 前处理模块提供了一个强大的实体建模及网格划分工具,用户可以方便地构造有限元模型; 分析计算模块包括结构分析(可进行线性分析、非线性分析和高度非线性分析)、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析,可模拟多种物理介质的相互作用,具有灵敏度分析及优化分析能力; 后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示(可看到结构内部)等图形方式显示出来,也可将计算结果以图表、曲线形式显示或输出。 软件提供了100种以上的单元类型,用来模拟工程中的各种结构和材料。该软件有多种不同版本,可以运行在从个人机到大型机的多种计算机设备上,如PC,SGI,HP,SUN,DEC,IBM,CRAY等。

悬臂梁有限元模拟分析步骤

Introduction to Simulation I-DEAS Tutorials: Simulation Projects Simulation involves three major steps: Pre-processing (modeling, applying boundary conditions, meshing); solving the model; and post-processing (displaying the results). Learn how to: ?create a finite element model ?apply boundary conditions ?mesh the FE model ?solve the FE model ?display the results

Before you begin... Prerequisite tutorials:?Introducing the I-DEAS Interface Quick Tips to Using I-DEAS –and– Creating Parts ?Extruding and Revolving Features

If you didn’t start I-DEAS with a new (empty) model file, open a new one now and give it a unique name. File Open Open Model File form Model File name: any unique name OK Simulation Master Modeler Set your units to mm. Options Units mm (milli newton)

专业课设,悬臂梁有限元分析

1 研究目的与问题阐述 1.1 基本研究目的 (1) 掌握ANSYS软件的基本几何形体构造、网格划分、边界条件施加等方法。 (2) 熟悉有限元建模、求解及结果分析步骤和方法。 (3) 利用ANSYS软件对梁结构进行有限元计算。 (4) 研究不同泊松比对同一位置应力的影响。 1.2 基本问题提出 图1.1 模型示意图 如图1.1所示,当EX=3.01e6,F=5000N,悬臂梁杆一端固定,另一端为自由端。当悬臂梁的泊松比u为:0.2、0.25、0.3、0.35、0.4时,确定同一位置的应力分布,得出分布云图。 二维模型,3*0.09m。 2 软件的介绍与使用 2.1 ANSYS 简介 ANSYS程序是一个功能强大的灵活的设计分析及优化、融结构、流体、电场、磁场、声场分析于一体的大型通用有限元商用分析软件,可广泛应用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、造船、生物医学、轻工、地矿、水利、日用家电等一般工业及科学研究。该软件提供了一个不断改进的功能清单,集体包括:结构高度非线性分析、电磁分析、计算流体动力分析、设计优化、接触分析、自适应网格划分、大应变/有限转动工功能一接利用ANSYS参数设计的扩展宏命令功能。 ANSYS由世界上最大的有限元分析软件公司之一的美国ANSYS开发,它能与多数

系统下生成的集合数据传入ANSYS,如Pro/Engineer, NASTRAN, Alogor, I-DEAS, AutoCAD等,并通过必要的修补可准确地在该模型上划分网格并求解。 2.2 ANSYS软件的功能介绍 ANSYS软件含有多种有限元分析的能力,包括从简单线性静态分析到复杂非线性动态分析。一个典型的ANSYS分析过程可分为以下三个步骤: 创建有限元模型; 施加载荷进行求解; 查看分析结果; 在有限元的分析过程中,程序通常使用以下三个部分:前处理模块,分析求解模块和后处理模块。 前处理模块提供了一个强大的实体建模及网格划分工具,通过这个模块用户可以建立自己想要的工程有限模型。 分析求解模块即是对已建立好的模型在一定的载荷和边界条件下进行有限元计算,求解平衡微分方程。包括结构分析(可进行线性分析、非线性分析和高度非线性分析)、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析(热-应力耦合、流-固耦合以及电-磁-热-应力耦合)等,可模拟多种物理介质的相互作用,具有灵敏度分析及优化分析能力; 后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示(可看到结构内部)等图形方式显示出来,也可将计算结果以图表、曲线形式显示或输出。 下面对ANSYS软件的三种模块的功能进行简要介绍: 1.前处理模块 ANSYS软件的前处理模块主要实现三种功能:参数定义、实体建模和网格划分。 (1)参数定义 ANSYS程序在进行结构建模的过程中,首先要对所有被建模型的材料进行参数定义。包括定义使用单位制,定义所有使用单元的类型,定义单元的实常数,定义材料的特性以及使用材料库文件等。 (2)实体建模

悬臂梁—有限元ABAQUS线性静力学分析实例

线性静力学分析实例——以悬臂梁为例 线性静力学问题是简单且常见的有限元分析类型,不涉及任何非线性(材料非线性、几何非线性、接触等),也不考虑惯性及时间相关的材料属性。在ABAQUS中,该类问题通常采用静态通用(Static,General)分析步或静态线性摄动(Static,Linear perturbation)分析步进行分析。 线性静力学问题很容易求解,往往用户更关系的是计算效率和求解效率,希望在获得较高精度的前提下尽量缩短计算时间,特别是大型模型。这主要取决于网格的划分,包括种子的设置、网格控制和单元类型的选取。在一般的分析中,应尽量选用精度和效率都较高的二次四边形/六面体单元,在主要的分析部位设置较密的种子;若主要分析部位的网格没有大的扭曲,使用非协调单元(如CPS4I、C3D8I)的性价比很高。对于复杂模型,可以采用分割模型的方法划分二次四边形/六面体单元;有时分割过程过于繁琐,用户可以采用精度较高的二次三角形/四面体单元进行网格划分。 悬臂梁的线性静力学分析 1.1 问题的描述 一悬臂梁左端受固定约束,右端自由,结构尺寸如图1-1所示,求梁受载后的Mises应力、位移分布。 ν 材料性质:弹性模量3 = E=,泊松比3.0 2e 均布载荷:F=103N 图1-1 悬臂梁受均布载荷图 1.2 启动ABAQUS 启动ABAQUS有两种方法,用户可以任选一种。 (1)在Windows操作系统中单击“开始”--“程序”--ABAQUS 6.10 --

ABAQUS/CAE。 (2)在操作系统的DOS窗口中输入命令:abaqus cae。 启动ABAQUS/CAE后,在出现的Start Section(开始任务)对话框中选择Create Model Database。 1.3 创建部件 在ABAQUS/CAE顶部的环境栏中,可以看到模块列表:Module:Part,这表示当前处在Part(部件)模块,在这个模块中可以定义模型各部分的几何形体。可以参照下面步骤创建悬臂梁的几何模型。 (1)创建部件。对于如图1-1所示的悬臂梁模型,可以先画出梁结构的二维截面(矩形),再通过拉伸得到。 单击左侧工具区中的(Create Part)按钮,或者在主菜单里面选择Part--Create,弹出如图1-2所示的Create Part对话框。 图1-2 Create Part对话框 在Name(部件名称)后面输入Beam,Modeling Space(模型所在空间)设

悬臂梁—有限元ABAQUS线性静力学分析实例

悬臂梁—有限元ABAQUS线性静力学分析实例

线性静力学分析实例——以悬臂梁为例 线性静力学问题是简单且常见的有限元分析类型,不涉及任何非线性(材料非线性、几何非线性、接触等),也不考虑惯性及时间相关的材料属性。在ABAQUS中,该类问题通常采用静态通用(Static,General)分析步或静态线性摄动(Static,Linear perturbation)分析步进行分析。 线性静力学问题很容易求解,往往用户更关系的是计算效率和求解效率,希望在获得较高精度的前提下尽量缩短计算时间,特别是大型模型。这主要取决于网格的划分,包括种子的设置、网格控制和单元类型的选取。在一般的分析中,应尽量选用精度和效率都较高的二次四边形/六面体单元,在主要的分析部位设置较密的种子;若主要分析部位的网格没有大的扭曲,使用非协调单元(如CPS4I、C3D8I)的性价比很高。对于复杂模型,可以采用分割模型的方法划分二次四边形/六面体单元;有时分割过程过于繁琐,用户可以采用精度较高的二次三角形/四面体单元进行网格划分。 悬臂梁的线性静力学分析 1.1 问题的描述 一悬臂梁左端受固定约束,右端自由,结构尺寸如图1-1所示,求梁受载后的Mises应力、位移分布。 ν 材料性质:弹性模量3 = E=,泊松比3.0 2e 均布载荷:F=103N 图1-1 悬臂梁受均布载荷图 1.2 启动ABAQUS 启动ABAQUS有两种方法,用户可以任选一种。 (1)在Windows操作系统中单击“开始”--“程序”--ABAQUS 6.10 --

ABAQUS/CAE。 (2)在操作系统的DOS窗口中输入命令:abaqus cae。 启动ABAQUS/CAE后,在出现的Start Section(开始任务)对话框中选择Create Model Database。 1.3 创建部件 在ABAQUS/CAE顶部的环境栏中,可以看到模块列表:Module:Part,这表示当前处在Part(部件)模块,在这个模块中可以定义模型各部分的几何形体。可以参照下面步骤创建悬臂梁的几何模型。 (1)创建部件。对于如图1-1所示的悬臂梁模型,可以先画出梁结构的二维截面(矩形),再通过拉伸得到。 单击左侧工具区中的(Create Part)按钮,或者在主菜单里面选择Part--Create,弹出如图1-2所示的Create Part对话框。 图1-2 Create Part对话框 在Name(部件名称)后面输入Beam,Modeling Space(模型所在空间)

悬臂梁应力分析有限元程序设计

题目:悬臂梁应力分析有限元程序设计

毕业设计(论文)外文摘要

本科毕业设计(论文)第Ⅰ页共Ⅰ页 目录 1 引言 (1) 2 有限元理论 (2) 2.1 有限元法产生的动因分析 (2) 2.2 有限元的发展历程 (3) 2.3 有限元分析的研究特点 (4) 2.4 有限元法的分析过程 (4) 2.5 有限元的发展趋势 (6) 3 悬臂梁应力分析有限元程序开发 (9) 3.1 Matlab语言指南 (9) 3.1.1 Matlab语言简介 (9) 3.1.2 Matlab的优点 (10) 3.2 悬臂梁应力分析程序设计 (11) 3.2.1平面问题的4节点矩形单元描述 (11) 3.2.2 平面问题4节点矩形单元的MATLAB程序 (16) 3.2.3 悬臂梁应用举例 (20) 结束语 (34) 致谢 (35) 参考文献 (36)

1 引言 悬臂梁在工程力学受力分析中,是一种比较典型的简化模型。在实际工程分析中,大部分实际工程受力部件都可以简化为悬臂梁。根据有限元法的基本原理和解决问题的基本思路,对悬臂梁所受的应力进行有限元分析有着重要的作用。尽管目前已有不少从国外引进的大型通用程序,但由于这些程序通用性很强,语句多,要求计算内存大,在计算具体问题时往往占用机时多,计算成本高,在PC微机广泛普及的今天,编制一些便于推广应用的专用程序无论对于工业设计还是教学实践都是具有一定意义的。 目前,悬臂梁结构在实际工程中被得到广泛的应用,是一种较为常用的结构,尤其在机械设计、建筑设计中更是常见。悬臂梁结构在实际的使用过程中,经常要承受各种集中载荷、分布载荷、弯矩和扭矩的作用,在梁的任意一处都有可能产生较大的应力和变形,从而使得悬臂梁结构破坏或失效。悬臂梁的强度及刚度是否满足要求将关系到整个设备的安全使用[1]。因此,在对悬臂梁结构设计的过程中,如何对悬臂梁的应力进行分析,具有工程实用价值和现实意义。有限元分析是用来决定复杂机械结构中的应力和变形的一种非常有效的方法,当前用计算机进行的应力分析几乎全部都是以有限元理论为基础的。反过来说,有限元方法的广泛应用也是以计算机技术的发展为其前提。机械结构的有限元模型可以看做是一个弹性系统,当该系统有载荷作用时,系统中所有单元发生变形直到全部力达到平衡为止。对于每一单元,可写出其节点位移和作用力的关系式,从而达到对悬臂梁进行应力分析的目的。悬臂梁应力分析有限元程序设计方法是计算机程序开发方法的一种变革,是利用计算机解决问题的一种新的思维方式,它使程序设计更加贴近现实。悬臂梁在计算机程序中表示为对象,其目的在于使实际问题中的悬臂梁与程序中的对象具有一一对应的关系,实现利用计算机解决实际问题的目的[2,3]。

悬臂梁的有限元分析

工程地质数值模拟成绩考核 ——昆明理工大学本科生课程 姓名:漆江 学院:国土资源工程学院 科系:地科系 专业:勘查111 学号:201110110117 2014年11 月8 日

悬臂梁的有限元分析 1.问题概述。 悬臂梁为矩形截面的钢梁,长10m宽1m、高2m,不计梁的自重,弹性模量为220GPa,泊松比为0.2,在悬臂端作用一集中荷载P=1200kN。试分析该悬臂梁的内力和变形情况。 2.启动ANSYS程序。 (1)在【开始】菜单中依次选取【所有程序】/【ANSYS8.0】/【ConfigureANSYSProducts】选项,打开【ANSYS8.0Launcher】对话框。 (2)选中【FileManagement】选项卡,输入目录名: “D:\ANSYSFX\zhang1\Exam01\ANSYSjs”,输入项目名:“Z101Beam”。 (3)单击按钮运行程序,进入ANSYS使用界面。 3.定义材料、实常数和单元类型。 (1)在【ANSYSMainMenu】菜单中依次选取【Preprocessor】(前处理)/【ElementType】/【Add/Edit/Delete】选项,打开单元类型对话框。单击按钮,打开单元类型库对话框,在右侧两个列表框中分别选取【Beam】选项和【2Delastic3】选项(简称为Beam3单元,以后叙述中记为【Beam】-【2Delastic3】单元,类似的情况记法相同),如图1-16所示。单击按钮,再单击【ElementType】对话框中的按钮。

图1-16【LibraryofElementTypes】对话框 (2)在【ANSYSMainMenu】菜单中依次选取【Preprocessor】/【RealConstants】/【Add/Edit/Delete】选项,打开实常数对话框, 如图1-17所示。单击按钮,打开Beam3实常数对话框,按照提示输入相应的面积、惯性矩和梁高参数,如图1-18所示。单击 按钮,单击【RealConstants】对话框中的按钮。 图1-17【RealConstants】对话框

悬臂梁在均布荷载作用下有限元分析

悬臂梁承受集中荷载作用问题的弹塑性分析 何方平邹里 (湘潭大学土木工程与力学学院,湖南湘潭411105) [摘要]本文针对曲杆在水平力作用下的受力性能,结合弹性力学基本方程和塑性力学中Mises屈服条件,得到了弹性阶段应力、位移之间的关系,以及材料发生塑性变形时,处于临界状态点的应力、应变值。同时,利用有限元分析软件ABAQUS,进行了数值模拟,分析结果与理论值吻合较好,证明所建立的有限元模型是合理的。 关键词:悬臂梁;集中荷载 THE ELASTIC-PLASTIC ANALYSIS OF THE CANTILEVER BEAM UNDER concentrated load He Fang-Ping Zhou Li (College of Civil Engineering & Mechanics, XiangTan University, Xiangtan 411105, China) 【Abstract】This article in view of the force performance of CANTILEVER BEAM UNDER concentrated load, combined with elastic mechanics basic equations and the plastic mechanics Mises yield conditions, obtained the elastic stage between stress and displacement, and the relationship between material happen plastic deformation, a critical state points of stress and strain value. At the same time, the finite element analysis software ABAQUS, the numerical simulation and analysis results and a good agreement with the theoretical value, show that the established finite element model is reasonable. Keywords: CANTILEVER BEAM concentrated load

悬臂梁分析

有限元法悬臂梁分析 算例: 如下图所示的悬臂梁,受均布载荷q =1N /mm 2作用。E =2.1×105N /mm 2, μ=0.3厚度h =10mm 。现用有限元法分析其位移及应力。 梁可视为平面应力状态,先按图示尺寸划分为均匀的三角形网格,共有8×10=80个单元,5×ll =55个节点,坐标轴以及单元与节点的编号如图。将均布载荷分配到各相应节点上,把有约束的节点5l 、52、53、54、55视作固定铰链,建立如图所示的离散化计算模型。 程序计算框图: (续左) (接右) 函数功能介绍 1. LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym)――该函数用于计算平面应力情况下弹性模量为E 、泊松比为NU 、厚度为t 、第一个节点坐标为(xi,yi)、第二个节点坐标为(xj,yj)、第三个节 开 始 输入材料参数 计算具有代表性的单元刚阵 K<=0 将各单元刚阵按整体编号集成到整体刚阵 处理根部约束,修改【K 】【Q 】 求解[K][δ]=[Q] 整理[δ] 并画图 计算单元应力,并输出 结束

点坐标为(xm,ym)时的线性三角形元的单元刚度矩阵.该函数返回6×6的单位刚度矩阵k. 2. LinearTriangleAssemble(K,k,i,j,m)――该函数将连接节点i,j,m的线性三角形元的单元刚度矩阵k集成到整体刚度矩阵K。每集成一个单元,该函数都将返回2N×2N的整体刚度矩阵K. 3. LinearTriangleElementStresses(E,NU,t,xi,yi,xj,yj,xm,ym,u)-- 该函数计算在平面应力情况下弹性模量为E、泊松比为NU、厚度为t、第一个节点坐标为(xi,yi)第二个节点坐标为(xj,yj)、第三个节点坐标为(xm,ym)以及单元位移矢量为u时的单元应力。该函数返回单元应力矢量。 函数源代码: function y = LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym) A = (xi*(yj-ym) + xj*(ym-yi) + xm*(yi-yj))/2;%三角形单元面积,单元节点应该按逆时针排序,保证每个三角形单元的面积都为正值(也可作为一个小函数:LinearTriangleElementArea) betai = yj-ym;betaj = ym-yi;betam = yi-yj;gammai = xm-xj; gammaj = xi-xm;gammam = xj-xi; B = [betai 0 betaj 0 betam 0 ; 0 gammai 0 gammaj 0 gammam ; gammai betai gammaj betaj gammam betam]/(2*A); %B为应变矩阵,其中betai=yi-ym,betaj=ym-yi,betam=yi-yj. gammai=xm-xj, gammaj=xi-xm, gammam=xj-xi. D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2]; %D为弹性矩阵,分为平面应力问题和平面应变问题 对于平面应力问题D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2]; 对于平面应变问题E1=E/(1-NU*NU),NU1=NU/(1-NU) y = t*A*B'*D*B;%单元刚度矩阵 function y = LinearTriangleAssemble(K,k,i,j,m) K(2*i-1,2*i-1) = K(2*i-1,2*i-1) + k(1,1); K(2*i-1,2*i) = K(2*i-1,2*i) + k(1,2); K(2*i-1,2*j-1) = K(2*i-1,2*j-1) + k(1,3); K(2*i-1,2*j) = K(2*i-1,2*j) + k(1,4); K(2*i-1,2*m-1) = K(2*i-1,2*m-1) + k(1,5); K(2*i-1,2*m) = K(2*i-1,2*m) + k(1,6); K(2*i,2*i-1) = K(2*i,2*i-1) + k(2,1); K(2*i,2*i) = K(2*i,2*i) + k(2,2); K(2*i,2*j-1) = K(2*i,2*j-1) + k(2,3); K(2*i,2*j) = K(2*i,2*j) + k(2,4); K(2*i,2*m-1) = K(2*i,2*m-1) + k(2,5); K(2*i,2*m) = K(2*i,2*m) + k(2,6); K(2*j-1,2*i-1) = K(2*j-1,2*i-1) + k(3,1); K(2*j-1,2*i) = K(2*j-1,2*i) + k(3,2); K(2*j-1,2*j-1) = K(2*j-1,2*j-1) + k(3,3); K(2*j-1,2*j) = K(2*j-1,2*j) + k(3,4); K(2*j-1,2*m-1) = K(2*j-1,2*m-1) + k(3,5); K(2*j-1,2*m) = K(2*j-1,2*m) + k(3,6); K(2*j,2*i-1) = K(2*j,2*i-1) + k(4,1); K(2*j,2*i) = K(2*j,2*i) + k(4,2); K(2*j,2*j-1) = K(2*j,2*j-1) + k(4,3); K(2*j,2*j) = K(2*j,2*j) + k(4,4); K(2*j,2*m-1) = K(2*j,2*m-1) + k(4,5); K(2*j,2*m) = K(2*j,2*m) + k(4,6); K(2*m-1,2*i-1) = K(2*m-1,2*i-1) + k(5,1); K(2*m-1,2*i) = K(2*m-1,2*i) + k(5,2); K(2*m-1,2*j-1) = K(2*m-1,2*j-1) + k(5,3); K(2*m-1,2*j) = K(2*m-1,2*j) + k(5,4); K(2*m-1,2*m-1) = K(2*m-1,2*m-1) + k(5,5); K(2*m-1,2*m) = K(2*m-1,2*m) + k(5,6); K(2*m,2*i-1) = K(2*m,2*i-1) + k(6,1); K(2*m,2*i) = K(2*m,2*i) + k(6,2); K(2*m,2*j-1) = K(2*m,2*j-1) + k(6,3); K(2*m,2*j) = K(2*m,2*j) + k(6,4); K(2*m,2*m-1) = K(2*m,2*m-1) + k(6,5); K(2*m,2*m) = K(2*m,2*m) + k(6,6); y = K;%对号入座,如前所述,每集成一次都将返回2N×2N的整体刚度矩阵K.此题为110×110 function y = LinearTriangleElementStresses(E,NU,t,xi,yi,xj,yj,xm,ym,u) A = (xi*(yj-ym) + xj*(ym-yi) + xm*(yi-yj))/2; betai = yj-ym;betaj = ym-yi;betam = yi-yj;gammai = xm-xj;gammaj = xi-xm;gammam = xj-xi; B = [betai 0 betaj 0 betam 0 ; 0 gammai 0 gammaj 0 gammam ; gammai betai gammaj betaj gammam betam]/(2*A); D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2];%平面应力和平面应变问题两种情况 y = D*B*u;%单元应力计算

悬臂梁MATLAB有限元算例注释

用有限元法对悬臂梁分析的算例 算例: 如下图所示的悬臂梁,受均布载荷q=1N/mm2作用。E=2.1×105N/mm2, μ=0.3厚度h=10mm。现用有限元法分析其位移及应力。 梁可视为平面应力状态,先按图示尺寸划分为均匀的三角形网格,共有8×10=80个单元,5×ll=55个节点,坐标轴以及单元与节点的编号如图。将均布载荷分配到各相应节点上,把有约束的节点5l、52、53、54、55视作固定铰链,建立如图所示的离散化计算模型。 程序计算框图: (续左)

程序中的函数功能介绍及源代码 1. LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym)――该函数用于计算平面应力情况下弹性模量为E、泊松比为NU、厚度为t、第一个节点坐标为(xi,yi)、第二个节点坐标为(xj,yj)、第三个节点坐标为(xm,ym)时的线性三角形元的单元刚度矩阵.该函数返回6×6的单位刚度矩阵k. 2. LinearTriangleAssemble(K,k,i,j,m)――该函数将连接节点i,j,m的线性三角形元的单元刚度矩阵k集成到整体刚度矩阵K。每集成一个单元,该函数都将返回2N×2N的整体刚度矩阵K. 3. LinearTriangleElementStresses(E,NU,t,xi,yi,xj,yj,xm,ym,u)-- 该函数计算在平面应力情况下弹性模量为E、泊松比为NU、厚度为t、第一个节点坐标为(xi,yi)第二个节点坐标为(xj,yj)、第三个节点坐标为(xm,ym)以及单元位移矢量为u时的单元应力。该函数返回单元应力矢量。 函数源代码: function y = LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym) A = (xi*(yj-ym) + xj*(ym-yi) + xm*(yi-yj))/2;%三角形单元面积,单元节点应该按逆时针排序,保证每个三角形单元的面积都为正值(也可作为一个小函数:LinearTriangleElementArea) betai = yj-ym;betaj = ym-yi;betam = yi-yj;gammai = xm-xj; gammaj = xi-xm;gammam = xj-xi; B = [betai 0 betaj 0 betam 0 ; 0 gammai 0 gammaj 0 gammam ; gammai betai gammaj betaj gammam betam]/(2*A); %B为应变矩阵,其中betai=yi-ym,betaj=ym-yi,betam=yi-yj. gammai=xm-xj, gammaj=xi-xm, gammam=xj-xi. D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2]; %D为弹性矩阵,分为平面应力问题和平面应变问题 对于平面应力问题D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2]; 对于平面应变问题E1=E/(1-NU*NU),NU1=NU/(1-NU) y = t*A*B'*D*B;%单元刚度矩阵 function y = LinearTriangleAssemble(K,k,i,j,m) K(2*i-1,2*i-1) = K(2*i-1,2*i-1) + k(1,1); K(2*i-1,2*i) = K(2*i-1,2*i) + k(1,2); K(2*i-1,2*j-1) = K(2*i-1,2*j-1) + k(1,3); K(2*i-1,2*j) = K(2*i-1,2*j) + k(1,4); K(2*i-1,2*m-1) = K(2*i-1,2*m-1) + k(1,5); K(2*i-1,2*m) = K(2*i-1,2*m) + k(1,6); K(2*i,2*i-1) = K(2*i,2*i-1) + k(2,1); K(2*i,2*i) = K(2*i,2*i) + k(2,2); K(2*i,2*j-1) = K(2*i,2*j-1) + k(2,3); K(2*i,2*j) = K(2*i,2*j) + k(2,4); K(2*i,2*m-1) = K(2*i,2*m-1) + k(2,5); K(2*i,2*m) = K(2*i,2*m) + k(2,6); K(2*j-1,2*i-1) = K(2*j-1,2*i-1) + k(3,1); K(2*j-1,2*i) = K(2*j-1,2*i) + k(3,2); K(2*j-1,2*j-1) = K(2*j-1,2*j-1) + k(3,3); K(2*j-1,2*j) = K(2*j-1,2*j) + k(3,4); K(2*j-1,2*m-1) = K(2*j-1,2*m-1) + k(3,5); K(2*j-1,2*m) = K(2*j-1,2*m) + k(3,6); K(2*j,2*i-1) = K(2*j,2*i-1) + k(4,1); K(2*j,2*i) = K(2*j,2*i) + k(4,2); K(2*j,2*j-1) = K(2*j,2*j-1) + k(4,3); K(2*j,2*j) = K(2*j,2*j) + k(4,4); K(2*j,2*m-1) = K(2*j,2*m-1) + k(4,5); K(2*j,2*m) = K(2*j,2*m) + k(4,6); K(2*m-1,2*i-1) = K(2*m-1,2*i-1) + k(5,1); K(2*m-1,2*i) = K(2*m-1,2*i) + k(5,2); K(2*m-1,2*j-1) = K(2*m-1,2*j-1) + k(5,3); K(2*m-1,2*j) = K(2*m-1,2*j) + k(5,4); K(2*m-1,2*m-1) = K(2*m-1,2*m-1) + k(5,5); K(2*m-1,2*m) = K(2*m-1,2*m) + k(5,6); K(2*m,2*i-1) = K(2*m,2*i-1) + k(6,1); K(2*m,2*i) = K(2*m,2*i) + k(6,2); K(2*m,2*j-1) = K(2*m,2*j-1) + k(6,3); K(2*m,2*j) = K(2*m,2*j) + k(6,4); K(2*m,2*m-1) = K(2*m,2*m-1) + k(6,5); K(2*m,2*m) = K(2*m,2*m) + k(6,6); y = K;%对号入座,如前所述,每集成一次都将返回2N×2N的整体刚度矩阵K.此题为110×110 function y = LinearTriangleElementStresses(E,NU,t,xi,yi,xj,yj,xm,ym,u) A = (xi*(yj-ym) + xj*(ym-yi) + xm*(yi-yj))/2; betai = yj-ym;betaj = ym-yi;betam = yi-yj;gammai = xm-xj;gammaj = xi-xm;gammam = xj-xi; B = [betai 0 betaj 0 betam 0 ;

相关文档
最新文档