有限元模拟卷解题思路
1. 用兹法求解下列边值问题：
2210 01
(0)0, (1)0
d y
y x dx y y ++=<<==， （14分）
该问题的泛函为：
1
22011[()]22I y x y y y dx ??
'=-- ???? （1）
设 ()()()1y x x x x ?α≈=- （2）
将（2）代入（1）， 然后 对[()]I y x 求导 令
[()]
0I x ?α
?=? 便可求出α 这样 ()()1y x x x α=-
2. 对于轴对称问题，推导矩形4节点单元的形函数。（14分） 解： 选用四边形环单元 设位移函数为
12345678(,)(,)u r z r z rz v r z r z rz αααααααα=+++=+++
这里18......αα为待定系数
然后代入4个节点的坐标及位移
12345678i i i i i i i i i i
u r z r z v r z r z αααααααα=+++=+++ i=1 ,2, 3, 4
这样就可以求得18......αα的表达式 这样型函数就出来了
四面体四节点形函数的求解方法与之类似 只不过这时候 位移表达式要写成：
123456789101112(,,)(,,)(,,)u x y z x y z v x y z x y z w x y z x y z
αααααααααααα=+++=+++=+++
这时候再代入4个节点的位移和坐标 就可以求得112......αα
1．平面应力问题有限元模型的网格划分如图所示：
1）试对模型的节点及单元进行编号，使得总刚度矩阵的带宽最小。 2）若n
ij k 表示第n 个单元的单元刚度矩阵第i 行第j 列的元素，试写出总刚度矩阵的表达式。 （12分）
这题的解法和作业题某题的解法类似 我们来看作业题
这一题考的是 从节点编号 到 总刚度矩阵合成 再到边界条件下刚度矩阵修改 的过程
这些都是我上课强调过的
我们来看 为了让总刚度矩阵半带宽最小 就是让全零元素次对角线最多
比如这道题目 如果让编号1和6在同一个三角形单元内 那么显然无非零次对角线 因此这种情况下半带宽显然最大，因为这时候 16 和61位置都有非零元素。
I II III IV 为单元编号
因此可以调节节点编号使得半带宽最小 如下：
11121322
23
2433
3435
36
44
4655
56
66000000I I I I II I II II I II III IV II III IV III IV II III III IV IV III IV
k k
k k
k
k k k k
k k
k k k k +++++++++???????????????????
?
对
称
那么这个考题的做法和作业题相同，只不过图形换了一下而已。
2．用里兹法求解区间[0,1]上的一维函数y(x)，且y(x)满足
()()()()
(0)0, (1)0
y x y x y x f x a x b y y '''-+=≤≤==， （10分）
这题的做法和第一题类似，只不过泛函的表达式写法不同，具体泛函求法可参看ppt 第一章具体公式。
2020年数值分析模拟试卷(三)
数值分析模拟试卷（三）班级学号姓名一、填空题（共2分，每题2分） 1、设x*=3149578…，取5位有效数字，则所得的近似值x=_______________ ； . 2、设一阶差商，，则二阶差商__________ ； 3、数值微分中，已知等距节点的函数值，则由三点的求导公式，有_______________ ； 4、求方程的近似根，用迭代公式，取初始值，那么x1=_________ ； 5、解初始值问题近似解的梯形公式是yk+1=_________ ； 6、，则A的谱半径______ ，cond (A)＝______ ； 7、设，则______ ， ______ ； 8、若线性代数方程组AX=b 的系数矩阵A为严格对角占优阵，则雅可比迭代和高斯-塞德尔迭代都_______ ； 9、解常微分方程初值问题的欧拉（Euler）方法的局部截断误差为_____ 1、设，当____________时，必有分解式A=LLT，其中L为下三角阵．二、计算题（共6分，每题15分） 1、（1）设试求f(x)在上的三次Hermite插值多项式使满足； （2）写出余项的表达式． 2、已知，满足，试问如何利用构造一个收敛的简单迭代函数，使…收敛？ 3、试确定常数A，B，C和a，使得数值积分公式有尽可能高的代数精度．所得的数值积分公式代数精度是多少？是否为Gauss型的？ 4、推导常微分方程的初值问题的数值解公式 三、证明题（共2分，每题1分） 1、设，（1）写出解 f(x)=的Newton迭代格式； （2）证明此迭代格式是线性收敛的． 2、设R=I－CA，如果，证明 （1）A、C都是非奇异的矩阵； （2）
数值分析试题及答案汇总
数值分析试题 一、 填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则x 有 2 位 有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ，则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商 公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所以当 系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收 敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的，其中的残差 r i ＝ (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ，(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时，若在迭代区间存在唯一解，且f (x )
螺栓连接的有限元分析(汇编)
1 概述 螺栓是机载设备设计中常用的联接件之一。其具有结构简单,拆装方便,调整容易等优点,被广泛应用于航空、航天、汽车以及各种工程结构之中。在航空机载环境下，由于振动冲击的影响，设备往往产生较大的过载，对作为紧固件的螺栓带来强度高要求。螺栓是否满足强度要求，关系到机载设备的稳定性和安全性。 传统力学的解析方法对螺栓进行强度校核，主要是运用力的分解和平移原理，解力学平衡方程，借助理论和经验公式，理想化和公式化。没有考虑到连接部件整体性、力的传递途径、部件的局部细节(如应力集中、应力分布)等等。通过有限元法，整体建模，局部细化，可以弥补传统力学解析的缺陷。用有限元分析软件MSC.Patran/MSC.Nastran提供的特殊单元来模拟螺栓连接，过程更方便，计算更精确，结果更可靠。因此，有限元在螺栓强度校核中的应用越来越广泛。 2 有限元模型的建立 对于螺栓的模拟，有多种模拟方法，如多点约束单元法和梁元法等。 多点约束单元法(MPC)即采用特殊单元RBE2来模拟螺栓连接。在螺栓连接处，设置其中一节点为从节点(Dependent)，另外一个节点为主节点(Independent)。主从节点之间位移约束关系使得从节点跟随主节点位移变化。比例因子选为1,使从节点和主节点位移变化协调一致，从而模拟实际工作状态下，螺栓对法兰的连接紧固作用。 梁元法模拟即采用两节点梁单元Beam，其能承受拉伸、剪切、扭转。通过参数设置，使梁元与螺栓几何属性一致。 本文分别用算例来说明这两种方法的可行性。 2.1 几何模型 如图1所示组合装配体，底部约束。两圆筒连接法兰通过8颗螺栓固定。端面受联合载荷作用。
基于有限元的电磁场仿真与数值计算介绍
鼠笼异步电动机磁场的有限元分析 摘要 鼠笼异步电动机具有结构简单、价格低廉、运行可靠、效率较高、维修方便等一系列的优点，在国民经济中得到广泛的应用。工业、农业、交通运输、国防工程以及日常生活中都大量使用鼠笼异步电动机。随着大功率电子技术的发展，异步电动机变频调速得到越来越广泛的应用，使得鼠笼异步电动机在一些高性能传动领域也得到使用。 鼠笼异步电动机可靠性高，但由于种种原因，其故障仍时有发生。由于电动机结构设计不合理，制造时存在缺陷，是造成故障的原因之一。对电机内部的电磁场进行正确的磁路分析，是电机设计不可或缺的步骤。利用有限元法对电机内部磁场进行数值分析，可以保证磁路分析的准确性。本文利用Ansys Maxwell软件，建立了鼠笼式异步电机的物理模型，并结合数学模型和边界条件，完成了对鼠笼式异步电动机的磁场仿真，得到了物理模型剖分图，磁力线和磁通分布图，为电机的进一步设计研究提供了依据。 关键词：Ansys Maxwell；鼠笼式异步电机；有限元分析
一、前言 当电机运行时，在它的内部空间，包括铜与铁所占的空间区域，存在着电磁场，这个电磁场是由定、转子电流所产生的。电机中电磁场在不同媒介中的分布、变化及与电流的交链情况，决定了电机的运行状态与性能。因此，研究电机中的电磁场对分析和设计电机具有重要的意义。 在对应用于交流传动的异步电机进行电磁场的分析计算时，传统的计算方法因建立在磁场简化和实验修正的经验参数的基础之上，其计算精度就往往不能满足要求。如果从电磁场的理论着手，研究场的分布，再根据课题的要求进行计算，就有可能得到满意的结果。电机电磁场的计算方法大致可以分为解析法、图解法、模拟法和数值计算法。数值解法是将所求电磁场的区域剖分成有限多的网格或单元，通过数学上的处理，建立以网格或单元上各节点的求解函数值为未知量的代数方程组。由于电子计算机的应用日益普遍，所以电机电磁场的数值解法得到了很大发展，它的适用范围超过了所有其它的解法，并能达到足够的精度。对于电机电磁场问题，常用的数值解法有差分法和有限元法两种。用有限元法时单元的剖分灵活性大，适用性强，解的精度高。因此我们采用有限元法对电机电磁场进行数值计算。 Maxwell2D 是一个功能强大、结果精确、易于使用的二维电磁场有限元分析软件。在这里，我们利用Ansys的Maxwell2D 有限元分析工具对一个三相四极电机进行有限元分析，构建鼠笼式异步电机电动机的物理模型，并结合电机的数学模型、边界条件进行磁场分析。
数值分析报告模拟试卷1,2,3
数值分析模拟试卷1 一、填空（共30分，每空3分） 1 设??? ? ??-=1511A ，则A 的谱半径=)(a ρ______，A 的条件数)(1A cond =________. 2 设 Λ ,2,1,0,,53)(2==+=k kh x x x f k ，则 ],,[21++n n n x x x f ＝________, ],,[321+++n n n n x x x x f ，＝________. 3 设?????≤≤-++≤≤+=2 1,121 0,)(232 3x cx bx x x x x x S ，是以0，1，2为节点的三次样条函数，则 b=________,c=________. 4 设∞ =0)]([k k x q 是区间［0，1］上权函数为x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x q ，则 ?=1 )(dx x xq k ________，=)(2x q ________. 5 设 ?? ?? ? ?????=11001a a a a A ，当 ∈ a ________时，必有分解式 ，其中L 为下三角阵，当其对角线元素 )3,2,1(=i L ii 满足条件________时，这种分解是唯一的.
二、（14分）设4 9,1,41,)(2102 3 === =x x x x x f , （1）试求)(x f 在]4 9,41[上的三次Hermite 插值多项式)(x H 使满足 2,1,0),()(==i x f x H i i ，)()(11x f x H '='. （2）写出余项)()()(x H x f x R -=的表达式. 三、（14分）设有解方程0cos 2312=+-x x 的迭代公式为n n x x cos 3 2 41+ =+， （1） 证明R x ∈?0均有?∞ →=x x n x lim （? x 为方程的根）； （2） 取40=x ，用此迭代法求方程根的近似值，误差不超过 ，列出各次迭代值； （3）此迭代的收敛阶是多少？证明你的结论.
数值分析模拟试题
数值分析模拟试题 一、填空题（每小题3分，共30分） 1、已知近似值* 2.4560x =是由真值x 经四舍五入得到，则相对误差限为 。 2 、为减少舍入误差的影响，应将10改写成 。 3、设(1,1,2,3)T x =-，则12_______,_______,_______x x x ∞===。 4、设1123A -??=????，则1________,________F A A ==，A 的谱半径()A ρ=。 5、用Gauss-Seidel 迭代法解方程组1212423 x ax ax x +=??+=-?，其中a 为实数，则该方法收敛的充要 条件是a 满足 。 6、迭代法12213k k k x x x +=+收敛于*x =，此迭代格式是 阶收敛的。 7、设01(),(),,()n l x l x l x 是以01,, ,n x x x 为节点的Lagrange 插值基函数，则0()n i i l x ==∑。 8、设3()321f x x x =++，则差商[0,1,2,3]_____,[0,1,2,3,4]_____f f ==。 9、数值积分的辛普森公式为()b a f x dx ≈?。 10、数值积分公式0()()n b k k a k f x dx A f x =≈∑?中，0n k k A ==∑。 二、设函数2()(3)x x a x ?=+-，由迭代公式1()k k x x ?+=产生的序列为{}k x ，试讨论 ⑴当a 为何值时，序列{}k x 收敛； ⑵当a 取何值时，收敛速度最快，并指出迭代法收敛的阶。（12分） 三、设4()[0,2]f x C ∈，且(0)2,(1)1,(2)0,'(1)0f f f f ==-==，试求函数()f x 的三次 插值多项式()P x ，并求余项表达式。（14分） 四、用矩阵的直接三角分解法（即LU 分解）解方程组Ax b =，其中
有限元与数值方法-讲稿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 σσ=-=
有限元分析软件比较分析
有限元分析软件 有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50 年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC 四个比较知名比较大的公司，其中ADINA、ABAQUS 在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件，ANSYS、MSC 进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析，但是除ADINA 以外其它三个必须与别的软件搭配进行迭代分析，唯一能做到真正流固耦合的软件只有ADINA。ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。 结构分析能力排名：ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名：ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名：ADINA、ANSYS、MSC、ABAQUS 性价比排名：最好的是ADINA，其次ABAQUS、再次ANSYS、最后MSC ABAQUS 软件与ANSYS 软件的对比分析： 1．在世界范围内的知名度：两种软件同为国际知名的有限元分析软件，在世界范围内具有各自广泛的用户群。ANSYS 软件在致力于线性分析的用户中具有很好的声誉，它在计算机资源的利用，用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题，其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。由于ANSYS 产品进入中国市场早于ABAQUS，并且在五年前ANSYS 的界面是当时最好的界面之一，所以在中国，ANSYS 软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北京办事处的成立，ABAQUS软件的用户数目和市场占有率正在大幅度和稳步提高，并可望在今后的几年内赶上和超过ANSYS。 2．应用领域：ANSYS 软件注重应用领域的拓展，目前已覆盖流体、电磁场和多物理场耦合等十分广泛的研究领域。ABAQUS 则集中于结构力学和相关领域研究，致力于解决该领域的深层次实际问题。 3．性价比：ANSYS 软件由于价格政策灵活，具有多种销售方案，在解决常规的
螺栓连接的有限元分析
1 概述螺栓是机载设备设计中常用的联接件之一。其具有结构简单, 拆装方便,调整容易等优点, 被广泛应用于航空、航天、汽车以及各种工程结构之中。在航空机载环境下，由于振动冲击的影响，设备往往产生较大的过载，对作为紧固件的螺栓带来强度高要求。螺栓是否满足强度要求，关系到机载设备的稳定性和安全性。 传统力学的解析方法对螺栓进行强度校核，主要是运用力的分解和平移原理，解力学平衡方程，借助理论和经验公式，理想化和公式化。没有考虑到连接部件整体性、力的传递途径、部件的局部细节( 如应力集中、应力分布) 等等。通过有限元法，整体建模，局部细化，可以弥补传统力学解析的缺陷。用有限元分析软件 MSC.Patran/MSC.Nastran 提供的特殊单元来模拟螺栓连接，过程更方便，计算更精确，结果更可靠。因此，有限元在螺栓强度校核中的应用越来越广泛。 2 有限元模型的建立 对于螺栓的模拟，有多种模拟方法，如多点约束单元法和梁元法等。 多点约束单元法(MPC)即采用特殊单元RBE2来模拟螺栓连接。在螺栓连接处，设置其中一节点为从节点(Dependent) ，另外一个节点为主节点(Independent) 。主从节点之间位移约束关系使得从节点跟随主节点位移变化。比例因子选为1, 使从节点和主节点位移变化协调一致，从而模拟实际工作状态下，螺栓对法兰的连接紧固作用。 梁元法模拟即采用两节点梁单元Beam其能承受拉伸、剪切、扭转。通过参数设置，使梁元与螺栓几何属性一致。 本文分别用算例来说明这两种方法的可行性。 2.1 几何模型 如图 1 所示组合装配体，底部约束。两圆筒连接法兰通过8 颗螺栓固定。端面受联合载荷作用。
工程数值方法与有限元分析
工程数值方法与有限元分析 （机械工程学院机械类专业） 课程号： 周学时：4 学分：3 课程类别： 预修课程：高等数学，线性代数，力学基础课 面向对象：机械类专业学生 教学方式：多媒体教学 教学目的和教学要求： 在科学研究与工程技术中，经常遇到数学模型的求解问题。然而在许多情况下，要获得模型问题的准确解往往是十分困难的，甚至是不可能的。因此，研究各种数学问题的近似解法非常必要。计算方法是一门与计算机应用密切结合的实用性很强的课程，它专门研究各种数学问题的一类近似解法，从一组原始数据出发，按照确定的运算规则进行有限步运算，最终获得问题的数值形式且满足精度要求的近似解。 通过对《计算方法》的学习，掌握数值计算的基本概念和基本理论，深入理解方法的设计原理与处理问题的技巧，重视误差分析与收敛性、数值稳定性，注重利用计算机进行科学计算能力的培养,并熟练掌握Matlab 软件,会用Matlab实现各种计算方法。 在此基础上进一步学习数值计算的集大成者-有限元方法, 了解有限元方法的基础知识及其在机械、机械电子领域中的应用，掌握有限元方法的基本原理与分析过程，包括静力学、动力学、非线性力学、热场、电磁场等的建模及分析。学生可使用有限元软件进行机械零件及系统的实例分析，并对分析结果进行评价，指导和优化机械零件及系统的设计。本课程面向机械电子专业及机械类相关专业的高年级本科生 课程简介： 内容主要包括：计算机上常用的数值计算方法以及有关的基本概念与理论，主要有误差、非线性方程求根、线性代数方程组的解法、插值与拟合、数值微分与数值积分、常微分方程初值问题的数值解法。并且算法面向计算机，注重培养学生运用计算机进行科学计算解决工程问题的能力。并熟练掌握Matlab 软件，会用Matlab实现各种计算方法。 有限元的分析与建模是一个机械工程师必须掌握的方法和技能。本课程为机械类专业的高年级学生核心课，使学生了解有限元方法的基本概念和基本理论，掌握有限元分析的基本处理方法，熟悉常用有限元分析软件在实际工程中的应用，最终培养学生在机械设计、机电系统设计中能有效的应用有限元方法。 主要内容及学时分配： 每周4学时，共16周 主要内容： ( O ) 绪论1学时 (一)误差2学时
有限元分析在结构分析和计算机仿真中的应用
第20卷增刊重庆交通学院学报2001年11月VoI.20Sup.JOURNAL OF CHONGOING JIAOTONG UNIVERSITY NoV.， ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 2001文章编号：1001-716 （2001）S0-0124-03 有限元分析在结构分析和计算机仿真中的应用" 韩西，钟厉，李博 （重庆交通学院结构工程部级重点实验室，重庆400074） 摘要：简要论述了有限元分析方法在结构分析和计算机仿真的发展趋势和应用情况. 关键词：有限元分析；结构分析；计算机仿真 中图分类号：TU311.41文献标识码：B 自1943年数学家Courant第一次尝试应用定义在三角形区域上的分片连续函数的最小位能原理求解St.Venant扭转问题以来，许多数学家、物理学家和工程师由于各种原因都涉及过有限单元的概念.但由于即使一个小规模的工程问题，用有限元分析都将产生较大的计算工作量.直到1960年后，随着计算机技术的发展，有限元分析这门特别依赖数值计算的学科才真正进入了飞速发展阶段.到目前为止，有限元法已成为最强有力的数值分析方法之一，在固体力学、流体力学、机械工程、土木工程、电气工程等领域得到了广泛的应用.由于其所涉及问题和算法基本上全部来源于工程实际、应用于工程中，其解决工程实际问题的能力愈来愈强.由于计算机技术作为有限元分析的计算平台和应用支撑工具，故有限元分析成为CAE（Computer Aided Engineering，计算机辅助工程）这一学科类的主要研究内容.与此同时，由于有限元分析所建模型具有和实际结构相对应的几何、材料、力学特性，对实际结构具有“真实”的模拟特性，和单纯的几何仿真有着本质的区别，所以可以称之为“真实的仿真”（ReaIity SimuIation），可以想象，其模型和计算的数据量将比单纯的几何仿真要大得多，当前，计算机并行多处理器技术正迅猛发展，如SGI OONU-MA体系使计算能力达到工程应用水平，极大地促进了有限元分析计算的发展［2］. 1现状与发展趋势 1.1现状 1956年，Tuner，CIough等人将刚架位移法推广应用于弹性力学平面问题，并用于分析飞机结构，这是现代有限元法第一次成功的尝试.他们第一次给出了用三角形单元求解平面应力问题的正确解答，其研究工作打开了利用计算机求解复杂平面问题的新局面.1963~1964年，BesseIing、MeIosh 和Jones等证明有限单元法是基于变分原理的Ritz 法的另一种形式，从而使Ritz分析的所有理论基础都适用于有限元法，确认了有限元法是处理连续介质问题的一种普遍方法. 几十年来，有限元法的应用已由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定问题、动力问题和波动问题；分析的对象从弹性材料扩展到塑性、粘塑性和复合材料等；从固体力学扩展到流体力学、传热学等连续介质力学领域.在工程分析中的作用已从分析比较扩展到优化设计并和CAD（计算机辅助设计）结合越来越紧密. 有限元分析理论的逐步成熟主要经历了60年代的探索发展时期，70~80年代的独立发展专家应用时期和90年代与CAD相辅相成的共同发展、推广使用时期. 有限元分析作为一种强有力的数值分析方法，在结构分析和仿真计算中有着极大的应用价值.目前，结构仿真中的静力分析、动力分析、稳定性计 "收稿日期：2001-03-19 作者简介：韩西（1964-），男，重庆人，工学博士，副教授，主要从事振动工程、结构损伤识别、结构动力及计算机仿真分析研究.
2014-2015数值分析考试试题卷
太原科技大学硕士研究生 2014/2015学年第1学期《数值分析》课程试卷 一、填空题（每空4分，共32分） 1、设?????≤≤-++<≤+=2 1,1321 0,)(2 323x x bx x x x x x s 是以0,1,2为节点三次样条函数，则b=__-2___ 2、解线性方程组12312312388 92688 x x x x x x x x x -++=-?? -+=??-+-=? 的Jacobi 迭代格式（分量形式）为 ?? ???+--=++-=++=+++)(2)(1)1(3) (3)(1)1(2) (3)(2)1(1882/)96(88k k k k k k k k k x x x x x x x x x ，其相应的迭代矩阵为??????????-0812/102/9810。 3、方程03 =-a x 的牛顿法的迭代格式为__3 12 3k k k k x a x x x +-=-__________，其收敛的阶为 2 。 4、已知数x 的近似值0.937具有三位有效数字，则x 的相对误差限是310534.0-? 解：x 1≈0.937， 31102 1 )(-?≤ x ε 3 31111 10(x )2 (x )0.53410x 0.937 r εε--?=≤=? 5、用列主元高斯消去法解线性方程组 ??? ??=--=++=++2333220221 321321x x x x x x x x 作第1次消元后的第2，3个方程分别为? ? ?=+--=-5.35.125 .15.03232x x x x 6、设???? ??-=3211A ，则=∞)(A Cond __4____.
滚压有限元数值分析
滚压模型有限元分析 姓名：黄肖 学号：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,进入如下界面，置材料属性。 分别设置工件和滚针的材料属性：弹性模量和泊松比，屈服强度和屈服应力，
ansys有限元计算
1.1 课程设计的意义、目的 1）ANSYS模态分析用于确定设计结构或机器部件的振动特性（固有频率和振型），即结构的固有频率和振型，他们是承受动态载荷的重要参数，也可作为其他动力学分析的起点，是进行谱分析或模态叠加法普响应分析或瞬态动力学分析所必需的前期分析过程。模态分析在动力学分析过程中必不可少的一个步骤，在谱响应分析、瞬态动力学分析的分析过程中均要求先进行模态分析才能进行其他步骤。 2）根据课堂讲授内容，学生做相应的自主练习，消化课堂所讲解的内容；通过调试典型例题或习题积累调试ANSYS程序的经验；通过完成课程设计中中的编程题，逐渐培养学生的编程能力、用ANSYS解决实际问题的能力。 1.2课程设计研究的内容 求解外受两端压力带孔薄板的系统或局部的位移、应变、应力。 ANSYS详细设计步骤 1问题分析 如图所示，E=30e6,两端压力100，中心孔内线压分布力500向外。中心孔直径为5。板厚为1。
基于ANSYS分析的简要步骤 (1)启动ANSYS，进入ANSYS界面。 (2)定义工作文件名 GUI : Utility Menu > File > Change Jobname 单击Utility Menu菜单下File中的Change Jobname按钮，会弹出Change Jobname对话框， 输入有限元分析作为工作文件名，单击Ok。 (3)定义分析标题 GUI:Utility Menu > File > Change Title 在弹出的对话框中，输入cui作为分析标题，单击OK。 (4)重新显示 GUI: Utility Menu>Plot>Replot 单击该按钮后，所命名的分析标题和工作文件名会出现在ANSYS中。 (5)选择分析类型 在弹出的对话框中，选择分析类型，由于此例属于结构分析，故选择Structural 这一项，单击OK。 (6)定义单元类型 GUI:Main Menu > Preprocessor > Element Type > Add/Edit/Delete 单击弹出对话框中的Add按钮，弹出单元库对话框，在材料的单元库中选择Plane82单元。即在左侧的窗口中选取Solid单元，在右侧选择8节点的82单元，然后单击OK。
华南理工大学数值分析试题-14年下-C
华南理工大学研究生课程考试 《数值分析》试卷C （2015年1月9日） 1. 考前请将密封线内各项信息填写清楚； 所有答案请按要求填写在本试卷上； 课程代码：S0003004； 4. 考试形式：闭卷； 5. 考生类别：硕士研究生； 本试卷共八大题，满分100分，考试时间为150分钟。 一、（12分）解答下列问题： 1）设近似值0x >，x 的相对误差为δ，试证明ln x 的绝对误差近似为δ。 2）利用秦九韶算法求多项式 542()681p x x x x x =-+-+ 在3x =时的值（须写出计算形式），并统计乘法次数。 （12分）解答下列问题： 1）设()235f x x =+，求[]0,1,2f 和[]0,1,2,3f 。 2）利用插值方法推导出恒等式： 33220,0[]j j i i x j i x i j =≠=-=-∑∏ 。
（1）设{}∞ =0)(k k x q 是区间[]1,0上带权1=ρ而最高次项系数为1的正交多项式族，其中1)(0=x q ，求1()q x 和2()q x 。 （2）求形如2y a bx =+的经验公式，使它与下列数据拟合： 四、（14分）对积分()10I f x dx = ?，试 （1）构造一个以012113,,424 x x x ===为节点的插值型求积公式； （2）指出所构造公式的代数精度； （3）用所得数值求积公式计算积分1 203x dx ?的精确值； （4）指出所得公式与一般的Newton-Cotes 型公式在形式上的重要区别。
（1）设?? ????=4321A ,计算1A 、()Cond A ∞和()A ρ。 （2）用列主元Gauss 消去法解方程组： 12312315410030.112x x x ????????????=????????????-?????? 六、（13分）对2阶线性方程组 11112212112222 a x a x b a x a x b +=??+=? （11220a a ≠ ） （1）证明求解此方程组的Jacobi 迭代与Gauss-Seidel 迭代同时收敛或同时发散； （2）当同时收敛时，试比较它们的收敛速度。
浅析拉格朗日有限元数值计算法
中国矿业大学2011级硕士研究生课程考试试卷 考试科目岩土工程数值计算法 考试时间2011.11.27 学生姓名夏明 学号ZS11020068 所在院系矿业工程学院 任课教师徐志伟 中国矿业大学研究生院培养管理处印制
对论文《基于拉格朗日差分法的露天边坡稳定性研 究》中数值计算法浅析 《基于拉格朗日差分法的露天边坡稳定性研究》一文以易门铜矿露天开采境界优化方案下的边坡为工程背景，通过岩体构造调查，质量分类，室内力学试验，力学参数的工程处理。建立铜厂露天开采边坡的三维地质模型，采用DIMINE 数字矿山软件的耦合集成技术—四面体网格化将地质模型转换为力学模型，应用基于拉格朗日法的有限差分(FLAC3D)大变形方法对铜厂露天开采边坡的稳定性进行了数值模拟，分析了基于强度折减理论计算出的边坡安全系数以及基于莫尔库伦屈服准则的边坡开采后的位移、应力等的变化状况，得出了易门铜厂露天矿露天境界优化方案下的边坡的稳定性状况。 岩质边坡稳定性评价的方法分主要有：极限平衡法、数值模拟计算、地质力学物理模拟试验和其它新方法。随着计算机技术和计算方法的发展，复杂的工程问题可以采用离散化的数值计算方法并借助计算机得到满足工程要求的数值解，数值模拟技术是现代工程学形成和发展的重要动力之一。通过计算模拟，可以模拟并得到模拟体内部的应力—应变关系，再现其变形甚至破坏过程及其机制。在岩土工程数值分析中最常用的数值方法有有限元法、离散元法、边界元法等。 拉格朗日差分法((FLAC法)源于流体力学。它首先是Cundail在80年代提出来的，其基本原理类似于离散元法，但它却能像有限元那样适用于多种材料模式与边界条件的非规则区域的连续问题求解。在求解的过程中，FLAC又采用了离散元的动态松弛法，不需要求解大型联立方程，便于在微机上实现。另一方面，同以往的差分分析相比，FLAC在以下几个方面做了较大的改进和发展:它不但能处理一般的大变形问题，而且能模拟岩体沿某一弱面产生的滑移变形。一般有限单元法可以用来解决材料非线性问题，但对于大变形的几何非线性问题，有限单元法和边界元法都无能为力。拉格朗日法是分析非线性大变形问题的数值方法，它依然遵循连续介质的假设，基于拖带坐标系的基本原理。用差分法或按时步显式迭代求解，不但可以解决几何非线性，也能解决材料非线性问题。 下面是原文的部分描述： 一、露天矿边坡三维地质模型的建立 1.1计算几何模型范围及模拟方案设计 a模拟计算方案设计 决定边坡稳定性的因素很多，其中包括岩体强度，结构特征，水理性质和边坡几何参数等。根据现场工程地质调查，力学试验参数，对露天境界边坡进行修改，修改方案同第五章简化实体模型的建立。
数值分析2007第二学期期末考试试题与答案(A)
期末考试试卷（A 卷） 2007学年第二学期 考试科目： 数值分析 考试时间：120 分钟 学号 姓名 年级专业 一、判断题（每小题2分，共10分） 1. 用计算机求 1000 1000 1 1 n n =∑时，应按照n 从小到大的顺序相加。 （ ） 2. 为了减少误差,进行计算。 （ ） 3. 用数值微分公式中求导数值时，步长越小计算就越精确。 （ ） 4. 采用龙格－库塔法求解常微分方程的初值问题时，公式阶数越高，数值解越精确。（ ） 5. 用迭代法解线性方程组时，迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有 关，与常数项无关。 （ ） 二、填空题（每空2分，共36分） 1. 已知数a 的有效数为0.01，则它的绝对误差限为________，相对误差限为_________. 2. 设1010021,5,1301A x -????????=-=-????????-???? 则1A =_____，2x =______，Ax ∞ =_____. 3. 已知5 3 ()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= . 4. 为使求积公式 1 1231 ()((0)f x dx A f A f A f -≈++? 的代数精度尽量高，应使1A = ，2A = ，3A = ，此时公式具有 次的代数精度。 5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 . 6. 用迭代法解线性方程组AX B =时，使迭代公式(1) ()(0,1,2,)k k X MX N k +=+= 产 生的向量序列{ }() k X 收敛的充分必要条件是 . 7. 使用消元法解线性方程组AX B =时，系数矩阵A 可以分解为下三角矩阵L 和上三角矩
有限元数值分析课程论文
《有限元数值分析》课程论文要求 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）
电磁仿真算中的有限元法
1电磁仿真算法中的有限元法 1.1常规的电磁计算方法简介 从上世纪50年代以来,伴随着计算机技术的进步,电磁仿真算法也蓬勃发展起来,这其中主要包括:单矩法、矩量法和有限元法等属于频域技术的算法; 传输线矩阵法、时域积分方程法以及时域有限差分法等属于时域技术的算法。除了这些以外, 还有属于高频技术的集合衍射理论等。本文根据国内外计算电磁学的发展状况,对日常生活中比较常用的电磁计算方法做了介绍,并对有限元法做了重点说明。 ⑴矩量法 矩量法属于电磁场的数值计算方法中频域技术的一种, 它的基本原理是利用把待解的微积分方程转化成的算子方程, 然后将由一组线性组合表示的待求函数代入第一步中的算子方程, 然后将算子方程转化成矩阵方程, 最后再通过计算机进行大量的数值计算从而得到数值结果。该方法在求解非均勻和不规则形状对象时,面很广,但会生成病态矩阵,所以会在一定程度上受到限制。矩量法的特点就是适用于求解微积分方程, 并且求解方法统一简单。但缺点就是会占用大量计算机内存,影响计算速度。 (2)单矩法 单矩法是一种解析方法和数值方法相结合的混合数值算法法,该方法的关键在于,如何合理的选择一个球面最小的半径,使得能够将分析对象的结构全部包含在内,以便将内外场进行隔离。外边的散射场单独使用其他函数表示,而包围的内部区域使用有限元法亥姆赫兹(Helmholtz)方程。此方法对于计算复杂形体乃至复杂埋入体内的电磁散射是种极为有效的手段。 (3)时域有限差分法 时域有限差分法(FDTD）近几年来越来越受到各方的重视, 因为一方面它处理庞大的电磁福射系统方面和复杂结构的散射体时很突出,另外一方面则在于它不是传统的频域算法, 它是种时域算法, 直接依靠时间变量求解麦克斯韦方程组,可以在有限的时间和体积内对场进行数据抽样, 这样同时也能够保证介质边界
数值分析模拟试题
1、 方程组中，，则求解方程组的Jacobi 迭代与Gauss-Seidel 迭代均收敛的a 的范围是___________。 2、，则A 的LDL T 分解中，。 3、，则__________，_______________. 4、已 知，则用复合梯形公式计算求 得，用三点式求得____________. 5、，则_________ ，三点高斯求积公式______________. 6设* 2.40315x =是真值 2.40194x =的近似值，则* x 有________位有效数字。 7 3()1,[0,1,2,3]f x x x f =+-=设 则差商(均差)_____________，[0,1,2,3,4]f =________________。 8 求方程()x f x =根的牛顿迭代格式是__________________。 9.梯形求积公式和复化梯形公式都是插值型求积公式_____（对或错）。 10.牛顿—柯特斯求积公式的系数和()0n n k k C ==∑__________________。 11.用二次拉格朗日插值多项式2()sin0.34L x 计算的值。插值节点和相应的函数值是（0，0），（0.30，0.2955），（0.40，0.3894）。 12.用二分法求方程3()10[1.0,1.5]f x x x =--=在 区间内的一个根，误差限 210ε-=。 13.用列主元消去法解线性方程组 1231231 232346,3525,433032.x x x x x x x x x ++=??++=??++=? 14. 确定求积公式
012()()(0)()h h f x dx A f h A f A f h -≈-++? 。 中待定参数i A 的值(0,1,2)i =，使求积公式的代数精度尽量高；并指出此时求积公式的代数精度。 15、 试求使求积公式的代数精度 尽量高，并求其代数精度。 16.证明区间[a,b]上带权()x ρ的正交多项式(),1,2,n P x n = 的n 个根都是单根，且位于区间(a,b)内。 17.设()()[,],max ()n n a x b f x C a b M f x ≤≤∈=，若取 21cos ,1,2,,222k a b a b k x k n n +--=+= 作节点，证明Lagrange 插值余项有估计式21()max ()!2n n n a x b M b a R x n -≤≤-≤ 18用n=10的复化梯形公式计算时， （1）试用余项估计其误差 （2）用n=10的复化梯形公式计算出该积分的近似值。 19已知方程组AX ＝f,其中 （1）列出Jacobi 迭代法和Gauss-Seidel 迭代法的分量形式。 （2）求出Jacobi 迭代矩阵的谱半径，SOR 迭代法的最佳松弛参数 和SOR 法 的谱半径(可直接用现有结论) 20试确定常数A ，B ，C 和，使得数值积分公式 有尽可能高的代数精度。试问所得的数值积分公式代数精度是多少？ 21证明方程=)(x f x 2-x -3=0在区间(2,3)内有且仅有一个根,并用迭代法求方程在区间(2,3)内的根,精确到小数点后4位。 22设f (1)=2，f (3)=4，f (4)=6，用拉格朗日插值法求f (x )的二次插值多项式P 2(x )，并求f (2)的近似值。