有限元作业
有限元大作业
250250试题 5：图示为带方孔（边长为 80mm ）的悬臂梁，其上受部分均布载荷（P=10KN/m ）作用，试采用一种平面单元，对图示两种结构进行有限元分析，并就方孔的布置（即方位）进行分析比较，如将方孔设计为圆孔，结果有何变化？（板厚为 1mm ，材料为钢）。
3001KN9003001KN图6-1一、几何建模与分析由图6-1及问题描述可知，板的长宽尺寸远远大于厚度，研究结构为一很薄的等厚度薄板，满足平面应力的几何条件；作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布，而在两板面上无外力作用，满足平面应力的载荷条件。
故该问题属于平面应力问题，薄板所受的载荷为面载荷，分布情况及方向如图6-1所示，建立几何模型，进行求解。
薄板的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3。
二、有限元分析及其计算结果选取PLANE182作为分析的单元，来分析薄板的位移和应力，由于此问题是平面应力问题，并在K3选择str w/thk ，设置THK 为1。
1）方孔竖直制，划分方式采用自由方式，划分后网格的模型如图6-2所示。
计算得到的位移和应力分布如图6-3所示。
图6-2 方孔竖直的网格划分图6-3 位移及应力分布云图2）方孔正直制，划分方式采用自由方式，划分后网格的模型如图6-4所示。
计算得到的位移和应力分布如图6-5所示。
图6-4 方孔正直的网格划分图6-5 位移及应力分布云图3）圆孔按图6-1所示模型进行建模。
并用PLANE182单元进行划分网格，网格大小采用全局网格控制，划分方式采用自由方式，划分后网格的模型如图6-6所示。
计算得到的位移和应力分布如图6-7所示。
图6-4 方孔正直的网格划分图6-5 位移及应力分布云图根据以上的模型分析的位移和应力图，可以得出方孔竖直、方孔正直、圆孔的最大最小位移应力的分布如表6-1所示。
三、比较与分析1）方孔竖直与方孔正直的比较，发现方孔正直的位移变形较小，应力相差不大2）圆孔与方孔比较，发现圆孔的位移变性最小，应力也最小，故可以得出圆孔的布置结构对整体布置的效果最好。
《有限元分析》课程作业
《有限元分析》课程作业任课教师：徐亚兰学生姓名：陈新杰学号：班级：1304012时间：2016-01-05一、问题描述及分析问题：如图1所示，有一矩形平板，在右侧受到P=10KN/m 的分布力，材料常数为：弹性模量Pa E 7101⨯=；泊松比3/1=μ；板的厚度为t=；试按平面应力问题利用三角形与矩形单元分别计算各个节点位移及支座反力。
图1 平面矩形结构的有限元分析分析：使用两种方案：一、基于3节点三角形单元的有限元建模，将矩形划分为两个3节点三角形单元；二、基于4节点矩形单元的有限元建模，使用一个4节点矩形单元。
利用MATLAB 软件计算出各要求量，再将两种方案的计算结果进行比较、分析、得出结论。
二、有限元建模及分析1、基于3节点三角形单元的有限元建模及分析 （1）结构的离散化与编号如图2所示，将平面矩形结构分为两个3节点三角形单P=10KN/m1m1m元。
单元①三个节点的编号为1，2，4，单元②三个节点的编号为3，4，2，各个节点的位置坐标为(),,1,2,3,4i i x y i =，各个节点的位移（分别沿x 方向和y 方向）为(),,1,2,3,4i i u v i =。
图2 方案一：使用两个3节点三角形单元（2）各单元的刚度矩阵及刚度方程 a.单元的几何和节点描述单元①有6个节点位移自由度（DOF ）。
将所有节点上的位移组成一个列阵，记作(1)q ；同样，将所有节点上的各个力也组成一个列阵，记作(1)F ，则有(1)112244,,,,,)q u v u v u v =（(1)112244(,,,,,)x y x y x y F F F F F F F =同理，对于单元②，有(2)334422,,,,,)q u v u v u v =（1234X y ①②(2)334422(,,,,,)x y x y x y F F F F F F F =b.单元的位移场描述对于单元①，设位移函数012012(,)(,)u x y a a x a y v x y b b x b y ⎫=++⎪⎬=++⎪⎭(1-1)由节点条件，在,i i x x y y ==处，有(,)(,)i i i i i i u x y u v x y v =⎫⎬=⎭1,2,4i = (1-2) 将式(1-1)代入节点条件式(1-2)中，可求出式（1-1）中待定系数，即011122211223444411()22u x y a u x y a u a u a u AAu x y ==++ (1-3) 11122112234441111()221u y a u y b u b u b u AAu y ==++ (1-4) 21122112234441111()221x u a x u c u c u c u AAx u ==++ (1-5) 01122341()2b a v a v a v A =++(1-6) 11122341()2b b v b v b v A =++(1-7) 21122341()2b c v c v c v A =++(1-8)在式(1-3)~式(1-8)中1122123441111()221x y A x y a a a x y ==++ (1-9)2212442442124421244(1,2,3)1111x y a x y x y x y y b y y y x c x x x ⎫==-⎪⎪⎪⎪=-=-⎬⎪⎪⎪==-+⎪⎭ (1-10) 上式中的符号(1,2,3)表示下标轮换，如12,23,31→→→同时更换。
有限元大作业
1．推导有限元计算格式，理解有限元原理：建立图示受拉直杆在自重（设单位长度重度为q ，截面积为A ）和外力P 作用下的拉伸问题的微分方程，并分别利用不同的原理（变分求极值（最小势能或虚功原理）、加权残值法）推导有限元计算格式（取两个单元）。
手工求出端点的位移（自己给定参数值）。
设杆长为L ，截面面积为A(x)，弹性模数为E,单位长重量q ，受拉杆x 处的位移为u(x)。
取微元dx 的力平衡，建立受拉杆位移所满足的微分方程()du x dx ε=，()du x E E dxσε== dx 上下截面内力与微元自重相等得()*()()*()A x dx x dx A x x dx qdx σσ++-+=-(()())dA x x q dxσ∴=- (())d duEA x q dx dx=- 0x L << ()0u x = 0x =()duEA x p dx= x L = 得解析解：2()2q x P u Lx x EA EA=-+将其分为两个单元，节点为1,2,3，得22382qL PL u EA EA=+232qL PL u EA EA=+有限元法：1)位移函数01u α= 2111u u l α-=得1211(1)x x u u u l l =-+ 令11(1)x N l =-21x N l = 11122122u u N u N u N N u⎧⎫⎪⎪⎡⎤=+=⎨⎬⎣⎦⎪⎪⎩⎭{}1u N d ⎡⎤=⎣⎦ 2)应变、应力表达{}{}111211du dN d d dx dx l l ε⎡⎤⎡⎤===-⎢⎥⎣⎦⎣⎦{}1B d ε⎡⎤=⎣⎦ {}1E E B d σε⎡⎤==⎣⎦ {}1S d σ⎡⎤=⎣⎦3)势能表示{}{}(){}{}(){}{}{}{}{}1111''112211''121112210111111111111111121221222T V ll T T T T T U W D dV F u F u qdx u u d B E d Adx F u F u ql EA EA ql l l d d d F d EA EA ql l l εε⎡⎤=-=-+-⎣⎦+⎡⎤=-+-⎣⎦⎡⎤⎡⎤-⎢⎥⎢⎥⎢⎥⎢⎥=--⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦⎢⎥⎣⎦∏⎰⎰⎰4)单元平衡方程 a)最小势能原理110u ∂=∂∏120u ∂=∂∏111111212112112ql F u AE l u ql F ⎧⎫-⎪⎪⎧⎫⎡⎤-⎪⎪⎪⎪=⎨⎬⎨⎬⎢⎥-⎣⎦⎪⎪⎪⎪⎩⎭+⎪⎪⎩⎭b)虚位移原理{}(){}(){}TeTdd F qdx d δδδεσΩ+=Ω⎰⎰{}{}1B d σεδ⎡⎤=⎣⎦ {}1E E B d σεδ⎡⎤==⎣⎦{}(){}{}(){}111111TTT l d F d B E B d Adxδδ⎡⎤=⎣⎦⎰ 由虚位移任意性得，{}{}1111T lF B E B Adxd ⎡⎤=⎣⎦⎰ 积分得111111212112112ql F u AE l u ql F ⎧⎫-⎪⎪⎧⎫⎡⎤-⎪⎪⎪⎪=⎨⎬⎨⎬⎢⎥-⎣⎦⎪⎪⎪⎪⎩⎭+⎪⎪⎩⎭ 记为{}{}111k d F ⎡⎤=⎣⎦ 同理222212323112112ql F u AE l u ql F ⎧⎫-⎪⎪⎧⎫⎡⎤-⎪⎪⎪⎪=⎨⎬⎨⎬⎢⎥-⎣⎦⎪⎪⎪⎪⎩⎭+⎪⎪⎩⎭{}{}222k d F ⎡⎤=⎣⎦ {}{}ei i eF R =∑ 12220F F += 23F P =11111112211223222022202EAEAql F l l u ql ql EA EA EA EA u l l l l u ql EAEA P l l ⎡⎤⎧⎫-⎢⎥+⎪⎪⎢⎥⎧⎫⎪⎪⎢⎥⎪⎪⎪⎪⎢-+-⎥=+⎨⎬⎨⎬⎢⎥⎪⎪⎪⎪⎢⎥⎩⎭⎪⎪+⎢⎥⎪⎪--⎢⎥⎩⎭⎣⎦可得：22382qL PLu EA EA=+232qL PL u EA EA=+与解析解结果一致。
有限元分析大作业报告
有限元分析大作业报告一、引言有限元分析是工程领域中常用的数值模拟方法，通过将连续的物理问题离散为有限个子区域，然后利用数学方法求解，最终得到数值解。
有限元分析的快速发展和广泛应用，为工程领域提供了一种强大的工具。
本报告将介绍在大作业中所进行的有限元分析工作及结果。
二、有限元模型建立本次大作业的研究对象是工程结构的应力分析。
首先，通过对结构进行几何建模，确定了结构的尺寸和形状。
然后，将结构离散为有限个单元，每个单元又可以看作一个小的子区域。
接下来，为了求解结构的应力分布，需要为每个单元确定适当的单元类型和单元属性。
最后，根据结构的边界条件，建立整个有限元模型。
三、材料属性和加载条件在建立有限元模型的过程中，需要为材料和加载条件确定适当的参数。
本次大作业中，通过实验获得了结构材料的弹性模量、泊松比等参数，并将其输入到有限元模型中。
对于加载条件，我们选取了其中一种常见的加载方式，并将其施加到有限元模型中。
四、数值计算和结果分析为了求解结构的应力分布，需要进行数值计算。
在本次大作业中，我们选用了一种常见的有限元求解器进行计算。
通过输入模型的几何形状、材料属性和加载条件，求解器可以根据有限元方法进行计算，并得到结构的应力分布。
最后，我们通过对计算结果进行分析，得出了结论。
五、结果讨论和改进方法根据计算结果，我们可以对结构的应力分布进行分析和讨论。
根据分析结果，我们可以得出结论是否满足设计要求以及结构的强度情况。
同时，根据分析结果，我们还可以提出改进方法，针对结构的特点和问题进行相应的优化设计。
六、结论通过对工程结构进行有限元分析，我们得到了结构的应力分布，并根据分析结果进行了讨论和改进方法的提出。
有限元分析为工程领域提供了一种有效的数值模拟方法，可以帮助工程师进行结构设计和分析工作，提高设计效率和设计质量。
【1】XXX，XXXX。
【2】XXX，XXXX。
以上是本次大作业的有限元分析报告，总结了在建立有限元模型、确定材料属性和加载条件、数值计算和结果分析等方面的工作，并对计算结果进行讨论和改进方法的提出。
有限元作业整理版
一、写出下图所示的三结点三角形单元的插值函数(形函数)Ni,Nj,Nm及插值函数矩阵[N],应变矩阵[B]。
二、如下图所示的三结点三角形单元,厚度为t,弹性模量是E,泊松比μ.试求:插值函数(形函数)矩阵[N],应变矩阵[B],应力矩阵[S],单元刚度矩阵[K]e.
三、下图所示的三结点三角形单元在jm边作用有线性分布的面载荷（x方向）,试求等效结点载荷向量.
四、如下图所示,一正方形平板,厚度为t,边长为a,弹性模量E,泊松比μ.划分为两个三角形单元.求:1,3点的位移.
五、下图所示矩形板，分成四个三结点三角形单元.要求:
(1)写出由单元刚度矩阵组装总体刚度矩阵的表达式;
(2)如1234就是一个矩形单元,试求此单元的单元刚度矩阵.。
有限元分析大作业
有限元大作业一题目要求：图1所示为一悬臂梁，在端部承受载荷，材料弹性模量为E，泊松比为1/3，悬臂梁的厚度（板厚）为t，若该粱被划分为两个单元，单元和节点编号如图所示，试按平面应力问题计算各个节点位移计支反力。
一、单元划分1.计算简图及单元划分如下所示：2.进行节点及单元编号节点i j m单元① 2 3 4② 3 2 13.节点坐标值节点号1 2 3 4坐标值X 2 2 0 0Y 1 0 1 0二、计算单元刚度矩阵1、计算每个单元面积△以及i b ，i c （m j i i ,,=） ①②单元的面积相等，即12121=⨯⨯=∆ 单元①的i b ，i c⎩⎨⎧=--==-=0)(1m j i m j i y x c y y b ⎩⎨⎧=--==-=2)(0i m ji m j x x c y y b ⎩⎨⎧-=--=-=-=2)(1j i mj i m y x c y y b 对平面应力问题，其表达式为[]⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-+-+-+-+∆-=s r s r sr s r s r s r s r s r b b uc c cb u b uc b c u c ub c c u b b u Et Krs 21212121)1(42 然后对单元①求解单元刚度子矩阵2==i r 2==i s []⎥⎦⎤⎢⎣⎡=3/1001329)1(22Et K 2==i r 3==j s []⎥⎦⎤⎢⎣⎡=03/23/20329)1(23Et K2==i r 4==m s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)1(24Et K 3==j r 3==j s []⎥⎦⎤⎢⎣⎡=4003/4329)1(33Et K 3==j r 2==i s []⎥⎦⎤⎢⎣⎡=03/23/20329)1(32Et K 3==j r 4==m s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)1(34Et K 4==m r 4==m s []⎥⎦⎤⎢⎣⎡=3/133/43/43/7329)1(44Et K 4==m r 2==i s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)1(42Et K 4==m r 3==j s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)1(43Et K由子矩阵[]e rs K 合成单元刚度矩阵[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=3/133/443/23/13/23/43/73/23/43/2143/24003/23/23/403/43/203/13/203/23/103/213/2001329)1(Et K将单元①的单元刚度矩阵补零升阶变为单元刚度矩阵，其在总体刚度矩阵中的位置为：节点号→单元②的i b ，i c⎩⎨⎧=--=-=-=0)(1m j im j i y x c y y b ⎩⎨⎧-=--==-=2)(0i m ji m j x x c y y b ⎩⎨⎧=--==-=2)(1j i mj i m y x c y y b 然后对单元 求解单元刚度子矩阵：3==i r 3==i s []⎥⎦⎤⎢⎣⎡=3/1001329)2(33Et K 3==i r 2==j s []⎥⎦⎤⎢⎣⎡=03/23/20329)2(32Et K 3==i r 1==m s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)2(31Et K 1 2 3 412[])1(22K[])1(23K[])1(24K3[])1(32K[])1(33K[])1(34K4[])1(42K[])1(43K[])1(44K2==j r 2==j s []⎥⎦⎤⎢⎣⎡=4003/4329)2(22Et K 2==j r 3==i s []⎥⎦⎤⎢⎣⎡=03/23/20329)2(23Et K 2==j r 1==m s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)2(21Et K 1==m r 1==m s []⎥⎦⎤⎢⎣⎡=3/133/43/43/7329)2(11Et K 1==m r 3==i s []⎥⎦⎤⎢⎣⎡----=3/13/23/21329)2(13Et K 1==m r 2==j s []⎥⎦⎤⎢⎣⎡----=43/23/23/4329)2(12Et K 由子矩阵[]e rs K 合成单元刚度矩阵[]⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡----------------=3/133/443/23/13/23/43/73/23/43/2143/24003/23/23/403/43/203/13/203/23/103/213/2001329)2(Et K将单元②的单元刚度矩阵补零升阶变为单元贡献矩阵，其在总体刚度矩阵中的位置为：节点号→1 2 3 41 [])2(11K[])2(12K[])2(13K2 [])2(21K[])2(22K[])2(23K3 [])2(31K [])2(32K [])2(33K 4三、计算总体刚度矩阵总体刚度矩阵是由各单元的贡献矩阵迭加而成)2()1(][][][][K K K K e +==∑四、进行节点约束处理根据节点约束情况，在总刚矩阵中可采用划行划列处理约束的方法，由题目易知，节点3和4的已知水平位移和垂直位移都为零，划去其相对应的行和列，则总刚矩阵由8阶变为4阶，矩阵如下：⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------2/02/03/13043/203/73/23/443/23/133/43/23/43/43/73292211p p v u v u Et329][Et K =1 2 3 413/133/43/43/743/23/23/4----3/13/23/21----000243/23/23/4----3/13003/73/43/403/13/23/21----33/13/23/21----3/43/403/13003/743/23/23/4----40003/13/23/21----43/23/23/4----3/133/43/43/7化简⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------Et p Et p v u v u 3/1603/160130122072412213424472211 五、求解线性方程组方法：采用LU 分解法 1.求解矩阵[]U 各元素⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------75/10775/640075/6475/353007/767/27/7502447~7/877/87/7607/87/337/207/767/27/7502447~13012207241221342447⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----353/44900075/6475/353007/767/27/7502447~ 得到的[]U 矩阵如下:[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=353/44900075/6475/353007/767/27/7502447U 2.求解矩阵[]L 各元素⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----13012207241221342447353/44900075/6475/353007/767/27/75024471353/6475/767/20175/27/40017/40001 得到的[]L 矩阵如下：[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--------=13012207241221342447L3.进行求解⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧---=⇒⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧--=Et p Et p Et p y Et p Et p Ly 79425/850800225/323/1603/1603/160⎪⎪⎭⎪⎪⎬⎫⎪⎪⎩⎪⎪⎨⎧---=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----⇒=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡Et p Et p Et p v u v u y v u v u U 79425/850800225/323/160353/44900075/6475/353007/7675/27/750244722112211 解得Et p v /422.82-= Et p u /497.12-= Et p v /028.91-= Et p u /897.11=于是求得各节点的位移为：⎩⎨⎧-==Etp v Etp u /028.9/897.111 ⎩⎨⎧-=-=Etp v Etp u /422.8/497.122 ⎩⎨⎧==033v u ⎩⎨⎧==044v u 六、求解相应的支反力（运用静力学的平衡方程进行求解）3号节点和4号节点的支反力如下图所示：。
有限元作业
Bierenzuode,kanbudong作业1： 有一个等截面两节点二力杆，杆长为L ，截面积为A ，材料弹性模量为E 。
每个节点只考虑一个水平位移，对于图 (a)、(b) 所示的坐标系统和位移插值函数，分别求相应的[B]矩阵和单元刚度矩阵[K]。
解：（a ）、212()u x x αα=+，由边界条件确定常数1α、2α：当0x =时，1i u α=；当x l =时，212j u l αα=+，可得2222()1i j x x u x u u l l ⎛⎫=-+ ⎪⎝⎭因每个节点只考虑一个水平位移故以矩阵形式表示的单元位移函数为：{}{}{}2211122222()1u u x x f x u N N u u ll ⎧⎫⎧⎫⎧⎫==-=⎨⎬⎨⎬⎨⎬⎩⎭⎩⎭⎩⎭单元的几何矩阵：[]''122222x x B N N l l ⎡⎤⎡⎤==-⎣⎦⎢⎥⎣⎦{}{}{}12x u E E B u σε⎧⎫=⨯=⎨⎬⎩⎭，即[][]D E = 对于矩形截面梁单元，积分：yzd dA =⎰⎰为单元横截面面积。
梁单元刚度矩阵[]0leT EA K B B dx ⎡⎤⎡⎤⎡⎤⎣⎦⎣⎦⎣⎦=⎰222202222l x x x l EA x ll l dx ⎡⎤⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦-=-⎰44334433EAEA l l EAEA ll --⎡⎤⎢⎥⎢⎥⎣⎦=（b ）、212()u x x x αα=+，由边界条件确定常数1α、2α：当2l x =-时，21224i a l a l u =-+；当2l x =时，21224j a l a l u =+可得222222()i j x lx x lxu x u u l l -+=+ 因每个节点只考虑一个水平位移故以矩阵形式表示的单元位移函数为：{}{}{}221112222222()u u x lxx lx f x u N N u u ll ⎧⎫⎧⎫⎧⎫-+===⎨⎬⎨⎬⎨⎬⎩⎭⎩⎭⎩⎭单元的几何矩阵：[]''122244x l x l B N N l l -+⎡⎤⎡⎤==⎣⎦⎢⎥⎣⎦{}{}{}12x u E E B u σε⎧⎫=⨯=⎨⎬⎩⎭，即[][]D E =对于矩形截面梁单元，积分：yzd dA =⎰⎰为单元横截面面积。
岩土工程有限元大作业
岩土工程有限元大作业题目：平面问题3结点三角形有限元法的Matlab程序实现专业：建筑与土木工程（岩土方向）姓名：学号：_________2014年12 月31 日1.概述对于二维平面问题，运用有限单元法求解时，单元的类型有3结点三角形单元、4结点矩形单元以及高精度的三角形单元（6结点、10结点）等。
三角形单元对复杂边界有较强的适应能力，因此很容易将一个二维域离散成有限个三角形单元，在边界上以若干段直线近似原来的曲线边界，随着单元增多，这种拟合将越精确。
本文即基于相对比较简单且运用广泛的3结点三角形单元的有限元法编制Matlab 程序，解决一个悬臂梁端点承受竖向集中荷载作用时的位移计算问题，并绘出变形前后的图形作对比。
在分析过程中，将此问题简化为平面应力问题，取单位厚度进行分析。
以下一共分为两部分介绍，第一部分是理论部分，介绍三角形有限元法的基本求解过程；第二部分是程序部分，介绍具体的Matlab 程序实现过程。
2.三角形有限元法的基本求解过程 2.1 假设单元的位移函数如图2.1所示为一端部承受竖向集中荷载作用的悬臂梁，将此问题简化为二维平面问题，梁取单位厚度，固定端简化为两个不动铰支座，已知悬臂梁的长度为8米，高度为1米，材料为钢材，它的弹性模量取2.06e11Pa ，泊松比取0.25，右端处作用的竖向荷载为100kN 。
二维域Ω被离散成如图所示的九个三角形单元，后续编程中的结点号和单元号都以此为基础，任取一三角形单元，设其结点编码为m j i ,,，以逆时针编码为正向(见图2.2)。
该三角形单元在x 、y 方向的位移函数表达式为：m m j j i i u N u N u N u ++= (2.1) m m j j i i v N v N v N v ++=(2.2)其中()m j i i N i ,, =称为单元的插值函数或形函数()y c x b a AN i i i i ++=21(2.3)图2.1 任意区域的三角形单元离散 j (x j , y j )xyi (x i , y i )m (x m , y m )图2.2 3结点三角形单元),,(m j i xx c y y b y x y x a m j im j i j m m j i ⎪⎩⎪⎨⎧+-=-=-= (2.4)上式),,(m j i 表示下标轮换，即i m m j j i →→→,,。
有限元课程大作业
金属坯料挤压过程有限元分析一、前言：金属挤压是将放在挤压模具内的金属锭坯从一端施加外力，强迫其从特定的模孔中流出，获得所需要的断面形状和尺寸的制品的技术。
冷挤压时由于材料是在冷态下成形，而且变形量一般都很大，挤压过程中作用在模具上的单位压力很大,此时模具有开裂破坏的可能，对压力机也构成威胁，金属坯料在通过模具过程中，坯料与模具之间产生相当大的应力，这就要求模具需要有相当大的强度、硬度、以及耐磨性，因此冷挤压时要进行挤压力的计算。
挤压力的计算是模具设计的重要依据，也是选择挤压设备的依据。
模具角度、接触表面的摩擦系数、坯料变形量都会影响应力变化，在保证加工要求的前提下，应当通过适当方式降低坯料及模具之间的应力。
通过有限元分析，得出应力分布图，分析变形区域、死区，对模具进行优化改进。
二、有限元介绍：ANSYS概述ANSYS软件是融结构、热、流体、电磁、声学于一体的大型通用有限元软件，可广泛地用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、生物医学、水利、日用家电等一般工业及科学研究。
该软件提供了不断改进的功能清单，具体包括：结构高度非线性分析、电磁分析、计算流体力学分析、设计优化、接触分析、自适应网格划分及利用ANSYS参数设计语言扩展宏命令功能。
ANSYS软件功能强大，主要特点有：实现多场及多场耦合分析；实现前后处理、求解及多场分析统一数据库的一体化；具有多物理场优化功能；强大的非线性分析功能;多种求解器分别适用于不同的问题及不同的硬件设备；支持异种、异构平台的网络浮动，在异种、异构平台上用户界面统一、数据文件全部兼容；强大的并行计算功能支持分布式并行及共享内存式并行；多种自动网格划分技术；良好的用户开发环境。
ANSYS不仅支持用户直接创建模型，也支持与其他CAD软件进行图形传递，其支持的图形传递有：SAT、Parasolid、STEP。
相应地，可以进行接口的常用CAD 软件有：Unigraphics、Pro/Engineer、I-Deas、Catia、CADDS、SolidEdge、SolidWorks等。
中科大工程中的有限元作业答案（最新）
、最小势能原理法
2
1 1 12 1 2 22 1 3 32 1 1 2 2
P = 2 k d + 2 k d + 2 k d - Fq - F q
2的位移。
d1 , d 2 , d 3分别代表弹簧 (1)，
( 2 )，( 3)的伸长量，q1和q2表示节点1，
= êêk21(1) k22(1) 0úú + êê0 k11(2) k12(2) úú
êë 0 0 0úû êë0 k21(2) k22(2) úû
k12(1)
0ù
ék11(1)
ê
(1)
(1)
(2)
= êk21 k22 + k11 k12(2) úú
êë 0
k21(2)
k22(2) úû
EA
é E1A1
ë
û
2
2
sinq cos
-cos q
-sinq cos
q
qù
é cos q
2
2
ê
AE sinq cosq
sin q
-sinq2cosq
-sin q úú
2
K=
ê
L ê -cos q
cos q
sinq cos
-sinq cos
q
qú
2
2
ê-sin cos
sinq cosq
cos q úû
ë q q -sin q
3
Þ RB = ql
8
3
R Bl
ql
, y BR =
= 3EI
8EI
4
3
ql
有限元分析题及大作业题答案
有限元分析及应用作业报告试题10一、问题描述确定图示扳手中的应力, E=210Gpa，μ=0.3, 假设厚度为10mm；并讨论采用何种处理可降低最大应力或改善应力分布。
图1为扳手的基本形状和基本尺寸图二、数学建模与分析由图1及问题描述可知，板手的长宽尺寸远远大于厚度，研究结构为一很薄的等厚度薄板，满足平面应力的几何条件；作用于薄板上的载荷平行于板平面且作用在沿厚度方向均匀分布在办手柄的左边缘线，而在两板面上无外力作用，满足平面应力的载荷条件。
故该问题属于平面应力问题，薄板所受的载荷为面载荷，分布情况及方向如图1所示，建立几何模型，并进行求解。
薄板的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3三、有限元建模1、单元选择：选取三节点常应变单元来计算分析薄板扳手的位移和应力。
由于此问题为平面应力问题，：三节点常应变单元选择的类型是PLANE42（Quad 4node42），该单元属于是四节点单元类型，在网格划分时可以对节点数目控制使其蜕化为三节点单元。
2、定义材料参数：ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK3、生成几何模型：a.创建关键点点：ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入16个点的坐标→OKb、将这16个关键点有直线依次连起来，成为线性模型4、生成实体模型：ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →连接特征点→生成两个area→Operate→Subtract→拾取整个扳手区域→OK→生成扳手模型5、结点布置及规模6、网格划分方案ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →Mesh: Areas, Shape: Tri,Free →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)7、载荷及边界条件处理8、求解控制A、模型施加约束给模型施加x方向约束ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On Lines →拾取模型左部的竖直边→OKB、给模型施加载荷ANSYS Main Menu: Solution →Define Loads →Apply →Structural →force→on keypoints→拾取上面左端关键点→700N/mm→okC、分析计算：ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load Step window) →OK6)结果显示：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 →displacement vector sum，von mises stress→OK四、计算结果及结果分析1、三节点常应变单元1）三节点单元的网格划分图2 常应变三节点单元的网格划分平面图图3 常应变三节点单元的网格划分立体图2）三节点单元的约束受载情况图4 常应变三节点单元的约束受载图3)三节点单元的位移分析图5 常应变三节点单元的位移分布图4)三节点单元的应力分析图6 常应变三节点单元的应力分布图2、六节点三角形单元1)六节点三角形单元网格划分图7 六节点三角形单元网格划分图2)六节点三角形单元约束和受载情况分析图8 六节点三角形单元约束受载图3）六节点三角形单元位移分析图9 六节点三角形单元的变形分布图4) 六节点三角形单元的应力分析图9 六节点三角形单元的应力分布图图10 六节点三角形单元的局部应力分布图根据以上位移和应力图，可以得出常应变三节点单元和六节点三角形单元的最小最大位移应力如表1-1所示。
有限元例题及答案
例 8-1：E ，A ，L ，s σ 杆I 弹塑性； 杆II 弹性。
求s AF σ3=下2点位移。
解：（1）理论解在荷载s A F σ3=作用下，杆I 屈服而有内力（拉力）S A N σ=1，杆II 内力（压力）为s II A N σ2=，中点2位移δ取决于杆II 的变形，即*===∆=δσσδ22)2(EL AE L A l S S II式中E Ls σδ=*（屈服位移）（2）直接迭代法杆I 和杆II 的刚度分别为⎩⎨⎧=**≤〉)()(δδδδδσL EAAI S k L EA k II =①迭I 迭代步迭代从*=δδ0开始，这时有L EAk k K II I 20=+=*-====δσσδ5.15.123101EL L EA A F K S S②第2迭代步杆I 进入塑性，有L EA A k s I 67.01==δσ杆Ⅱ完全弹性，刚度不变。
因此，总刚为L EAk k K II I 67.11=+=*-====δσσδ8.18.167.13112E L LEA A F k S s 整个迭代过程见表8-1。
表8-1 直接迭代法各次迭代结果(3)切线刚度法杆Ⅰ和杆Ⅱ的切线刚度分别为⎩⎨⎧=**≤〉)()(0δδδδLEAI k L EA k II =①第1迭代步初始状态时，00=δ，杆Ⅰ，Ⅱ中应力、应变均匀为零。
总刚为：L EAk k K T TI T 21=+=由F K T -=δψ，得S A σψ30-=由n Tn n K ψδ1--=∆得，*=--=∆δσδ5.1)3(10S A L由式n n n δδδ∆+=+1得，s δδ5.11=杆中应力：S SI σσσσ5.111-==杆中内力：S SI A N A N σσ5.111-==②第2迭代步由于杆I 已进入塑性，杆Ⅱ仍处弹性，总刚：L EAk k K TIITI T =+=2由F K T -=δψ，得S S S A A A σσσψ5.035.21-=-=由n Tn n K ψδ1--=∆得，*=--=∆δσδ5.0)5.0(11S A LEA由式n n n δδδ∆+=+1得，*=∆+=σδδδ0.2112杆中应力：S II SI A N A N σσ0.222-==检验F K T -=δψ，有030.32=-=S S A A σσψ迭代平衡。
有限元作业-内容十分详细
第一题1.题目概况矩形板尺寸如下图1，板厚为5mm。
材料弹性模量为52E=⨯泊松210N/mm μ。
选择以下一种工况讨论：比27=.0本次分析选取的是1和2两种情况。
由于1,2种情况十分类似，所以这里主要分析第一种情况的步骤。
2．模型建立2.1 单元选择及其分析本次问题中的矩形薄板的应力分析属于平面应力分析，是结构静力学问题。
定义单元类型为二维四边形单元。
（1）图2，首先在Preference菜单中定义分析类型为Structural。
图2（2）在Preprocessor/Elementtype/Add/Edit/Delet中定义单元属性为二维四边形单元，如图3所示。
图3相应的选项设置如图4所示：（3）定义材料特性：EX=200000，PRXY=0.27。
如图5，图6所示：图5（4）定义平板厚度为5，如图7所示：图72.2 模型建立及网格划分（1）图8在XY面内建立矩形，输入如图中所示数据，完成后创建的模型如图8所示。
图8图9（2）划分网格。
点击Preprocessor/Meshing/Meshingtool后，设置网格的属性。
定义四边形网格的边长为5如图10所示,点击Mesh后，开始网格划分如图10所示。
图10图112.3 载荷处理（1）定义分析类型。
点击Solution/Analysis Type/New Analysis,设为static，即结构静态分析。
如图12所示。
（2）施加约束。
点击Solution/Define Loads/Apply/Structural/Displacement/on Nodes后，点击c,d两个节点后，设置选项如图13所示，约束后的模型如图14所示。
图14（3）施加载荷。
点击Solution/Define Loads/Operate/Apply/Structual/Pressure/On lines,选择a,b边后出现选项卡后，点击设置如图15所示参数。
设置完成后载荷如图16所示。
有限元大作业
1、某齿轮弯曲应力分析
某直齿圆柱齿轮相关参数为18，模数4.1，压力角22.5°，变位系数为0.36168。
齿轮材料为40Cr ，假设其径向力2364N ，切向力为6400N ，试求其弯曲应力分布情况。
2、动力学作业：
某发动机外壳模型Motor_Cover_5.x_t ，材料为钢材，壁厚为0.05in ，材料参数为：弹性模量E=2E11Pa ，泊松比u=0.28
，密度。
边界条件说明如图所示。
试计算其固有频率。
要求：
1、A4纸，边距均为2.5mm ；标题4号字，其它为小4号；文内图形宽不大于8mm 、高不大于6mm ，有图名和标号；表格用3线格形式，有编号和名称。
2、写成文章形式，包括问题、分析过程、结果及分析、结论。
给出必要的几何建模（或外部模型导入）、网格划分、求解及后处理等关键步骤。
必须有必要的结果分析讨论。
3、提交完整分析报告：封面，目录，正文。
正文为小四宋体，封面参考附只约束法
向位移
固定约束
件。
有限元作业
1 有限元的基础理论有限单元法是20世纪50年代以来随着电子计算机的广泛应用而发展起来的有一种数值解法。
有限元分析（FEA，Finite Element Analysis）的基本概念是用较简单的问题代替复杂问题有限元分析后再求解。
它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。
这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。
有限元法的基本思路和基本原则以结构力学中的位移法为基础，把复杂的结构或连续体看成有限个单元的组合，各单元彼此在节点处连接而组成整体。
把连续体分成有限个单元和节点，称为离散化。
先对单元进行特性分析，然后根据各节点处的平衡和协调条件建立方程，综合后作整体分析。
这样一分一合，先离散再综合的过程，就是把复杂结构或连续体的计算问题转化为简单单元的分析与综合的问题。
因此，一般的有限元解法包括三个主要步骤：离散化、单元分析、整体分析。
(1)离散化一个复杂的弹性体可以看作由无限个质点组成的连续体。
为了进行解算，可以将此弹性体简化为有限个单元组成的集合体，这些单元只在有限个节点上铰接，因此，这集合体只具有有限个自由度，这就为解算提供了可能。
有无限个质点的连续体转化为有限个单元的集合体，就称为离散化。
(2)单元分析单元分析首先要进行单元划分。
在工程结构中，一般采用四种类型的基本单元，即标量单元、线单元（杆、梁单元）、面单元和体单元。
四种基本单元的若干例子及各单元节点自由度（节点位移）表示在图（1-1）中。
而单元划分一般注意下面几点：一、从有限元本身来看，单元划分的越细，节点布置得越多，计算的结果越精确。
但计算时间和计算费用的增加。
所以在划分单元时对应兼顾这两个方面。
二、在边界比较曲折，应力比较集中，应力变化较大的地方，单元应划分的细点，而在应力变化平缓处单元划分的大些。
张年梅有限元方法作业2016
张年梅有限元方法作业2016有限元方法是一种常用的工程数值分析方法，可以用于求解各种工程问题的数值解。
张年梅是一位工程师，他在2016年的有限元方法作业中，面临着以下几个问题：1. 简要介绍有限元方法的原理和应用领域；2. 描述有限元方法的基本步骤和流程；3. 以张年梅研究的某个工程问题为例，说明如何利用有限元方法进行数值分析；4. 讨论有限元方法的优点和局限性。
有限元方法是一种将复杂问题离散化为一系列简单子问题的方法，然后通过数值计算的方式求解这些子问题，最终得到原问题的数值解。
它的应用领域非常广泛，包括结构力学、流体力学、热传导、电磁场等等。
有限元方法的基本步骤和流程可以概括为以下几个步骤：1. 建立几何模型：首先，需要根据实际问题建立一个几何模型，包括问题的几何形状、边界条件等。
2. 离散化：将几何模型离散化为有限个有限元，每个有限元由若干个节点和相应的单元组成。
离散化的过程需要根据具体问题选择合适的网格划分方法。
3. 建立方程：根据所研究问题的物理学原理，建立相应的数学模型和方程。
这些方程通常是偏微分方程。
4. 组装方程：将各个有限元的局部方程组合成一个整体方程，得到整个系统的方程。
5. 解方程：通过数值计算方法，求解得到整个系统的方程的数值解。
常用的求解方法包括直接法、迭代法等。
6. 后处理：对数值解进行后处理，包括计算各个节点的位移、应力等。
以张年梅研究的某个工程问题为例，假设他研究的是一根梁的弯曲问题。
他首先需要建立梁的几何模型，包括梁的长度、截面形状等。
然后，他将梁离散化为若干个有限元，每个有限元由两个节点和一个单元组成。
接下来，他根据梁的力学原理，建立相应的偏微分方程。
然后，他将各个有限元的局部方程组装成一个整体方程。
最后，他通过数值计算方法，求解得到整个系统的方程的数值解。
通过后处理，他可以得到梁的位移、应力等。
有限元方法具有以下几个优点：1. 适用性广：有限元方法适用于各种不同的工程问题，包括结构、流体、热传导、电磁场等。
有限元作业
第一章工程实例描述所分析的是汽车连杆的模型，连杆的厚度是0.5in，模型如下图所示，在小头孔的内侧90。
范围内承受P=1000psi的面载荷作用，利用有限元分析该连杆的受力状态。
连杆的材料属性为杨氏模量E=30X106psi，泊松比为0.3。
图1 连杆的几何模型第二章有限元软件ansys分析2.1 定义单元类型及材料属性图2-1 设置单元类型图2-2 设置单元选项图2-3 设置材料属性2.2 建立二维模型图2-4 连杆二维模型2.3生成2D网格图2-5 2D网格2.4 拖动生成3D网格图2-6 生成3D网格2.5 施加约束、载荷和求解图2-7 施加约束和载荷2.6 浏览分析结果图2-8 变形形状结果图2-9 节点位移结果图2-10 节点的von mises应力MAXIMUM ABSOLUTE VALUES（最大位移）NODE 1036 2092 19 335VALUE 0.43595E-04 0.33341E-04 0.13754E-02 0.13760E-02MINIMUM VALUES（最小应力）NODE 1557 1557 2100 1791 1791 VALUE -16.528 -61.165 -1087.2 0.25937E-06 0.24537E-06 MAXIMUM VALUES（最大应力）NODE 1023 24 1103 2098 2098VALUE 2273.0 222.25 43.989 3192.8 2813.72.7 生成源程序/BATCH/COM,ANSYS RELEASE 11.0SP1 UP2007083009:29:56 06/26/2010 /TITLE,ansys of c-rod/REPLOT/PREP7!*ET,1,MESH200!*KEYOPT,1,1,7 KEYOPT,1,2,0!*ET,2,SOLID95!*!* MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,30e6 MPDATA,PRXY,1,,0.3 ESIZE,0.2,0, MSHAPE,0,2D MSHKEY,0!*CM,_Y,AREAASEL, , , , 2 CM,_Y1,AREA CHKMSH,'AREA' CMSEL,S,_Y!*AMESH,_Y1!*CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2!*TYPE, 2MAT, 1 REAL,ESYS, 0 SECNUM,!*TYPE, 2EXTOPT,ESIZE,3,0,EXTOPT,ACLEAR,0!*EXTOPT,ATTR,0,0,0MAT,1REAL,_Z4ESYS,0!*!*VOFFST,2,0.5, ,/VIEW,1,1,1,1/ANG,1/REP,FAST/VIEW,1,1,2,3/ANG,1/REP,FAST/VIEW,1,,,1/ANG,1/REP,FASTFINISH/SOL/VIEW,1,1,1,1/ANG,1/REP,FAST/VIEW,1,1,2,3/ANG,1/REP,FAST/VIEW,1,1,2,3/ANG,1/REP,FAST/VIEW,1,,,1/ANG,1/REP,FAST/VIEW,1,1,1,1/ANG,1/REP,FAST/VIEW,1,,,1/ANG,1/REP,FAST/VIEW,1,1,1,1/ANG,1/REP,FASTFLST,2,4,5,ORDE,2FITEM,2,27FITEM,2,-30DA,P51X,SYMM!*/PSF,DEFA, ,2,0,1/PBF,DEFA, ,1/PIC,DEFA, ,1/PSYMB,CS,0/PSYMB,NDIR,0/PSYMB,ESYS,0/PSYMB,LDIV,0/PSYMB,LDIR,0/PSYMB,ADIR,0/PSYMB,ECON,0/PSYMB,XNODE,0/PSYMB,DOT,1/PSYMB,PCONV,/PSYMB,LAYR,0/PSYMB,FBCS,0!*/PBC,ALL, ,1/REP!*/VIEW,1,1,1,1/ANG,1/REP,FAST/VIEW,1,1,2,3/ANG,1/REP,FASTFLST,2,4,5,ORDE,2FITEM,2,31FITEM,2,-34/GO!*SFA,P51X,1,PRES,1000EQSLV,PCG,1E-8/STATUS,SOLU SOLVE FINISH/POST1 PLDISP,1/VIEW,1,1,2,3 /ANG,1/REP,FAST/VIEW,1,1,1,1 /ANG,1/REP,FAST/VIEW,1,,,1/ANG,1/REP,FAST/VIEW,1,1,2,3 /ANG,1/REP,FAST/VIEW,1,1,1,1 /ANG,1/REP,FAST/VIEW,1,1,2,3 /ANG,1/REP,FAST/VIEW,1,1,2,3 /ANG,1/REP,FAST/VIEW,1,1,1,1 /ANG,1/REP,FAST/VIEW,1,1,2,3 /ANG,1/REP,FAST/VIEW,1,1,2,3 /ANG,1/REP,FAST/VIEW,1,1,1,1 /ANG,1/REP,FAST/VIEW,1,1,2,3 /ANG,1/REP,FAST/VIEW,1,,,1/ANG,1/REP,FAST /VIEW,1,1/ANG,1/REP,FAST/VIEW,1,,1/ANG,1/REP,FAST/VIEW,1,,,-1/ANG,1/REP,FAST/VIEW,1,-1/ANG,1/REP,FAST/VIEW,1,,-1/ANG,1/REP,FAST/VIEW,1,1,1,1/ANG,1/REP,FAST/VIEW,1,1,2,3/ANG,1/REP,FAST/ZOOM,1,BACK/ANG,1,30,XS,1/REP,FAST/ANG,1,-30,XS,1/REP,FAST/ANG,1,-30,XS,1/REP,FAST/ANG,1,30,XS,1/REP,FAST/ANG,1,-30,YS,1/REP,FAST/ANG,1,30,YS,1/REP,FAST/VIEW,1,1,2,3/ANG,1/REP,FAST/VIEW,1,1,1,1/ANG,1/REP,FAST/VIEW,1,1,2,3/ANG,1/REP,FAST!*/EFACET,1PLNSOL, U,SUM, 0,1.0!*/EFACET,1PLNSOL, U,SUM, 1,1.0!*/EFACET,1PLNSOL, S,EQV, 1,1.0!*PRNSOL,S,PRIN!*PRNSOL,U,COMPFINISH! /EXIT,AL参考文献[1] 邓凡平. ANSYS 10.0有限元分析自学手册. 北京:人民邮电出版社,2007[2] 博弈创作室. ANSYS 9.0经典产品高级分析技术与实例详解. 北京:中国水利水电出版社，2005[3] 李又村，程赫明，周友坤. 仿真在薄板小孔应力集中现象实验中的应用. 昆明冶金高等专科学校学报，2008，(1)[4] 张洪信，赵清海.ansys有限元分析完全自学手册.机械工业出版社，2008。
有限元作业
有限元作业一、名称:平面结构的静力学分析二、问题描述如图所示不规则板，板厚0.5，孔直径为1，倒角半径为0.4。
约束条件：在左边的孔约束全部自由度，载荷施加在右边圆的下半部分，两端施加最小值（50）逐渐均匀变大到底部的最大值（500）。
三、过程:1.选择单元格类型：新建材料类型Main Menu> Preprocessor > Material Props> Material Models> Structural > Linear > Elastic >Isotropic.并设置属性四、建立模型1.运行Main Menu> Preprocessor> Modeling> Create> Areas> Rectangle> By Dimensions分别输入：X1 = 0 X2 = 6 Y1 = -1 Y2 = 1X1 = 4 X2 = 6 Y1 = -1 Y2 = -3（1）在两个矩形的基础上依次建立两个圆形1运行Utility Menu> PlotCtrls> Pan, Zoom, Rotate2单击small dot3运行Utility Menu> WorkPlane> Display Working Plane4运行Utility Menu> WorkPlane> WP Settings点Polar .Grid. Triad5运行Main Menu> Preprocessor> Modeling> Create> Areas> Circle> Solid Circle6输入WP X = 0 WP Y = 07建立半径为1的圆.8同样的方法建立第二个圆（2）进行布尔运算运行Main Menu> Preprocessor> Modeling> Operate> Booleans> Add>Areas（3）倒圆角1.运行Utility Menu> PlotCtrls> Numbering点line numbering.2.运行Utility Menu> WorkPlane> Display Working Plane Main Menu>Preprocessor> Modeling> Create> Lines> Line Fillet选择17 和 8.3.输入半径0.44.运行Utility Menu> Plot> Lines5.运行Utility Menu> PlotCtrls> Pan, Zoom, Rotate点Zoom button.放大该区域6.运行Main Menu> Preprocessor> Modeling> Create> Areas> Arbitrary> ByLines选线4.5和17.运行Utility Menu> Plot> Areas分别建立两个孔1.运行Utility Menu> WorkPlane> Display Working Plane2.运行Main Menu> Preprocessor> Modeling> Create> Areas> Circle> SolidCircle输入WP X = 0WP Y = 03.建立半径为0.4的孔4.运行Utility Menu> WorkPlane> Offset WP to> Global Origin5.运行Main Menu> Preprocessor> Modeling> Create> Areas> Circle> SolidCircle选取点WP X = 0WP Y = 0建立第二个孔6.运行Utility Menu> WorkPlane> Display Working Plane7.运行Utility Menu> Plot> Replot和Utility Menu> Plot> Lines进行布尔运算运行Main Menu> Preprocessor> Modeling> Operate> Booleans> Add> Areas【一】设置材料属性运行Main Menu> Preprocessor> Material Props> Material Models如下图设置【二】定义实常数运行Main Menu> Preprocessor> Real Constants> Add/Edit/Delete输入【三】划分网格1.运行Main Menu> Preprocessor> Meshing> Mesh Tool Set Global Sizecontrol. 输入0.5。
有限元上机作业
有限元上机作业目录1分布载荷作用下的悬臂梁应力计算 (1)1.1问题描述 (1)1.2计算模型 (1)1.2.1有限元模型及网格划分 (1)1.2.2边界条件 (1)1.2.3材料模型 (2)1.3计算结果 (2)1.4结果检验 (2)1.5小结 (3)2悬臂式连接环的应力与变形分析 (4)2.1问题描述 (4)2.2计算模型 (4)2.2.1有限元模型及网格划分 (4)2.2.2边界条件 (5)2.2.3材料模型 (5)2.3计算结果 (5)2.4结果检验 (6)2.5小结 (7)3具有中心孔的薄壁圆筒受均匀拉伸分析 (7)3.1问题描述 (7)3.2计算模型 (8)3.2.1有限元模型及网格划分 (8)3.2.2边界条件 (8)3.2.3材料模型 (8)3.3计算结果 (9)3.4结果检验 (10)3.5小结 (11)4刚架与弹簧混合建模与分析 (12)4.1问题描述 (12)4.2计算模型 (12)4.2.1有限元模型及网格划分 (12)4.2.2边界条件 (13)4.2.3材料模型 (13)4.3计算结果 (14)4.4结果检验 (15)4.5小结 (17)5过盈配合圆环受力分析 (17)5.1问题描述 (17)5.2计算模型 (18)5.2.1有限元模型及网格划分 (18)5.2.2边界条件 (18)5.2.3材料模型 (18)5.3计算结果 (19)5.4结果检验 (19)5.5小结 (20)6带中心圆孔板的热应力分析 (20)6.1问题描述 (20)6.2计算模型 (21)6.2.1有限元模型及网格划分 (21)6.2.2边界条件 (21)6.2.3材料模型 (22)6.3计算结果 (22)6.4小结 (23)7壳—梁组合结构建模计算 (24)7.1问题描述 (24)7.2计算模型 (24)7.2.1有限元模型及网格划分 (24)7.2.2边界条件 (25)7.2.3材料模型 (25)7.3计算结果 (26)7.4小结 (27)8复杂3D实体建模及受力分析 (27)8.1问题描述 (27)8.2计算模型 (28)8.2.1有限元模型及网格划分 (28)8.2.2边界条件 (28)8.2.3材料模型 (28)8.3计算结果 (29)8.4小结 (29)9含裂缝结构体的应力强度因子 (30)9.1问题描述 (30)9.2计算模型 (30)9.2.1有限元模型及网格划分 (30)9.2.2边界条件 (31)9.2.3材料模型 (31)9.3计算结果 (31)9.4小结 (32)10受热载荷的正方形烟囱建模与温度场求解 (32)10.1问题描述 (32)10.2计算模型 (33)10.2.1有限元模型及网格划分 (33)10.2.2边界条件 (33)10.2.3材料模型 (33)10.3计算结果 (34)11旋转车轮的建模与应力计算 (34)11.1问题描述 (34)11.2计算模型 (34)11.2.1有限元模型及网格划分 (34)11.2.2边界条件 (35)11.2.3材料模型 (35)11.3计算结果 (35)12U型支架的模态分析 (36)12.1问题描述 (36)12.2计算模型 (36)12.2.1有限元模型及网格划分 (36)12.2.2边界条件 (37)12.2.3材料模型 (37)12.3计算结果 (37)13异形截面梁的几何特性和扭转切应力分布 (38)13.1问题描述 (38)13.2计算模型 (39)13.2.1有限元模型及网格划分 (39)13.2.2边界条件 (40)13.2.3材料模型 (40)13.3计算结果 (40)1分布载荷作用下的悬臂梁应力计算1.1 问题描述分析模型如图1-1 所示, 梁的横截面为矩形（长x宽x高= 10x1x2 m2），受到沿长度方向线性分布载荷作用，q=(10000-1000x)N/m。
现代设计方法(关于有限元)作业
《现代设计方法》作业关于有限元法的研究学院：机械工程学院专业：机械制造及其自动化0.有限元法有限元法分析起源于50年代初杆系结构矩阵的分析。
随后，Clough于1960年第一次提出了“有限元法”的概念。
其基本思想是利用结构离散化的概念，将连续介质体或复杂结构体划分成许多有限大小的子区域的集合体，每一个子区域称为单元（或元素），单元的集合称为网格，实际的连续介质体（或结构体）可以看成是这些单元在它们的节点上相互连接而组成的等效集合体；通过对每个单元力学特性的分析，再将各个单元的特性矩阵组集成可以建立整体结构的力学方程式，即力学计算模型；按照所选用计算程序的要求，输入所需的数据和信息，运用计算机进行求解。
当前，有限元方法/理论已经发展的相当成熟和完善，而计算机技术的不断革新，又在很大程度上推进了有限元法分析在工程技术领域的应用。
然而，如此快速地推广和应用使得人们很容易忽视一个前提，即有限元分析软件提供的计算结果是否可靠、满足使用精度的前提，是合理地使用软件和专业的工程分析。
有限元法分析一般包括四个步骤：物理模型的简化、数学模型的程序化、计算模型的数值化和计算结果的分析。
每一个步骤在操作过程中都或多或少地引入了误差，这些误差的累积最终可能会对计算结果造成灾难性的影响，进而蒙蔽我们的认识和判断。
1.受内压空心圆筒的轴对称有限元分析例图1.1所示为一无限长的受内压的轴对称圆筒，该圆筒置于内径为120mm的刚性圆孔中，试求圆筒内径处的位移。
结构的材料参数为：200=，0.3E GPaμ=。
图1 结构图对该问题进行有限元分析的过程如下。
（1）结构的离散化与编号由于该圆筒为无限长，取出中间一段（20mm高），采用两个三角形轴对称单元，如图1.2所示。
对该系统进行离散，单元编号及结点编号如图1.3所示，有关结点和单元的信息见表1.1。
图1.2 有限元模型图1.3 节点位移编号及单元编号表1.1 单元编号及结点编号单元编号结 点 编 号 ①②1 2 3 2 3 4结构的结点位移列阵为11223344[]T r r r r u w u w u w u w δ= (1.1) 结构的结点外载列阵12[000000]T r r F F F = （1.2）1r F 和2r F 为由内压作用而等效在结点1和结点2上的载荷，其大小为1122240202//502622r r r h p F N F N ππ-⨯⨯⨯==== （1. 3） 约束的支反力矩阵123344[00T z z r z r z R R R R R R R = ] (1.4)其中1z R 和2z R 为结点1和结点2在Z 方向的约束支反力，（3r R ，3z R ）和（4r R ，4z R ）为结点3和结点4在r 方向和Z 方向的约束支反力。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1.练习题目内容【2以后为逐步落实练习的过程或结果，正反面为一个练习】
对简单的空间钢架桥梁进行有限元分析。
该钢架桥由方钢组成，顶梁及侧梁、桥身弦杆和底梁分别采用3种不同型号的方钢。
桥长L=40m，桥高H=6m。
桥身由8段桁架组成，每段长5m。
该桥梁可以通行卡车，这里仅考虑卡车位于中间的位置，假设卡车的质量为4800kg，
可以将卡车对桥梁的作用力简化为f1，f1，f3，其中f1=f3=6000N，f2=12000N。
2具体过程
（1）进入ansys程序,设置计算类型
anasys main menu:preferences>structural>ok
（2）定义单元类型
anasys main menu:preprocessor>element type>add/edit/delete>add>beam:188>ok>close
（3）定义截面形状及参数
anasys main menu:preprocessor>sections>beam>common sections>beam tool: 1 , 0.1，0.1>apply,
2, 0.05，0.05>apply,3, 0.15，0.15>apply。
（4）定义材料参数
anasys main menu:preprocessor>material models>structual>linear>elastic>isotropic>input ex: 2.06e11,
Prxy:0.3（泊松比）>ok>density>input dens：7800（材料密度），>ok>close
（5）构造桁架桥模型
生成桥体几何模型
anasys main menu:preprocessor>modeling>create>keypoints>in active cs>npt kepoint number:1,
x,y,z location in active cs:0,0,0>apply>同样输入其他31个特征点坐标(5,0,0),(10,0,0),(15,0,0),(20,0,0)(25,0,0) (30,0,0),(35,0,0),(40,0,0),(5,6,0),(10,6,0),(15,6,0),(20,6,0),(25,6,0),(30,6,0),(35,6,0),(0,0,5),(5,0,5)(10,0,5),(15,0,5) (20,0,5),(25,0,5),(30,0,5),(35,0,5),(40,0,5),(5,6,5),(10,6,5),(15,6,5),(20,6,5),(25,6,5),(30,6,5),(35,6,5) >line>lines
>straight line>依次分别连接特征点>ok
网格划分
anasys main menu:preprocessor>meshing>mesh attributes>picked lines>选择桥顶梁及侧梁>ok>
select sections:1>apply>选择桥体弦杆>ok>select sections:2>apply选择桥底梁>ok>select sections:3>ok
anasys main menu:preprocessor>meshing>mesh tool>位于size controls下的lines:set>element size on
Picked>pick all>apply>ndiv:1>ok>mesh>lines>pick all>ok
（6）模型加约束
anasys main menu:solution>define loads>apply>structual>displacement>on nodes>选取桥身左端节点
>ok>select lab2:all dof>apply>选取桥身右端节点>ok>select lab2:uy >ok
（7）施加载荷
anasys main menu:solution>define loads>apply>strudtual>force/moment>on keypoints>选取底梁上卡车
两侧关键点（x坐标为15,25）>ok>select lab:fy,value:-6000>apply>选取底梁上卡车中部关键点
（x坐标为20）>ok>select lab:fy,value:-12000>ok
anasys utility menu:>select>everything
（8）计算分析
anasys main menu:solution>solve>current ls>ok
（9）结果显示
（1）y方向位移云图
anasys main menu:general postproc>plot results>deformed shape>def shape only>ok>contour plot
>nodal solu>dof solution,y-component of displacement>ok
2.结果显示。
（1）有限元模型
（2）加载后的图形
（3）变形图
（4）y方向位移图。