钢结构有限元分析

2 受料仓与给料机的钢结构有限元分析

2.1建立有限元模型

如图2.1破碎站主视图和图2.2破碎机布置图，它的工作过程是：卸料卡车间歇把最大入料粒度为1500mm的煤块倒入受料仓，受料仓存储大粒度煤块。刮板给料机把受料仓的大粒度的煤块连续的刮给破碎平台的破碎机。破碎机把最大入料粒度为1500mm 的煤块破碎成最大排料粒度为300mm的煤块，煤块由底部的传送带传出。

图2.1 破碎站主视图

图2.2 破碎机布置图

破碎站钢结构的弹性模量E＝200000MPa，泊松比μ=0.3，质量密度ρ=7.8×10-3kg/cm3。破碎站由支撑件H型钢和斜支撑(角钢)组成。在结构离散化时，由于角钢和其它部位铰接，铰接是具有相同的线位移，而其角位移不同。承受轴向力，不承受在其它方向的弯矩，相当于二力杆，所以H型钢用梁单元模拟，角钢用杆单元模拟。破碎站是由受料仓与给料机和破碎平台与控制室两部分组成，故计算时是分别对这两部分进行的。离散后，受料仓和给料机共686个单元，其中梁单元598 个，杆单元88个，节点总数为597个，有限元模型如图2.3和图2.4所示。

图2.3 受料仓与给料机有限元模型

图2.4 受料仓与给料机有限元模型俯视图

2.2载荷等效计算

2.2.1主要结构截面几何参数

破碎站主要结构采用H型钢梁，截面尺寸如图2.5所示，各截面横截面积A，截面

惯性矩I

y ，I

和极惯性矩I如下。

图2.5 截面尺寸

料仓及给料机支撑结构

料仓及给料机六根支撑立柱(H500×400×12×20)

A= 215.2mm2，I y＝101947×104mm4，I z＝21340×104mm4，I＝240×104mm4料仓B－B面横梁和给料机E－E、F－F面横梁(H400×300×12×20)

A=16320mm2，I y＝48026×104mm4，I z＝9005×104mm4，I＝181×104mm4料仓C－C面和D－D面横梁(H400×400×12×20)

A=20320mm2，I y＝62479×104mm4，I z＝21339×104mm4，I＝234×104mm4给料机两根纵梁(H550×400×12×20)

A=22120mm2，I y＝125678×104mm4，I z＝21341×104mm4，I＝243×104mm4给料机六根横梁(H400×400×12×20)

A=20320mm2，I y＝62479×104mm4，I z＝21339×104mm4，I＝234×104mm4其它横梁(H400×300×12×20)

A=16320mm2，I y＝48026×104mm4，I z＝9005×104mm4，I＝181×104mm4

斜支撑的横截面积

∠125×12：A＝2856mm2

∠75× 6：A＝864mm2

2.2.1实际载荷情况

给料机自

重载荷：

65000kg

相对应立柱梁单元局部坐标z轴的弯矩如图2.18所示，最大弯矩位于B－B面中风载作用面的立柱，底部最大弯矩0.126×108N·mm，在B－B面拉筋连结点处弯矩数值为-0.196×108N·mm，D－D面第一根水平横梁处弯矩为0.25563×108N·mm，

图2.18 局部坐标z轴的弯矩图

对应的应力如图2.19所示，B－B面中风载作用面立柱底部应力为19.625MPa，底部连结点处应力为30.605MPa，D－D面第一根水平横梁处应力为24.024MPa。

图2.19 局部坐标z轴的弯曲应力图

2.2.3两根纵梁结果

两根纵梁轴向力如图2.20，可见两根纵梁轴向应力很小，最大轴向力192337N，位于F面和纵梁与斜支撑的接触之间。这里的变形也是最大的。轴的方向与大部分的载荷

方向近似垂直。在斜支撑与纵梁连接到纵梁的前端只承受弯矩，不承受轴向力。

图2.20 两根纵梁轴向力图

两根纵梁轴向应力如图2.21所示，可见两根纵梁轴向应力很小，最大轴向应力＝13.066MPa，位于F面和纵梁前端之间。

图2.21 两根纵梁轴向应力图

相对于纵梁梁单元局部坐标y轴的弯矩如图2.22所示，其最小弯矩M

＝-0.176×

109N·mm，位于斜支撑与纵梁连接处，这里的弯矩最大驱动扭矩作用在悬臂梁顶端。

图2.22 局部坐标y轴的弯矩图

对应的应力如图2.23所示，最小应力＝-97.755MPa，位于F－F面处，。

图2.23 局部坐标y轴的弯曲应力图

＝0.174×相对于纵梁梁单元局部坐标Z轴的弯矩如图2.24所示，其最大弯矩M

108N·mm，这里是承受给料机尾部受仓压载荷，位置是纵梁和C_C截面相连接。

图2.24 局部坐标Z轴的弯矩图

对应的应力如图2.25所示，最小应力＝-47.911MPa。位于C－C面处。

图2.25 局部坐标Z轴的弯曲应力图

专业ABAQUS有限元建模经验笔记

基于ABAQUS的有限元分析和应用第一章绪论 1.有限元分析包括下列步骤： 2.为了将试验数据转换为输入文件，分析者必须清楚在程序中所应用的和由实验人员提供的材料数据的应力和应变的度量。 3.ABAQUS建模需注意以下内容： 4.对于许多包含过程仿真的大变形问题和破坏分析，选择合适的网格描述是非常重要的，需要认识网格畸变的影响，在选择网格时必须牢牢记住不同类型网格描述的优点。第二章ABAQUS基础 1.一个分析模型至少要包含如下的信息：离散化的几何形体、单元截面属性、材料数据、载荷和边界条件、分析类型和输出要求。 ①离散化的几何形体：模型中所有的单元和节点的集合称为网格。 ②载荷和边界条件： 2.功能模块：（1）Assembly（装配）：一个ABAQUS模型只能包含一个装配件。（2）Interaction（相互作用）：相互作用与分析步有关，这意味着用户必须规定相互作用是在哪些分析步中起作用。（3）Load（载荷）：载荷和边界条件与分析步有关，这意味着用户指定载荷和边界条件是在哪些分析步中起作用。（4）Job（作业）：多个模型和运算可以同时被提交并进行监控。 3.量纲系统 ABAQUS没有固定的量纲系统，所有的输入数据必须指定一致性的量纲系统，常用的一致性量纲系统如下：

4.建模要点（1）创建部件：设定新部件的大致尺寸的原则必须是与最终模型的最大尺寸同一量级。（2）用户应当总是以一定的时间间隔保存模型数据（例如，在每次切换功能模块时）。（3）定义装配：在模型视区左下角的三向坐标系标出了观察模型的方位。在视区中的第2个三向坐标系标出了坐标原点和整体坐标系的方向（X，Y和Z轴）。（4）设置分析过程：（5）在模型上施加边界条件和荷载：用户必须指定载荷和边界条件是在哪个或哪些分析步中起作用。所有指定在初始步中的力学边界条件必须赋值为零，该条件是在ABAQUS/CAE中自动强加的。在许多情况下，需要的约束方向并不一定与整体坐标方向对齐，此时用户可定义一个局部坐标系以施加边界条件。在ABAQUS中，术语载荷通常代表从初始状态开始引起结构响应发生变化的各种因素，包括：集中力、压力、非零边界条件、体力、温度（与材料热膨胀同时定义）。

matlab有限元分析实例

MATLAB: MATLAB是美国MathWorks公司出品的商业数学软件，用于数据分析、无线通信、深度学习、图像处理与计算机视觉、信号处理、量化金融与风险管理、机器人，控制系统等领域。 MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室），软件主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。 MATLAB有限元分析与应用:

《MATLAB有限元分析与应用》是2004年4月清华大学出版社出版的图书，作者是卡坦，译者是韩来彬。内容简介：《MATLAB有限元分析与应用》特别强调对MATLAB的交互应用，书中的每个示例都以交互的方式求解，使读者很容易就能把MATLAB用于有限分析和应用。另外，《MATLAB有限元分析与应用》还提供了大量免费资源。《MATLAB有限元分析与应用》采用当今在工程和工程教育方面非常流行的数学软件MATLAB来进行有限元的分析和应用。《MATLAB有限元分析与应用》由简单到复杂，循序渐进地介绍了各种有限元及其分析与应用方法。书中提供了大量取自机械工程、土木工程、航空航天工程和材料科学的示例和习题，具有很高的工程应用价值。

基于有限元法和极限平衡法的边坡稳定性分析

目录摘要 (1) 1引言 (1) 2 简要介绍有限元和极限平衡方法 (1) 3影响边坡稳定性的因素 (2) 3.1水位下降速度的影响 (2) 3.2 不排水粘性土对边坡失稳的影响 (5) 3.3 裂缝位置的影响 (9) 4 总结和结论 (12)

基于有限元法和极限平衡法的边坡稳定性分析摘要：相较于有限元分析法，极限平衡法是一种常用的更为简单的边坡稳定性分析方法。这两种方法都可用于分析均质和不均质的边坡，同时考虑了水位骤降，饱和粘土和存在张力裂缝的条件。使用PLAXIS8.0(有限元法)和SAS-MCT4.0(极限平衡方法)进行了分析，并对两种方法获得的临界滑动面的安全系数和位置进行了比较。关键词：边坡稳定；极限平衡法；有限元法；PLAXIS；SAS-MCT 1.引言近年来，计算方法，软件设计和高速低耗硬件领域都得到快速发展，特别是相关的边坡稳定性分析的极限平衡法和有限元方法。但是，使用极限平衡方法来分析边坡，可能会在定位临界滑动面(取决于地质)时出现几个计算困难和前后数值不一致，因此要建立一个安全系数。尽管极限平衡法存在这些固有的局限性，但由于其简单，它仍然是最常用的方法。然而，由于个人电脑变得更容易获得，有限元方法已越来越多地应用于边坡稳定性分析。有限元法的优势之一是，不需要假设临界破坏面的形状或位置。此外，该方法可以很容易地用于计算压力，位移，路堤空隙压力，渗水引起的故障，以及监测渐进破坏。邓肯(1996年)介绍了一个综合观点，用极限平衡和有限元两种方法对边坡进行分析。他比较了实地测量和有限元分析的结果，并且发现一种倾向，即计算变形大于实测变形。Yu 等人(1998年)比较了极限平衡法和严格的上、下界限法对于简单土质边坡的稳定性分析的结果，同时，他们也将采用毕肖普法和利用塑性力学上、下限原理的界限法得到的结果进行了比较。Kim等人(1999年)同时使用极限平衡法和极限分析法对边坡进行分析，发现对于均质土边坡，得自两种方法的结果大体是一致的，但是对于非均质土边坡还需要进行进一步分析工作。Zaki(1999年)认为有限元相对于极限平衡法更显优势。Lane和Griffiths (2000年) 提出一个看法，用有限元方法在水位骤降条件下评价边坡的稳定性，应绘制出适用于实际结构的操作图表。Rocscience有限公司(2001年)提出了一个文件，概述了有限元分析方法的能力，并通过与各种极限平衡方法的结果比较，提出了有限元方法更为实用。Kim等人(2002年)用上、下界限法和极限平衡法分析了几处非均质土体且几何不规则边坡的剖面。这两种方法给出了类似有限元分析法产生的安全系数，临界滑动面位置。 2.简要介绍有限元和极限平衡方法有限元法(FEM)是一个应用于科学和工程中，求解微分方程和边值问题的数值方法。进一步的细节，读者可参考Clough和Woodward(1967年)，Strang和Fix(1973年)，Hughes(1987年)，Zienkiewicz和Taylor(1989年)所做的研究工作。 PLAXIS 8版(Brinkgreve 2002年)是一个有限元软件包，应用于岩土工程二维的变形和折稳定性分析。该程序可以分析自然成型或人为制造的斜坡问题。安全系数的确定使用c

有限元分析及应用大课后复习

有限元分析及应用作业报告

目录有限元分析及应用作业报告....................................... I 目录 ........................................................ II 试题1 . (1) 一、问题描述 (1) 二、几何建模与分析 (2) 三、第1问的有限元建模及计算结果 (2) 四、第2问的有限元建模及计算结果 (7) 五、第3问的有限元建模及计算结果 (13) 六、总结和建议 (16) 试题5 (17) 一、问题的描述 (17) 二、几何建模与分析 (18) 三、有限元建模及计算结果分析 (18) 四、总结和建议 (26) 试题6 (27) 一、问题的描述 (27) 二、几何建模与分析 (27) 三、有限元建模及计算结果分析 (27) 五、总结和建议 (35)

试题1 一、问题描述图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； 2)分别采用不同数量的三节点常应变单元计算； 3)当选常应变三角单元时，分别采用不同划分方案计算。图1-1模型示意图及划分方案

二、几何建模与分析图1-2力学模型由于大坝长度>>横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看作无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作用于大坝上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力，满足平面应变问题的载荷条件。因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况及方向如图1-2所示，建立几何模型，进行求解。假设大坝的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3 三、第1问的有限元建模本题将分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算。1）设置计算类型：两者因几何条件和载荷条件均满足平面应变问题，故均取Preferences为Structural 2）选择单元类型：三节点常应变单元选择的类型是PLANE42（Quad 4node42），该单元属于是四节点单元类型，在网格划分时可以对节点数目控制使其蜕化为三节点单元；六节点三角形单元选择的类型是PLANE183（Quad 8node183），该单元属于是八节点单元类型，在网格划分时可以对节点数目控制使其蜕化为六节点单元。因研究的问题为平面应变问题，故对Element behavior（K3）设置为plane strain。 3）定义材料参数：按以上假设大坝材料为钢，设定：ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 → OK 4）生成几何模型： a. 生成特征点：ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS→依次输入三个点的坐标：

abaqus有限元分析过程

一、有限单元法的基本原理有限单元法（The Finite Element Method）简称有限元（FEM），它是利用电子计算机进行的一种数值分析方法。它在工程技术领域中的应用十分广泛，几乎所有的弹塑性结构静力学和动力学问题都可用它求得满意的数值结果。有限元方法的基本思路是：化整为零，积零为整。即应用有限元法求解任意连续体时，应把连续的求解区域分割成有限个单元，并在每个单元上指定有限个结点，假设一个简单的函数（称插值函数）近似地表示其位移分布规律，再利用弹塑性理论中的变分原理或其他方法，建立单元结点的力和位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程组，从而求解结点的位移分量. 进而利用插值函数确定单元集合体上的场函数。由位移求出应变, 由应变求出应力二、ABAQUS有限元分析过程有限元分析过程可以分为以下几个阶段 1.建模阶段: 建模阶段是根据结构实际形状和实际工况条件建立有限元分析的计算模型――有限元模型，从而为有限元数值计算提供必要的输入数据。有限元建模的中心任务是结构离散，即划分网格。但是还是要处理许多与之相关的工作：如结构形式处理、集合模型建立、单元特性定义、单元质量检查、编号顺序以及模型边界条件的定义等。

2.计算阶段:计算阶段的任务是完成有限元方法有关的数值计算。由于这一步运算量非常大，所以这部分工作由有限元分析软件控制并在计算机上自动完成 3.后处理阶段: 它的任务是对计算输出的结果惊醒必要的处理，并按一定方式显示或打印出来，以便对结构性能的好坏或设计的合理性进行评估，并作为相应的改进或优化，这是惊醒结构有限元分析的目的所在。下列的功能模块在ABAQUS/CAE操作整个过程中常常见到，这个表简明地描述了建立模型过程中要调用的每个功能模块。 “Part（部件）用户在Part模块里生成单个部件，可以直接在ABAQUS/CAE环境下用图形工具生成部件的几何形状，也可以从其它的图形软件输入部件。 Property（特性）截面（Section）的定义包括了部件特性或部件区域类信息，如区域的相关材料定义和横截面形状信息。在Property模块中，用户生成截面和材料定义，并把它们赋于(Assign)部件。 Assembly（装配件）所生成的部件存在于自己的坐标系里，独立于模型中的其它部件。用户可使用Assembly模块生成部件的副本（instance），并且在整体坐标里把各部件的副本相互定位，从而生成一个装配件。一个ABAQUS模型只包含一个装配件。

Matlab有限元分析操作基础

Matlab 有限元分析20140226 为了用Matlab 进行有限元分析，首先要学会Matlab 基本操作，还要学会使用Matlab 进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵 11221 21200k k k k k k k k -????-????--+??

2. Matlab有限元分析的基本操作（1）单元划分（选择何种单元，分成多少个单元，标号）（2）构造单元刚度矩阵（列出…）（3）组装系统刚度矩阵（集成整体刚度矩阵）（4）引入边界条件（消除冗余方程）（5）解方程（6）后处理（扩展计算）

3. Matlab有限元分析实战【实例1】

分析：步骤一：单元划分

步骤二：构造单元刚度矩阵 >>k1=SpringElementStiffness(100) >>…?

步骤三：构造系统刚度矩阵 a) 分析SpringAssemble库函数function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？今天的系统刚度矩阵是什么？因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ?? ?? -- ?? ？

有限元分析及优化设计

《有限元分析及优化设计》实验指导书桂林电子科技大学机电工程学院庄未编 2012年05月

实验一:平面问题的结构分析计算 1.实验目的 ?了解ANSYS软件的基本功能与应用范围； ?熟悉在计算机上运用ANSYS软件的基本步骤和方法； ?结合具体平面问题实例,利用ANSYS软件进行计算分析； ?时间许可,可对上述实例利用有限元方法进行计算，并与ANSYS计算结果进行分析比较. 2.实验内容 1. 结合具体平面问题实例,利用ANSYS软件进行计算分析； 2. 利用ANSYS软件进行建模，并施加约束和载荷； 3 对计算结果进行比较分析与讨论； 4. 时间许可,可对上述实例利用ANSYS的非交互模式（Batch Mode/命令流的方式）再进行一次计算，并与用ANSYS交互模式的计算结果进行分析比较. 3.实验预习报告内容要求实验预习报告在实验前写好，其主要内容应包括：复习有限元法基本原理、解题方法与步骤等,建立有限元模型应包含的内容；提供具体平面问题的结构简图,画出计算模型；对给定的平面问题实例的结果进行预估,以供计算后进行比较讨论用； 4.上机实践举例一）如图1所示的6结点4单元平面应力平板问题.各三角形单元的直角边的长度为α=10mm,假设平板的厚度t=5mm,材料均匀,其弹性模量E=200GPa, 泊松比μ=0.3.今在结点1处,竖直向下作用一个力P=1,若不计平板重量( 即设容重γ=0 ).利用ANSYS软件进行分析。

图1 二）、求解下图所示的平面问题。图2 实验二：轴对称实体结构静力有限元分析 1. 实验目的 ? 了解ANSYS 软件的基本功能与应用范围； ? 熟悉在计算机上运用ANSYS 软件的基本步骤和方法； ? 结合具体实体问题实例,利用ANSYS 软件进行计算分析； ? 时间许可,可对上述实例利用有限元方法进行计算，并与ANSYS 计算结果进行分析比较.

基于ansys的连杆机构的有限元分析

目录摘要 ............................................................................................ 错误！未定义书签。Abstract (2) 第一章分析方法和研究对象 ........................................... 错误！未定义书签。 1.1 有限单元法的概述....................................................... 错误！未定义书签。 1.1.1 有限单元法的历史 (4) 1.1.2 有限单元法的基本概念 (4) 1.2 ANSYS软件简介 (4) 1.2.1 ANSYS主要应用领域 (4) 1.2.2 ANSYS操作界面 (5) 1.2.3 ANSYS的主要功能 (6) 1.2.4 ANSYS主要特点 (7) 1.3 曲柄滑块机构简介 (7) 1.3.1 曲柄滑块定义 (8) 1.3.2 曲柄滑块机构特性应用以及分类 (8) 第二章曲柄滑块机构的求解 (10) 2.1 曲柄滑块机构的问题描述 (10) 2.2 曲柄滑块机构问题的图解法 (10) 2.2.1 图解法准备工作 (11) 2.2.2 图解法操作步骤 (11) 第三章有限元瞬态动力学概述 (14) 3.1 有限元瞬态动力学定义 (14) 3.2 瞬态动力学问题求解方法........................................... 错误！未定义书签。 3.2.1 完全法 (14) 3.2.2 模态分析法 (14) 3.2.2 缩减法 (15) 3.1 有限元结构静力学分析基本概念 (15) 3.1 有限元结构静力学分析步骤 (16) 第四章曲柄滑块的有限元瞬态动力学分析 (17) 4.1 曲柄滑块机构瞬态简要概述 (17) 4.2曲柄滑块有限元瞬态动力学分析步骤 (18)

abaqus 有限元分析(齿轮轴)

Abaqus分析报告（齿轮轴）名称：Abaqus齿轮轴姓名：班级：学号：指导教师：

一、简介所分析齿轮轴来自一种齿轮泵，通过用abaqus软件对齿轮轴进行有限元分析和优化。齿轮轴装配结构图如图1，分析图1中较长的齿轮轴。图1.齿轮轴装配结构图二、模型建立与分析通过part、property、Assembly、step、Load、Mesh、Job等步骤建立齿轮轴模型，并对其进行分析。 1.part 针对该齿轮轴，拟定使用可变型的3D实体单元，挤压成型方式。 2.材料属性材料为钢材，弹性模量210Gpa，泊松比0.3。

3.截面属性截面类型定义为solid，homogeneous。 4.组装组装时选择dependent方式。 5.建立分析步本例用通用分析中的静态通用分析（Static，General）。 6.施加边界条件与载荷对于齿轮轴，因为采用静力学分析，考虑到前端盖、轴套约束，而且根据理论，对受力部分和轴径突变的部分进行重点分析。边界条件：分别在三个轴径突变处采用固定约束，如图2。载荷：在Abaqus中约束类型为pressure，载荷类型为均布载荷,分别施加到齿轮接触面和键槽面，根据实际平衡情况，两力所产生的绕轴线的力矩方向相反，大小按比例分配。均布载荷比计算：矩形键槽数据：长度：8mm、宽度：5mm、高度：3mm、键槽所在轴半径：7mm 键槽压力面积：S1 = 8x3=24mm2 平均受力半径：R1=6.5mm 齿轮数据：= 齿轮分度圆半径：R2 =14.7mm、压力角：20°、单个齿轮受力面积：S2 ≈72mm2 通过理论计算分析，S1xR1xP1=S2xR2xP2，其中，P1为键槽均布载荷

优化设计有限元分析总结

目录目录 (1) 1. 优化设计基础 (2) 1.1 优化设计概述 (2) 1.2 优化设计作用 (3) 1.3 优化设计流程 (3) 2. 问题描述 (4) 3. 问题分析 (5) 4. 结构静力学分析 (6) 4.1 创建有限元模型 (6) 4.2 创建仿真模型并修改理想化模型 (7) 4.3 定义约束及载荷 (7) 4.4 求解 (8) 5. 结构优化分析 (9) 5.1 建立优化解算方案 (9) 5.2 优化求解及其结果查看 (11) 6. 结果分析 (13) 7. 案例小结 (14)

1.优化设计基础 1.1优化设计概述优化设计是将产品/零部件设计问题的物理模型转化为数学模型，运用最优化数学规划理论，采用适当的优化算法，并借助计算机和运用软件求解该数学

模型，从而得出最佳设计方案的一种先进设计方法，有限元被广泛应用于结构设计中，采用这种方法任意复杂工程问题，都可以通过它们的响应进行分析。如何将实际的工程问题转化为数学模型，这是优化设计首先要解决的关键问题，解决这个问题必须要考虑哪些是设计变量，这些设计变量是否受到约束，这个问题所追求的结果是在优化设计过程要确定目标函数或者设计目标，因此，设计变量、约束条件和目标函数是优化设计的3个基本要素。因此概括来说，优化设计就是：在满足设计要求的前提下，自动修正被分析模型的有关参数，以到达期望的目标。 1.2优化设计作用以有限元法为基础的结构优化设计方法在产品设计和开发中的主要作用如下： 1）对结构设计进行改进，包括尺寸优化、形状优化和几何拓扑优化。2）从不合理的设计方案中产生出优化、合理的设计方案，包括静力响应优化、正则模态优化、屈曲响应优化和其他动力响应优化等。 3）进行模型匹配，产生相似的结构响应。 4）对系统参数进行设别，还可以保证分析模型与试验结果相关联。 5）灵敏度分析，求解设计目标对每个设计变量的灵敏度大小。 1.3优化设计流程不同的优化软件其操作要求及操作步骤大同小异。一般为开始、创建有限元模型、创建仿真模型、定义约束及载荷，然后进行结构分析，判断是否收

支架的有限元分析ABAQUS

支架的线性静力学分析实例：建模和分析计算在此实例中读者将学习ABAQUS/CAE的以下功能。 1) Sketch功能模块:导人CAD二维图形，绘制线段、圆弧和倒角，添加尺寸，修改平面图，输出平面图。 2) Part功能模块:通过拉伸来创建几何部件，通过切割和倒角未定义几何形状。 3) Property功能模块:定义材料和截面属性。 4) Mesh功能模块:布置种子，分割实体和面，选择单元形状、单元类型、网格划分技术和算法，生成网格，检验网格质量，通过分割来定义承受载荷的面。 5) Assembly功能模块:创建非独立实体。 6) Step功能模块:创建分析步，设置时间增量步和场变量输出结果。 7) Interaction功能模块:定义分布榈合约束(distributing coupling constraint)。 8) Load功能模块:定义幅值，在不同的分析步中分别施加面载荷和随时间变化的集中力，定义边界条件。 9) Job功能模块:创建分析作业，设置分析作业的参数，提交和运行分析作业，监控运行状态。 10) Visualization功能模块:后处理的各种常用功能。结构静力学分析（static analysis）是有限元法的基本应用领域，适用于求解惯性及阻尼对结构响应不显著的问题。主要用来分析由于稳态外载荷引起的位移，应力和应变等。本章的静力学分析实例按照ABAQUS工程分析的流程对支架进行线性静力学分析，通过实例基本掌握了分析的流程，同时了解接触的定义。 1.问题描述所示的支架，一端牢固地焊接在一个大型结构上，支架的圆孔中穿过一个相对较软的杆件，圆孔和杆件用螺纹连接。材料的弹性模量E=2100000MPa,泊松比为0.3。

《有限元基础教程》_【MATLAB算例】3.3.7(2)__三梁平面框架结构的有限元分析(Beam2D2Node)

【MA TLAB 算例】3.3.7(2) 三梁平面框架结构的有限元分析 (Beam2D2Node) 如图3-19所示的框架结构，其顶端受均布力作用，结构中各个截面的参数都为：113.010Pa E =?，746.510I m -=?，426.810A m -=?。试基于MA TLAB 平台求解该结构的节点位移以及支反力。图3-19 框架结构受一均布力作用解答：对该问题进行有限元分析的过程如下。（1）结构的离散化与编号将该结构离散为3个单元，节点位移及单元编号如图3-20所示，有关节点和单元的信息见表3-5。（a ）节点位移及单元编号

（b）等效在节点上的外力图3-20 单元划分、节点位移及节点上的外载（2）各个单元的描述首先在MA TLAB环境下，输入弹性模量E、横截面积A、惯性矩I和长度L，然后针对单元1，单元2和单元3，分别二次调用函数Beam2D2Node_ElementStiffness，就可以得到单元的刚度矩阵k1(6×6)和k2(6×6)，且单元2和单元3的刚度矩阵相同。 >> E=3E11; >> I=6.5E-7; >> A=6.8E-4; >> L1=1.44; >> L2=0.96; >> k1=Beam2D2Node_Stiffness(E,I,A,L1); >> k2=Beam2D2Node_Stiffness(E,I,A,L2); （3）建立整体刚度方程将单元2和单元3的刚度矩阵转换成整体坐标下的形式。由于该结构共有4个节点，则总共的自由度数为12，因此，结构总的刚度矩阵为KK(12×12)，对KK清零，然后两次调用函数Beam2D2Node_Assemble进行刚度矩阵的组装。 >> T=[0,1,0,0,0,0;-1,0,0,0,0,0;0,0,1,0,0,0;0,0,0,0,1,0;0,0,0,-1,0,0;0,0,0,0,0,1] ; >> k3=T'*k2*T; >> KK=zeros(12,12); >> KK=Beam2D2Node_Assemble(KK,k1,1,2);

MD Nastran突破有限元分析的极限

MD Nastran突破有限元分析的极限作者：MSC.Software公司来源：汽车制造业有限元法FEM分析变得日益复杂，同时有限元分析模型的大小和细节设计要求也在不断增加。尤其是在汽车行业，这一趋势尤其明显。项目背景由数百万个单元和数百万的自由度组成的有限元网格的模型已经变得司空见惯，然而模型的尺寸仍在不断地增加。由于数学方法和软件工程学技术的改进，有限元法程序的工作效率和计算能力也在不断提升，同时构建模型和网格划分软件技术的飞速进步使模型的生成变得更加方便快捷。数年前，发动机引擎气缸体的网格划分需要几个月的时间，而现在只是几个小时的问题。德国汽车制造商宝马公司是大范围使用虚拟仿真技术的公司之一。在宝马公司和其他一些制造商中，为了缩短研发周期，减少物理样机和物理试验的次数，完整的汽车模型得到了最优化的使用，其基础便是日益复杂的有限元仿真模型，包括对噪音和舒适度的刚性评定、乘客安全性和空气动力学仿真等。在数值计算方法方面，使用了隐式线性分析和显式非线性瞬态分析。图1 “后天之模型”的基础是宝马X3汽车的车体早在2007年初，宝马公司便对计算机辅助工程CAE的流程重新进行了检测，以便发现将来可能由仿真模型尺寸增加引起的瓶颈问题。宝马公司的车体和零部件设计小组开发了迄今为止最大的有限元法模型作为基准测试的考题模型，被冠以“后天之模型（Model of the

Day After Tomorrow）”的名称。小组成员丹尼尔·海泽尔博士表示，“对我们来说，在标准的硬件和软件设备上进行此次基准测试是非常重要的，使用当前的基础设施解决基准模型问题的目的，并不是为了要减少计算时间，而是为了识别理论极限和当前方法的瓶颈。” 基准考题的目的是为了寻找标准分析（双载荷工况条件下的线性静态分析）中进行有限元法分析基本步骤的极限和时间： 1. 读取输入数据，对它们进行分类、制成表格，并进行一致性检查； 2. 计算单元刚体矩阵，并集成一个整体刚体矩阵； 3. 计算位移和应力数据； 4. 输出结果。宝马公司提出的问题是有限元分析还能应对这一增长趋势多长时间？用“后天之模型”作为考题的目的是如何突破近10年间所要面临的硬件和软件极限问题。MSC.Software公司同美国国际商用机器IBM公司合作，能够在短短的几个月的时间内解决这一问题。在一份用该模型分析的详细报告中，项目成员彼得·沙尔茨和杰拉德·希姆莱（MSC.Software公司），丹尼尔·海泽尔（宝马汽车制造公司）和D·皮特施（IBM公司）详细介绍了他们实现宝马公司苛刻要求的方法。图2 BMW X3减振器支座外壳模型（蓝色），MODAW部分描绘图（黄色）软、硬件的发展大多数有限元法分析程序都存在计算能力不在最佳状态的情形。1957年，雷W克拉夫和他的学生在一台内存只有16位的IBM701计算机上开发出了后来成为有限元法的程序。方程式大约在40个以上的问题需要out of core（即数据不全部存储在内存中，而是存储在硬盘的临时文件夹中）求解逻辑，这意味着要借助二级存储介质。10年之后，Nastran软件被开发出来之后，要求条件也非常类似。软件客户美国国家航空航天局（NASA）要求开

车架的有限元分析及优化

车架的有限元分析及优化作者：马迅盛…文章来源：湖北汽车工业学院点击数：1687 更新时间：2008-8-5 有限元法将设计人员丰富的实践经验与计算机高速精确的计算完美地结合在一起，大大提高了设计计算精度，缩短了产品开发时间。概念设计阶段车架的结构方案参考某一同类型车架，考虑到车身安装和其他总成的布置，将概念设计阶段的车架大致结构拟定如下：选用框架式平行梯形车架结构，由2根左右分开的纵梁和8根横梁组成，全长6.3m，宽0.8m，轴距3.65m。各梁的大致形状尺寸及板材厚度如表1所示。除第3、4根横梁外，其他各横梁的尺寸与参考的同类型车架几乎相同。由于参考车架的第3、4根横梁为上下两片形状复杂的钢板组合而成，无法用梁单元模拟，在概念车架中将之改用两根方型截面的等直梁代替。第1、6横梁为非等截面梁，其宽和高分别由两个尺寸表示。参考车架纵梁的前后两段和中间段的连接采用的是线性渐变的截面，在概念车架中用一等直梁来代替，等直梁的高度等于渐变梁的中间高度。纵横梁上所有的孔及连接板都不予以考虑。车架的有限元模型为了后续的优化设计，必须对车架进行参数化建模。选择表1中车架纵横梁的截面尺寸为模型参数，先建立左半个车架的几何模型，选用ANSYS中的二节点12自由度梁单元BEAM188号单元采用不同的梁单

元截面形式对其进行网格剖分；再将左边的几何模型和网格模型进行映射得到右边车架模型，最终合并对称面上的节点使左右车架模型“牢固的”“粘结起来”。在ANSYS中用BEAM188单元实施网格剖分时，为了保证单元的正确方向，应事先定义该单元的方向点并检查所要剖分的线的法向。单元截面形状和偏置量需用命令SECTYPE、SECOFFSET和SECDATA设定。单元总数为312，节点总数为626。网格剖分并映射后车架模型如图1所示。图中显示出了梁单元的截面形状。图1 车架的有限元模型边界条件车架刚度有多种，其中最重要的是车架的弯曲刚度和扭转刚度。参照车架的刚度试验方法确定车架弯扭刚度的边界条件。 1.弯曲工况的边界条件计算时约束前后桥在车架纵梁上的竖直投影点的垂直位移，让车架形成一简支梁结构，并在前后支承点中点处加一垂直向下的力，让车架产生纯弯曲变形，如图2所示。

用ANSYS进行四连杆机构的有限元分析

用ANSYS进行四连杆机构的有限元分析作者：谭辉日期：08年3月6日分析目的 1、利用ANSYS对典型的四连杆机构进行分析，主要包含各点的轨迹分析，例如X和Y方向的位移等。 2、为五连杆和六连杆机构的分析提供可行的分析方法以及原型代码。问题简述分析主动杆1绕节点1旋转一周时节点4的运动轨迹，杆2和杆3为从动杆，具体问题见下图：

分析思路 1、根据分析目的，在ANSYS选用link1单元进行单元建模，主要考虑是link1单元具有X和Y方向的自由度，可以获得各个节点的位移轨迹。之后可以用梁单元等实现更高级的分析目的，例如获得杆上的力，位移，加速度等相关信息。 2、该模型结构简单，可以利用直接建模方法进行有限元系统建模，主要命令：N，E。 3、利用自由度耦合对重合节点进行建模，例如节点2和节点3、节点4 和节点5进行建模，主要命令：cpintf，利用该命令可以一次性将重合节点生成自由度耦合。 4、利用表数组对于杆1（主动杆）的节点2进行瞬态边界条件的载荷施加，分析类型为瞬态分析，主要命令：*dim，d等。 5、生成节点位移的对应变量，从而获得节点4的随时间的位移曲线，主要命令：nsol，plvar等。命令流如下行号命令符号注释结束上一次的分析 1finish ！清除数据库，并读取启动配置文件2/clear,start ！ 3 ！设置图形显示的背景颜色 4/color,pbak,on,1,5 ！ 5 ！

6/units,si ！设置单位制：国际单位制 7*afun,deg ！设置三角函数运算采用度为单位 8 ！ 9/prep7 ！进入前处理模块 10et,1,link1 ！设置单元类型：link1 11mp,ex,1,2.07e11 ！设置材料的弹性模量 12r,1,1 ！设置单元的实常数，面积为1 13n,1,0,0,0 ！在（0，0，0）处建立节点1 14n,2,3,0,0 ！在（3，0，0）处建立节点2 15n,3,3,0,0 ！在（3，0，0）处建立节点3，和节点2重合 16n,4,8,7,0 ！在（8，7，0）处建立节点4 17n,5,8,7,0 ！在（8，7，0）处建立节点4，和节点4重合 18n,6,10,0,0 ！在（10，0，0）处建立节点6 19e,1,2 ！建立单元1（连接节点1和2） 20e,3,4 ！建立单元2（连接节点3和4） 21e,5,6 ！建立单元3（连接节点5和6） 22 ！ 23cpintf,all,1e-3 ！对于重合节点一次性的建立耦合自由度，容差1e－3 24 ！ 25/pnum,node,1 ！显示节点编号 26/pnum,elem,1 ！显示单元编号 27eplot ！显示单元

风力发电机组轮毂极限强度的有限元分析

风力发电机组轮毂极限强度的有限元分析文章是基于有限元理论，对兆瓦级风力发电机组的轮毂进行强度及疲劳计算。轮毂是风力发电机中的重要组成部分，铸造而成，是将机械能转换为电能的核心部件，其形状复杂，轮毂的设计质量会直接影响到整个机组的正常运行及使用寿命，在其受复杂风载荷的作用下，其强度和疲劳耐久性成为此行业关注的焦点。此分析利用大型有限元分析软件Ansys对轮毂模型分析。模型中包含轮毂、主轴及叶片，从轮毂的应力分布情况，从中找出最危险的部位，为轮毂的设计提供可靠依据。标签：风力发电机；轮毂；有限元分析；极限强度 1 绪论 1.1 课题研究背景经济发展过程中，我国作为世界上人口最多的发展中国家，能源消耗量不断增加，传统化石能源无以为继，面临的能源开发利用的资源约束越来越多，环境压力也越来越大。如今，生态环境承载能力弱、资源相对紧张。传统能源利用导致的环境问题越来越严重，以及全国范围内的雾霾天气都在提醒我们要努力做到全面、协调、可持续发展，以符合当今国情。在众多的可再生能源中，风能以其巨大的优越性和发展潜力受到人们的瞩目。 1.2 轮毂在大型风力发电机组的重要性在大型风力发电机组中，轮毂是核心构件，其不仅承担着与驱动连的链接，而且将叶片所受的风载荷通过主轴传递给齿轮箱，承担着风力发电机组容量增大而带来的更大的负荷。它需要有足够的强度和刚度，以保证机组在各种工况下能正常运行。由此可看出轮毂在风力发电机组的设计和制造过程中的重要性。 2 轮毂的强度校核计算 2.1 轮毂模型介绍轮毂模型结构见图1 此机组风轮由三片叶片对称安装在轮毂上构成，叶片间的夹角为120°。利用CAD绘图软件Solidworks，绘制了轮毂的三维实体几何简化模型。在保证计算精度的前提下，由于小的孔类、圆角及小凸台类结构对计算结果影响很小并且不是关键部位，已经略去。叶片产生的气动载荷以及由于风轮旋转和机舱对风轮转动引起的离心力、惯性力和重力通过三片叶片连接点传递到轮毂上，这些载荷和轮毂自身的重力构成了轮毂载荷。最终，轮毂简化后的几何模型如图1所示。

基于ANSYS Workbench的定位卡锁机构有限元分析

基于ANSYS Workbench的定位卡锁机构有限元分析摘要本文首先在Pro/E中建立了定位卡锁机构受最大外力时的简化模型，然后将该模型导入到ANSYS Workbench 13平台中进行了有限元模型的分析求解，最后结合求解结果用第四强度理论对定位卡锁机构各零件进行了强度校核，同时对该定位卡锁机构的改进提出了建议。关键词定位卡锁机构；有限元分析在某工程项目中应用的定位卡锁机构承担着为某输送设备准确定位的作用。由于该输送设备运行一个周期位就要启停一次，启停工作由定位卡锁机构配合实现。定位卡锁机构收回，输送设备开始运转，一个周期位后电机停转，定位卡锁机构伸出，进入与之配合的凹槽使输送设备完全停位。因此，定位卡锁机构成为该输送设备的关键部件，是保证输送设备正常工作的必备条件。所以，对定位卡锁机构的研究与分析有着重要的意义。定位卡锁机构在伸出状态受最大外力时，其所受最大应力不应超过材料的许用应力是保证定位卡锁机构实现其功能的充分条件。为了保证定位卡锁机构的工作可靠性，本文利用ANSYS Workbench对该机构进行有限元分析，研究在定位卡锁机构受最大外力时的受力及变形情况，并依据理论知识对其强度进行校核。 1 定位卡锁机构模型的建立与导入在对定位卡锁机构进行有限元分析之前，首先应建好定位卡锁机构的三维模型。一般在整个有限元分析的过程中，几何建模的工作量占据了非常多的时间，同时也是非常重要的过程[2]。ANSYS Workbench 13中，建模工作主要由ANSYS Workbench 自带的几何建模工具Design Modeler模块完成。对于小型或简单模型的建立可以直接在Design Modeler模块中建模，这样避免了从CAD系统中导入ANSYS的模型可能不能直接进行网格划分，需进行大量修补完善工作的麻烦。对于零部件较多的装配体的建模，通常先利用专业的三维建模软件完成模型的建立，然后再把它导入到ANSYS中进行分析。这样，工程技术人员就可以使用自己擅长的CAD软件建好模型，从而避免了重复现有CAD模型的劳动。本文采用PTC公司的Pro/Engineer对定位卡锁机构进行三维建模。定位卡锁机构简化模型由液压缸、卡锁活塞杆、端盖、螺塞、螺钉组成，建好的三维模型如图1所示。建好后的三维模型可以在Pro/E中直接导入到ANSYS Workbench 13 中进行有限元分析。图1 定位卡锁机构的三维模型 2 定位卡锁机构的有限元分析 2.1 定义模型材料属性

Matlab有限元分析操作基础共11页

Matlab有限元分析20140226 为了用Matlab进行有限元分析，首先要学会Matlab基本操作，还要学会使用Matlab进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵

3. Matlab有限元分析实战【实例1】

分析：步骤一：单元划分

>>k1=SpringElementStiffness(100)

a) 分析SpringAssemble库函数 function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？今天的系统刚度矩阵是什么？因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ????-- ??？