有限元分析及应用大作业

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

目录

有限元分析及应用作业报告....................................... 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→依次输入三个点的坐标：

input:1(0,0),2(6,0),3(0,10)→OK

b.生成坝体截面：ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接三个特征点，

1(0,0),2(6,0),3(0,10) →OK

5）网格化分：划分网格时，拾取lineAB和lineBC进行Size Conrotls，设定input NDIV 为15；拾取lineAC，设定input NDIV 为20，选择网格划分方式为

Tri+Mapped，最后得到600个单元。

6)模型施加约束：

约束采用的是对底面BC全约束。

大坝所受载荷形式为Pressure，作用在AB面上，分析时施加在L AB上，方向水平向右，载荷大小沿L AB由小到大均匀分布（见图1-2）。以B为坐标原点，BA方向为纵轴y，则沿着y方向的受力大小可表示为：

ρ

P-

=

g

=

=ρ

-

gh

y

{*

}

98000

98000

)

(Y

10

（1）

其中ρ为水的密度，取g为9.8m/s2，可知P max为98000N，P min为0。施加载荷时只需对L AB插入预先设置的载荷函数（1）即可。

网格划分及约束受载情况如图1-3(a)和1-4(a)所示。

7）分析计算

ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load Step window) →OK

8）结果显示

ANSYS Main Menu: General Postproc →Plot Results →Deformed Shape… →select Def + Undeformed →OK (back to Plot Results window)→Contour Plot →Nodal Solu… →select: DOF solution, UX,UY, Def + Undeformed , Stress ,SX,SY,SZ, Def + Undeformed→OK

四、计算结果及结果分析

4.1计算结果

（1）三节点常应变单元

图1-3(a)常应变三节点单元的网格划分及约束受载图

图1-3(b)常应变三节点单元的位移分布图

图1-3(c) 常应变三节点单元的应力分布图（2）六节点三角形单元

图1-4(a)六节点三角形单元网格划分及约束受载图

图1-4(b) 六节点三角形单元的变形分布图

图1-4(c)六节点三角形单元的应力分布图

根据以上位移和应力图，可以得出常应变三节点单元和六节点三角形单元的最小最大位移应力如表1-1所示。

4.2 结果分析

由以上各图和数据表可知，采用三节点和六节点的三角形单元分析计算：（1）最大位移都发生在A点，即大坝顶端，最大应力发生在B点附近，即坝底和水的交界处，且整体应力和位移变化分布趋势相似，符合实际情况；

（2）结果显示三节点和六节点单元分析出来的最大应力值相差较大，原因可能是B点产生了虚假应力，造成了最大应力值的不准确性。

（3）根据结果显示，最小三节点和六节点单元分析出来的最小应力值相差极为悬殊，结合理论分析，实际上A点不承受载荷，最小应力接近于零，显然六节点三角形单元分析在这一点上更准确。

（4）六节点的应力范围较大，所以可判断在单元数目相同的前提下，节点数目越多，分析精度就越大；但是节点数目的增多会不可避免地带来计算工作量增加和计算效率降低的问题。

五、第2问的有限元建模及计算结果

此次分析选择的单元类型为常应变三节点单元。选用三种不同单元数目情况进行比较分析。具体做法如下：有限元建模步骤与第1小题类似，只是在划分网格时，对L BC的AERA ：Size Controls不同：依次设置单元边长度Element edge length为1.6、0.4、0.1，所获得的单元数目依次为96（图1-9(a)）、1536（图1-10(a)）、20880（图1-11(a)）；分别计算并得到位移变化图如图1-9(b)、1-10(b)、1-11(c)所示；分别计算并得到应力变化云图如图1-9(c)、1-10(c)、1-11(c)所示。

（1）单元数目为96的常应变三节点单元

图1-9(a)单元数目为96的网格划分及约束受载图

图1-9(b)单元数目为96的位移分布图

图1-9(c)单元数目为96的应力分布图

（2）单元数目为1536的常应变三节点单元

图1-10(a)单元数目为1536的网格划分及约束受载图

图1-10(b)单元数目为1536的位移分布图

图1-10(c)单元数目为1536的应力分布图

（3）单元数目为20880的常应变三节点单元

图1-11(a)单元数目为20880的网格划分及约束受载图

图1-11(b) 单元数目为20880的位移分布图

图1-11(c)单元数目为20880的应力分布图

由以上不同单元数目的位移应力分布图可以看出，大坝截面所受位移和应力的变化趋势是相同的，最大应力都发生在坝底和水的交界点附近，最小应力发生在大坝顶端；最大变形位移也是发生在坝顶。不同单元数目下计算的数据如表1-2所示。

序号单元数最大位移（mm）最小应力（Pa）最大应力（Pa）

1 96 0.026

2 13656 275480

2 1536 0.0289 2730 499624

3 20880 0.0292 925 716492 （4）结果分析

由以上分析结果可知：

（1）随着单元数目的增加，最大位移变化不大，应力变化范围逐步增大；

（2）随着单元数目的增加，即网格划分越密，分析的结果准确度将会提高；但是单元数目的增加和节点数目的增加都会造成计算量的增加和计算速度的下降的问题。

（3）对于本次计算结果，仍可能存在虚假应力，应力的准确值无法准确得出，只是网格划分越密，计算结果越精确。所以减少虚假应力影响的措施之一就是增加单元的数目，提高网格划分的密度。

五、第3问的有限元建模及计算结果

由图1-1所示的划分方案可知，需采用手动划分网格：首先创建6个节点，然后采用不同的方式连接节点创建单元，从而分别得到两种不同的网格划分方式，见下图1-12所示。

对底边的三个节点施加全约束；载荷建立方程式并创建table；其他的处理方式与第1小题相同。

图1-12方案一和二的划分方案图

有限元模型建立完成后进行求解，则可得到方案一和方案二的的位移图和应力图，如图1-13(a)、1-13(b)、1-14（a）、1-14（b）所示。

图1-13（a）方案一网格划分方式下的位移图

图1-13（b）方案一网格划分方式下的应力图

图1-14（a）方案二网格划分方式下的位移图

图1-14（b）方案二网格划分方式下的应力图

由以上两种方案的位移和应力图可得出的最大位移和最小最大应力如表1-3所示：

表1-3 方案一和方案二计算数据表

由以上分析结果可知，由于方案一和二都只有四个单元，所以在计算应力和位移的时结果的准确度较低。分析应力图可知，方案二得出的最大应力不在坝底和水的交界处，不符合实际情况，而方案一的最大应力所在位置符合实际情况，所以总体来说，方案一的分析结果优于方案二。原因是方案一具有整体几何保形性的单元数目多于方案二的数目。

六、总结和建议

通过以上分析情况可以看出，如果要使分析结果较为精确，首先是单元的类型选择要恰当，由第1小题计算结果可知，不同的单元类型会造成结果的不同，节点较多可以保证计算精度较高；由第二小题的计算结果可知，划分网格时，单元数目也不能太少，单元数目的增加也可以提高计算的精度；但是对于实际工程而言，采用较多节点的单元反而会影响计算的工作量，这是不经济不必要的。因此在保证网格划分大小适当和均匀的前提下，使应力集中处划的密集些，这样也能得到较为精确的结果。

试题4

一、问题的描述

图示为带方孔（边长为80mm）的悬臂梁，其上受部分均布载荷

（p=10Kn/m）作用，试采用一种平面单元，对图示两种结构进行有限元分析，并就方孔的布置（即方位）进行分析比较，如将方孔设计为圆孔，结果有何变化？（板厚为1mm，材料为钢）。

图2-1

二、几何建模与分析

由图2-1及问题描述可知，板的长宽尺寸远远大于厚度，研究结构为一很薄的等厚度薄板，满足平面应力的几何条件；作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面上无外力作用，满足平面应力的载荷条件。故该问题属于平面应力问题，薄板所受的载荷为面载荷，分布情况及方向如图2-1所示，建立几何模型，进行求解。

薄板的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3

三、有限元建模及计算结果分析

选取三节点常应变单元plane42，来计算分析薄板的位移和应力。由于此问题为平面应力问题，所以分析时对每个单元类型的Element behavior（K3）都设置为plane str w/thk。

定义材料参数：ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK

定义实常数：ANSYS Main Menu: Preprocessor →Real Constant s… →Add… →select Type 1→OK→input THK:1→OK →Close (the Real Constants Window)

2个单元的分析按照

1）生成特征点

ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：

三节点常应变单元：input:1(0,0),2(0,500),3(450,500),4(900,500),5(900,250),

6(300,113.4),7(243.4,250),8(300,306.6),9(356.6,250)→OK

2）生成平板

ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →连接特征点→生成两个area→Operate→Subtract→拾取梯形A1→OK→拾取矩形A2→OK→生成平板

3）网格划分

ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →Mesh: Areas, Shape: Tri, Free →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)

4)模型施加约束

给模型施加x，y方向约束

有限元分析大作业报告

有限元分析大作业报告 试题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） 三节点常应变单元 三节点常应变单元的位移分布图 三节点常应变单元的应力分布图

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

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

目录 有限元分析及应用作业报告....................................... 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→依次输入三个点的坐标：

华科大有限元分析题及大作业题答案——船海专业(DOC)

姓名：学号：班级：

有限元分析及应用作业报告 一、问题描述 图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 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

有限元分析大作业试题

有限元分析习题及大作业试题 要求：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)

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

有限元分析报告大作业

有限元分析》大作业基本要求： 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 几何建模

重庆大学研究生有限元大作业教学内容

重庆大学研究生有限 元大作业

课程研究报告 科目：有限元分析技术教师：阎春平姓名：色学号： 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

有限元分析作业报告 英文版

有限元分析及应用上机实验报告 学院：机电工程学院 专业：机械工程 班级：硕士1606班 姓名：钱树生 学号：163712160 指导教师：李毅波 日期：2016.12.02

1.Question Fig.1.Schematic diagram of herringbone roof truss. Question:The geometric dimensions of the chevron roof is shown in Fig.1,you should analyze it by statics，as a result you should give the displacement and axial force and axial force diagram of the deformation diagram. Conditions:The ends of the roof truss were fixed,the sectional area of the truss is 0.01m2,elastic modulus is2.0×1011 N/m2,poisson's ratio is 0.3. 2.The software used ANSYS Finite element software(APDL 15.0) 3.Solving process Point1 was choosed as the Coordinate point, horizontal to the right was the X axis，the upright direction is choosed as the Yaxis to create a coordinate system. The nodes was numberedas shown in Figure 1，node 1 and node 5 was fixed，and the force on node 6,7,8was is 1k N，the direction is opposite to the Y-axis 3.1 The preparatory work before analysis （1）Specify the new file name. Select Utility>Menu> File>Change Jobname, then pop-up the dialog box Change Jobname,inputthe the working file name ‘2D-sp’ in the Enter New Jobname, click OK to finish the difinition, as shown in Fig.2. Fig.2.The difinition of working file name. （2）Specify a new title.Select Utility>Menu>File>Change Title，then pop-up the dialog box Change Title,inputthe the file name ‘2D-sp pro’ in the Enter New Title,click OK to finish the difinition, as shown in Fig.3.

有限元分析软件及应用

3.5 ANSYS软件加载、求解、后处理技术 3.5.1 ANSYS 3.5.1 ANSYS 荷载概述荷载概述 在这一节中将讨论: 有限元分析软件及应用 8 有限元分析软件及应用 8 A. 载荷分类 3.5 ANSYS 软件加载、求解、后处理技术 3.5 ANSYS 软件加载、求解、后处理技术 B. 加载 C. 节点坐标系 D. 校验载荷 孙瑛 孙瑛 E. 删除载荷 哈哈尔尔滨滨工工业业大学空大学空间结间结构研构研究中心究中心 2010秋 2010秋 SSRC SSRC 1/ 76 S Space pace S Stru truc ctu ture re R Res esear earc ch h C Center enter, H , HI IT, T, CH CHIN INA A

理技术 A. 载荷分类 B. 加载 A. 载荷分类 B. 加载 ANSYS中的载荷可分为: 可在实体模型或 FEA 模型节点和单元上加载自由度DOF - 定义节点的自由度( DOF )值结构分析_ 沿单元边界均布的压力 沿线均布的压力 位移集中载荷 - 点载荷结构分析_力面载荷 - 作用在表面的分布载荷结构分析_压力 在关键点处 在节点处约 约束体积载荷 - 作用在体积或场域内热分析_ 体积膨胀、内生 束 成热、电磁分析_ magnetic current density等实体模型 FEA 模型惯性载荷 - 结构质量或惯性引起的载荷重力、角速度等 在关键点加集中力在节点加集中力 SSR SSRC C SSR SSRC C 2/ 76 3/ 76 S Space pace S Stru truc ctu ture re R Res esear earc ch h C Center enter, H , HI IT, T, CH CHIN INA A S Space pace S Stru truc ctu ture re R Res esear earc ch h C Center enter, H , HI IT, T, CH CHIN INA A

有限元分析及应用大作业

有限元分析及应用大作业 作业要求： 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

有限元法大作业

有限元法大作业 一平面刚架的程序 用Visual C++编制的平面刚架的源程序如下： ///////////////////////////////////////////////////////程序开始////////////////////////////////////////////////////////////////// #include"iostream.h" #include"math.h" #include"stdlib.h" #include"conio.h" //***************** //声明必要变量 //***************** #define PI 3.141592654 int NE; //单元数 int NJ; //节点数 int NZ; //支承数 int NPJ; //有节点载荷作用的节点数 int NPF; //非节点载荷数 int HZ; //载荷码 int E; //单元码 int fangchengshu; double F[303]; //各节点等效总载荷数值 int dym_jdm[100][2]; //单元码对应的节点码:dym_jdm[][0], dym_jdm[][1]分为前后节点总码 int zhichengweizhi[300]; //记录支持节点作用点的数组 int fjzhzuoyongdanyuan[100]; //非节点载荷作用单元 int fjzhleixing[100]; //非节点载荷类型:1-均布,2-垂直集中,3-平行集中,4-力偶,5-角度集中 double fjzhzhi[100]; //非节点载荷的值 double fjzhzuoyongdian[100]; //非节点载荷在各竿的作用点 double fjzhjiaodu[100]; //非节点载荷作用角度 int jdzhzuoyongdian[100]; //节点载荷作用的节点数组 double jiedianzaihe[101][3];//节点载荷值,其jiedianzaihe[][0]-- jiedianzaihe[][2]分别为U, V, M double zhengtigangdu[303][303]; //整体刚度数组 double changdu[100]; //各单元竿长数组 double jiaodu[100]; //各单元角度数组 double tanxingmoliang[100]; //各单元弹性模量数组 double J_moliang[100]; //各单元J模量数组 double mianji[100]; //各单元面积数组 double weiyi[303]; //记录各个节点位移的数组 double dy_weiyi[100][6]; //各个单元在局部坐标系中的位移数组dy_weiyi[i][0]-dyweiyi[i][6]分别为第i+1单元的u1,v1,@1,u2,v2,@2 double dy_neili[100][6];//各个单元在局部坐标系中的固端内力dy_weiyi[i][0]-dyweiyi[i][6]分别为第i+1单元的U1,V1,M1,U2,V2,M2 double gan_neili[100][6];//各个单元的竿端内力数组,gan_neili[i][6]表示第i+1单元的6内力. //*******************

有限元分析及应用例子FEM14

第9章受内外压筒体的有限元建模与应力变形分析(Project 2) 计算分析模型如图9-1 所示, 习题文件名: cylinder。 X (a) σO=100N/mm2 σI =200N/mm2 γ =7.85g/cm3 μ =0.3 E =210000N/mm2 (b) 图9-1 计算分析模型 9.1进入ANSYS 程序→ANSYSED 6.1ed →Interactive →change the working directory into yours→input Initial jobname: cylinder→Run 9.2 设置计算类型 ANSYS Main Menu: Preferences…→select Structural →OK 9.3 选择单元类型 ANSYS Main Menu: Preprocessor → Element Type →Add/Edit/Delete… → Add… →select Solid Quad 4node 42 →Apply →select Solid Brick 8node 45 → OK → Close (the Element

Types window) 9.4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Materials Models →Structural→Lineal →Elastic→Isotropic…→input EX:2.1e5, PRXY:0.3→ OK 关闭材料定义窗口 9.5构造筒体模型 ?生成模型截平面 ANSYS Main Menu: Preprocessor →Modeling→Create →Keypoints →In Active CS… →按次序输入横截平面的十个特征点和旋转对称轴上两点坐标（十个特征点：(300,0,0), (480,0,0), (480,100,0), (400,100,0), (400,700,0), (480,700,0), (480,800,0), (300,800,0), (300,650,0), (300,150,0)，对称轴上两点：(0,0,0), (0,800,0)）（每次输入完毕，用Apply结束，0可以不输入） →Cancel (back to Create window) →-Areas- Arbitrary → Through KPs →依次连接截面边线上的十个特征点(注意在选完第10点后结束，不要再选第1点)→ OK ?对平面进行网格划分 ANSYS Main Menu: Preprocessor →Meshing→Mesh Tool →(Size Controls) Globl: Set →input SIZE (element edge length): 50 →OK (back to MeshTool window)→Mesh → Pick All (in Picking Menu) → Close( the MeshTool window) ?用旋转法生成筒体模型 ANSYS Main Menu: Preprocessor →Modeling→Operate →Extrude→Elem Ext Opts→select TYPE:SOLID 45→Element sizing options for extrusion No. Elem divs: 1→OK (back to Extrude window)→Areas →About Axis →Pick All(in Picking Menu)→OK→Pick the two keypoints (11,12) of the Symmetrical Axis → OK→input ARC: 90; NSEG: 3→ OK 9.6 模型加位移约束 ANSYS Main Menu: Solution→Define Loads →Apply→Structural→Displacement ?两截面分别加Z, X方向的约束 ANSYS Utility Menu: Select → Entities…→Nodes → By Location →select X coordinates →input 0→ OK (back to Displacement window)→On Nodes → Pick All(in Picking Menu) → select Lab2:UX →OK →ANSYS Utility Menu: Select → Everything ANSYS Utility Menu: Select → Entities…→ Nodes → By Location →select Z coordinates →input 0→ OK (back to Displacement window)→On Nodes →Pick All(in Picking Menu) → select Lab2:UZ →OK →ANSYS Utility Menu: Select →Everything ?底面加Y方向的约束 ANSYS Utility Menu: Select → Entities… → Nodes → By Location →select Y coordinates →input 0→ OK (back to Displacement window)→On Nodes →Pick All(in Picking Menu) →

有限元分析报告大作业

基于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。

《有限元分析及应用》配书盘说明

《有限元分析及应用》配书盘 曾攀 (清华大学机械工程系) 说明 该配书盘针对《有限元分析及应用》一书中有关有限元分析的自主程序开发、与ANSYS平台的衔接、基于ANSYS的有限元建模、基于MARC的有限元建模的章节，提供相应的电子材料及文档，以便在进行实际编程和应用国际著名商业软件进行建模和分析时参考。电子文档材料包括三大部分：(1)有限元分析源程序(f,c,ANSYS衔接)；(2) 四类问题有限元分析的操作指南(ANSYS,MARC)；(3) ANSYS一般性帮助文件。具体的文件目录和清单如下。 在目录/有限元分析源程序(f,c,ANSYS衔接)/中有以下内容 (1) 使用说明文件 自主程序开发使用说明(fortran,C,ANSYS平台衔接).pdf (2 ) 在子目录/fortran源程序及与ANSYS衔接(FEM2D)/中有以下文件 源程序文件： FEM2D.FOR 程序需读入的数据文件： BASIC.IN(模型的基本信息文件，需手工生成) NODE_ANSYS.IN (节点信息文件，可由ANSYS前处理导出，或手工生成) ELEMENT_ANSYS.IN(单元信息文件，可由ANSYS前处理导出，或手工生成)程序输出的数据文件： DATA.OUT (一般结果文件) FOR_POST.DAT(专供ANSYS进行后处理的结果数据文件) 与ANSYS后处理衔接的接口程序： USER_POST.LOG(在ANSYS中进行后处理的命令流文件) (3 ) 在子目录/c源程序及与ANSYS衔接(JIEKOU)/中有以下文件 源程序文件： JIEKOU.CPP 程序需读入的数据文件： NODE_ANSYS.IN(从ANSYS前处理导出的节点信息文件) ELEMENT_ANSYS.IN(从ANSYS前处理导出的单元信息文件) INPUT.DAT(包含除网格划分信息之外的所有前处理信息) 程序输出的数据文件：

有限元分析的典型应用领域

有限元分析的典型应用领域 6-4：对该问题进行有限元分析的过程如下。 (1)进入ANSYS（设定工作目录和工作文件） 程序→ANSYS→ANSYS Interactive →Working directory（设置工作目录）→Initial jobname（设置工作文件名）：Press →Run →OK (2)设置分析特性 ANSYSMain Menu：Preferences…→Structural →OK (3)定义单元类型 ANSYSMain Menu：Preprocessor →Element Type →Add/Edit/Delete... →Add…→Solid: Quad 4node 42 →OK（返回到Element Types窗口）→Options…→K3：Plane Strs w/thk（带厚度的平面应力问题）→OK →Close (4)定义材料参数 ANSYSMain Menu：Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic：EX：2.1e11（弹性模量），PRXY：0.3（泊松比）→OK →点击该窗口右上角的“×”来关闭该窗口 (5)定义实常数以确定平面问题的厚度 ANSYSMain Menu：Preprocessor →Real Constants…→Add/Edit/Delete →Add →Type 1 PLANE42 →OK →Real Constant Set No：1（第1号实常数），THK：3.4（平面问题的厚度）→OK →Close （6）生成几何模型 生成上拱形梁 ANSYSMain Menu：Preprocessor →Modeling →Create →Keypoints→In Active CS →NPT Keypoint number：1，X，Y，Z Location in active CS：-4.5，8.5→Apply →同样输入后5个特征点坐标（坐标分别为（-2.25，8.5），（2.25，8.5），（4.5，8.5），（0，13），（0，10.75））→OK →Lines →Lines →Straight Line 用鼠标分别连接特征点1，2和3，4生成直线→OK→Arcs →By End KPs & Rad →用鼠标点击特征点2，3 →OK →用鼠标点击特征点6 →OK →RAD Radius of the arc：2.25→Apply （出现Warning对话框，点Close关闭）→用鼠标点击特征点1，4 →OK →用鼠标点击特征点5 →OK →RAD Radius of the arc：4.5→OK（出现Warning对话框，点Close关闭）→Areas →Arbitrary →By Lines →用鼠标点击刚生成的线→OK 生成下拱形梁 ANSYSMain Menu：Preprocessor →Modeling →Create →Keypoints→In Active CS →NPT Keypoint number：7，X，Y，Z Location in active CS：-4.5，-8.5→Apply →同样输入后5个特征点坐标（坐标分别为（-2.25，-8.5），（2.25，-8.5），（4.5，-8.5），（0，-13），（0，-10.75）→OK →Lines→Lines →Straight Line →用鼠标分别连接特征点7，8和9，10生成直线→OK →Arcs →By End KPs & Rad →用鼠标点击特征点8，9 →OK用鼠标点击特征点12 →OK →RAD Radius of the arc：2.25→Apply （出现Warning对话框，点Close关闭）→用鼠标点击特征点7，10 →OK →用鼠标点击特征点11 →OK →RAD Radius of the arc：4.5→OK（出现Warning对话框，点Close关闭）→Areas →Arbitrary →By Lines →用鼠标点击刚生成的线→OK 生成两根立柱 ANSYSMain Menu：Preprocessor →Modeling →Create →Areas →Rectangle →By 2 Corners →WP X：-4.5，WP Y：-8.5，Width：2.25，Height：17→Apply →WP X：2.25，WP Y：-8.5，Width：2.25，Height：17→OK 粘结所有面