. 一、题目 如图1所示，一个厚度均匀的三角形薄板，在顶点作用沿板厚方向均匀分布的竖向载荷。已知：P=150N/m ，E=200GPa ，=0.25，t=0.1m ，忽略自重。试计算薄板的位移及应力分布。 要求： 1. 编写有限元计算机程序，计算节点位移及单元应力。（划分三角形 单元，单元数不得少于30个）； 2. 采用有限元软件分析该问题（有限元软件网格与程序设计网格必 须一致），详细给出有限元软件每一步的操作过程，并将结果与程序计算结果进行对比（任选取三个点，对比位移值）； 3. 提交程序编写过程的详细报告及计算机程序； 4. 所有同学参加答辩，并演示有限元计算程序。 有限元法中三节点三角形分析结构的步骤如下： 1）整理原始数据，如材料性质、荷载条件、约束条件等，离散结构并进行单元编码、结点编码、结点位移编码、选取坐标系。 2）单元分析，建立单元刚度矩阵。 3）整体分析，建立总刚矩阵。 4）建立整体结构的等效节点荷载和总荷载矩阵 5）边界条件处理。 6）解方程，求出节点位移。 7）求出各单元的单元应力。 8）计算结果整理。 一、程序设计 网格划分 如图，将薄板如图划分为6行，并建立坐标系，则
刚度矩阵的集成 建立与总刚度矩阵等维数的空矩阵，已变单元刚度矩阵的集成。 由单元分析已知节点、单元的排布规律，继而通过循环计算求得每个单元对应的节点序号。 通过循环逐个计算：（1）每个单元对应2种单元刚度矩阵中的哪一种； （2）该单元对应总刚度矩阵的那几行哪几列 （3）将该单元的单元刚度矩阵加入总刚度矩阵的对应行列 循环又分为3层循环：（1）最外层：逐行计算 （2）中间层：该行逐个计算 （3）最里层：区分为第 奇/偶 数个计算 单元刚度的集成：[ ][][][][][]' '''''215656665656266256561661e Z e e e Z e Z e e e e k k k K k k k k k k +?++=? =?==?==?=?????? 边界约束的处理：划0置1法 X Y P X Y P
一．问题描述 如图所示的平面矩形结构，设E=1，NU=0.25，h=1，考虑以下约束和外载： 位移边界条件BC（u）：U A=0，V A=0，U D=0， 力边界条件BC（p）：在CD边上有均布载荷q=1， 建模情形：使用四个四节点矩形单元， 试在该建模情形下，求各节点的位移以及各个单元的应力分布。
二．Matlab程序 (1).函数定义： function k= Quad2D4Node_Stiffness(E,NU,h,xi,yi,xj,yj,xm,ym,xp,yp,ID) syms s t; a = (yi*(s-1)+yj*(-1-s)+ym*(1+s)+yp*(1-s))/4; b = (yi*(t-1)+yj*(1-t)+ym*(1+t)+yp*(-1-t))/4; c = (xi*(t-1)+xj*(1-t)+xm*(1+t)+xp*(-1-t))/4; d = (xi*(s-1)+xj*(-1-s)+xm*(1+s)+xp*(1-s))/4; B1 = [a*(t-1)/4-b*(s-1)/4 0 ; 0 c*(s-1)/4-d*(t-1)/4 ; c*(s-1)/4-d*(t-1)/4 a*(t-1)/4-b*(s-1)/4]; B2 = [a*(1-t)/4-b*(-1-s)/4 0 ; 0 c*(-1-s)/4-d*(1-t)/4 ; c*(-1-s)/4-d*(1-t)/4 a*(1-t)/4-b*(-1-s)/4]; B3 = [a*(t+1)/4-b*(s+1)/4 0 ; 0 c*(s+1)/4-d*(t+1)/4 ; c*(s+1)/4-d*(t+1)/4 a*(t+1)/4-b*(s+1)/4]; B4 = [a*(-1-t)/4-b*(1-s)/4 0 ; 0 c*(1-s)/4-d*(-1-t)/4 ; c*(1-s)/4-d*(-1-t)/4 a*(-1-t)/4-b*(1-s)/4]; Bfirst = [B1 B2 B3 B4]; Jfirst = [0 1-t t-s s-1 ; t-1 0 s+1 -s-t ; s-t -s-1 0 t+1 ; 1-s s+t -t-1 0]; J = [xi xjxmxp]*Jfirst*[yi ;yj ; ym ; yp]/8; B = Bfirst/J; if ID == 1 D = (E/(1-NU*NU))*[1 NU 0 ; NU 1 0 ; 0 0 (1-NU)/2]; elseif ID == 2 D = (E/(1+NU)/(1-2*NU))*[1-NU NU 0 ; NU 1-NU 0 ; 0 0 (1-2*NU)/2]; end BD = J*transpose(B)*D*B; r = int(int(BD, t, -1, 1), s, -1, 1); z = h*r; k = double(z); end function z = Quad2D4Node_Assembly(KK,k,i,j,m,p) DOF(1)=2*i-1; DOF(2)=2*i; DOF(3)=2*j-1; DOF(4)=2*j; DOF(5)=2*m-1; DOF(6)=2*m; DOF(7)=2*p-1; DOF(8)=2*p; for n1=1:8
有限元分析大作业报告 试题1： 一、问题描述及数学建模 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： （1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； （2）分别采用不同数量的三节点常应变单元计算； （3）当选常应变三角单元时，分别采用不同划分方案计算。 该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况及方向如图所示。 二、采用相同单元数目的三节点常应变单元和六节点三角形单元计算 1、有限元建模 （1）设置计算类型：两者因几何条件和载荷条件均满足平面应变问题，故均取Preferences 为Structural （2）选择单元类型：三节点常应变单元选择的类型是Solid Quad 4 node182；六节点三角形单元选择的类型是Solid Quad 8 node183。因研究的问题为平面应变问题，故对Element behavior（K3）设置为plane strain。 （3）定义材料参数：弹性模量E=2.1e11，泊松比σ=0.3 （4）建几何模型：生成特征点；生成坝体截面 （5）网格化分：划分网格时，拾取lineAB和lineBC，设定input NDIV 为15；拾取lineAC，设定input NDIV 为20，选择网格划分方式为Tri+Mapped，最后得到600个单元。
（6)模型施加约束：约束采用的是对底面BC 全约束。大坝所受载荷形式为Pressure ，作用在AB 面上，分析时施加在L AB 上，方向水平向右，载荷大小沿L AB 由小到大均匀分布。以B 为坐标原点，BA 方向为纵轴y ，则沿着y 方向的受力大小可表示为： }{*980098000)10(Y y g gh P -=-==ρρ 2、 计算结果及结果分析 （1） 三节点常应变单元 三节点常应变单元的位移分布图 三节点常应变单元的应力分布图
结构有限元分析程序设计 绪论 §0.1 开设“有限元程序设计”课程的意义和目的 §0.2 课程特点 §0.3 课程安排 §0.4 课程要求 §0.5 基本方法复习 $0.1 意义和目的 1.有限元数值分析技术本身要求工程设计研究人员掌握 1). 有限元数值分析技术的完善标志着现代计算力学的真正成熟和实用化,已在各种 力学中得到了广泛的应用。比如：，已杨为工程结构分析中最得以收敛的技术手段，现代功用大致有： a). 现代结构论证。对结构设计从内力，位移等方面进行优劣评定，从而进 行结构优化设计。 b)可取代部份实验，局部实验+有限元分析，是现代工程设计研究方法的一大 特点。 c)结构的各种功能分析（疲劳断裂，可靠性分析等）都以有限元分析工具作为 核心的计算工具。 2). 有限元数值分析本身包括着理论+技术实现（本身功用所绝定的） 有限元数值分析本身包括着泛函理论+分片插值函数+程序设计 2. 有限元分析的技术实现（近十佘年的事）更依赖于计算机程序设计 有限元分析的技术取得的巨大的成就，从某种意义上说，得益于计算机硬件技术的发展和程序设计技术的发展，这两者的依赖性在当代表现得更加突出。（如可视化技术） 3.从学习的角度，不仅要学习理论，而且要从程序设计设计角度对这些理论的技术实现有 一个深入的了解，应当致力于掌握这些技术实现能力，从而开发它，发展它。（理论本身还有待于进一步完美相应的程序设计必须去开发） 4.程序设计不仅是实现有限元数值分析的工具和桥梁，而且在以下诸方面也有意义： 1). 精通基本概念，深化理论认识； 2). 锻炼实际工程分析，实际动手的能力； 3). 获得以后工作中必备的工具。（作业+老师给元素库） 目的：通过讲述有限元程序设计的技术与技巧，便能达到自编自读的能力。 §0.2 课程特点 总描述：理论+算法+数据结构（程序设计的意义） 理论：有限元算法，构造，步骤，解的等外性，收敛性，稳定性，误差分析 算法；指求解过程的技术方法，含两方面的含义；a. 有限元数值分析算法，b, 与数据结构有关的算法（总刚稀疏存贮，提取，节点优化编号等） 数据结构：指各向量矩阵存贮管理与实现，辅助管理结构（指针，数据记录等） 具体特点： 理论性强：能量泛函理论+有限元构造算法+数据结构构造算法 内容繁杂：理论方法+技术方法+技术技巧 技巧性强：排序，管理结构（指针生成，整型运算等）
重庆大学本科学生课程设计任务书 课程设计题目有限元程序设计 学院资源及环境科学学院专业工程力学年级2010级 已知参数和设计要求： 1.独立完成有限元程序设计。 2.独立选择计算算例，并能通过算例判断程序的正确性。 3.独立完成程序设计报告，报告内容包括理论公式、程序框图、程序本 体、计算算例，算例结果分析、结论等。 学生应完成的工作： 1.复习掌握有限单元法的基本原理。 2.掌握弹性力学平面问题3节点三角形单元或4节点等参单元有限元方法 的计算流程，以及单元刚度矩阵、等效节点载荷、节点应变、节点应力 和高斯积分等的计算公式。 3.用Fortran语言编写弹性力学平面问题3节点三角形单元或4节点等参 单元的有限元程序。 4.在Visual Fortran 程序集成开发环境中完成有限元程序的编辑和调试 工作。 5.利用编写的有限元程序，计算算例，分析计算结果。 6.撰写课程设计报告。 目前资料收集情况（含指定参考资料）： 1.王勖成，有限单元法，北京：高等教育出版社，2002。 2.O.C. Zienkiewicz, R. L. Taylor, Finite Element Method, 5th Eition, McGraw-Hall Book Company Limited, 2000。 3.张汝清，董明，结构计算程序设计，重庆：重庆大学出版社，1988。 课程设计的工作计划： 1.第1周星期一上午：教师讲解程序设计方法，程序设计要求和任务安 排。 2.第1周星期一至星期二完成程序框图设计。 3.第1周星期三至第2周星期四完成程序设计。 4.第2周星期五完成课程设计报告。 任务下达日期 2013 年 6 月 6 日完成日期 2013 年 07 月 03 日 指导教师（签名） 学生（签名）
姓名：学号：班级：
有限元分析及应用作业报告 一、问题描述 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； 2)分别采用不同数量的三节点常应变单元计算； 3)当选常应变三角单元时，分别采用不同划分方案计算。
二、几何建模与分析 图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）定义材料参数 4）生成几何模 a. 生成特征点 b.生成坝体截面 5）网格化分：划分网格时，拾取所有线段设定input NDIV 为10，选择网格划分方式为Tri+Mapped，最后得到200个单元。 6)模型施加约束： 约束采用的是对底面BC全约束。 大坝所受载荷形式为Pressure，作用在AB面上，分析时施加在L AB上，方向水平向右，载荷大小沿L AB由小到大均匀分布（见图1-2）。以B为坐标原点，BA方向为纵轴y，则沿着y方向的受力大小可表示为： ρ（1） = gh P- =ρ g = - 10 {* } 98000 98000 (Y ) y
湖南大学 《有限单元法》编程大作业 专业：土木工程 姓名： 学号： 2013年12月
目录 程序作业题目： (3) 1、程序编制总说明 (3) 2、Matlab程序编制流程图 (3) 3、程序主要标示符及变量说明 (4) 4、理论基础和求解过程 (5) 4.1、构造插值函数 (5) 4.2位移插值函数及应变应力求解 (5) 5.程序的验证 (6) 附录：程序代码 (15)
程序作业题目： 完成一个包含以下所列部分的完整的有限元程序( Project) 须提供如下内容的文字材料（1500字以上）： ①程序编制说明； ②方法的基本理论和基本公式； ③程序功能说明； ④程序所用主要标识符说明及主要流程框图； ⑤ 1～3 个考题：考题来源、输出结果、与他人成果的对比结果（误差百分比）； ⑥对程序的评价和结论（包括正确性、适用范围、优缺点及其他心得等）。 须提供源程序、可执行程序和算例的电子文档或文字材料。选题可根据各自的论文选题等决定。 1、程序编制总说明 a.该程序采用平面三角形等参单元，能解决弹性力学的平面应力、平面应变问题。 b.能计算单元受集中力的作用。 c.能计算结点的位移和单元应力。 d.考题计算结果与理论计算结果比较，并给出误差分析。 e.程序采用MATLAB R2008a编制而成。 2、Matlab程序编制流程图
图1 整个程序流程图 3、程序主要标示符及变量说明 1、变量说明： Node ------- 节点定义 gElement ---- 单元定义 gMaterial --- 材料定义，包括弹性模量，泊松比和厚度 gBC1 -------- 约束条件 gNF --------- 集中力 gk------------总刚 gDelta-------结点位移 输入结构控制参数 输入其它数据 形成整体刚度阵 引入支承条件 解方程，输出位移 求应力，输出应力 形成节点荷载向量 开始 结束 1 单元面积 求弹性矩阵 单元刚度矩阵 位移-应变矩阵 6 7 8 9 10 2 3 4 5
有限元大作业程序设计 学校：天津大学 院系：建筑工程与力学学院 专业：01级工程力学 姓名：刘秀 学号：\\\\\\\\\\\ 指导老师：
连续体平面问题的有限元程序分析 [题目]： 如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界 上受正向分布压力， m kN p 1=，同时在沿对角线y 轴上受一对集中压 力，载荷为2KN ，若取板厚1=t ，泊松比0=v 。 [分析过程]： 由于连续平板的对称性，只需要取其在第一象限的四分之一部分参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。
[程序原理及实现]： 用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN，经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。 该程序的特点如下： 问题类型：可用于计算弹性力学平面问题和平面应变问题 单元类型：采用常应变三角形单元 位移模式：用用线性位移模式 载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷 材料性质：弹性体由单一的均匀材料组成 约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独立约束 方程求解：针对半带宽刚度方程的Gauss消元法
输入文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT.IN 结果文件：输出一般的结果文件DATA.OUT 程序的原理如框图：
有限元分析习题及大作业试题 要求：1）个人按上机指南步骤至少选择习题中3个习题独立完成，并将计算结果上交； 2）以小组为单位完成有限元分析计算； 3）以小组为单位编写计算分析报告； 4）计算分析报告应包括以下部分： A、问题描述及数学建模； B、有限元建模（单元选择、结点布置及规模、网格划分方 案、载荷及边界条件处理、求解控制） C、计算结果及结果分析（位移分析、应力分析、正确性分 析评判） D、多方案计算比较（结点规模增减对精度的影响分析、单 元改变对精度的影响分析、不同网格划分方案对结果的 影响分析等） E、建议与体会 4）11月1日前必须完成，并递交计算分析报告（报告要求打印）。
习题及上机指南：（试题见上机指南） 例题1 坝体的有限元建模与受力分析 例题2 平板的有限元建模与变形分析 例题1：平板的有限元建模与变形分析 计算分析模型如图1-1 所示, 习题文件名: plane 0.5 m 0.5 m 0.5 m 0.5 m 板承受均布载荷：1.0e 5 P a 图1－1 受均布载荷作用的平板计算分析模型 1.1 进入ANSYS 程序 →ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: plane →Run 1.2设置计算类型 ANSYS Main Menu : Preferences →select Structural → OK 1.3选择单元类型 ANSYS Main Menu : Preprocessor →Element T ype →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK (back to Element T ypes window) → Options… →select K3: Plane stress w/thk →OK →Close (the Element T ype window) 1.4定义材料参数 ANSYS Main Menu : Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY :0.3 → OK 1.5定义实常数 ANSYS Main Menu: Preprocessor →Real Constant s… →Add … →select T ype 1→ OK →input THK:1 →OK →Close (the Real Constants Window)
有限元法分析与建模课程设计报告 学院：机械电子工程学院 专业：机械电子工程 指导教师：杜平安 学生：乔林 学号：201221080212 2012-12-10
摘要 摘要 连杆的作用是将活塞的往复运动变成曲轴的旋转运动, 并把活塞上的力传 给曲轴连杆工作的小端做往复运动, 大端作旋转运动, 杆身做复杂的平面运动。本文用Pro/E建立连杆的三维模型，并运用ANSYS强大的有限元分析和优化功能来实现连杆的分析ANSYS 是一款极其强大的有限元分析软件。通过数据接口，ANSYS 可以方便的实现从CAD 软件中导入实体模型。因此，将Pro/E强大的 建模功能与ANSYS 优越的有限元分析功能结合在一起可以极大地满足设计者 在设计过程中对建模与分析的需求。 关键词：连杆，有限元，Pro/E，ANSYS
ABSTRACT ABSTRACT The role oftheconnecting rodisthesmall end ofthereciprocation of the pistonintoarotational movementofthecrankshaft, and to transmittheforceon the pistontothecrankshaft connecting rodreciprocates, thebig endfor pivotal movement, Shaftdo complexplanar motion. The establishment ofalinkageof thethree-dimensionalmodelusingPro / E, thepowerfulANSYSfinite elementanalysis andoptimization capabilitiestoachievetheconnecting rodfatigueanalysisANSYSisan extremelypowerfulfinite element analysis software. Throughthedata interface, ANSYS canfacilitate the realization ofsolid modelsimportedfromCAD software. Therefore,thesuperiorpowerfulmodeling capabilitiesofPro / Eand ANSYSfinite elementanalysis capabilitiestogethercanmeetthedesignersin the design processmodelingand analysis. Keywords:rod, finite element, Pro / E, ANSYS
ansys有限元分析大作业
有限元大作业 设计题目: 单车的设计及ansys有限元分析 专业班级: 姓名: 学号: 指导老师: 完成日期: 2016.11.23
单车的设计及ansys模拟分析 一、单车实体设计与建模 1、总体设计 单车的总体设计三维图如下，采用pro-e进行实体建模。 在建模时修改proe默认单位为国际主单位(米千克秒 mks) Proe》文件》属性》修改
2、车架 车架是构成单车的基体，联接着单车的其余各个部件并承受骑者的体重及单车在行驶时经受各种震动和冲击力量，因此除了强度以外还应有足够的刚度，这是为了在各种行驶条件下，使固定在车架上的各机构的相对位置应保持不变，充分发挥各部位的功能。车架分为前部和后部，前部为转向部分，后部为驱动部分，由于受力较大，所有要对后半部分进行加固。
二、单车有限元模型 1、材料的选择 单车的车身选用铝合金（6061-T6）T6标志表示经过热处理、时效。 其属性如下： 弹性模量：） .6+ 90E （2 N/m 10 泊松比：0.33 质量密度：） 3 2.70E+ N/m （2 抗剪模量：） 60E .2+ N/m （2 10 屈服强度：） .2+ （2 75E 8 N/m 2、单车模型的简化 为了方便单车的模拟分析，提高电脑的运算
效率，可对单车进行初步的简化；单车受到的力的主要由车架承受，因此必须保证车架能够有足够的强度、刚度，抗振的能力，故分析的时候主要对车架进行分析。简化后的车架如下图所示。 3、单元体的选择 单车车架为实体故定义车架的单元类型为实体单元（solid）。查资料可以知道3D实体常用结构实体单元有下表。 单元名称说明 Solid45 三维结构实体单元，单元由8个节点定义，具有塑性、蠕变、应力刚化、 大变形、大应变功能，其高阶单元是 solid95
有限元分析软件 有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50 年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC 四个比较知名比较大的公司，其中ADINA、ABAQUS 在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件，ANSYS、MSC 进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析，但是除ADINA 以外其它三个必须与别的软件搭配进行迭代分析，唯一能做到真正流固耦合的软件只有ADINA。ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。 结构分析能力排名：ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名：ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名：ADINA、ANSYS、MSC、ABAQUS 性价比排名：最好的是ADINA，其次ABAQUS、再次ANSYS、最后MSC ABAQUS 软件与ANSYS 软件的对比分析： 1．在世界范围内的知名度：两种软件同为国际知名的有限元分析软件，在世界范围内具有各自广泛的用户群。ANSYS 软件在致力于线性分析的用户中具有很好的声誉，它在计算机资源的利用，用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题，其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。由于ANSYS 产品进入中国市场早于ABAQUS，并且在五年前ANSYS 的界面是当时最好的界面之一，所以在中国，ANSYS 软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北京办事处的成立，ABAQUS软件的用户数目和市场占有率正在大幅度和稳步提高，并可望在今后的几年内赶上和超过ANSYS。 2．应用领域：ANSYS 软件注重应用领域的拓展，目前已覆盖流体、电磁场和多物理场耦合等十分广泛的研究领域。ABAQUS 则集中于结构力学和相关领域研究，致力于解决该领域的深层次实际问题。 3．性价比：ANSYS 软件由于价格政策灵活，具有多种销售方案，在解决常规的
目录一. 前言 二．有限元设计部分 1 问题阐述 2 解析法求解 3 模型简化 4 ANSYS软件应用说明 5 结果分析 三．机械优化设计部分 1 问题阐述 2 解析算法 3 黄金分割法顺序流程图 4 C语言源程序代码 5 结果分析 四．设计心得 五. 参考文 一．前言
二．有限元设计部分 1、问题阐述 外伸梁上均布载荷的集中度为q=3kN/m，集中力偶矩M e=3kN·m列出剪力方程和弯矩方程，并绘制剪力图。材料力学Ι（刘鸿文第四版）P121
图2-1 外伸梁简化图 2、解析法求解 由梁的平衡方程，求出支反力为 F RA=14.5kN，F RB=3.5kN 梁的C A、ＡＤ、ＤＢ等三段内，剪力和弯矩都不能有同一个方程来表示，所以应分为三段考虑。对每一段都可以用同一个方法计算，列出剪力方程和弯矩方程，方程中ｘ以ｍ为单位，Ｆｓ（ｘ）以ｋＮ为单位，Ｍ（ｘ）以ｋＮ为单位。 在ＣＡ段内： Ｆｓ（ｘ）＝－ｑｘ＝－３ｘ（０＜=ｘ＜２ｍ） （g） Ｍ（ｘ）＝－（3／２）Ｘ２（０＜ｘ＜=2ｍ） （h）
在AD段内： Ｆｓ（ｘ）＝F RA－ｑｘ＝14.5－３ｘ（2m＜ｘ＜=6ｍ） （i） Ｍ（ｘ）＝F RA（x-2）-（1/2）X２=14.5（x-2）- （3/2）X 2 (j) (2m＜ ｘ6m) M（x）是x的二次函数,根据极值条件dM（x）/d(x)=0,得 14.5-3x=0 由此解出x=4.83m，亦即在那这一截面上，弯矩为极值。 代入（j） 式得AD段内的最大弯矩为 M=6.04kN·m 当截面取在DB段，用截面右侧的外力计算剪力和弯矩比较方
重庆大学研究生有限 元大作业
课程研究报告 科目：有限元分析技术教师：阎春平姓名：色学号： 2 专业：机械工程类别：学术 上课时间： 2015 年 11 月至 2016 年 1 月 考生成绩： 阅卷评语： 阅卷教师 (签名)
有限元分析技术作业 姓名: 色序号: 是学号: 2 一、题目描述及要求 钢结构的主梁为高160宽100厚14的方钢管，次梁为直径60厚10的圆钢管（单位为毫米），材料均为碳素结构钢Q235；该结构固定支撑点位于左右两端主梁和最中间。主梁和次梁之间是固接。试对在垂直于玻璃平面方向的2kPa 的面载荷（包括玻璃自重、钢结构自重、活载荷（人员与演出器械载荷）、风载荷等）作用下的舞台进行有限元分析。 二、题目分析 根据序号为069，换算得钢结构框架为11列13行。由于每个格子的大小为1×1（单位米），因此框架的外边框应为11000×13000（单位毫米）。 三、具体操作及分析求解 1、准备工作 执行Utility Menu：File → Clear＆start new 清除当前数据库并开始新的分析，更改文件名和文件标题，如图1.1。选择GUI filter，执行 Main Menu: Preferences → Structural → OK，如图1.2所示
图1.1清除当前数据库并开始新的分析 图1.2 设置GUI filter 2、选择单元类型。 执行Main Menu: Preprocessor →Element Type →Add/Edit/Delete →Add→ select→ BEAM188，如图2.1。之后点击OK(回到Element Types window) →Close
有限元分析软件的比较（购买必看）－转贴 随着现代科学技术的发展，人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能，需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响，看看是否会发生破坏性事故；分析计算核反应堆的温度场，确定传热和冷却系统是否合理；分析涡轮机叶片内的流体动力学参数，以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式，这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析（FEA，Finite Element A nalysis）方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中，有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃，主要表现在以下几个方面： 增加设计功能，减少设计成本； 缩短设计和分析的循环周期； 增加产品和工程的可靠性； 采用优化设计，降低材料的消耗或成本； 在产品制造或工程施工前预先发现潜在的问题； 模拟各种试验方案，减少试验时间和经费； 进行机械事故分析，查找事故原因。 在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分析计算，FEA在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在，世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PA FEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。 以下对一些常用的软件进行一些比较分析： 1. LSTC公司的LS-DYNA系列软件
有限元法分析与建模课程设计报告 学院：机电学院 专业：机械制造及其自动化指导教师：**** 学生：* *** 学号：2012011**** 2015-12-31
摘要 本文通过ANSYS10.0建立了标准光盘的离心力分析模型，采用有限元方法对高速旋转的光盘引起的应力及其应变进行分析，同时运用经典弹性力学知识来介绍ANSYS10.0中关于平面应力问题分析的基本过程和注意事项。力求较为真实地反映光盘在光驱中实际应力和应变分布情况，为人们进行合理的标准光盘结构设计和制造工艺提供理论依据。 关键词：ANSYS10.0；光盘；应力；应变。
目录 第一章引言 (3) 1.1 引言 (3) 第二章问题描述 (4) 2.1有限元法及其基本思想 (4) 2.2 问题描述 (4) 第三章力学模型的建立和求解 (5) 3.1设定分析作业名和标题 (5) 3.2定义单元类型 (6) 3.3定义实常数 (9) 3.4定义材料属性 (12) 3.5建立盘面模型 (14) 3.6对盘面划分网格 (22) 3.7施加位移边界 (27) 3.8施加转速惯性载荷并求解 (30) 第四章结果分析 (32) 4.1 旋转结果坐标系 (32) 4.2查看变形 (33) 4.3查看应力 (35) 总结 (38) 参考文献 (39)
第一章引言 1.1 引言 光盘业是我国信息化建设中发展迅速的产业之一，认真研究光盘产业的规律和发展趋势，是一件非常迫切的工作。光盘产业发展的整体性强，宏观调控要求高，因此，对于光盘产业的总体部署、合理布局和有序发展等问题，包括节目制作、软件开发、硬件制造、节目生产、技术标准等。 在高速光盘驱动器中，光盘片会产生应力和应变，在用ANSYS分析时，要施加盘片高速旋转引起的惯性载荷，即可以施加角速度。需要注意的是，利用ANSYS施加边界条件时，要将内孔边缘节点的周向位移固定，为施加周向位移，而且还需要将节点坐标系旋转到柱坐标系下。 本文通过ANSYS10.0建立了标准光盘的离心力分析模型，采用有限元方法对高速旋转的光盘引起的应力及其应变进行分析，同时运用经典弹性力学知识来介绍ANSYS10.0中关于平面应力问题分析的基本过程和注意事项。
有限元分析及应用大作业 作业要求： 1）个人按上机指南步骤至少选择习题中3个习题独立完成，并将计算结果上交； 也可根据自己科研工作给出计算实例。 2）以小组为单位完成有限元分析计算； 3）以小组为单位编写计算分析报告； 4）计算分析报告应包括以下部分： A、问题描述及数学建模； B、有限元建模（单元选择、结点布置及规模、网格划分方案、载荷及边界 条件处理、求解控制） C、计算结果及结果分析（位移分析、应力分析、正确性分析评判） D、多方案计算比较（结点规模增减对精度的影响分析、单元改变对精度的 影响分析、不同网格划分方案对结果的影响分析等） 题一：图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；（注意ANSYS中用四边形单元退化为三节点三角形单元） 2）分别采用不同数量的三节点常应变单元计算； 3）当选常应变三角单元时，分别采用不同划分方案计算。 解：1.建模： 由于大坝长度>>横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看作无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作
用于大坝上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力，满足平面应变问题的载荷条件。因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况P=98000-9800*Y；建立几何模型，进行求解；假设大坝的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3； 2：有限元建模过程： 2.1 进入ANSYS : 程序→ANSYS APDL 15.0 2.2设置计算类型: ANSYS Main Menu: Preferences →select Structural →OK 2.3选择单元类型: ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182，六节点三角形单元选择Solid Quad 8node 183)→OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window) 2.4定义材料参数： ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK 2.5生成几何模型： 生成特征点： ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS →依次输入四个点的坐标：input:1(0,0),2(10,0),3(1,5),4(0.45,5) →OK 生成坝体截面： ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点，1(0,0),2(6,0),3(0,10) →OK 2.6 网格划分： ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window) 2.7 模型施加约束： 给底边施加x和y方向的约束： ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK 给竖直边施加y方向的分布载荷： ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量；2) 在Result窗口中出现{X},写入所施加的载荷函数： 98000-9800*{Y}；3) File>Save(文件扩展名：func) →返回：Parameters →Functions →Read from file：将需要的.func文件打开，参数名取meng，它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边；OK →在下拉列表框中，选择：Existing table →OK →选择需要的载荷为meng参数名→OK 2.8 分析计算： ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load
基于ANSYS软件的有限元分析报告 机制1205班杜星宇U201210671 一、概述 本次大作业主要利用ANSYS软件对桌子的应力和应变进行分析,计算出桌子的最大应力和应变。然后与实际情况进行比较,证明分析的正确性,从而为桌子的优化分析提供了充分的理论依据，并且通过对ANSYS软件的实际操作深刻体会有限元分析方法的基本思想，对有限元分析方法的实际应用有一个大致的认识。 二、问题分析 已知：桌子几何尺寸如图所示，单位为mm。假设桌子的四只脚同地面完全固定，桌子上存放物品，物品产生的均匀分布压力作用在桌面，压力大小等于300Pa，其中弹性模量E=9.3GPa，泊松比μ=0.35，密度ρ=560kg/m3，分析桌子的变形和应力。
将桌脚固定在地面，然后在桌面施加均匀分布的压力，可以看作对进行平面应力分析，桌脚类似于梁单元。由于所分析的结构比较规整且为实体，所以可以将单元类型设为八节点六面体单元。 操作步骤如下： 1、定义工作文件名和工作标题 （1）定义工作文件名：执行Utility Menu/ File/Change Jobname，在弹出Change Jobname 对话框修改文件名为Table。选择New log and error files复选框。 （2）定义工作标题：Utility Menu/File/ Change Title，将弹出Change Title对话框修改工作标题名为The analysis of table。 （3）点击：Plot/Replot。 2、设置计算类型 （1）点击：Main Menu/Preferences,选择Structural,点击OK。