有限元分析大作业 展示
ansys有限元分析工程实例大作业
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。
图1钢桥的形式图2桥梁的简化平面模型(取桥梁的一半)图3刚桁架桥简图所用的桁架杆件有三种规格,见表1表1 钢桁架杆件规格杆件截面号形状规格顶梁及侧梁 1 工字形400X400X12X12桥身弦杆 2 工字形400X300X12X12底梁 3 工字形400X400X16X16所用的材料属性,见表2表2 材料属性参数钢材弹性模量EX泊松比PRXY 0.3密度DENS 78002 模型构建将下承式钢桁梁桥的各部分杆件,包括顶梁及侧梁,桥身弦杆,底梁均采用BEAM188单元,此空间梁单元可以考虑所模拟杆件的轴向变形; 定义了一套材料属性,各类杆件为钢材,其对应的参数如表2所示;根据表1中的杆件规格定义了三种梁单元截面,根据表1分别定义在相应的梁上;建模时直接建立节点和单元,在后续按照先建节点再建杆的次序一次建模。
机械结构有限元分析大作业
ห้องสมุดไป่ตู้
254 50.8 25.4 25.4
1334.5
260 48.2 24.1 24.1
1320.4
248 46.8 23.4 23.4
1330.8
240 45.8 22.9 22.9
1348.3
245 52.2 26.1 26.1
1318.4
一、建立几何模型 1、建立4个关键点 Main menu >Preprocessor >Modeling >Create >Keypoints >in Active CS,弹 出下面对话框,输入坐标值
2、连直线 Main Menu >Preprocessor >Modeling >Create >Lines >Lines >in Active CS,弹出下示对话框。用光标点屏幕上点(直线段)
3、组合成平面图形 Main Menu >Preprocessor >Modeling >Create >Areas >Arbitrary >By lines,弹出对话框。用光标点各线段。
四、加载和求解 1、定义分析类型 Main Menu >Solution >Analysis Type >New Analysis,选择Static 2、定义位移约束 Main Menu >Solution >Define Loads >Apply >Structral >Displacement > On Keypionts,弹出下面的对话框,点 关键点,点OK,弹出对话框后,选All DOF,值输入:0 3、加载荷 Main Menu >Solution >Define Loads >Apply >Structral >Force/Moment >On Keypionts,弹出下面的对话框,点右边 加载点,点OK,弹出对话框,选FY,输入 1320.4
《有限元分析》课程作业
《有限元分析》课程作业任课教师:徐亚兰学生姓名:陈新杰学号:班级: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. 以小组为单位完成有限元分析计算,并将计算结果上交;2. 以小组为单位撰写计算分析报告;3. 按下列模板格式完成分析报告;4. 计算结果要求提交电子版,一个算例对应一个文件夹,报告要求提交电子版和纸质版。
《有限元分析》大作业小组成员:储成峰李凡张晓东朱臻极高彬月Job name :banshou完成日期:2016-11-22一、问题描述(要求:应结合图对问题进行详细描述,同时应清楚阐述所研究问题的受力状况和约束情况。
图应清楚、明晰,且有必要的尺寸数据。
)如图所示,为一内六角螺栓扳手,其轴线形状和尺寸如图,横截面为一外接圆半径为0.01m的正六边形,拧紧力F为600N,计算扳手拧紧时的应力分布<图1扳手的几何结构(要求:针对问题描述给出相应的数学模型,应包含示意图,示意图中应有必要的尺寸数据;数学模型图2数学模型如图二所示,扳手结构简单,直接按其结构进行有限元分析。
三、有限元建模3.1单元选择(要求:给出单元类型,并结合图对单元类型进行必要阐述, 包括节点、自由度、 实常数等。
)图3单元类型扳手截面为六边形,采用4节点182单元,182单元可用来对固体结构进行如进行了简化等处理,此处还应给出文字说明。
)Figure 1B2.1: PLANE1S2 Geometry二维建模。
182单元可以当作一个平面单元,或者一个轴对称单元。
它由4个结点组成,每个结点有2个自由度,分别在x,y方向。
扳手为规则三维实体,选择8节点185单元,它由8个节点组成,每个节点有3个自由度,分别在x,y,z方向。
3.2实常数(要求:给出实常数的具体数值,如无需定义实常数,需明确指出对于本问题选择的单元类型,无需定义实常数。
)因为该单元类型无实常数,所以无需定义实常数3.3材料模型(要求:指出选择的材料模型,包括必要的参数数据。
)对于三维结构静力学,应力主要满足广义虎克定律,因此对应ANSYS^的线性,弹性,各项同性,弹性模量EX 2e11 Pa,泊松比PRXY=0.33.4几何建模由于扳手结构比较简单,所以可以直接在ANSYS软件上直接建模,在ANSYS建立正六边形,再创立直线,面沿线挤出体,得到扳手几何模型VULUHES 1TYPE NUMEZ图4几何建模3.5网格划分方案(要求:指出网格划分方法,网格控制参数,最终生成的单元总数和节点总数,此外还应附上最终划分好的网格截图。
有限元分析大作业
有限元大作业一题目要求:图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号节点的支反力如下图所示:。
有限元分析大作业
《有限元分析及应用》大作业——齿根弯曲应力计算报告班级:无可奉告姓名:无可奉告学号:无可奉告指导老师:无可奉告目录目录 (2)1.概述 (3)1.1工程问题描述 (3)1.2问题分析 (3)2.建模过程 (4)2.1几何建模 (4)2.2CAE网格划分与计算 (5)2.3后处理 (8)3.多方案比较与结果分析 (9)3.1多方案比较 (9)3.2结果分析 (11)1.概述1.1工程问题描述我在本次作业中的选题为齿根弯曲应力的计算与校核。
通过对机械设计的学习,我们可以知道,齿轮的失效形式主要是齿面接触疲劳和齿根弯曲断裂,而闭式传动硬齿面齿轮的失效形式以齿根弯曲断裂,这个时候进行齿根弯曲应力的校核才比较有意义,在设计问题的时候应当选取这种类型的算例。
设计计算的另一个主要思路是将有限元计算的结果与传统机械设计的结算结果进行对比,以从多方面验证计算结果的准确性。
综上,我们最终选取了《机械原理》(第三版)P50例3-1中的问题进行校核计算。
已知起重机械用的一对闭式直齿圆柱齿轮,传动,输入转速n1=730r/min,输入功率P1=35kW,每天工作16小时,使用寿命5年,齿轮为非对称布置,轴的刚性较大,原动机为电动机,工作机载荷为中等冲击。
z1=29,z2=129,m=2.5mm,b1=48mm,b2=42mm,大、小齿轮均为20CrMnTi,渗碳淬火,齿面硬度为58~62HRC,齿轮精度为7级,试验算齿轮强度。
齿面为硬齿面,传动方式为闭式传动。
根据设计手册查出的许用接触应力为1363.6Mpa,计算结果为1260Mpa,强度合格。
根据设计手册查出的许用弯曲应力为613.3MPa,计算结果为619Mpa,强度略显不够。
1.2问题分析大小齿轮啮合,小齿轮受载荷情况较为严峻,故分析对象应当为小齿轮。
可以看出,由于齿轮单侧受载荷,传动过程中每个齿上载荷的变化过程是相同的,故问题可被简化为反对称问题,仅需研究单个齿。
ansys有限元分析作业经典案例
a n s y s有限元分析作业经典案例本页仅作为文档封面,使用时可以删除This document is for reference only-rar21year.March有限元分析作业作业名称输气管道有限元建模分析姓名陈腾飞学号 3070611062班级 07机制(2)班宁波理工学院题目描述:输气管道的有限元建模与分析计算分析模型如图1所示承受内压:1.0e8 PaR1=0.3R2=0.5管道材料参数:弹性模量E=200Gpa;泊松比v=0.26。
图1受均匀内压的输气管道计算分析模型(截面图)题目分析:由于管道沿长度方向的尺寸远远大于管道的直径,在计算过程中忽略管道的断面效应,认为在其方向上无应变产生。
然后根据结构的对称性,只要分析其中1/4即可。
此外,需注意分析过程中的单位统一。
操作步骤1.定义工作文件名和工作标题1.定义工作文件名。
执行Utility Menu-File→Chang Jobname-3070611062,单击OK按钮。
2.定义工作标题。
执行Utility Menu-File→Change Tile-chentengfei3070611062,单击OK按钮。
3.更改目录。
执行Utility Menu-File→ change the working directory –D/chen 2.定义单元类型和材料属性1. 设置计算类型ANSYS Main Menu: Preferences →select Structural → OK2.选择单元类型。
执行ANSYS Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 8node 82 →applyAdd/Edit/Delete →Add →select Solid Brick 8node 185 →OKOptions…→select K3: Plane strain →OK→Close如图2所示,选择OK接受单元类型并关闭对话框。
有限元分析与应用大作业
有限元分析及应用大作业课程名称: 有限元分析及应用班级:姓名:试题2:图示薄板左边固定,右边受均布压力P=100Kn/m作用,板厚度为0.3cm;试采用如下方案,对其进行有限元分析,并对结果进行比较。
1)三节点常应变单元;(2个和200个单元)2)四节点矩形单元;(1个和50个单元)3)八节点等参单元。
(1个和20个单元)图2-1 薄板结构及受力图一、建模由图2-1可知,此薄板长和宽分别为2m和1.5m,厚度仅为0.3cm,本题所研究问题为平面应力问题。
经计算,平板右边受均匀载荷P=33.33MPa,而左边被固定,所以要完全约束个方向的自由度,如图2-2所示。
取弹性模量E=2.1×11Pa,泊松比μ=0.3。
P=33.33MPa图2-2 数学模型二、第一问三节点常应变单元(2个和200个单元)三节点单元类型为PLANE42,设置好单元类型后,实常数设置板厚为0.3M。
采用2个单元的网格划分后的结果如图2-3,200个单元的网格划分图如图2-6所示。
约束的施加方式和载荷分布如图2-2中所示。
约束右边线上节点全部自由度。
计算得到的位移云图分别如图2-4、7所示,应力云图如图2-5、8所示。
图2-3 2个三角形单元的网格划分图图2-4 2个三角形单元的位移云图图2-5 2个三角形单元的应力云图图2-7 200个三角形单元的位移云图三、第二问四节点矩形单元的计算四节点单元类型为PLANE42,设置好单元类型后,实常数设置板厚为0.3M。
采用1个单元的网格划分后的结果如图2-9,50个单元的网格划分图如图2-12所示。
约束的施加方式和载荷分布如图2-2中所示。
约束右边线上节点全部自由度。
计算得到的位移云图分别如图2-10、11所示,应力云图如图2-13、14所示。
图2-9 1个四边形单元的网格划分图图2-11 1个四边形单元的应力云图图2-12 50个四边形单元的网格划分图图2-13 50个四边形单元的位移云图图2-14 50个四边形单元的应力云图四、第三问八节点等参单元的计算四节点单元类型为PLANE82,设置好单元类型后,实常数设置板厚为0.3M。
ABAQUS 有限元分析大作业 凹槽成型实例
有限元分析大作业-凹槽成型一:前处理-利用ABAQUS/CAE创建模型。
1、定义并创建四个部件如下图:图1:可变形毛坯Blank图2:刚性冲头Punch图3:刚性夹具Holder2、定义材料及截面特性图5:分别根据提供的数据定义Elastic和Plastic两个材料特性图6:创建均匀实体截面提交给材料Steel并赋予Blank此截面属性图7:建立一个局部数据坐标系(在随着毛坯运动的共旋坐标系下显示应力和应变的值): 2、装配部件图8:装配图(根据相互关系进行装配)3、创建几何集合创建6个几何集合:每个刚性体参考点各一个,毛坯对称面一个,毛坯中面的每段各一个。
图9:创建六个几何集合4、定义分析步和输出要求图10:创建分析步1 Establish contact 1图11:创建分析步2,Remove right constraint图12:创建分析步3 Holder force图13:创建分析步4 Establish contact 2图14:创建分析步5 Move punch图15:编辑场输出图16:编辑历史输出5、监控自由度的值图17:定义RefPunch的监控自由度6、定义接触相互作用图18:首先在Interaction模块中定义以上5个表面图19:定义一个无摩擦接触相互作用属性图20:再定义一个有摩擦相互作用属性,摩擦系数取为0.1 最后定义三个表面间的相互作用:图21:定义三个表面之间的相互作用(具体见下面三个图)图22:定义Die-Blank 相互作用 图23:定义Holder-Blank 相互作用图24:定义Punch-Blank 相互作用7、各分析步的边界条件首先在STEP1总添加每个初始边界条件:图25:Center、MidLeft、MidRight边界条件图26:RefDie、RefHolder、RefPunch边界条件图27:分析步1边界条件图28:Step2边界条件图29:Step3边界条件图30:Step4边界条件图31:Step5边界条件图:32集中力与负压力施加8、划分网格和定义作业图33:首先在毛坯上下与左右表面分别撒种图34:选择单元类型图35:划分网格二:求解1、创建作业图36:创建作业2、提交求解并监视求解过程图:37:求解完成图:38监视点U2向位移图三:后处理1、成型过程图39:开始图40:开始变形图41:变形扩大图42:最后成型结果2、绘制塑性应变等值线图3、绘制冲头上的反作用力4、绘制接触压力等值线图图45:接触压力等值线图。
有限元分析-案例
有限元分析作业作业名称横臂梁有限元建模分析姓名学号班级一、问题描述图25所示为一工字钢梁,两端均为固定端,其截面尺寸为,16,2.0,.0===。
试建立该工字钢梁的三维实体模型,并.0=0.1=,l03cmmdm02b,.0mma在考虑重力的情况下对其进行结构静力分析。
其他已知参数如下:弹性模量(也称杨式模量) E= 206GPa ;泊松比3.0=u ;材料密度3/7800m kg =ρ;重力加速度2/8.9s m g =;作用力Fy 作用于梁的上表面沿长度方向中线处,为分布力,其大小Fy=-5000N二、实训目的本实训的目的是使学生学会掌握ANSYS 在三维实体建模方面的一些技术,并深刻体会ANSYS 软件在网格划分方面的强大功能。
三、结果演示使用ASSY S 8。
0软件对该工字钢梁进行结构静图26单元类型库对话框图25 工字钢结构示意图力分析,显示其节点位移云图。
四、实训步骤(一)ASSYS8.0的启动与设置与实训1第一步骤完全相同,请参考。
(二)单元类型、几何特性及材料特性定义1定义单元类型。
点击主菜单中的“Preprocessor>Element Type >Add/Edit/Delete ”,弹出对话框,点击对话框中的“Add…”按钮,又弹出一对话框(图26),选中该对话框中的“Solid ”和“Brick 8node 45”选项,点击“OK ”,关闭图26对话框,返回至上一级对话框,此时,对话框中出现刚才选中的单元类型:Solid45,如图27所示。
点击“Close ”,关闭图27所示对话框。
注:Solid45单元用于建立三维实体结构的有限元分析模型,该单元由8个节点组成,每个节点具有X 、Y 、Z 方向的三个移动自由度。
2.定义材料特性。
点击主菜单中的 “Preprocessor>Material Props >Material Models ”,弹出窗口如图28所示,逐级双击右框中“Structural\ Linear\ Elastic\ Isotropic ”前图标,弹出下一级对话框,在“弹性模量”(EX )文本框中输入:2.06e11,在“泊松比”(PRXY )文本框中输入:0.3,如图29所示,点击“O K ”图28 材料特性参数对话框按钮,回到上一级对话框,然后,双击右框中的“Density ”选项,在弹出对话框的“DENS ”一栏中输入材料密度:7800,点击“OK ”按钮关闭对话框。
有限元分析大作业
一、有限元方法的手工计算结果与ansys分析结果的对比1分析的问题描述如图1所示,桁架的杆截面面积为8,由钢制成(E=200GPa)。
用有限元法计算出每个节点的位移以及反作用力。
(1)(2)(3)图1对于上述问题,本文将用手工计算和ansys软件分别计算出结果,对计算出来的结果进行对比。
2手工计算2.1桁架结构的有限元计算方法对于桁架结构,每个单元的刚度矩阵为,(2-1)YX图2其中,为桁架单元在整体坐标系中与X轴的夹角;,A为桁架的截面积,E 为弹性模量,L为桁架长度。
在固体力学问题中,有限元公式通常由如下的一般形式,Ku=F(2-2)其中,K为刚度矩阵,u为位移矩阵,F为载荷矩阵。
运用公式(2-3),就能求出反作用力,R=Ku-F(2-3)其中,R为反作用力矩阵。
2.2计算过程计算每个桁架单元的刚度,用公式(2-1)计算每个每个桁架单元的刚度矩阵,将每个单元放入总刚度矩阵,他们的位置分别为:10-100000 00000000 -10100000 00000000 00000000 00000000 00000000 0000000000000000 00000000 0010-1000 0000000000-101000 00000000 00000000 000000003.9-4.90000-3.9 4.9-4.9 6.10000 4.9-6.1 00000000 00000000 00000000 00000000-3.9 4.90000 3.9-4.9 4.9-6.10000-4.9 6.100000000 00000000 00000000 000 1.28000-1.28 00000000 00000000 00000000 000-1.28000 1.2800000000 00000000 00000000 00000000 0000 3.9 4.9-3.9-4.9 0000 4.9 6.1-4.9-6.1 0000-3.9-4.9 3.9 4.9 0000-4.9-6.1 4.9 6.1将个刚度矩阵相加得到总刚度矩阵为,19.9-4.90-16000-3.9 4.9 -4.9 6.100000 4.9-6.1 -160320-16000 00012.8000-12.8 00-16019.9 4.9-3.9-4.9 0000 4.9 6.1-4.9-6.1 -3.9 4.900-3.9-4.97.80 4.9-6.10-12.8-4.9-6.1025应用边界条件施加载荷,将总刚度矩阵带入式(2-2)得:19.9-4.90-16000-3.9 4.9Ux1 -4.9 6.100000 4.9-6.1Uy1 -160320-16000Ux2 00012.8000-12.8Uy200-16019.9 4.9-3.9-4.9Ux3 0000 4.9 6.1-4.9-6.1Uy3 -3.9 4.900-3.9-4.97.80Ux4 4.9-6.10-12.8-4.9-6.1025Uy4带入边界条件解得:将结果带入(2-3)得:=Fx1Fy1Fx2Fy2Fx3Fy3Fx4Fy43用ansys软件求解(单位统一N,mm,Mpa)(1)选择单元(图3)图3(2)附材料属性(图4)图4(3)创建模型(图5)图5(4)施加载荷(图6)图6(5)求解每个节点的位移(图7)图7节点的反力(图8)图8(6)模型变形图(7)位移等值线分布图4结果对比及分析手算结果ansys 计算结果位移(mm)Ux100Uy100Ux2-0.0016-0.0016Uy2-0.0468-0.0468Ux300Uy300Ux4-0.0066-0.0066Uy4-0.0317-0.0317表1手算结果ansys计算结果节点反力(N)Fx1-1027.8-1027.8 Fy11608.31608.3 Fx2 5.60 Fy2-100 Fx32066.72063.1 Fy32257.12255.4 Fx48.80 Fy4-2.70表2由表1和表2可以看出,手工计算的结果与ansys计算的结果基本一致。
有限元受力分析大作业
1. 定义分析类型
Main Menu >Solution > Analysis Type > New Analysis, 选择Static 2. 定义位移约束
Main Menu >Solution> Define Loads >Apply >Structural>Displacement > On Areas, 在 弹出对话框中选Circle,点击OK后,在弹出的对话框中填入:0
后处理
1. 查看总体变形
Main Menu >General Postproc >Plot Results > Deformed shape
(2)
(2)
四. 加载和求解
1. 定义分析类型
Main Menu >Solution > Analysis Type > New Analysis, 选择Static 2. 定义位移约束
Main Menu >Solution> Define Loads >Apply >Structural>Displacement > On Areas, 在 弹出对话框中选Circle,点击OK后,在弹出的对话框中填入:0
七. 后处理
1. 查看总体变形
Main Menu >General Postproc >Plot Results > Deformed shape
S
作业三:轴承座盈利分析
1、 题目
2、建模
3、ansys前处理
导入
ansys前处理
1、定义单元类型 Main Menu>Proprocessor>Element Type >Add/Edit/Delete 弹出对话框 中后,点“Add” 。双弹出对话框,选“Solid”和“10node 92”,点 “OK”,退回到前一个对话框。
有限元分析技术课程大作业
有限元分析技术课程大作业1 工程介绍现需要对某露天大型玻璃平面舞台的钢结构进行分析,该钢结构布置在xy 平面内。
学生序号为079,分格的列数(x向分格)=0×10+7+5=12,分格的行数(y向分格)=9+4=13,共有156个分格。
每个分格x方向尺寸为1m,y方向尺寸为1m。
钢结构的主梁为高160宽100厚14的方钢管;次梁为直径60厚10的圆钢管(单位为毫米),材料均为碳素结构钢Q235;该结构固定支撑点位于左右两端主梁和最中间的次梁的两端。
玻璃采用四点支撑与钢结构连接(采用四点支撑表明垂直作用于玻璃平面的面载荷将传递作用于玻璃所在钢结构分格四周的节点处,表现为点载荷;试对在垂直于玻璃平面方向的22/KN m的面载荷(包括玻璃自重、钢结构自重、活载荷(人员与演出器械载荷)、风载荷等)作用下的舞台进行有限元分析(每分格面载荷对于每一支撑点的载荷可等效于0.5KN的点载荷)。
作业提交的内容:(1)屏幕截图显示该结构的平面布置结构;(2)该结构每个支座的支座反力;(3)该结构节点的最大位移及其所在位置;(4)对该结构中最危险单元(杆件)进行强度校核。
2有限元模型的建立该钢结构中每一分格x方向尺寸为1m,y方向尺寸为1m,x方向分格数量为12,y方向分格数量为13。
该钢结构由主梁和次梁构成,其中主梁为高160mm、宽100mm、厚14mm的方钢管,次梁为直径60mm、厚10mm的圆钢管。
由于在该结构中所有构件均为梁单元,而Ansys程序中提供了多种梁单元,以模拟不同场合的应用,且对于每种梁单元类型都有特定的算法。
在本次建模过程中,考虑到需要对该结构中的危险单元进行强度校核,因此,选择了BEAM188单元类型来建立本钢架结构,进而对其进行有限元分析。
BEAM188为三维线性有限应变梁单元,该单元基于铁木辛哥的梁结构理论,考虑了剪切变形的影响,能够满足本次分析的需求。
以下为基于ANSYS图形界面(Graphic User Interface , GUI)的菜单操作流程(1) 进入ANSYS(设定工作目录和工作文件)程序→ ANSYS → ANSYS Interactive → Working directory(设置工作目录)→Initial Jobname(设置工作文件名):Analysis → Run → OK(2) 设置计算类型ANSYS Main Menu:Preferences → Structural → OK(3) 定义单元类型ANSYS Main Menu:Preprocessor → Element Type → Add/Edit/Delete... → Add → Beam: 3D 2node 188 → OK(返回到Element Types窗口)→ Close(4) 定义材料参数ANSYS Main Menu: Preprocessor → Material Props → Material Models → Structural → Linear → Elastic → Isotropic → input EX: 2.0E5, PRXY: 0.3(定义泊松比及弹性模量) → OK → Close(关闭材料定义窗口)(5)定义梁单元截面ANSYS Main Menu:Preprocessor →Sections→Beam→Common Sections→Beam Tool(6) 构造梁模型生成舞台几何模型ANSYS Main Menu:Preprocessor → Modeling → Create → Keypoints → In Active CS → NPT Keypoint number:1,X,Y,Z Location in active CS:0,0,0 → Apply 通过复制关键点操作,形成14行13列的关键点。
solidworks有限元分析16例
注意:本文件内容只是一个简短的分析报告样板,其内相关的分析条件、设置和结果不一定是正确的,您还是要按本书正文所教的自行来做。
一、范例名: (Gas Valve气压阀)1 设计要求:(1)输入转速1500rpm。
(2)额定输出压力5Mpa,最大压力10Mpa。
2 分析零件该气压泵装置中,推杆活塞、凸轮轴和箱体三个零件是主要的受力零件,因此对这三个零件进行结构分析。
3 分析目的(1)验证零件在给定的载荷下静强度是否满足要求。
(2)分析凸轮轴零件和推杆活塞零件的模态,在工作过程中避开共振频率。
(3)计算凸轮轴零件的工作寿命。
4 分析结果1.。
推杆活塞零件材料:普通碳钢。
在模型上直接测量得活塞推杆的受力面积S为:162mm2,由F=PS计算得该零件端面的力F为:1620N。
所得结果包括:1 静力计算:(1)应力。
如图1-1所示,由应力云图可知,最大应力为21Mpa,静强度设计符合要求。
(2)位移。
如图1-2所示,零件变形导致的最大静位移为2.2e-6m。
(3)应变。
如图1-3所示,应变云图与应力云图的对应的,二者之间存在一转换关系。
图1-1 应力云图图1-2 位移云图图1-3 应变云图图1-4 模态分析2 模态分析:图1-4的“列举模式”对话框中列出了“推杆活塞”零件在工作载荷下,其前三阶的模态的频率远远大于输入转速的频率,因此在启动及工作过程中,该零件不会发生共振情况。
模态验证符合设计要求。
2。
凸轮轴零件材料:45钢,屈服强度355MPa。
根据活塞推杆的受力情况,换算至该零件上的扭矩约为10.5N·m。
1 静力分析:如图1-5所示为“凸轮轴”零件的应力云图,零件上的最大应力为212Mpa,平均应力约为120MPa,零件的安全系数约为1.7,符合设计要求。
图1-5 应力云图图1-6 模态分析2 模态分析图1-6的“列举模式”对话框中列出了“推杆活塞”零件在工作载荷下的模态参数,“模式1”的结果为其自由度内的模态,不作为校核参考。
有限元分析大作业
基于ANSYS软件的有限元分析报告机制1205班杜星宇U201210671一、概述本次大作业主要利用ANSYS软件对桌子的应力和应变进行分析,计算出桌子的最大应力和应变.然后与实际情况进行比较,证明分析的正确性,从而为桌子的优化分析提供了充分的理论依据,并且通过对ANSYS软件的实际操作深刻体会有限元分析方法的基本思想,对有限元分析方法的实际应用有一个大致的认识。
二、问题分析已知:桌子几何尺寸如图所示,单位为mm。
假设桌子的四只脚同地面完全固定,桌子上存放物品,物品产生的均匀分布压力作用在桌面,压力大小等于300Pa,其中弹性模量E=9。
3GPa,泊松比μ=0。
35,密度ρ=560kg/m3,分析桌子的变形和应力.将桌脚固定在地面,然后在桌面施加均匀分布的压力,可以看作对进行平面应力分析,桌脚类似于梁单元。
由于所分析的结构比较规整且为实体,所以可以将单元类型设为八节点六面体单元。
操作步骤如下:1、定义工作文件名和工作标题(1)定义工作文件名:执行UtilityMenu/ File/ChangeJobname,在弹出Change Jobname 对话框修改文件名为Table。
选择New log anderrorfiles复选框。
(2)定义工作标题:Utility Menu/File/Change Title,将弹出ChangeTit le对话框修改工作标题名为The analysis of table。
(3)点击:Plot/Replot。
2、设置计算类型(1)点击:Main Menu/Preferences,选择Structural,点击OK。
3、定义单元类型和材料属性(1)点击:Main Menu/Preprocessor/Element Type/Add/Edit/Delete,点击Add,选择Solid〉Brick 8node 185,点击OK,点击Close。
(2)点击Main menu/preprocessor/Material Props/Material Models / Structural/ Linear/ Elastic/Isotropic,设置EX为9.3e9,PRXY为0。
有限元分析大作业报告
有限元分析大作业报告试题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) 三节点常应变单元三节点常应变单元的位移分布图三节点常应变单元的应力分布图(2)六节点三角形单元六节点三角形单元的变形分布图六节点三角形单元的应力分布图单元类型最小位移(mm)最大位移(mm)最小应力(Pa)最大应力(Pa)三节点0 0.0284 5460.7 392364六节点0 0.0292 0.001385 607043①最大位移都发生在A点,即大坝顶端,最大应力发生在B点附近,即坝底和水的交界处,且整体应力和位移变化分布趋势相似,符合实际情况;②结果显示三节点和六节点单元分析出来的最大应力值相差较大,原因可能是B点产生了虚假应力,造成了最大应力值的不准确性。
有限元分析大作业
机电工程学院有限元分析及应用直齿圆柱齿轮的模态分析学号:S314070064专业:机械工程学生姓名:***任课教师:*** 教授2014年12月一 研究目的齿轮传动是机械工程领域应用最广泛的传动之一,模态分析技术已经成为振动系统分析与设计中广泛使用的重要手段,它是研究结构动力特性的一种近代方法,是系统辨别方法在工程振动领域中的应用。
模态是机械结构的固有振动特性,每一个模态具有特定的固有频率、阻尼比和模态振型。
齿轮在传递运动和动力时,传动系统通过各种外部激励和内部激励传递给齿轮系统,从而使齿轮在传动过程中产生振动。
齿轮的固有频率是齿轮的动态特性之一,对动载荷的产生与传递以及系统的振动形式有很重要的影响,因此分析齿轮的动态特性,对齿轮的设计和改进以及整个传递系统的动态性能的改进都有非常重要的实际意义。
二 齿轮模态求解分析齿轮副在啮合过程中,因加工误差、齿侧间隙和轮齿受载弹性变形及热变形,会产生啮合合成基节误差,使轮齿啮合时产生转速差异与突变,引起振动,也就是固有频率,从传统的静力学分析,固有频率可有下式近似计算mk f π210= (1) 式中:m 和k 分别为齿轮的等效质量和刚度系数,其大小根据查阅手册选取或据经验而定。
传统的模态分析技术无法有效地处理含有接触关系的非线性系统的装配体模态分析问题,为处理此问题,人们采取了一些线性化的近似处理方法,例如将装配体视为单一实体零件,或在将零件间的联接简化成线性弹簧等。
这种线性化的简化分析方法,难以对含有非线性接触联接的装配体进行准确分析。
而且往往要多次计算,消耗大量人力物力,为此在材料力学基础上产生了弹性力学的有限元法。
其中,齿轮系统的运动微分方程为()t F KX X C X M =++ (2)式中:M,C,K 分别是齿轮系统质量矩阵、阻尼矩阵和刚度矩阵,F 为收到外界激振力向量。
若无外力作用,即F(t)=0,则是系统自由振动方程,刚度矩阵与约束有关,但由于啮合部位的接触面积不断变化,K 也会发生相应变化,用传统的线性分析方法不易分析,有限元采用分段逼近方法,模拟连续体的约束条件是求解问题的关键。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
二 不同单元对计算结果的影响
实例:一简支梁,各项数据如表所示:
q / kN l / m b / m h / m E / GPa
0.25
Hale Waihona Puke 100161
3
200
分析其变形和 受力情况。 (分析均取梁 结构的一半进 行计算)
对本例分别采用3节点三角形、4节点矩形、8节点六 面体、20节点六面体进行计算分析,且用平面单元分析时, 还考虑网格的划分不同对结果的影响,分别对梁高划分15 和30等份进行计算, 比较结果。 用平面plane 42 单元建立模型,建模图如下:
谱分析整体结果列表
之后读取结果,在File>Read Input from 里面选择 “.MCOM” 文件,查看各个结果。
谱分析6个节点的反作用力表
读取第一组结果的变形和应力图,如下图 所示:
第一组结果的变形图
第一组结果的应力图
3. 时程分析 实例:上题的三层框架结构,施加了两个x方向大小随时 间变化的力(如下图所示),进行时程分析,最后选择节 点,绘制时程图。
梁高划分30等份三角形3节点建模图
梁高划分15等份三角形3节点建模图
梁高划分30等份4节点矩形建模 图
梁高划分15等份4节点矩形建模图
用六面体8节点单元solid 45和六面体20节点单元solid 95进行分析,建模图如下:
六面体8节点建模图
六面体20节点建模图
各个分析数据均是结构的x方向的应力图,最后进过 汇总,所得的结果如下:
计算分以下几歩完成: 1. 设置为模态分析,进行模态分析求解; 2. 设置为谱分析,输入单点响应谱的频率和 谱值之后(阻尼为0),获得谱解; 3. 再设置为模态分析,扩展模态,扩展5阶模 态,进行扩展模态求解; 4. 最后设置成谱分析,进行模态叠加,选择 SRSS方法,输出类型为Displacement , 最后求解,得出结果,整体结果如下图:
将梁高划 六面体8 六面体20 分15等份 节点单元 节点单元 4节点矩 形单元 2.1521 MPa 0.8% 2.1519 MPa 2.1521 MPa
X方向应 力/MPa 误差
0.7916% 0.8%
通过分析,得出如下结论: 1. 对比平面三角形单元的数据,发现单元划分越细,计算越精确。 2. 对比平面三角形和矩形单元的数据,发现对于此题,选择三角 形单元,计算的结果更加精确。 3. 对比平面矩形单元和六面体单元的数据,发现划分为平面单元 和立体单元,计算结果差别不大,均可采用。
四 结构的动力分析
1. 模态分析
实例:模型仍然为上题 三层框架结构,各个参 数均不变,去掉500N的 荷载,对其进行模态分 析。 求解类型为Modal,模 态提取方法为BLOCK Lanczos,提取前5阶模 态,之后进行求解。
模态分析建模图
模态信息列表
第一阶模态振型图
第二阶模态振型图
第三阶模态振型图
三 空间结构的受力变形分析
实例:三层钢材梁板结构,6个点承受x方向500N的点荷载, 受重力作用,重力加速度为9.8N/kg 。计算x方向变形和应力。 材料为A3钢,板壳厚度为2e-3,其他参数如下表:
2e11Pa
0.3
7.8e3
kg / m3
板壳选用shell 181单元,立柱 选用 beam 188单元
建好的模型图如下图所示:(选用beam188单元)
超静定梁的模型图
求解之后的y方向的位移图如下图所示:
超静定梁y方向位移图
超静定梁y方向位移图
用列表显示y方向的位移,其结果如下:
其结果显示为,最大 位移为SMN=-0.00204 列表图显示为2节点 位移最大,为 -0.20399E-02 ,两者 接近。
第四阶模态振型图
第五阶模态振型图
2. 单点响应谱分析 实例:同样是前题的三层框架结构,计算在Y方向的地震位 移响应谱左右下整个结构的响应情况。地震响应谱如下表:
频率/Hz 0.5 1.0 2.4 3.8 17 18 20 32
位移/mm 1.0 0.5 0.9 0.8 1.2 0.75 0.86 0.2
理论解为2.135MPa,将所有结果进行统一列表,并计算 出各种不同方式分析的误差,从而从中发现问题,列表详细如 下图:
项目
将梁高划 分30等份 3节点三 角形单元 2.1349 MPa 0.00468 %
将梁高划 分15等份 3节点三 角形单元 2.1034 MPa 1.48 %
将梁高划 分30等份 4节点矩 形单元 2.153 MPa 0.843%
荷载F1随时间变化图
荷载F2随时间变化图
施加动荷载通过Parameters>Array Parameters>Define / Edit 来实现。
建模图及加载情况如下图所示:
荷载F2施加 在90号节点 上 荷载F1 施加在 33号节 点上
时程分析建模图
建立好模型后,在进行求解的选项设置(Solution>Sol’n Controls里面,Analysis Options 选择Large Displacement Transient; Time at end of loadstep 设置成1; Number of substeps 设置成300;Frequency选择write every substep),设置好之后,再进行求解,求解结果如 下图:
时程分析求解结果图
时程分析求解整体结果列表图
可以查看上图列表显示任意时间点的变形和应力图,选择 71时间点的变形和应力图,如下图所示:
变形图
应力图
最后进行时域后处理,查看71号节点的x方向变形随时间 的响应关系图,如下图所示:
71号 节点 时程分析71号节点x方向变形随时 间响应关系图 71号节点位置图
还可根据需要显示任意一点的变形随时间的响应关系图, 在此不一一列举。
ANSYS的结构分析实例
目录
结构的静力分析 不同单元对计算结果的影响 空间结构的受力变形分析 结构的动力分析 总结
一 结构的静力分析
实例:一次超静定梁,长为l的梁左端X向和Y向以 及旋转受限,右端Y向受限。梁中点受竖直向下集 中力F。取l=400mm,F=1000N,截面为b=50, h=150的矩形截面。( 2e5, 0.3 )
用材料力学的方法,对梁中点进行截面弯曲应力简单计算:
1 2 1 M ql 105 162 3200kN m 8 8 1 1 3 I z bh 1 33 2.25m 4 12 12 3 M 320010 c y 1.5 2.135MPa Iz 2.25
梁几何性质:截面面积 1.6e 5m 2 , 惯性矩为64 / 3e 12m 4 , 宽度为4e 3m, 高度为4e 3m。
X方向 500N的 荷载
三层框架结构建模图
求解完成后,得出结果:
X方向变形前后对比图
X方向应力前后对比图
框架结构x方向变形最大值
框架结构x方向应力最大值与最小值
划分梁高30等份三角形单元x 方向应力图
划分梁高15等份三角形单元x方向 应力图
划分梁高30等份矩形单元x方向 应力图
划分梁高15等份矩形单元x方向应 力图
六面体8节点单元x方向应力图
六面体20节点单元x方向应力图
将每次的分析结果通过列表显示,提取关键点的数据,也 就是梁中点的最下方的x方向的应力,在ansys列表里面为节 点1,每次记录下节点1的值,最后在通过比较分析,得出结 论。