传热学中的有限元法数值分析
有限元论文
机械结构有限元分析 作业名称：基于ANSYS的机械结构仿真学生姓名：陆宁 学号： 班级：机械电子工程103班 指导教师：谢占山老师 作业时间： 2013.05.28 二零一二----二零一三第二学习期
基于ANSYS的机械结构仿真 摘要：介绍了ANSYS优化设计模块，并针对机械结构优化设计给出了具体设计步骤，利用实例分析介绍ANSYS在机械结构优化设计中的应用。证明了ANSYS优化设计模块在机械结构优化设计上的方便性和可行性，为从事机械优化设计人员提供了新的方法和思路。 关键词：机械结构；ANSYS；优化设计；悬臂梁 前言：有现场合，比如，在研究桥梁的受迫振动时，由于激振载荷和和桥梁自重比较接近，所以桥梁自重是必须考虑的因素。激振载荷是正弦载荷，桥梁自重是静载荷，此时桥梁同时受静载荷和正弦载荷的作用。当结构只作用于静载荷时，可以用静力学分析计算其应力、应变等；当结构只作用于正弦载荷时，可以对其进行谐响分析。但是当结构同时作用于静载荷和正弦载荷时，却无法单独用静力学分析或谐响应分析来求解问题，因静力学分析要求载荷恒定，谐响应分析施加的载荷都是正弦载荷。如果用瞬态分析，则载荷就不能是从负无穷时刻到正无穷时刻的周期函数，即施加载荷要对正弦载荷进行加窗处理，势必存在误差，此时就应用有限元法进行分析。
一、基于ANSYS参数化语言的机械结构优化设计概述 机械最优化设计是在现代计算机广泛应用的基础上发展起来的一门新学科，是根据最优化原理和方法综合各方面的因素，以人机配合方式或/自动探索0方式在计算机上进行的半自动或自动设计，以选出在现有工程条件下最佳设计方案的一种现代设计方法.人机连接的传媒是靠一些编程语言来实现，例如C、C十十、VC、FOR-TRAM 等等，这些语言要求用户必须有深厚的理论知识，对于普通用户实现起来就显得很困难。 ANSYS软件是容结构、热、流体、电磁、声学于一体的大型通用有限元分析软件，其内嵌的参数化设计语言(APDL)用建立智能分析的手段为用户提供了自动完成循环的功能，即程序的输入可设定为根据指定的函数、变量以及选出的分析标准作决定.这样的功能扩展完全满足优化设计的要求，而且其强大的前处理建模、可视化界面也是其他优化语言所无法比拟的，更重要的是ANSYSAPDL编程语句简单，更具人性化即使是普通用户也能够掌握。 目前，关于利用ANSYS进行机械优化设计的文献鲜有报道[C17，本文具体剖析了ANSYS优化设计模块，并运用ANSYS12.0的参数化语言求解机械工程设计中的优化问题，给出了在机械优化设计方面的实现方法和具体实例，旨在为从事机械优化设计的人员提供一种新的方法和思路。
《结构分析中的有限元法》2015-有限元习题-参考答案
本科有限元习题参考答案
2015年3月10日作业 1、简述力学课程中介绍的各种力学模型的简化条件、基本假设和适用范围（包括有拉压杆模型、弯曲梁模型、平面应力和平面应变模型、轴对称模型、板模型、壳模型等） 2、给出弹性力学问题中平衡方程、几何方程、物理方程的表达式及其意义。 （1）平衡方程：
zy yz xz zx yx xy z yz xz z y xy zy y x zx yx x f y x z f x z y f z y x ττττττττσττσττσ====+??+??+??=+??+??+??=+??+??+??,000， 物理意义：应力分量与体力分量之间的关系。 （2）几何方程： z u x w y w z v x v y u z w y v x u zx yz xy z y x ??+??=??+??=??+??=??=??=??=γγγεεε,,,, 物理意义：应变分量与位移分量之间的关系。 （3）物理方程： [] [] [] zx zx yz yz xy xy y x z z z x y y z y x x G G G E E E τγτγτγσσμσεσσμσεσσμσε1,1,1) (1 ) (1 )(1 ===+-=+-=+-= 物理意义：应变分量与应力分量之间的关系。 3、简述最小势能原理的主要内容和主要公式。 根据虚功原理得到：??=-Γ T Ω T T 0Td Γδu d Ω）F δu -σδε（，由 )(21εδσεδδεU T T =?? ? ??=则0)21((=Γ-Ω-=∏??ΩΓ）Td u d F u T T T p σεδδ 其中，??ΩΓ Γ-Ω-=∏Td u d F u T T T p )21 (σε即为系统的总势能，它是弹性体变 形势能和外力势能之和。上面变分为零式表明：在所有区域内满足几何关系，在边界上满足给定位移条件的可能位移中，真实位移使系统的总势能取驻值(可证
基于有限元的电磁场仿真与数值计算介绍
鼠笼异步电动机磁场的有限元分析 摘要 鼠笼异步电动机具有结构简单、价格低廉、运行可靠、效率较高、维修方便等一系列的优点，在国民经济中得到广泛的应用。工业、农业、交通运输、国防工程以及日常生活中都大量使用鼠笼异步电动机。随着大功率电子技术的发展，异步电动机变频调速得到越来越广泛的应用，使得鼠笼异步电动机在一些高性能传动领域也得到使用。 鼠笼异步电动机可靠性高，但由于种种原因，其故障仍时有发生。由于电动机结构设计不合理，制造时存在缺陷，是造成故障的原因之一。对电机内部的电磁场进行正确的磁路分析，是电机设计不可或缺的步骤。利用有限元法对电机内部磁场进行数值分析，可以保证磁路分析的准确性。本文利用Ansys Maxwell软件，建立了鼠笼式异步电机的物理模型，并结合数学模型和边界条件，完成了对鼠笼式异步电动机的磁场仿真，得到了物理模型剖分图，磁力线和磁通分布图，为电机的进一步设计研究提供了依据。 关键词：Ansys Maxwell；鼠笼式异步电机；有限元分析
一、前言 当电机运行时，在它的内部空间，包括铜与铁所占的空间区域，存在着电磁场，这个电磁场是由定、转子电流所产生的。电机中电磁场在不同媒介中的分布、变化及与电流的交链情况，决定了电机的运行状态与性能。因此，研究电机中的电磁场对分析和设计电机具有重要的意义。 在对应用于交流传动的异步电机进行电磁场的分析计算时，传统的计算方法因建立在磁场简化和实验修正的经验参数的基础之上，其计算精度就往往不能满足要求。如果从电磁场的理论着手，研究场的分布，再根据课题的要求进行计算，就有可能得到满意的结果。电机电磁场的计算方法大致可以分为解析法、图解法、模拟法和数值计算法。数值解法是将所求电磁场的区域剖分成有限多的网格或单元，通过数学上的处理，建立以网格或单元上各节点的求解函数值为未知量的代数方程组。由于电子计算机的应用日益普遍，所以电机电磁场的数值解法得到了很大发展，它的适用范围超过了所有其它的解法，并能达到足够的精度。对于电机电磁场问题，常用的数值解法有差分法和有限元法两种。用有限元法时单元的剖分灵活性大，适用性强，解的精度高。因此我们采用有限元法对电机电磁场进行数值计算。 Maxwell2D 是一个功能强大、结果精确、易于使用的二维电磁场有限元分析软件。在这里，我们利用Ansys的Maxwell2D 有限元分析工具对一个三相四极电机进行有限元分析，构建鼠笼式异步电机电动机的物理模型，并结合电机的数学模型、边界条件进行磁场分析。
数值分析_迭代法
华北科技学院上机报告 系（部） 专业、班级 学号 课程名称数值分析 上机题目实验六，实验七 任课教师 指导教师 成绩（优、良、中、及格、不及格） 华北科技学院基础部
实验六 解线性方程组的迭代法 1.目的与要求： 1） 熟悉求解线性方程组的有关理论哈方法。 2） 会编制雅可比迭代和高斯—塞得尔迭代法。 3） 通过实际计算，进一步了解各算法的优缺点，选择合适的数值方法。 2．雅可比迭代法 算法 设方程组AX=b 的系数矩阵的对角元素0(1,2, ,),ii i n a ≠=M 为迭代次数容许的最大值,ε 为容许误差. ① 取初始向量(0)(0)(0),(,,,)12T x x x x n =令k=0; ② 对1,2,,i n =计算 (1) ()11();n k k i i ij j j ii j i x b a x a +=≠= -∑ ③ 如果 (1)() 1 ,n k k i i i x x ε+=-<∑则输出(1) k x +,结束;否则执行④, ④ 如果,k M ≥则不收敛,终止程序;否则1,k k ←+转②. 1.分别用雅可比迭代法与高斯-塞德尔迭代法解下列方程组: 2),311300010000 151335901100002709311000000 230010793000090,0 00305770502000007473000120000030410070000500272700 2 2910RI V R V =---????????---????????---????---??????==----???-??????--???--??????--? ???其中???? ??? ? 1.用雅可比迭代法计算： #include "stdafx.h" #include "iostream.h"
有限元与数值方法-讲稿19 弹塑性增量有限元分析课件
材料非线性问题有限元方法 教学要求和内容 1.掌握弹塑性本构关系和塑性力学的基本法则； 2.掌握弹塑性增量分析的有限元格式； 3.学习常用非线性方程组的求解方法： （1）直接迭代法； （2） Newton-Raphson 方法，修正的N－R 方法； （3）增量法等。 请大家预习，争取对相关内容有大概的了解和把握。
弹塑性增量有限元分析 一．材料弹塑性行为的描述 弹塑性材料进入塑性的特点：存在 不可恢复的塑性变形； 卸载时：非线性弹性材料按原路径 卸载； 弹塑性材料按不同的路径卸载，并 且有残余应变，称为塑性应变。
1．单向加载 1） 弹性阶段: 卸载时不留下残余变形; 2） 初始屈服：s σσ= 3） 强化阶段：超过初始屈服之后,按弹性规律卸载，再加载弹性范 围扩大：ss σσ'>，s σ'为相继屈服应力。
4） 鲍氏现象（Bauschinger ）： 二．塑性力学的基本法则 1．初始屈服准则： 00(,)0ij F k σ= 已经建立了多种屈服准则： （1） V . Mises 准则：000(,)()0ij ij F k f k σσ=-= 2 2 001 1 ()(),()2 3ij ij ij s f s s J k σσ===第二应力不变量1122221 ,() 3 ij ij ij m m s σδσσσσσ=-=++偏应力张量：平均应力： （2） Tresca 准则（最大剪应力准则）： 0max ()0ij s F S ττ=-=
2．流动法则 V . Mises 流动法则： 0(,)()ij ij p ij ij ij F k f d d d σσελ λ σσ??==??， 0d λ> 待定有限量 塑性应变增量 p ij d ε 沿屈服面当前应力点的法线方向增加。 因此，称为法向流动法则。 3．硬化法则： （1）各向同性硬化：(,)()0ij ij F k f k σσ=-=
梁结构应力分布ANSYS分析汇总
J I A N G S U U N I V E R S I T Y 先进制造及模具设计制造实验 梁结构应力分布ANSYS分析 学院名称：机械工程学院 专业班级：研1402 学生姓名：XX 学生学号：S1403062 2015年5 月
梁结构应力分布ANSYS分析 （XX，S1403062，江苏大学） 摘要：本文比较典型地介绍了如何用有限元分析工具分析梁结构受到静力时的应力的分布状态。我们遵循对梁结构进行有限元分析的方法，建立了一个完整的有限元分析过程。首先是建立梁结构模型，然后进行网格划分，接着进行约束和加载，最后计算得出结论，输出各种图像供设计时参考。通过本论文，我们对有限元法在现代工程结构设计中的作用、使用方法有个初步的认识。 关键词：梁结构；应力状态；有限元分析；梁结构模型。 Beam structure stress distribution of ANSYS analysis （Dingrui, S1403062, Jiangsu university） Abstract: This article is typically introduced how to use the finite element analysis tool to analyze the stress of beam structure under static state distribution. We follow the beam structure finite element analysis method, established the finite element analysis of a complete process. Is good beam structure model is established first, and then to carry on the grid, then for constraint and load, calculated the final conclusion, the output of images for design reference. In this article, we have the role of the finite element method in modern engineering structural design, use method has a preliminary understanding. Key words: beam structure; Stress state; The finite element analysis; Beam structure model. 1引言 在现代机械工程设计中，梁是运用得比较多的一种结构。梁结构简单，当是受到复杂外力、力矩作用时，可以手动计算应力情况。手动计算虽然方法简单，但计算量大，不容易保证准确性。相比而言，有限元分析方法借助计算机，计算精度高，
数值计算迭代法
习题二 3、证明：当X 0=1.5时，迭代法X k+1=Xk +410和X k+1=21k X 310-都收敛于方程f(x)=x 3+4x 2-10=0在区间[1,2]内唯一实根x *，并分别用上述迭代法求满足于精度要求︱X k+1-X k ︱≤10-5的近似根。 解：证明：{先用迭代法求f(x)=x 3+4x 2-10=0的根。 （a ）对x 3+4x 2-10=0变形有：4x 2=10-x 3 所以：X=21310X - 则相应的迭代公式为：X k+1=21k X 310- 取：X 0=1.5，根据计算可以看出看，我们认为得到的迭代序列是 收敛的。}（此行可忽略） { 由 f(x)=x 3+4x 2-10=0得迭代方程：X=21310X -=g （x ） 先证明在区间【1,2】上x=g （x ）有实根。由于[1,2]上g ‘（x ）存在，所以g （x ）连续。作Q （x ）=x-g(x),则Q(x)在[1,2]上也连续。由定理1条件2有：Q （1）=1-g （1）≤0，Q （,2）=1-g （2）≥0 故存在x *∈[1,2]使Q *（x ）=0，即x *= Q *（x ） 又因为，x *是方程f(x)=x 3+4x 2-10=0在区间[1,2]内的唯一实根，（由定理一条件 2）对任意的x 0∈[1,2]时，X k ∈[1,2]（k=0,1,2,3…） 因为：x *- X k+1=g （x *）-g （X k ）=g ‘（h k ）（x *- X k ）故由条件1知： ︱X *-X k+1︱≤L ︱X *-X k ︱（k=0,1,2,3…）于是有：0≤︱X *-X k ︱≤L k ︱X *-X 0︱,0＜L ＜1，立即可知：lim （k 趋于无穷）︱X *-X k ︱=0,从而lim （k 趋于无穷）X k= X *。所以当X 0=1.5时，迭代法X k+1=Xk +410和X k+1=21k X 310-都是由迭代法X k+1=g （X k ）产生的迭代序列{ X k }收敛于方程f(x)=x 3+4x 2-10=0在区间[1,2]内唯一实 根x *。 正解如下： (1) （牛顿迭代法）： 证明：对方程f(x)=x 3+4x 2-10=0在区间[1,2]内， （a ） f ‘(x)=3x 2+8x ，f ’‘(x)=6x+8，f ’‘(x)在区间[1,2]内连续； （b ） f （1）=-5，f （2）=14，f （1）f （2）＜0； （c ） 对于任意的x ∈[1,2]，都有f ‘(x)=/(不等于)0； （d ） f ’‘(x)在[1,2]上保号； 综上所述，当X 0=1.5时，迭代法X k+1=Xk +410和X k+1=21k X 310-都收敛于方程f(x)=x 3+4x 2-10=0在区间[1,2]内唯一实根x *。 （2）用牛顿迭代法求近似根。 方程f(x)=x 3+4x 2-10=0有唯一实根x *∈[1,2]，容易验证，f(x)=x 3+4x 2-10在[1,2]
ANSYS实体建模有限元分析-课程设计报告
南京理工大学 课程设计说明书(论文) 作者:学号： 学院(系):理学院 专业:工程力学 题目:ANSYS实体建模有限元分析 指导者： (姓名) (专业技术职务) 评阅者： (姓名) (专业技术职务) 20 年月日
练习题一 要求： 照图利用ANSYS软件建立实体模型和有限元离散模型，说明所用单元种类、单元总数和节点数。 操作步骤： 拟采用自底向上建模方式建模。 1.定义工作文件名和工作标题 1)选择Utility Menu>File>Change Jobname命令，出现Change Jobname对话框，在[/FILNAM ] Enter new jobname文本框中输入工作文件名learning1，单击OK按钮关闭该对话框。 2)选择Utility Menu>File>Change Title命令，出现Change Title对话框，在[/TITLE] Enter new title文本框中输入08dp,单击OK按钮关闭该对话框。 2.定义单元类型 1)选择Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令，出现Element Types对话框，单击Add按钮,出现 Library of Element Types 对话框。在Library of Element Types 列表框中选择 Structural Solid, Tet 10node 92，在Element type reference number文本框中输入1，单击OK按钮关闭该对话框。 2)单击Element Types对话框上的Close按钮，关闭该对话框。 3.创建几何模型 1)选择Utility Menu>P1otCtrls>Style>Colors>Reverse Video命令，设置显示颜色。 2)选择Utility Menu>P1otCtrls>View Settings>Viewing Direction命令，出现Viewing Direction对话框，在XV,YV,ZV Coords of view point文本框中分别输入1, 1, 1，其余选项采用默认设置，单击OK按钮关闭该对话框。 3)建立支座底块 选择Main Menu>Preprocessor> Modeling>Create>volumes>Block>By Demensios 命令，出现Create Block by Demensios对话框，在X1,X2 X-coor dinates文本框
工程数值方法与有限元分析
工程数值方法与有限元分析 （机械工程学院机械类专业） 课程号： 周学时：4 学分：3 课程类别： 预修课程：高等数学，线性代数，力学基础课 面向对象：机械类专业学生 教学方式：多媒体教学 教学目的和教学要求： 在科学研究与工程技术中，经常遇到数学模型的求解问题。然而在许多情况下，要获得模型问题的准确解往往是十分困难的，甚至是不可能的。因此，研究各种数学问题的近似解法非常必要。计算方法是一门与计算机应用密切结合的实用性很强的课程，它专门研究各种数学问题的一类近似解法，从一组原始数据出发，按照确定的运算规则进行有限步运算，最终获得问题的数值形式且满足精度要求的近似解。 通过对《计算方法》的学习，掌握数值计算的基本概念和基本理论，深入理解方法的设计原理与处理问题的技巧，重视误差分析与收敛性、数值稳定性，注重利用计算机进行科学计算能力的培养,并熟练掌握Matlab 软件,会用Matlab实现各种计算方法。 在此基础上进一步学习数值计算的集大成者-有限元方法, 了解有限元方法的基础知识及其在机械、机械电子领域中的应用，掌握有限元方法的基本原理与分析过程，包括静力学、动力学、非线性力学、热场、电磁场等的建模及分析。学生可使用有限元软件进行机械零件及系统的实例分析，并对分析结果进行评价，指导和优化机械零件及系统的设计。本课程面向机械电子专业及机械类相关专业的高年级本科生 课程简介： 内容主要包括：计算机上常用的数值计算方法以及有关的基本概念与理论，主要有误差、非线性方程求根、线性代数方程组的解法、插值与拟合、数值微分与数值积分、常微分方程初值问题的数值解法。并且算法面向计算机，注重培养学生运用计算机进行科学计算解决工程问题的能力。并熟练掌握Matlab 软件，会用Matlab实现各种计算方法。 有限元的分析与建模是一个机械工程师必须掌握的方法和技能。本课程为机械类专业的高年级学生核心课，使学生了解有限元方法的基本概念和基本理论，掌握有限元分析的基本处理方法，熟悉常用有限元分析软件在实际工程中的应用，最终培养学生在机械设计、机电系统设计中能有效的应用有限元方法。 主要内容及学时分配： 每周4学时，共16周 主要内容： ( O ) 绪论1学时 (一)误差2学时
数值分析拉格朗日插值法上机实验报告
课题一：拉格朗日插值法 1.实验目的 1．学习和掌握拉格朗日插值多项式。 2.运用拉格朗日插值多项式进行计算。 2.实验过程 作出插值点（1.00，0.00），（-1.00，-3.00），（2.00，4.00）二、算法步骤 已知：某些点的坐标以及点数。 输入：条件点数以及这些点的坐标。 输出：根据给定的点求出其对应的拉格朗日插值多项式的值。 3.程序流程： （1）输入已知点的个数； （2）分别输入已知点的X坐标； （3）分别输入已知点的Y坐标； 程序如下： #include
插值算法*/ { int i,j; float *a,yy=0.0; /*a作为临时变量，记录拉格朗日插值多项*/ a=(float*)malloc(n*sizeof(float)); for(i=0;i<=n-1;i++) { a[i]=y[i]; for(j=0;j<=n-1;j++) if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]); yy+=a[i]; } free(a); return yy; } int main() { int i; int n; float x[20],y[20],xx,yy; printf("Input n:");
scanf("%d",&n); if(n<=0) { printf("Error! The value of n must in (0,20)."); getch();return 1; } for(i=0;i<=n-1;i++) { printf("x[%d]:",i); scanf("%f",&x[i]); } printf("\n"); for(i=0;i<=n-1;i++) { printf("y[%d]:",i);scanf("%f",&y[i]); } printf("\n"); printf("Input xx:"); scanf("%f",&xx); yy=lagrange(x,y,xx,n); printf("x=%f,y=%f\n",xx,yy); getch(); } 举例如下：已知当x=1，-1,2时f(x)=0，-3,4，求f(1.5)的值。
数值分析5-用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组
作业六：分别编写用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组Ax=B的标准程序，并求下列方程组的解。 可取初始向量 X(0) =(0，0，0)’； 迭代终止条件||x(k+1)-x(k)||<=10e-6 （1） = （2） = Jacobi迭代法： 流程图 开 始 判断b中的最大值 有没有比误差大 给x赋初值 进行迭代 求出x，弱到100次还没到，警告不收 结束
程序 clear;clc; A=[8,-1,1;2,10,01;1,1,-5]; b=[1;4;3]; e=1e-6; x0=[0;0;0]'; n=length(A); x=zeros(n,1); k=0; r=max(abs(b)); while r>e for i=1:n d=A(i,i); if abs(d)
程序结果（1）
（2）
Gauss-Seidel迭代法： 程序 clear;clc; %A=[8,-1,1;2,10,01;1,1,-5]; %b=[1;4;3]; A=[5,2,1;-1,4,2;2,-3,10]; b=[-12;20;3]; m=size(A); if m(1)~=m(2) error('矩阵A不是方阵'); end n=length(b); %初始化 N=0;%迭代次数 L=zeros(n);%分解A=D+L+U,D是对角阵，L是下三角阵，U是上三角阵U=zeros(n); D=zeros(n); G=zeros(n);%G=-inv(D+L)*U d=zeros(n,1);%d=inv(D+L)*b x=zeros(n,1); for i=1:n%初始化L和U for j=1:n if i
有限元分析课程论文2011
《ANSYS10.0基础及工程应用》考查要求 一、课程考核方式 撰写课程结课论文。 二、论文撰写范围 在掌握有限元基本理论及方法的基础上，运用《ANSYS10.0基础及工程应用》课程所学的建模，分网，加载，求解及后处理知识，针对某一你所熟悉的产品、设备或零件进行有限元计算分析。 三、论文撰写要求 １．论文按科技论文的标准格式撰写，包括有题目、作者、单位（班级、学号、联系方式）、摘要（200字左右）、关键词（3—4个）、正文及参考文献（包括作者姓名、文献名、出版社所在地、出版社名、出版时间等），正文引用文献要标出，严禁抄袭。２．全文字数不少于3000字。 ３．参考文献至少５篇。 4．统一以武汉理工大学华夏学院论文纸。
有限元分析课程要求 要求：1）个人至少分析3种方案并独立完成（可选择一个模型三种不同方案或三个不同模型的有限元分析；题目可从上机指南，有限元分析大作业试题中选择或自行选择算例），并将计算 结果分析在论文中较详细分析说明（包括几何模型视图、单元模型视图、结果云图，矢量 分布图，列表，命令流等及结果分析说明。） 2）课程论文应包括以下部分：（正文5号字体） A、引言； B、问题描述及几何建模； C、有限元建模（单元选择、节点布置及规模、网格划分方案、载荷及边界条件 处理、求解控制） D、计算结果及结果分析（位移分析、应力分析、正确性分析评判，如同一模型 则必须进行多方案计算比较，需讨论节点规模增减对精度的影响分析、单元 改变对精度的影响分析、不同网格划分方案、不同结构对结果的影响分析等） E、结论 F、参考文献 3）12月1日前必须完成，并递交课程论文报告（报告要求打印）。 4）学生的课程总评成绩由平时成绩（占30%）和期末考查成绩（占70%）两部分构成。平时成绩中包括出勤、作业、上机操作、学习主动性等。
结构分析及有限元分析基础知识
第一章结构分析及有限元分析基础知识 注：摘自《NX知识工程应用技术——CAD/CAE篇》 洪如瑾编译 清华大学出版社 [目标] 本章将简述结构分析及有限元分析的基础知识，为学习与应用结构分析做好准备，包括： ※ 结构与结构分析定义 ※ 结构的线性静态分析 ※ 材料行为与故障 ※ 有限元分析的基本概念 ※ 有限元模型 1.1结构分析基础知识 1.1.1结构基本概念 1.结构定义 结构可以定义为一个正承受作用的载荷处于平衡中的系统。平衡条件意味着结构是不移动的。一个自由的支架不是一个结构，它未被连接到任一物体上并无载荷作用与它。仅当它附着到外部世界，并且有作用力、压力或力矩时，支架成为一个结构。 例如横跨江面的大桥就是一个普通的结构，一个支架通过它的支撑连接到地面上，桥的重量是在结构上的一种载荷（力）。当汽车通过桥时，附加的力作用于桥的不同位置。 一个好的结构必须满足以下标准： （1） 当预期的载荷作用时，结构必须不出现故障。这个似乎是显而易见的，并意味着结构必须是“强度足够的”。故障意味着结构破裂、分离、弯曲，以及支撑作用载荷失败。 注意：考虑到意外的载荷，通常在设计中提供安全余量。余量常常利用安全因素来描述。例如，如果在结构上期待载荷是10 000磅，规定安全因素是2.0，则结构将设计成能经受住20 000磅载荷。 （2） 当载荷作用时，结构必须不产生过分变形。这意味着结构必须“刚度足够”。 变形可接受的极限（弯曲度、挠度、拉伸等）取决于特定情况。例如，在通常住宅中的地板由足够的吊带支撑，以防止当人在地板岸上行走时有“柔软”的感觉。 （3） 在它的服务生命周期，结构的行为应不会恶化。这意味着结构必须“足够耐用”，必须考虑环境影响和“磨损与破裂”。如果一座桥假定维持50年，则桥的设计必须提供整个50年寿命的结构完整性与充分的安全余量。2.结构分析 结构分析是用于决定一个结构是否将正确完成任务的工程分析过程。结构将在某些方式中进行模拟和求解描述它的行为的数学方程。分析可以人工方法或用计算机方法来完成。 结构分析的结果（答案）用于评估性能，摘要如下： （1）“强度足够吗？”：应力必须是在一可接受的范围内。 （2）“刚度足够吗？”：位移必须是在一可接受的范围内。 （3）“耐用度足够？”：对一个长的疲劳周期应力必须足够低。
滚压有限元数值分析
滚压模型有限元分析 姓名：黄肖 学号：201721000362 指导老师：陈勇
目录 1.模型的参数 (1) 2.模型的建立 (1) 3.部件的属性模块 (4) 4.网格的划分 (6) 5.创建分析步 (8) 6.创建接触 (8) 7.创建相关 (10) 8.施加载荷 (11) 9..创建作业及结果 (11) 10．结果分析 (14)
滚压有限元数值模拟是制定滚压工艺，预测滚压工作表面残余应力分布，以及判定工件疲劳性的重要工具。目前滚压数值模拟主要集中在对曲周以及回转体的分析，少有的对平面滚压数值模拟，大多数也只分析了单圈或者不到一圈的滚程，而且与实际滚压工艺存在较大区别。本文采用有限元商业软件ABAQUS模拟研究了滚压力的大小对残余应力分布规律的影响。 1.模型的参数 圆柱状滚针滚压平面工件，其数值模拟工件材料为Q235，物性参数如表1所示。滚针相对于工件变形较小，可以忽略不计。 计算模型材料工艺参数 滚针与工件的接触表面其相比远小于滚针的直径，取工件的局部尺寸为2mm×15mm滚针的直径为1mm,依次选取滚压力为120N,160N,180N.摩擦系数为0.3进行计算。 2.模型的建立 (1)打开ABAQU软件，点击创建Great Part。
(2)开始绘制模型 输入（0，0），然后enter，再输入（2，15），然后enter。
再点击Great Part。建立滚针的模型图。 建立滚针的模型图 （3）将两个部件装配到一起形成一个整体。
装配后的效果图如下： 3.部件的属性模块 （1）点击Module中的Propety,进入如下界面，置材料属性。 分别设置工件和滚针的材料属性：弹性模量和泊松比，屈服强度和屈服应力，
有限元分析论文
用有限元分析Hyperworks结构 机制1091 19号何志强 论文关键词:拓扑优化形状优化精密铸造后悬置支架有限元分析 论文摘要: 本文主要阐述借助于Alatir公司的Hyperworks结构优化软件，对精密铸造产品进行结构优化设计，且以对某汽车驾驶室后悬置支架的结构优化为例，着重介绍了拓扑优化和形状优化在精密铸造产品结构设计上的应用方法及功能。事实表明拓扑优化和形状优化的联合应用，对精密铸造产品的结构设计起到非常关键的帮助作用，最后通过此软件对优化后的产品结构进行有限元分析，验证优化后产品结构的强度和刚度。 HyperWorks在精密铸造产品优化设计中的应用 一、引言 在当前的汽车工业中，减轻设计重量和缩短设计周期是两个突出的问题，在传统的设计中，由于机械产品机构的复杂性，长期以来主要应用经验类比设计，对产品结构作定性分析和经验类比估算，在决定实际结构时，一般都取较大的安全系数，结果使得产品都是“傻”、“大”、“粗”，使材料的潜力得不到充分发挥，产品的性能也得不到充分的把握。所以传统的汽车设计思路已经不能满足当前设计的需要。汽车轻量化设计开始占据了汽车发展中的主要地位，它既可以提高车辆的动力性，降低成本，减少能源消耗又能减少污染。但是，简单的汽车轻量化设计却是一把双刃剑，它在减轻汽车重量的同时，也牺牲了车辆的强度和刚度，甚至对产品的结构寿命也产生影响，在此情况下，有限元分析方法在汽车设计中的合理应用就得到了充分体现，经过近几年的实践证明，Altair公司的有限元分析技术以及拓扑优化技术在汽车行业获得了非常成功的应用。特别是对于一些结构复杂的汽车铸造结构件，Hyperworks 的有限元分析技术、拓扑优化和形状优化技术的推广使得材料的潜能及铸造的优势得到了充分的发挥。 本文将详细介绍利用Hyperworks的拓扑优化和形状优化技术对东风商用车驾驶室后悬置支架进行减重优化设计的应用过程。以及如何应用Hyperworks验证改进结构后的应力和应变情况，使该后悬置支架减重优化后的结构能够满足产品的使用性能和铸造工艺性要求。 二、有限元法的概念和优化设计流程确立 2.1有限元法和有限单元的概念 有限元法又称有限单元法，是结构分析的一种数值计算方法，它随着计算机的发展而应运而生，并得到了广泛应用，目前已成为工程数值分析的有力工具。在实际工程应用中，我们首先把CAD模型分割成有限个实体或者壳单元。一般作为实体单元所适合的结构，是具有三维形状变化的物体，不太适合棒状、平板状的物体。实体单元是利用3D-CAD所作
浅析拉格朗日有限元数值计算法
中国矿业大学2011级硕士研究生课程考试试卷 考试科目岩土工程数值计算法 考试时间2011.11.27 学生姓名夏明 学号ZS11020068 所在院系矿业工程学院 任课教师徐志伟 中国矿业大学研究生院培养管理处印制
对论文《基于拉格朗日差分法的露天边坡稳定性研 究》中数值计算法浅析 《基于拉格朗日差分法的露天边坡稳定性研究》一文以易门铜矿露天开采境界优化方案下的边坡为工程背景，通过岩体构造调查，质量分类，室内力学试验，力学参数的工程处理。建立铜厂露天开采边坡的三维地质模型，采用DIMINE 数字矿山软件的耦合集成技术—四面体网格化将地质模型转换为力学模型，应用基于拉格朗日法的有限差分(FLAC3D)大变形方法对铜厂露天开采边坡的稳定性进行了数值模拟，分析了基于强度折减理论计算出的边坡安全系数以及基于莫尔库伦屈服准则的边坡开采后的位移、应力等的变化状况，得出了易门铜厂露天矿露天境界优化方案下的边坡的稳定性状况。 岩质边坡稳定性评价的方法分主要有：极限平衡法、数值模拟计算、地质力学物理模拟试验和其它新方法。随着计算机技术和计算方法的发展，复杂的工程问题可以采用离散化的数值计算方法并借助计算机得到满足工程要求的数值解，数值模拟技术是现代工程学形成和发展的重要动力之一。通过计算模拟，可以模拟并得到模拟体内部的应力—应变关系，再现其变形甚至破坏过程及其机制。在岩土工程数值分析中最常用的数值方法有有限元法、离散元法、边界元法等。 拉格朗日差分法((FLAC法)源于流体力学。它首先是Cundail在80年代提出来的，其基本原理类似于离散元法，但它却能像有限元那样适用于多种材料模式与边界条件的非规则区域的连续问题求解。在求解的过程中，FLAC又采用了离散元的动态松弛法，不需要求解大型联立方程，便于在微机上实现。另一方面，同以往的差分分析相比，FLAC在以下几个方面做了较大的改进和发展:它不但能处理一般的大变形问题，而且能模拟岩体沿某一弱面产生的滑移变形。一般有限单元法可以用来解决材料非线性问题，但对于大变形的几何非线性问题，有限单元法和边界元法都无能为力。拉格朗日法是分析非线性大变形问题的数值方法，它依然遵循连续介质的假设，基于拖带坐标系的基本原理。用差分法或按时步显式迭代求解，不但可以解决几何非线性，也能解决材料非线性问题。 下面是原文的部分描述： 一、露天矿边坡三维地质模型的建立 1.1计算几何模型范围及模拟方案设计 a模拟计算方案设计 决定边坡稳定性的因素很多，其中包括岩体强度，结构特征，水理性质和边坡几何参数等。根据现场工程地质调查，力学试验参数，对露天境界边坡进行修改，修改方案同第五章简化实体模型的建立。
基于abaqus的ujoint有限元分析有限元分析论文大学论文
有限元分析课程论文 课程名称：有限元分析 论文题目：ujoint有限元分析学生班级； 学生姓名： 任课教师： 学位类别： 评分标准及分值选题与参阅资料 （分值） 论文内容 （分值） 论文表述 （分值） 创新性 （分值） 评分 论文评语： 总评分评阅教师: 评阅时间 年月日 注：此表为每个学生的论文封面，请任课教师填写分项分值
基于abaqus的ujoint有限元分析 摘要：万向传动装置在汽车中起到了传递扭矩的关键作用,在abaqus中导入ujoint实体模型，之后对其进行坐标系建立，wire 建立，以及各部件之间的连接关系的建立，最后对该模型施加边界条件，令其运动。 关键词：abaqus、有限元、ujoint 一问题的描述 对导入的ujoint在所有步骤完成后，施加力：在step initial：均设为0；step SPIN：doundary1：限制除 UR2的所有，且把UR2值设为：pi。在boundary2 中，限制UR1和UR3自由度。 二在abaqus中导入ujoint实体模型 启动abaqus CAE，在文件下拉菜单中选择：import ， 选择最终文件位置or 输入ws_connector_ujoint.py.inp 打开文件ujoint。（如下图所示）
2.1 创建坐标系 单机操作界面中的tool，从下拉菜单中选择datum，再出来的窗口中选择coordinate，3points。首先选择origin，在选择x正方向，Y正方向、z正方向。创建完成。 2.2创建VERT和CROSS之间的2坐标系。 根据 2.1所述操作步骤创建坐标系V-C 和V-G （VERT和GROUND）。 Notice：1、创建过程中为了清晰分辨，可将IN的suppress，创建完成后再将其resume。其他同样 2、在V-C和I-C中，x轴与cross转动所绕轴平行。
有限元数值分析课程论文
《有限元数值分析》课程论文要求 1．简述平面问题有限元法（自选单元类型）。 2．写出计算程序框图。 3．自编一份平面问题有限元弹性应力计算程序。 4．以图1为算例，，分析计算结果。 图1
1.简述平面问题有限元法 1.1有限元法的基本思想： 1.1.1假想把连续系统分割成数目有限的单元，单元之间只在数目有限的指定点处（称为节点）相互连接，构成一个单元集合体来代替原来的连续系统。在节点上引进等效载荷，代替实际作用于系统上的外载荷。 1.1.2对每个单元由分块挖的思想，按一定的规则建立求解未知量与节点相互作用之间的关系。 1.1.3把所有单元的这种特性关系按一定的条件集合起来，引入边界条件，构成一组以节点变量为未知量的代数方程组，求解之就得到有限个节点处的待求变量。 1.2平面问题包括平面应力和平面应变问题 1.2.1平面应力问题 平面应力问题研究等厚度薄板状弹性体，厚度尺寸远远小于截面尺寸， 受力方向沿板面方向，不沿厚度变化。 1.2.2平面应变问题 平面应变问题处理面内受力但垂直于平面方向上不产生变形的二维受力问题，沿截面方向的截面形状和大小相同且厚度尺寸远远大于截面尺寸，载荷垂直于厚度方向且沿厚度均匀分布。 1.3结合所给模型说明平面问题的三角形单元求解 1.3.1 结构的离散化 将所给实际模型简化，并划分为如图所示的2个三节点三角形单元，并对其进行编码，节点编号和坐标、单元编号以及相关节点如表1、表2所示。 图2
表1 节点编码及其坐标 表2 单元编码及其相关节点 表3 节点及节点载荷 表4 节点及节点位移 节点载荷列阵(1) 节点位移列阵：（2）1.3.2 分别计算单元的刚度矩阵 将每个单元局部编码，都可看做如图4的形式，每个单元都有其对应的i、j、m。 图3 单元刚度矩阵（t为板料厚度）（3）
数值分析常用的插值方法
数值分析报告 班级： 专业： 流水号： 学号： 姓名：
常用的插值方法 序言 在离散数据的基础上补插连续函数，使得这条连续曲线通过全部给定的离散数据点。插值是离散函数逼近的重要方法，利用它可通过函数在有限个点处的取值状况，估算出函数在其他点处的近似值。 早在6世纪，中国的刘焯已将等距二次插值用于天文计算。17世纪之后，牛顿、拉格朗日分别讨论了等距和非等距的一般插值公式。在近代，插值法仍然是数据处理和编制函数表的常用工具，又是数值积分、数值微分、非线性方程求根和微分方程数值解法的重要基础，许多求解计算公式都是以插值为基础导出的。 插值问题的提法是：假定区间[a，b〕上的实值函数f（x）在该区间上n＋1个互不相同点x0，x1……x n处的值是f（x0），……f（x n），要求估算f（x）在[a，b〕中某点的值。其做法是：在事先选定的一个由简单函数构成的有n＋1个参数C0， C1，……C n的函数类Φ(C0，C1，……C n)中求出满足条件P（x i）＝f（x i）（i＝0，1，……n）的函数P(x)，并以P(x)作为f(x)的估值。此处f（x）称为被插值函数，x0，x1，……xn 称为插值结(节)点，Φ(C0，C1，……C n)称为插值函数类，上面等式称为插值条件，Φ（C0，……C n）中满足上式的函数称为插值函数，R（x）＝f（x）－P（x）称为插值余项。
求解这类问题，它有很多种插值法，其中以拉格朗日(Lagrange)插值和牛顿(Newton)插值为代表的多项式插值最有特点，常用的插值还有Hermit 插值，分段插值和样条插值。 一．拉格朗日插值 1.问题提出： 已知函数()y f x =在n+1个点01,,,n x x x L 上的函数值01,,,n y y y L ，求任意一点 x '的函数值()f x '。 说明：函数()y f x =可能是未知的；也可能是已知的，但它比较复杂，很难计算其函数值()f x '。 2.解决方法： 构造一个n 次代数多项式函数()n P x 来替代未知（或复杂）函数()y f x =，则 用()n P x '作为函数值()f x '的近似值。 设()2012n n n P x a a x a x a x =++++L ，构造()n P x 即是确定n+1个多项式的系数 012,,,,n a a a a L 。 3.构造()n P x 的依据： 当多项式函数()n P x 也同时过已知的n+1个点时，我们可以认为多项式函数 ()n P x 逼近于原来的函数()f x 。根据这个条件，可以写出非齐次线性方程组： 20102000 20112111 2012n n n n n n n n n n a a x a x a x y a a x a x a x y a a x a x a x y ?++++=?++++=?? ? ?++++=?L L L L L 其系数矩阵的行列式D 为范德萌行列式： ()20 0021110 2111n n i j n i j n n n n x x x x x x D x x x x x ≥>≥= = -∏L L M M M M L