有限元及程序设计
1.下列关于高精度单元描述正确的是（）。
A.等参元的位移模式和坐标变换采用不同的形函数
B.矩形单元形状规则，因而使用范围较广
C.6结点三角形单元、10结点三角形单元、8结点矩形单元和12结点矩形单元的单元刚度矩阵的建立过程是不一样的
D.6结点三角形单元较容易模拟物体的边界形状
【参考答案】：D
2.φ=cxy能解决矩形板（）问题。
A.左右均布拉压
B.上下均布拉压
C.纯剪切
D.纯弯曲
【参考答案】：C
3.下列关于等参元的叙述不正确的是（）。
A.精度较高
B.能较好的模拟边界条件
C.输入的信息量较少
D.输入的信息量较多
【参考答案】：D
4.薄板的边界不包括（）。
A.简支边界
B.固定边界
C.自由边界和荷载边界
D.非固定边界
【参考答案】：D
5.下列属于平面应力问题的是（）。
A.平板坝的平板支墩
B.挡土墙
C.重力水坝
D.受内水压力作用的圆管
【参考答案】：A
6.在应力函数上任意增减一个（），对应力分量无影响。
A.线性项
B.二次项
C.三次项
D.常数项
【参考答案】：A
7.下列不属于提高单元精度的方法是（）。
A.增加单元结点数目
B.在单元内增设结点
C.减少单元结点数目
D.设等参元
【参考答案】：C
8.空间问题的基本平衡微分方程有（）个。
A.2
B.3
C.4
D.5
【参考答案】：C
9.φ=by2能解决矩形板（）问题。
A.左右均布拉压
B.上下均布拉压
C.纯剪切
D.纯弯曲
【参考答案】：B
10.下列属于不规则单元的有（）。
A.正四面体单元
B.正三棱体单元
C.任意四面体单元
D.正六面体单元
【参考答案】：C
11.空间问题的基本未知位移分量有（）个。
A.2
B.3
C.4
D.5
【参考答案】：B
1.薄板小挠度弯曲理论的基本假定是（）。
A.直法线假定
B.法向位移假定
C.中面位移假定
D.板内无挤压假定
【参考答案】：A|C|D
2.弹性力学平面问题按应力求解具体可分为（）两种。
A.逆解法
B.顺解法
C.半逆解法
D.半顺解法
【参考答案】：A|C
3.弹性力学的边界条件有（）。
A.位移边界条件
B.应力边界条件
C.应变边界条件
D.混合边界条件
【参考答案】：A|B|D
4.弹性力学的基本方程包括（）。
A.平衡方程
B.应力方程
C.几何方程
D.物理方程
【参考答案】：A|C|D
5.有限元收敛性准则是（）。
A.完备性准则
B.对称性准则
C.平衡性准则
D.协调性准则
【参考答案】：A|D
6.总体刚度矩阵具有（）性质。
A.对称性
B.稀疏性
C.带状分布性
D.奇异性
【参考答案】：A|B|C|D
7.边界条件的处理方法有（）。
A.划0置1法
B.置大数法
C.划1置0法
D.置小数法
【参考答案】：A|B
8.弹性力学平面问题包括（）两种。
A.平面物理问题
B.平面几何问题
C.平面应力问题
D.平面应变问题
【参考答案】：C|D
9.用有限元法处理实际问题主要分为（）三部分。
A.离散化
B.单元分析
C.局部分析
D.整体分析
【参考答案】：A|B|D
10.弹性力学问题的求解方法有（）。
A.按应变求解
B.按应力求解
C.按体力求解
D.按位移求解
【参考答案】：B|D
11.弹性力学的基本假定有（）。
A.假设物体是连续的
B.假设物体是均匀的和各向同性的
C.假设物体是完全弹性的
D.假设物体内无初应力
E.假设物体的变形是很小的
【参考答案】：A|B|C|D|E
12.离散化过程应遵循（）的原则。
A.位移近似
B.几何近似
C.物理近似
D.形状近似
【参考答案】：B|C
13.下面说法错误的是（）。
A.应力以正面正向，负面负向为正
B.正应变以伸长为正
C.剪应变以夹角小为正
D.体力、面应、位移以沿坐标轴正向为正
【参考答案】：A|B|C|D
14.总体刚度矩阵的形成方法有（）。
A.结点平衡法
B.直接刚度法
C.间接刚度法
D.位移法
【参考答案】：A|B
15.下列属于高精度空间单元的有（）。
A.10结点30自由度四面体单元
B.20结点60自由度六面体单元
C.6结点三角形单元
D.4结点48自由度四面体单元
【参考答案】：A|B|D
16.边界内应力的处理最常用的方法有（）。
A.结点法
B.绕结点平均法
C.二单元平均法
D.应力平均法
【参考答案】：B|C
1.如果均采用三角形单元，则平面应力问题和平面应变问题单元刚度矩阵相同。
N.错误Y.正确
【参考答案】：N
2.解平面应力和平面应变问题采用的应力矩阵相同。
N.错误Y.正确
【参考答案】：Y
3.解平面应力和平面应变采用的单元刚度矩阵不相同。
N.错误Y.正确
【参考答案】：Y
4.在体力为常量时，平面问题应力分布与材料的弹性常数无关。
N.错误Y.正确
【参考答案】：Y
5.位移模式的选取必须遵循几何近似和物理近似两条原则。
N.错误Y.正确
【参考答案】：N
6.有限元离散结构时，在应力变化较大的区域应采用较密网格。
N.错误Y.正确
【参考答案】：Y
7.重力水坝属于弹性力学平面应力问题。
N.错误Y.正确
【参考答案】：N
8.按位移法求解弹性力学问题时需满足平衡方程和应变协调方程。
N.错误Y.正确
【参考答案】：N
9.应力函数线性项可解决矩形板拉伸问题。
N.错误Y.正确
【参考答案】：N
10.在体力为常量时，平衡方程、相容方程及应力边界条件中均不含弹性常数E和μ，故我们可以由一种材料替代另一种材料，用平面应力问题替代平面应变问题作实验，得到的应力是完全一样的。
N.错误Y.正确
【参考答案】：Y
11.剪力墙墙体受自重作用属于平面应变问题，天然气管道属于平面应力问题。
N.错误Y.正确
【参考答案】：N
12.有限元解薄板问题的基本单元有矩形薄板单元和三角形薄板单元，其中矩形薄板单元是完全协调单元。
N.错误Y.正确
【参考答案】：N
13.三角形单元是常应力单元，但不是常应变单元。
N.错误Y.正确
【参考答案】：N
14.在应力函数中任意减少一个线性项对应力无影响。
N.错误Y.正确
【参考答案】：Y
15.受自重作用的门属于平面应力问题，挡土墙受土压力属于平面应变问题。N.错误Y.正确
【参考答案】：Y
16.应力函数φ=ax2能解决矩形板纯弯曲问题。
N.错误Y.正确
【参考答案】：N
17.矩形单元的公式推导简单，但精度比三结点三角形单元低。
N.错误Y.正确
【参考答案】：N
18.挡土墙受土压力属于弹性力学平面应变问题。
N.错误Y.正确
【参考答案】：Y
19.按应力求解要同时满足平衡方程、相容方程和应力边界条件。
N.错误Y.正确
【参考答案】：Y
20.有限元的单元局部编码，应按逆时针规律编排。
N.错误Y.正确
【参考答案】：Y
21.三结点三角形单元是常应力单元，但不是常应变单元。
N.错误Y.正确
【参考答案】：N
22.直角坐标与极坐标系下基本方程完全类似，相应未知量也相同。
N.错误Y.正确
【参考答案】：N
23.应力函数φ=ay3能解决矩形板受均布拉（压）力的问题。
N.错误Y.正确
【参考答案】：N
. 一、题目 如图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
受均匀内压作用的厚壁圆筒： 问题描述： 受均匀内压p=12.5N/mm 2作用的厚壁圆筒。其几何参数为：内径R i =100mm ， 外径R e =200mm ，桶壁后h=100mm ，材料参数为：E=8666.67Mpa ，v=0.3， s σ=17.32Mpa ，材料符合Mise 屈服条件。 (a)求理想塑性材料的解，给出应力r σ和θσ沿径向r 的分布曲线，并求完全卸载 后圆筒内的残余应力分布。 (b)求线性强化材料(E 1=0.6E 或E 1=0.6E)的解，即应力r σ和θσ沿径向r 的分布 曲线。 (c)求幂硬化材料的解并绘出当弹塑性比例系数为m=0,1/4,1/2,2/3和m=1.0时， 即应力r σ和θσ沿径向r 的分布曲线。 求解分析： 由于该厚壁筒模型是轴对称模型，所以在求解过程中，我们选取了1/4模型进行了进行建模分析，具体如下图： 建模时取了柱坐标系下厚壁筒从0。~90。范围内的部分，高度取为100mm ，模型完成后进行网格的划分，这里利用了Patran 的Mesh Seed 功能，通过在径向、周向，高度方向撒种生成Mesh 网格，网格划分如上图。 考虑到实体的变形情况，关于模型的边界条件，定义如下： (1)模型的上、下表面为两个平面，在该两平面上限制z 方向的位移为0； (2)对于模型的内外两圆弧面，为了方便定义边界条件，建立了柱坐标，该两平
是延径向变形的，所以ρ坐标是放开的，为了限制模型的刚体移动，这里限制角坐标θ为0。 (3)对于模型两个侧平面，是属于模型的对称面，所以该两平面的单元在垂直于平面的方向上位移为零，这里利用柱坐标，即沿周向的位移为零，所以同样要限制角坐标θ为0。 由于厚壁筒受到均匀内压，所以在施加载荷时选择均布载荷Pressure，大小为p=12.5N/mm2，作用在内圆弧表面上。 对于材料塑性的定义，首先定义样式模量和泊松比，然后在弹塑性对话框里定义屈服载荷和硬化系数或通过在Stress/Strain Curve栏中添加事先定义的材料属性场来表征弹塑性比例系数m。 对于求解分析，求解器选择Nastran进行计算分析，单元属性选择3D Solid 属性，分析类型定义为非线性并设置大变形和跟随力及载荷增量步等，以此来进行弹塑性的非线性求解。 结果分析： (a)对于理性塑性材料，即硬化系数为0，求解结果如下： 该图为100%载荷作用下模型的应力云图及变形情况。观察可知，筒内壁应力较高且首先达到屈服应力发生塑性变形，沿径向方向向外，各层应力逐渐递减，且外层部分属于弹性变形的范畴，模型某一层为弹塑性变形的分界面。
有限元分析大作业报告 试题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, 与数据结构有关的算法（总刚稀疏存贮，提取，节点优化编号等） 数据结构：指各向量矩阵存贮管理与实现，辅助管理结构（指针，数据记录等） 具体特点： 理论性强：能量泛函理论+有限元构造算法+数据结构构造算法 内容繁杂：理论方法+技术方法+技术技巧 技巧性强：排序，管理结构（指针生成，整型运算等）
姓名：学号：班级：
有限元分析及应用作业报告 一、问题描述 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 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
有限元分析》大作业基本要求： 1.以小组为单位完成有限元分析计算，并将计算结果上交； 2.以小组为单位撰写计算分析报告； 3.按下列模板格式完成分析报告； 4.计算结果要求提交电子版，一个算例对应一个文件夹，报告要求提交电子版和纸质版。 有限元分析》大作业 小组成 员： 储成峰李凡张晓东朱臻极高彬月 Job name ：banshou 完成日 期： 2016-11-22 一、问题描述 （要求：应结合图对问题进行详细描述，同时应清楚阐述所研究问题的受力状况 和约束情况。图应清楚、明晰，且有必要的尺寸数据。）如图所示，为一内六角螺栓扳手，其轴线形状和尺寸如图，横截面为一外 接圆半径为0.01m的正六边形，拧紧力F为600N，计算扳手拧紧时的应力分布 图1 扳手的几何结构 数学模型
要求：针对问题描述给出相应的数学模型，应包含示意图，示意图中应有必要的尺寸数据；
图 2 数学模型 如图二所示，扳手结构简单，直接按其结构进行有限元分析。 三、有限元建模 3.1 单元选择 要求：给出单元类型， 并结合图对单元类型进行必要阐述， 包括节点、自由度、 实常数等。） 图 3 单元类型 如进行了简化等处理，此处还应给出文字说
扳手截面为六边形，采用4 节点182单元，182 单元可用来对固体结构进行
二维建模。182单元可以当作一个平面单元，或者一个轴对称单元。它由4 个结点组成，每个结点有2 个自由度，分别在x,y 方向。 扳手为规则三维实体，选择8 节点185单元，它由8 个节点组成，每个节点有3 个自由度，分别在x，y，z 方向。 3.2 实常数 （要求：给出实常数的具体数值，如无需定义实常数，需明确指出对于本问题选择的单元类型，无需定义实常数。） 因为该单元类型无实常数，所以无需定义实常数 3.3材料模型 （要求：指出选择的材料模型，包括必要的参数数据。） 对于三维结构静力学，应力主要满足广义虎克定律，因此对应ANSYS中的线性，弹性，各项同性，弹性模量EX：2e11 Pa, 泊松比PRXY=0.3 3.4几何建模由于扳手结构比较简单，所以可以直接在ANSYS软件上直接建模，在ANSYS建 立正六 边形，再创立直线，面沿线挤出体，得到扳手几何模型 图4 几何建模
有限元分析习题及大作业试题 要求：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)
有限元大作业程序设计 学校：天津大学 院系：建筑工程与力学学院 专业：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 程序的原理如框图：
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 软件由于价格政策灵活，具有多种销售方案，在解决常规的
重庆大学本科学生课程设计任务书 课程设计题目有限元程序设计 学院资源及环境科学学院专业工程力学年级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 日 指导教师（签名） 学生（签名）
ansys有限元分析工程实例大作业
————————————————————————————————作者：————————————————————————————————日期：
辽宁工程技术大学 有限元软件工程实例分析 题目基于ANSYS钢桁架桥的静力分析专业班级建工研16-1班（结构工程）学号 471620445 姓名 日期 2017年4月15日
基于ANSYS钢桁架桥的静力分析 摘要：本文采用ANSYS分析程序，对下承式钢桁架桥进行了有限元建模;对桁架桥进行了静力分析，作出了桁架桥在静载下的结构变形图、位移云图、以及各个节点处的结构内力图（轴力图、弯矩图、剪切力图），找出了结构的危险截面。 关键词：ANSYS；钢桁架桥；静力分析；结构分析。 引言：随着现代交通运输的快速发展，桥梁兴建的规模在不断的扩大，尤其是现代铁路行业的快速发展更加促进了铁路桥梁的建设，一些新建的高速铁路桥梁可以达到四线甚至是六线，由于桥面和桥身的材料不同导致其受力情况变得复杂，这就需要桥梁需要有足够的承载力，足够的竖向侧向和扭转刚度，同时还应具有良好的稳定性以及较高的减震降噪性，因此对其应用计算机和求解软件快速进行力学分析了解其受力特性具有重要的意义。 1、工程简介 某一下承式简支钢桁架桥由型钢组成，顶梁及侧梁，桥身弦杆，底梁分别采用3种不同型号的型钢，结构参数见表1，材料属性见表2。桥长32米，桥高5.5米，桥身由8段桁架组成，每个节段4米。该桥梁可以通行卡车，若只考虑卡车位于桥梁中间位置，假设卡车的质量为4000kg，若取一半的模型，可以将卡车对桥梁的作用力简化为P1，P2，和P3，其中P1=P3=5000N，P2=10000N,见图2,钢桥的形式见图1，其结构简图见图3。
重庆大学研究生有限 元大作业
课程研究报告 科目：有限元分析技术教师：阎春平姓名：色学号： 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系列软件
有 限 元 大 作 业 程 序 设 计 学校：天津大学 院系：建筑工程与力学学院 专业：01级工程力学 姓名：刘秀 学号：\\\\\\\\\\\ 指导老师： 连续体平面问题的有限元程序分析 [题目]： 如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界 上受正向分布压力， m kN p 1=，同时在沿对角线y 轴上受一对集中压 力，载荷为2KN ，若取板厚1=t ，泊松比0=v 。 [分析过程]： 由于连续平板的对称性， 只需要取其在第一象限的四分之一部分
参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。 [ 用和单元信息文件DATA.OUT。 位移模式：用用线性位移模式 载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷 材料性质：弹性体由单一的均匀材料组成 约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独立约束 方程求解：针对半带宽刚度方程的Gauss消元法 输入文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT.IN 结果文件：输出一般的结果文件DATA.OUT 程序的原理如框图：
（1） ID ： ID=2时为平面应变问题 （平面问题） ，LJK_ELE(I,1)，LJK_ELE(I,2)， X(I),Y(I)分别存放节点I 的x ，y 表示第I 个作用有节点载荷的节点x,y 方向的节点载荷数值 存放节点载荷向量，解方程后该矩 （2 READ_IN ： 读入数据 BAND_K ： 形成半带宽的整体刚度矩阵 FORM_KE ： 计算单元刚度矩阵 FORM_P ： 计算节点载荷 CAL_AREA ：计算单元面积 DO_BC ： 处理边界条件 CLA_DD ： 计算单元弹性矩阵 SOLVE ： 计算节点位移 CLA_BB ： 计算单元位移……应变关系矩阵 CAL_STS ：计算单元和节点应力 （3）文件管理： 源程序文件： chengxu.for 程序需读入的数据文件：
有限元分析及应用大作业 作业要求： 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。
有限元编程的c++实现算例 1. #include<> 2. #include<> 3. 4. 5. #define ne 3 #define nj 4 #define nz 6 #define npj 0 #define npf 1 #define nj3 1 2 #define dd 6 #define e0 #define a0 #define i0 #define pi 16. 17. 18. int jm[ne+1][3]={{0,0,0},{0,1,2},{0,2,3},{0,4,3}}; /*gghjghg*/ 19. double gc[ne+1]={,,,}; 20. double gj[ne+1]={,,,}; 21. double mj[ne+1]={,a0,a0,a0}; 22. double gx[ne+1]={,i0,i0,i0}; 23. int zc[nz+1]={0,1,2,3,10,11,12}; 24. double pj[npj+1][3]={{,,}}; 25. double pf[npf+1][5]={{0,0,0,0,0},{0,-20,,,}}; 26. double kz[nj3+1][dd+1],p[nj3+1]; 27. double pe[7],f[7],f0[7],t[7][7]; 28. double ke[7][7],kd[7][7]; 29. 30. 31. 36. void jdugd(int); 38. void zb(int); 39. void gdnl(int); 40. void dugd(int);