有限元第二次作业教学提纲
2-2图示悬臂板，属于平面应力问题，其网格图及单元、节点编号见图2-1 , E=2.1X 1011 ,
u=0.28，演算其单刚阵到总刚阵的组集过程，并用MATLAB软件计算总刚阵。
答：根据图2-1所示列出单元节点列表:
节点
i j k
单儿
1 3 5 4
2 2 5 3
3 2 6 5
4 1 6 2
(1 )计算单元刚度阵
%
单元1的刚度矩阵：k10 0 0 0 0 0 0 0 0 0 0 0 0 0 比kJ k3,5 0 0 0 k；3 k；4 k4,5 0
；
0 0 比k5,4 k5,4 0 0 0 0 0 0 0
单元2的刚度矩阵：k2
0 0 0 0 0 0
0 k 2
k2,2
k2
k2,3
0 k 2
k2,5
0 k 2
k3,2
k2
k3,3
0 k 2
k3,5
0 0 0 0 0 0
；
0 k
2
k5,2
k；,3
0 k
2
k5,4
0 0 0 0 0 0 k5,4
k4,4
Matlab 程序语言的编写：
fun ctio n Idex
global gNode gEleme nt gMaterial gNode=[0.0 0.01
0.5 0.01 1.0 0.01 1.0 0.0 0.5 0.0
4
K k e
e 1
k 1
k 2
k 3
k 4
总刚度矩阵:
0 0 0 0 0 0 0 k ；,2 0 0 k ；5
k ；,6 0 0 0 0 0 0 0
0 0
0 ;
0 k ；2
0 0
k ；,5 k ；6 0 k ；2 0 0 k ；,5 k ：6
k ：1 k ：2
单元4的刚度矩阵：
k 4
kJ
k ；,2
k ；,6 ,
k ：1 k ：,2
k ；6
k ：1
k ：1 k 4
k 1,2
k 4
k 2,2
0 0 0 0 0 0
k ；6
0 0 0 0
0 k 4
0 0 0 0 0 0
0 0 0 0 0
k ：
1
0 0 0 k 6,6
k 4 k 1,1
k 4
k 2,1
k ；,2 ki\
k ；,2 k/
k ；,2
■ 1 .2 k 3,3
k 3,3
k i 4 ,6
k ；,6
4 2,6
k 6,i 4,3 1 2 k 5,3 k 5,3
k 3,4 k ；
4
k ；,4
k 5,5
i 4,5
k ：5 k ；5
3
6,5
k 6,6
4 6,6
单元3的刚度矩阵：k 3
k 3,5
2 3,5
0.0 0.0]
%gNode 同样是一个矩阵，每一行表示一个结点，第点的y 坐标
gElement=[3 4 5
2 3 5
2 5 6
1 2 6 ];
%gElement 是一个矩阵，每一行表示一个单元，第单元的第2 个结点号。
Return
function k=StiffnessMatrix(ie) %计算单元刚度矩阵函数global gNode gElement k=zeros(6,6);
E=2.1*10A11;
u=0.28 ;
t=0.01;
xi=gNode(gElement(ie,1),1);
yi=gNode(gElement(ie,1),2);
xj=gNode(gElement(ie,2),1);
yj=gNode(gElement(ie,2),2);
xm=gNode(gElement(ie,3),1);
ym=gNode(gElement(ie,3),2);
分量
ai=xj*ym-xm*yj;
1 列是结点的x 坐标，第
2 列是结
1 行是单元的第1 个结点号，第
2 行是
%6x6 单元刚阵
%材料特性
%材料特性
%材料特性
%计算节点坐标
aj=xm*yi-xi*ym;
am=xi*yj-xj*yi;
bi=yj-ym;
bj=ym-yi;
bm=yi-yj; ci=-
(xj-xm); cj=-
(xm-xi); cm=-
(xi-xj);
d=[1,xi,yi;1,xj,yj;1,xm,ym];
area=det(d);
B=[bi 0 bj 0 bm 0 ;0 ci 0 cj 0 cm;ci bi cj bj cm bm]; B=B/2/area;
D=[1 u 0;u 1 0;0 0 (1-u)/ 2];
D=D*E/(1-u A2);
k=transpose(B)*D*B*t*abs(area); 阵
Return
%计算单元面积%计算单元刚度矩
function gK=AssembleStiffnessMatrix
% 计算总刚阵
global gElement gK ie
gK=zeros(12,12);
for ie =1:1:4 %单元循环k=StiffnessMatrix(ie);
for i=1:1:3 %节点循环for j=1:1:3 %节点循环
for p=1:1:2 %自由度循环for q=1:1:2 %自由度循环
m=(i-1)*2+p; %每个节点有2个自由度，i 节点的第p 个自由度为
(i-1)*2+p
n=(j-1)*2+q; %每个节点有2个自由度，i节点的第p个自由度为
(i-1)*2+p
M=(gElement(ie,i)-1)*2+p;
N=(gElement(ie,j)-1)*2+q;
gK(M,N)=gK(M,N)+k(m,n);
end
end
end
end
end
Return
则单元 1 的刚度矩阵为
>> StiffnessMatrix(1)
ans =
1.0e+010 *
2.0508 0 -2.0508 0.0410 0 -0.0410
0 5.6966 0.0319 -5.6966 -0.0319 0
2.0508 0.0319 2.0531 -0.0729 -0.0023 0.0410
0.0410 -5.6966 -0.0729 5.6974 0.0319 -0.0008
0 -0.0319 -0.0023 0.0319 0.0023 0
-0.0410 0 0.0410 -0.0008 0 0.0008
单元 2 的刚度矩阵
>> StiffnessMatrix(2)
ans =
1.0e+010
2.0531 -0.0729 -2.0508 0.0319 -0.0023 0.0410
0.0729 5.6974 0.0410 -5.6966 0.0319 -0.0008
2.0508 0.0410 2.0508 0 0 -0.0410
0.0319 -5.6966 0 5.6966 -0.0319 0
0.0023 0.0319 0 -0.0319 0.0023 0
0.0410 -0.0008 -0.0410 0 0 0.0008
单元 3 的刚度矩阵为
>> StiffnessMatrix(3) ans =
1.0e+010
0.0023 0 -0.0023 0.0319 0 -0.0319
0 0.0008 0.0410 -0.0008 -0.0410 0
0.0023 0.0410 2.0531 -0.0729 -2.0508 0.0319
0.0319 -0.0008 -0.0729 5.6974 0.0410 -5.6966
0 -0.0410 -2.0508 0.0410 2.0508 0
-0.0319
单元 4 的刚度矩阵
>> StiffnessMatrix(4) ans =
1.0e+010
0 0.0319 -5.69660 5.6966
2.0531 -0.0729 -2.0508 0.0319 -0.0023 0.0410
0.0729 5.6974 0.0410 -5.6966 0.0319 -0.0008
2.0508 0.0410 2.0508 0 0 -0.0410
0.0319 -5.6966 0 5.6966 -0.0319 0
0.0023 0.0319 0 -0.0319 0.0023 0
0.0410 -0.0008 -0.0410 0 0 0.0008 总刚度矩阵为ans = 精品文档
华科大有限元分析题及大作业题答案——船海专业(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
有限元上机实验指导书
弹性力学及有限元实验上机指导书 土木工程教学部 20１5年6月 一、ANSYS软件安装及有限元建模方法演示１、实验目的 掌握Ansys商用有限元软件的安装及了解有限元建模方法２、实验任务 １）Ansys商用软件8.1安装过程详解。
２）采用有限元直接建模法创建杆系模型演示。 3）采用有限元间接建模法创实体模型演示。 模型一：平面桁架如下图所示，长度单位为m，求支座反力和各杆内力。弹性模量2.1e+11，泊松比0.3，杆件截面面积均为0.01m2。 1/6 1/3 1/2 图1 平面桁架模型 模型二：正方形带孔平板，边长1m，小孔直径0.1m，板厚0.05m。弹性模量2.1e+11，泊松比0.3。上下边受均匀压力1000N。 图2 带孔正方形平板 3、实验方法 实验方法同课堂操作演示。 命令流(模型一) /PREP7 ET,1,LINK1 R,1,0.01, , MP,EX,1,,2.1e+11
MP,PRXY,1,,0.3 n,1,0,0 n,2,1,0 n,3,1,0.5 n,4,2,0 n,5,2,0.833 n,6,3,0 n,7,3,1 n,8,4,0 n,9,4,0.833 n,10,5,0 n,11,5,0.5 n,12,6,0 e,1,2 e,1,3 e,2,3 e,2,4 e,3,5 e,4,5 e,3,4 e,4,6 e,5,7 e,6,7 e,4,7 e,6,8 e,7,9 e,8,9 e,7,8 e,8,10
e,9,11 e,10,11 e,8,11 e,10,12 e,11,12 save 命令流(模型二) 二、利用ANSYS创建杆系或实体结构有限元模型 １、实验目的 掌握有限元建模的基本方法并能创建简单的杆系结构和实体结构有限元模型。 ２、实验任务 （1）采用直接建模法创建上节平面桁架结构模型。 （2）采用间接建模法创建上节带孔平板实体模型。 3、实验方法 实验方法同操作演示。 三、有限元求解及结果后处理演示 １、实验目的 掌握基本参数的设置、荷载施加方法及后处理操作。 ２、实验任务 （1）读入数据文件（命令流）的形式生成杆系结构有限元模型 （2）实常数、材料参数、求解参数设置演示 （3）位移约束、集中荷载施加方法演示 （4）计算求解与后处理操作演示。 3、实验方法 实验方法同课堂操作演示 附：后处理GUI操作及命令流操作 A、通过后处理提取节点计算结果 三种后处理操作： 1、plot results
重庆大学有限元第二次作业(刘静老师)
【有限元分析技术】第二次作业 科 目： 有限元分析技术 教 师： 姓 名： 学 号： 班 级： 类 别： 学术型 上课时间： 2016 年 11 月至 2017 年 1 月 考 生 成 绩： 卷面成绩 平时成绩 课程综合成绩 阅卷评语： 阅卷教师 (签名) 大学研究生院
第一章 题目概况 1.1 原始数据 矩形板尺寸如下图，板厚为5mm ，弹性模量为522.010/E N mm =? ，泊松比为0.27μ= 图1.1 原始计算简图 1.2工况选择 （1）试按下表的载荷约束组合，任选2种进行计算，并分析其位移、应力分布的异同。 表1 两种不同工况的载荷及约束 序号 载荷 约束 备注 1 向下均布载荷P=5N/mm,作用于ab 边 c ，d 点固定 2 向下均布载荷P=5N/mm,作用于ab 边 a ，b 点固定 3 向下均布载荷P=5N/mm,作用于ab 边 a ，c 边固定 还可讨论a ，c 点固定 4 向下均布载荷P=5N/mm,作用于cd 边 c ，d 点简支 5 向下均布载荷P=5N/mm,作用于cd 边 a ，b 点简支 6 向下均布载荷P=5N/mm,作用于cd 边 a ，c 边固定 还可讨论a ，c 点固定 7 向下集中载荷F=1000N,作用于ab 边中点 c ，d 点简支 8 向下集中载荷F=1000N,作用于ab 边中点 a ，b 点简支 9 向下集中载荷F=1000N,作用于ab 边中点 a ，c 边固定 还可讨论a ，c 点固定 10 向下集中载荷F=1000N,作用于cd 边中点 c ，d 点简支 11 向下集中载荷F=1000N,作用于cd 边中点 a ，b 点简支 12 向下集中载荷F=1000N,作用于cd 边中点 a ，c 边固定 还可讨论a ，c 点固定 1.3 工况选择结果及分析任务 （1）工况选择结果 根据表1的工况，选取工况1，2，8进行对比分析，选取结果如表2所示，为了方便下文中分别将序号1、2、8的工况称为工况一、工况二、工况三。 表2 分析工况的载荷及约束 序号 载荷 约束 备注 1 向下均布载荷P=5N/mm,作用于ab 边 c ，d 点固定 工况一 2 向下均布载荷P=5N/mm,作用于ab 边 a ，b 点固定 工况二 8 向下集中载荷F=1000N,作用于ab 边中点 a ， b 点简支 工况三
有限元答案
1.1有限单元法中“离散”的含义是什么？有限单元法是如何将具有无限自由度的连续介质问题转变成有限自由度问题的？位移有限元法的标准化程式是怎样的？ （1）离散的含义即将结构离散化，即用假想的线或面将连续体分割成数目有限的单元，并在其上设定有限个节点；用这些单元组成的单元集合体代替原来的连续体，而场函数的节点值将成为问题的基本未知量。 （2）给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律，即通过插值以单元节点位移表示单元内任意点的位移。因节点位移个数是有限的，故无限自由度问题被转变成了有限自由度问题。 （3）有限元法的标准化程式：结构或区域离散，单元分析，整体分析，数值求解。 1.3单元刚度矩阵和整体刚度矩阵各有哪些性质？各自的物理意义是什么?两者有何区别? 单元刚度矩阵的性质：对称性、奇异性（单元刚度矩阵的行列式为零）。 整体刚度矩阵的性质：对称性、奇异性、稀疏性。 单元Kij物理意义Kij即单元节点位移向量中第j个自由度发生单位位移而其他位移分量为零时，在第j个自由度方向引起的节点力。 整体刚度矩阵K中每一列元素的物理意义是：要迫使结构的某节点位移自由度发生单位位移，而其他节点位移都保持为零的变形状态，在所有个节点上需要施加的节点荷载。 2.2什么叫应变能？什么叫外力势能？试叙述势能变分原理和最小势能原理，并回答下述问题：势能变分原理代表什么控制方程和边界条件？其中附加了哪些条件？ (1)在外力作用下，物体内部将产生应力ζ和应变ε，外力所做的功将以变形能的形式储存起来，这种能量称为应变能。 (2)外力势能就是外力功的负值。 (3)势能变分原理可叙述如下：在所有满足边界条件的协调位移中，那些满足静力平衡条件的位移使物体势能泛函取驻值，即势能的变分为零 δΠp=δUε+δV=0 此即变分方程。对于线性弹性体，势能取最小值，即 δ2ΠP=δ2Uε+δ2V≧0 此时的势能变分原理就是著名的最小势能原理。 势能变分原理代表平衡方程、本构方程和应力边界条件，其中附加了几何方程和位移边界条件。 2.3什么是强形式？什么是弱形式？两者有何区别？建立弱形式的关键步骤是什么？ 等效积分形式通过分部积分，称式 ∫ΩC T(v)D(u)dΩ+∫ΓE T(v)F(u)dΓ为微分方程的弱形式，相对而言，定解问题的微分方程称为强形式。 区别：弱形式得不到解析解。 建立弱形式的关键步骤：对场函数要求较低阶的连续性。 2.4为了使计算结果能够收敛于精确解，位移函数需要满足哪些条件？为什么？ 只要位移函数满足两个基本要求，即完备性和协调性，计算结果便收敛于精确解。 2.6为什么采用变分法求解通常只能得到近似解？变分法的应用常遇到什么困难？Ritz法收敛的条件是什么？ （1）在Ritz 法中，N决定了试探函数的基本形态，待定参数使得场函数具有一定的任意性。如果真实场函数包含在试探函数之内，则变分法得到的解答是精确的；如果试探函数取自完全的函数序列，则当项数不断增加时，近似解将趋近于精确解。然而，通常情况下试探函数不会将真实场函数完全包含在内，实际计算时也不可能取无穷多项。因此，试探函数只能是真实场函数的近似。可见，变分法就是在某个假定的范围内找出最佳解答，近似性就源于此。 （2）采用变分法近似求解，要求在整个求解区域内预先给出满足边界条件的场函数。通常情况下这是不可能的，因而变分法的应用受到了限制。 （3）Ritz 法的收敛条件是要求试探函数具有完备性和连续性，也就是说，如果试探函数满足完备性和连续性的要求，当试探函数的项数趋近于无穷时，则Ritz 法的近似解将趋近于数学微分方程的精确解。 3.1构造单元形函数有哪些基本原则？形函数是定义于单元内坐标的连续函数。单元位移函数通常采用多项式，其中的待定常数应该与单元节点自由度数相等。为满足完备性要求，位移函数中必须包括常函数和一次式，即完全一次多项式。多项式的选取应由低阶到高阶，尽量选择完全多项式以提高单元的精度。若由于项数限制而不能选取完全多项式时，也应使完全多项式具有坐标的对称性，并且一个坐标方向的次数不应超过完全多项式的次数。有时为了使位移函数保持一定阶次的完全多项式，可在单元内部配置节点。然而，这种节点的存在将增加有限元格式和计算上的复杂性，除非不得已才加以采用。形函数应保证用它定义的位移函数满足收敛要求，即满足完备性要求和协调性条件。 3.1构造单元形函数有哪些基本原则？试采用构造单元的几何方法，构造T10 单元的形函数，并对其收敛性进行讨论。 通常单元位移函数采用多项式，其中的待定常数由节点位移参数确定，因此其个数应与单元节点自由度数相等。根据实体结构的几何方程，单元的应变是位移的一次导数。为了反映单元刚体位移和常应变即满足完备性要求，位移函数中必须包含常数项和一次项，即完全一次多项式。 3.3何谓面积坐标？其特点是什么？为什么称其为自然坐标或局部坐标？ （1）三角形单元中，任一点P(x,y)与其3个角点相连形成3个子三角形,其位置可以用下述称为面积坐标的三个比值来确定： L1=A1/A L2=A2/A L3=A3/A 其中A1,A2,A3分别为P23,P31,P12的面积。 （2）面积坐标的特点： a T3单元的形函数Ni就是面积坐标Li b面积坐标与三角形在整体坐标系中的位置无关。 c三个节点的面积坐标分别为节点1（1, 0, 0）、节点2（0, 1, 0）、节点3（0, 0, 1），形心的面积坐标为（1/3, 1/3, 1/3）。 d单元边界方程为Li=0(i=1,2,3) e在平行于23边的一条直线上，所有点都有相同的面积坐标L1（L1对应的三角形具有相同的高和底边），而且L1就等于此直线至23边的距离与节点1至23边的距离之比值。
大作业报告参考2有限元学习心得
有限元学习心得 吴清鸽车辆工程 50110802411 短短八周的有限元课已经结束。关于有限元，我一直停留在一个很模糊的概念。我知道这是一个各个领域都必须涉及的点，只要有关于CAE分析的，几乎都要涉及有限元。总体来说，这是一门非常重要又有点难度的课程。 有限元方法(finite element method) 或有限元分析(finite element analysis)，是 求取复杂微分方程近似解的一种非常有效的工具，是现代数字化科技的一种重要 基础性原理。将它用于在科学研究中，可成为探究物质客观规律的先进手段。将 它应用于工程技术中，可成为工程设计和分析的可靠工具。本课程教学基本内容 有固体力学和结构力学简介；有限元法基础；桁架、梁、刚架、二维固体、板和 壳、三维固体的有限元法；建模技术；热传导问题的有限元分析；PATRAN软件 的使用. 通过有限元分析课程学习使我了解和掌握了一些有限元知识： 1．简要了解二维和三维固体以及桁架、梁和板结构的三组基本力学方程，即表示位移-应变关系的几何方程，表示应力-应变关系的本构方程和表示内力-外力关系的平衡方程。 2．了解利用能量法形成有限元离散系统方程的基本原理，即哈密尔顿原理。掌握有限元分 析的基本方法及步骤，包括域的离散、位移插值、构造形函数、单元有限元方程 的建立、坐标变换、整体有限元方程的组装、整体有限元方程的求解技术。 3．具体深入的了解并掌握桁架结构、梁结构、刚架结构、二维固体、板和壳结构、三维固体的有限元法分析技术，包括他们具体的形函数构造，应变矩阵，局部坐标系和整体坐标系中的单元矩阵。各种结构的实例研究。 4．了解并掌握建立高质量建模所涉及的各种关键技术。包括单元类型的选择，单元畸形的限制，不同阶数单元混用时网格的协调性问题，对称性的应用（平面对称、轴对称、旋转对称、重复对称），由多点约束方程形成刚域及应用（模拟偏移、不同自由度单元的连接、网格协调性的施加）等，以及多点约束方程的求解。以PATRAN有限元通用软件为例了解一般商业有限元软件的组成及结构。掌握PATRAN软件的基本使用。利用PATRAN软件上机实践完成两个上机练习：刚架结构有限元分析和三维固体有限元分析。 课程的具体学习内容： 内容： 1.三节点三角形单元：单元分析、总刚度矩阵组装、引入约束条件修正总刚度 矩阵、载荷移置、方程求解； 2.四边形单元分析、四节点四面体单元分析、八节点六面体单元分析；
有限元 第二次作业教学提纲
2-2 图示悬臂板，属于平面应力问题，其网格图及单元、节点编号见图2-1，E=2.1×1011，u=0.28，演算其单刚阵到总刚阵的组集过程，并用MATLAB 软件计算总刚阵。 图2-1 答：根据图2-1所示列出单元节点列表： i j k 1 3 5 4 2 2 5 3 3 2 6 5 4 1 6 2 （1）计算单元刚度阵 单元1的刚度矩阵：[] ????? ?????=15,514 ,513 ,515,414,41 3 ,415,314,313,3 1k k k k k k k k k k ，[] ?????????? ??????????=000000 000 00000000 000000 14 ,514 ,513 ,515,414,41 3,41 5,31 4,31 3 ,31 k k k k k k k k k k ； 单元2的刚度矩阵：[] ??? ? ????? ?=25 ,523 ,522 ,525,323,322,325,223,222,22 k k k k k k k k k k ，[] ????????? ?????? ?? ???=00 000000000000000000 0024,523,522,525,323 ,322,32 5 ,22 3 ,22 2,22k k k k k k k k k k ； 节点 单元
单元3的刚度矩阵：[] ??? ? ????? ?=36 ,635 ,632 ,636,535,532,536,235,232,23 k k k k k k k k k k ，[] ????????? ? ????? ?? ???=36,635 ,632,636,535,532 ,536,23 5 .23 2,2300000000 00 0000000000000000k k k k k k k k k k ； 单元4的刚度矩阵：[] ???? ? ???? ?=46,642 ,641,646,242,241 ,246,142,141 ,14 k k k k k k k k k k ，[] ?? ? ??? ? ?? ? ??????????=46,641 ,646,242,241.246,142 ,141,140000 000000000000000000000000k k k k k k k k k ； 总刚度矩阵：[][][][][][]4 3 2 1 4 1 k k k k k K e e +++=∑== []??????? ?? ? ????? ?????++++++++++++=4 6,636,635 ,642 ,632,641 ,636 ,535 ,525,515,514,523 ,513,53 2,522,515,414,413,425 ,315,314,323 ,313,322 ,346,236,235 ,225,223 ,242 ,232,222,241,246,142 ,141 ,10 00000 00 000k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k K Matlab 程序语言的编写： function Idex global gNode gElement gMaterial gNode=[0.0 0.01 0.5 0.01 1.0 0.01 1.0 0.0 0.5 0.0
重庆大学有限元第一次作业
有限元分析技术课程大作业 科 目：有限元分析技术 教 师： 姓 名： 学 号： 专 业： 机械设计及理论 类 别： 学 术 上课时间： 2016 年 11 月至 2017 年 1 月 考 生 成 绩： 阅卷评语： 阅卷教师 (签名) 重庆大学研究生院
第一章 问题提出 1.1工程介绍 某露天大型玻璃平面舞台的钢结构如图1所示，每个分格（图2中每个最小的矩形即为一个分格）x 方向尺寸为1m ，y 方向尺寸为1m ；分格的列数（x 向分格）=学生序号的百位数值×10+十位数值+5，分格的行数（y 向分格）=学生序号的个位数值+4，如序号为041的同学分格的列数为9，行数为5，111号同学分格的列数为16，行数为5。 钢结构的主梁（图1中黄色标记单元）为高160宽100厚14的方钢管，其空间摆放形式如图3所示；次梁（图1中紫色标记单元）为直径60厚10的圆钢管（单位为毫米），材料均为碳素结构钢Q235；该结构固定支撑点位于左右两端主梁和最中间（如不是正处于X 方向正中间，偏X 坐标小处布置）的次梁的两端，如图2中标记为UxyzRxyz 处。 玻璃采用四点支撑与钢结构连接（采用四点支撑表明垂直作用于玻璃平面的面载荷将传递作用于玻璃所在钢结构分格四周的节点处，表现为点载荷，如图4所示）；试对在垂直于玻璃平面方向的22 /KN m 的面载荷（包括玻璃自重、钢结构自重、活载荷（人员与演出器械载荷）、风载荷等）作用下的舞台进行有限元分析.（每分格面载荷对于每一支撑点的载荷可等效于0.5KN 的点载荷）。 1.2 作业内容 （1）屏幕截图显示该结构的平面布置结构，图形中应反映所使用软件的部分界面，如图1-2； （2）该结构每个支座的支座反力； （3）该结构节点的最大位移及其所在位置； （4）对该结构中最危险单元（杆件）进行强度校核。 图1-1
有限单元法部分课后题答案
1.1 有限单元法中“离散”的含义是什么？有限单元法是如何将具有无限自由度的连续介质问题转变成有限自由度问题的？位移有限元法的标准化程式是怎样的？ （1）离散的含义即将结构离散化，即用假想的线或面将连续体分割成数目有限的单元，并在其上设定有限个节点；用这些单元组成的单元集合体代替原来的连续体，而场函数的节点值将成为问题的基本未知量。 （2）给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律，即通过插值以单元节点位移表示单元内任意点的位移。因节点位移个数是有限的，故无限自由度问题被转变成了有限自由度问题。 （3）有限元法的标准化程式：结构或区域离散，单元分析，整体分析，数值求解。 1.3 单元刚度矩阵和整体刚度矩阵各有哪些性质？各自的物理意义是什么?两者有何区别?单元刚度矩阵的性质：对称性、奇异性（单元刚度矩阵的行列式为零）。整体刚度矩阵的性质：对称性、奇异性、稀疏性。单元 Kij 物理意义 Kij 即单元节点位移向量中第 j 个自由度发生单位位移而其他位移分量为零时，在第 j 个自由度方向引起的节点力。整体刚度矩阵 K 中每一列元素的物理意义是：要迫使结构的某节点位移自由度发生单位位移，而其他节点位移都保持为零的变形状态，在所有个节点上需要施加的节点荷载。 2.2 什么叫应变能？什么叫外力势能？试叙述势能变分原理和最小势能原理，并回答下述问题：势能变分原理代表什么控制方程和边界条件？其中附加了哪些条件？ (1)在外力作用下，物体内部将产生应力σ和应变ε，外力所做的功将以变形能的形式储存起来，这种能量称为应变能。 (2)外力势能就是外力功的负值。 (3)势能变分原理可叙述如下：在所有满足边界条件的协调位移中，那些满足静力平衡条件的位移使物体势能泛函取驻值，即势能的变分为零 δ∏p=δ Uε+δV=0 此即变分方程。对于线性弹性体，势能取最小值，即 δ2∏P=δ2Uε+δ2V≥0 此时的势能变分原理就是著名的最小势能原理。 势能变分原理代表平衡方程、本构方程和应力边界条件，其中附加了几何方程和位移边界条件。 2.3 什么是强形式？什么是弱形式？两者有何区别？建立弱形式的关键步骤是什么？ 等效积分形式通过分部积分，称式 ∫ΩCT(v)D(u)dΩ+∫ΓET(v)F(u)dΓ 为微分方程的弱形式，相对而言，定解问题的微分方程称为强形式。 区别：弱形式得不到解析解。建立弱形式的关键步骤：对场函数要求较低阶的连续性。2.4 为了使计算结果能够收敛于精确解，位移函数需要满足哪些条件？为什么？ 只要位移函数满足两个基本要求，即完备性和协调性，计算结果便收敛于精确解。 2.6 为什么采用变分法求解通常只能得到近似解？变分法的应用常遇到什么困难？Ritz 法收敛的条件是什么？ （1）在 Ritz 法中，N 决定了试探函数的基本形态，待定参数使得场函数具有一定的任意性。如果真实场函数包含在试探函数之内，则变分法得到的解答是精确的；如果试探函数取自完全的函数序列，则当项数不断增加时，近似解将趋近于精确解。然而，通常情况下试探函数不会将真实场函数完全包含在内，实际计算时也不可能取无穷多项。因此，试探函数只能是真实场函数的近似。可见，变分法就是在某个假定的范围内找出最佳解答，近似性就源于此。 （2）采用变分法近似求解，要求在整个求解区域内预先给出满足边界条件的场函数。通常情况下这是不可能的，因而变分法的应用受到了限制。 （3）Ritz 法的收敛条件是要求试探函数具有完备性和连续性，也就是说，如果试探函数满足完备性和连续性的要求，当试探函数的项数趋近于无穷时，则 Ritz 法的近似解将趋近于数学微分方程的精确解。 3.1 构造单元形函数有哪些基本原则？ 形函数是定义于单元内坐标的连续函数。单元位移函数通常采用多项式，其中的待定常数应该与单元节点自由度数相等。为满足完备性要求，位移函数中必须包括常函数和一次式，即完全一次多项式。多项式的选取应由低阶到高阶，尽量选择完全多项式以提高单元的精度。若由于项数限制而不能选取完全多项式时，也应使完全多项式具有坐标的对称性，并且一
有限元作业第二次作业
土木工程专业 有限元第二次作业 姓名： 班级： 学号： 指导教师： 二〇一五年6月12日
习 题：平面应力问题的八节点等参元，已给定8个节点 的坐标。试查资料并论述： 1、单元中位移函数u （ξ，η），v （ξ，η）和单元节点位 移{δe }的关系式； 2、[ B ]矩阵的计算步骤和计算式； 3、单元刚度矩阵[ k e ]的一般计算方法和计算步骤； 4、论述相邻单元间公共边界上位移的连续性； 5、如果给定母单元中点A ， （ξ，η），怎样求实际单元中与 A ， 相对应的点A （x ，y ）；反之，如果给定实际单元中的点A （x ，y ），怎样求其在母单元中对应点A ， （ξ，η）？ 6、如果已经求解得到单元8个节点的位移值{δe }怎样求单 元中某一点B （x ，y ）的应力？ 实际单元 1 2 6 7 Y 1 2 43 67 8η= 1η=﹣1 母单元 ξ= 1 ξ=﹣1
解： 1、此题分两步进行： 单元位移场的表达： 如图1所示,在任意四边形的每边中间设一附加节点，则单元边界就变成二次曲线的了。如果直接在整体坐标系(),x y 下，像八节点矩形元那样，构造双二次多项式的位移插值函数，则因曲边四边形单元边界是二次曲线，故边界上的位移是()x y 或的五次多项式， 它不能由曲边上三个节点的位移分量唯一地决定，从而不能保证相邻两个单元在公共边上位移的协调条件，所以在整体坐标系(),x y 下构造完全协调的位移插值 函数是很困难的，利用坐标变换，可将曲边四边形单元变换成基本单元，如图2所示的在自然坐标(),ξη下具有边长为2的八节点正方形单元，自然坐标系(),ξη是外节点坐标值为±1的局部坐标系。在自然坐标系的单元上构造 协调的位移插值函数，其形状函数是较普通的，取位移分量为,ξη的双二次多项式, 即： 2222 12345678222 2910111213141516u a a a a a a a a v a a a a a a a a ξηξξηηξηξηξηξξηηξηξη?=+++++++??=+++++++?? （1-1） 利用8 个节点的16 个位移分量可唯一确定16 个待定常数1216,,a a a …,，若代入8个节点的局部坐标值，得： 图1：在总坐标系中具有二 次曲边的四边形单元 图2：在自然坐标系中的 曲边四边形的基本单元
现代设计方法(关于有限元)作业
《现代设计方法》作业关于有限元法的研究 学院：机械工程学院 专业：机械制造及其自动化
0.有限元法 有限元法分析起源于50年代初杆系结构矩阵的分析。随后，Clough于1960年第一次提出了“有限元法”的概念。其基本思想是利用结构离散化的概念，将连续介质体或复杂结构体划分成许多有限大小的子区域的集合体，每一个子区域称为单元（或元素），单元的集合称为网格，实际的连续介质体（或结构体）可以看成是这些单元在它们的节点上相互连接而组成的等效集合体；通过对每个单元力学特性的分析，再将各个单元的特性矩阵组集成可以建立整体结构的力学方程式，即力学计算模型；按照所选用计算程序的要求，输入所需的数据和信息，运用计算机进行求解。 当前，有限元方法/理论已经发展的相当成熟和完善，而计算机技术的不断革新，又在很大程度上推进了有限元法分析在工程技术领域的应用。然而，如此快速地推广和应用使得人们很容易忽视一个前提，即有限元分析软件提供的计算结果是否可靠、满足使用精度的前提，是合理地使用软件和专业的工程分析。有限元法分析一般包括四个步骤：物理模型的简化、数学模型的程序化、计算模型的数值化和计算结果的分析。每一个步骤在操作过程中都或多或少地引入了误差，这些误差的累积最终可能会对计算结果造成灾难性的影响，进而蒙蔽我们的认识和判断。 1.受内压空心圆筒的轴对称有限元分析 例图1.1所示为一无限长的受内压的轴对称圆筒，该圆筒置于内径为120mm的刚性圆孔中，试求圆筒内径处的位移。结构的材料参数
为：200 =，0.3 E GPa μ=。 图1 结构图 对该问题进行有限元分析的过程如下。 （1）结构的离散化与编号 由于该圆筒为无限长，取出中间一段（20mm高），采用两个三角形轴对称单元，如图1.2所示。对该系统进行离散，单元编号及结点编号如图1.3所示，有关结点和单元的信息见表1.1。 图1.2 有限元模型
有限元分析报告大作业
有限元分析》大作业基本要求： 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-2 图示悬臂板,属于平面应力问题,其网格图及单元、节点编号见图2-1,E=2、1×1011,u=0、28,演算其单刚阵到总刚阵得组集过程,并用MATLAB 软件计算总刚阵。 图2-1 答:根据图2-1所示列出单元节点列表: i j k 1 3 5 4 2 2 5 3 3 2 6 5 4 1 6 2 (1)计算单元刚度阵 单元1得刚度矩阵: ,; 单元2得刚度矩阵:,; 单元3得刚度矩阵:,; 单元4得刚度矩阵:,; 总刚度矩阵: []??????? ?? ? ????? ?????++++++++++++=4 6,636,635 ,642 ,632,641 ,636 ,535 ,525,515,514,523 ,513,53 2,522,515,414,413,425 ,315,314,323 ,313,322 ,346,236,235 ,225,223 ,242 ,232,222,241,246,142 ,141 ,10 00000 00 000k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k K 节点 单元
Matlab 程序语言得编写: function Idex global gNode gElement gMaterial gNode=[0、0 0、01 0、5 0、01 1、0 0、01 1、0 0、0 0、5 0、0 0、0 0、0] %gNode 同样就是一个矩阵,每一行表示一个结点,第1 列就是结点得x 坐标,第2 列就是结点得y坐标 gElement=[3 4 5 2 3 5 2 5 6 1 2 6 ]; %gElement 就是一个矩阵,每一行表示一个单元,第1 行就是单元得第1 个结点号,第2 行就是单元得第2个结点号。 Return function k=StiffnessMatrix(ie) %计算单元刚度矩阵函数 global gNode gElement k=zeros(6,6); %6x6单元刚阵 E=2、1*10^11; %材料特性 u=0、28 ; %材料特性 t=0、01; %材料特性 xi=gNode(gElement(ie,1),1);
有限元复习题答案
1、何为有限元法？其基本思想是什么？ 有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法，该方法以计算机为手段，采用分片近似，进而逼近整体的研究思想求解物理问题。 基本思想是化整为零集零为整。 2、为什么说有限元法是近似的方法，体现在哪里？ 有两点：用离散单元的组合体来逼近原始结构，体现了几何上的近似；而用近似函数逼近未知变量在单元内的真实解，体现了数学上的近似。 3、单元、节点的概念？ 节点：表达实际结构几何对象之间相互连接方式的概念 单元：网格划分中的每一个小部分称为单元，网格间相互联结点称为节点 4、有限元法分析过程可归纳为几个步骤？ 结构离散化、单元分析、整体分析 5、有限元方法分几种？本课程讲授的是哪一种？ 位移法、力法、混合法本课程讲授位移法 6、弹性力学的基本变量是什么？何为几何方程、物理方程及虚功方程？弹性矩阵的特点？ 弹性力学变量:外力、应力、应变和位移。 描述弹性体应变分量与位移分量之间的方程称为几何方程；物理方程描述应力分量与应变分量之间的关系；弹性体上外力在虚位移发生过程中所做的虚功与储存在弹性体内的需应变能相等。 弹性矩阵由材料的弹性模量和泊松比确定，与坐标位置无关。 7、何为平面应力问题和平面应变问题？ 平面应力问题：在结构上满足a几何条件：研究对象是等厚度薄板。b载荷条件：作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面无外力作用。 平面应变问题：满足a几何条件：长柱体，即长度方向的尺寸远远大于横截面的尺寸，且横截面沿长度方向不变。b载荷条件：作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力两条件的弹性力学问题。 1、何为结构的离散化？离散化的目的？何为有限元模型？ ①离散化:把连续的结构看成由有限个单元组成的集合体。②目的：建立有限元计算模型③通常把由节点,单元及相应的节点载荷和节点约束构成的模型称为有限元模型2、结构离散化时，划分单元数目的多少以及疏密分布，将直接影响到什么？确定单元数量的原则？通常如何设置节点？
有限元08上机作业
调试书本26到30页程序 开列数组维数： DIMENSION LOC(4,3),CX(6),CY(6),IFIX(6),F(12), 1GK(12,12),STRES(4,3),BAK(4,3,6) 结点集中力输入： DO 10 I=1,ND 10 F(I)=0.0 F(2)=-1.0 数据文件输入： 6,4,12,6,1.0E0,0.0,1.0,0.0,1 3,1,2 5,2,4 2,5,3 6,3,5 0.0,2.0 0.0,1.0 1.0,1.0 0.0,0.0 1.0,0.0 2.0,0.0 1,3,7,8,10,12 数据文件输出： NN NE ND NFIX E ANU T GM NTYPE 6 4 12 6 0.1000E+01 0.000 1.000 0.0000E+00 1 NODE X-LOAD Y-LOAD 1 0.000000E+00 -0.100000E+01 2 0.000000E+00 0.000000E+00 3 0.000000E+00 0.000000E+00 4 0.000000E+00 0.000000E+00 5 0.000000E+00 0.000000E+00 6 0.000000E+00 0.000000E+00 NODE X-DISP Y-DISP 1 -0.879121E-15 -0.325275E+01 2 0.879121E-16 -0.125275E+01 3 -0.879121E-01 -0.373626E+00 4 0.117216E-1 5 -0.835165E-15 5 0.175824E+00 -0.293040E-15 6 0.175824E+00 0.263736E-15 ELEMENT X-STR Y-STR XY-STR 1 -0.879121E-01 -0.200000E+01 0.439560E+00
有限元大作业
风电主轴承有限元分析 XXX 摘要：基于有限元法在接触问题中的应用，对风电主轴承进行非线性分析。以轴承外圈的内表面和内圈的外表面为目标面，以滚子为接触面创建接触对分析滚子的接触应力情况。最大应力值出现在滚子边缘出，对最大承载滚子环向接触应力分析表明，有限元分析结果与理论计算结果相近，验证了利用有限元法分析风电主轴承应力状态的可行性。 关键词：风电主轴承；接触应力；有限元分析 0 引言 随着传统能源的日益枯竭以及环境污染问题愈发严重，风能作为一种清洁的的可再生能源近些年受到越来越多的关注。风力发电技术已广泛运用于世界各地。一些发达国家风力发电产业已得到了迅猛发展，技术日趋成熟，并开始走向产业化规模化发展阶段[1-3]。 风电主轴承是风力发电机重要的组成部分。其结构形式图下图1所示。据统计，如今安装的所有风力发电机中，采用主轴轴承支撑原理的占总数的75-80％[4]，这种支撑是轴承内圈安装在旋转的主轴上，外圈固定在单独的轴承座上，相对于圆锥滚子轴承或圆柱滚子轴承来说，主轴轴承位置处轴产生变形，需要轴承具有一定的调心作用，所以都采用了调心滚子轴承。近年来由于计算机技术的飞速发展，轴承的受力分析计算已经普遍采用有限元分析的方法，能够准确合理地解决轴承复杂的非线性接触问题，为轴承的分析和计算提供了一种新的方法，成为未来的一个发展方向。在机械设备的设计过程中,对受力较大且复杂的零件进行受力分析,校核其整体和局部强度并进行合理的布局设计,是为了防止因应力过大而导致在实际工作中损坏或寿命降低[5]。本文主要运用ANSYS Workbench有限元软件对风电主轴承进行静力学计算，分析轴承内部结构参数对轴承载荷分布和最大接触应力的影响规律。 图1 风电主轴承结构及安装图 1 有限元分析过程 1.1 风电轴承有限元分析基本步骤 不同的物理性质和数学模型的问题,有限元法求解的基本步骤是相同的,只不过 具体公式推导和运算求解不尽相同。有限元分析求解问题的基本计算步骤[6]: 1.问题及求解域定义； 2.求解域离散化； 3.确定状态变量及控制方法； 4.单元推导；
有限元分析及其应用思考题附答案2012
有限元分析及其应用-2010 思考题： 1、有限元法的基本思想是什么？有限元法的基本步骤有那些？其中“离散”的含义是什 么？是如何将无限自由度问题转化为有限自由度问题的？ 答：基本思想：几何离散和分片插值。 基本步骤：结构离散、单元分析和整体分析。 离散的含义：用假想的线或面将连续物体分割成由有限个单元组成的集合，且单元之间仅在节点处连接，单元之间的作用仅由节点传递。当单元趋近无限小，节点无限多，则这种离散结构将趋近于实际的连续结构。 2、有限元法与经典的差分法、里兹法有何区别？ 区别：差分法：均匀离散求解域，差分代替微分，要求规则边界，几何形状复杂精度较低； 里兹法：根据描述问题的微分方程和相应的定解构造等价的泛函表达式，求得近似解； 有限元：基于变分法，采用分片近似进而逼近总体的求解微分方程的数值计算方法。 3、一根单位长度重量为q的悬挂直杆，上端固定，下端受垂直向下的外力P，试 1）建立其受拉伸的微分方程及边界条件； 2）构造其泛函形式； 3）基于有限元基本思想和泛函求极值构造其有限元的计算格式（即最小势能原理）。4、以简单实例为对象，分别按虚功原理和变分原理导出有限元法的基本格式（单元刚度矩 阵）。 5、什么是节点力和节点载荷？两者有何区别？ 答：节点力：单元与单元之间通过节点相互作用 节点载荷：作用于节点上的外载 6、单元刚度矩阵和整体刚度矩阵各有何特点？其中每个矩阵元素的物理意义是什么（按自 由度和节点解释）？ 答：单元刚度矩阵：对称性、奇异性、主对角线恒为正 整体刚度矩阵：对称性、奇异性、主对角线恒为正、稀疏性、带状性。 Kij，表示j节点产生单位位移、其他节点位移为零时作用i节点的力，节点力等于节点位移与单元刚度元素乘积之和。 7、单元的形函数具有什么特点？有哪些性质？ 答：形函数的特点：Ni为x，y的坐标函数，与位移函数有相同的阶次。 形函数Ni在i节点的值为1，而在其他节点上的值为0； 单元内任一点的形函数之和恒等于1； 形函数的值在0~1间变化。 8、描述弹性体的基本变量是什么？基本方程有哪些组成？ 答：基本变量：外力、应力、应变、位移 基本方程：平衡方程、几何方程、物理方程、几何条件 9、何谓应力、应变、位移的概念？应力与强度是什么关系？ 答：应力：lim△Q/△A=S △A→0 应变：物体形状的改变 位移：弹性体内质点位置的变化 10、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系？何谓“强形 式”？何谓“弱形式”，两者有何区别？建立弱形式的关键步骤是什么？
ansys上机作业
实验一坝体的有限元建模及应力应变分析 一、实验目的： 1、掌握ANSYS软件基本的几何形体构造方法、网格划分方法、边界条件施加方法及各 种载荷施加方法。 2、熟悉有限元建模、求解及结果分析步骤和方法。 3、能利用ANSYS软件对结构进行有限元分析。 二、实验设备： 微机，ANSYS软件。 三、实验内容： 计算分析模型如图所示，分析坝体的应力、应变。 四、实验步骤： 1 进入ANSYS 程序→ANSYS →change the working directory into yours →input Initial jobname: dam 2设置计算类型 ANSYS Main Menu: Preferences →select Structural → OK 3选择单元类型
ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK (back to Element Types window)→Options… →select K3: Plane Strain →OK→Close (the Element Type window) 4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear→Elastic→Isotropic→input EX:2.1e11, PRXY:0.3→ OK 5生成几何模型 生成特征点 ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input:1(0,0),2(1,0),3(1,5),4(0.45,5)→OK 生成坝体截面 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS→依次连接四个特征点，1(0,0),2(1,0),3(1,5),4(0.45,5) →OK 生成坝体截面如图一 图一 6 网格划分 ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条横边:OK→input NDIV: 15 →Apply→依次拾取两条纵边:OK →input NDIV: 40 →OK →(back to the mesh tool window)Mesh: Areas, Shape: Quad, Mapped →