(完整版)北航有限元分析与应用试题库
1． 证明3结点三角形单元的插值函数满足ij j i i y x N δ=),(，及1=++m j i N N N 。 2． 图示3三结点三角形单元，厚度为t ，弹性模量为E ，泊桑比ν＝0。试求：插
值函数矩阵N ，应变矩阵B ，应力矩阵S ，单位刚度矩阵K e
。
3． 以平面问题常应变三角形单元为例，证明单元刚度矩阵的任何一行（或列）元
素的总和为零。
4． 试证明面积坐标与直角坐标满足下列转换关系。
m m j j i i l x l x l x x ++= m m j j i i l y l y l y y ++=
5． 写出题5图所示三角形单元的插值函数Ni ，Nj ，Nm 以及应变矩阵B 。
6． 题5图中单元在jm 边作用有线性分布的面载荷（x 方向），试求结点载荷问题。 7． 证明常应变三角形单元发生在刚体位移时，单元中将不产生应力。
8． 求图示二次三角形单元在1 4 2边作用有均布侧压g 时的等效结点载荷，假设
结点坐标已知，单元厚度为t 。
9． 验证用面积坐标给出二次（三角形）单元的插值函数的N 1～N 6满足∑==6
~11i i N
10． 二维单元在xy 坐标平面内平移到不同位置，单元刚度矩阵相同吗?在平面内旋
转时怎样？单元旋转180o 后单元刚度矩阵与原来的相同吗？单元作上述变化时，应力矩阵S 如阿变化？
11． 图中两个三角形单元组成平行四边形，已知单元○
1按局部编码i ，j ，m 的单元刚度矩阵K ○
1和应力矩阵S ○1为 K ○1=?????????
???????????---------5.55.15.95.15.15.135.45.75.45.134012016626608 S ○1=????
??????-----5.15.05.15.102103000030300 按图中单元○
2的局部编码写出K ○2，S ○2。 12． 图示为二次四边形单元，试计算x N ??1和y N ??2在自然坐标为(1/2,1/2)的点
Q 的数值（因为单元的边是直线，可用4个结点定义单元的几何形状）。
13． 图示为一次三角形单元，试计算x N ??4和y N ??4在点P(1.5 ,2.0)的数值。
14． 垂直悬挂的等截面直杆受自重作用，截面积为A ，长度为l ，质量密度为ρ。
如图一维杆单元求解杆内的应力分布，问应采用多少结点的单元？在什么位置有限元结果可以达到解析解的精度？给出它们的数值。
15．有中心椭球孔的矩形板，两侧边受线性分布的侧压P，如图所示如何利用对称面条件减少求解的工作量，并画出计算模型，列出计算步骤。
有限元分析与应用详细例题
《有限元分析与应用》详细例题 试题1：图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比 较： 1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； 2）分别采用不同数量的三节点常应变单元计算； 3）当选常应变三角单元时，分别采用不同划分方案计算。 一．问题描述及数学建模 无限长的刚性地基上的三角形大坝受齐顶的水压作用可看作一个平面问题，简化为平面三角形受力问题，把无限长的地基看着平面三角形的底边受固定支座约束的作用，受力面的受力简化为受均布载荷的作用。 二．建模及计算过程 1. 分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算 下面简述三节点常应变单元有限元建模过程（其他类型的建模过程类似）： 1.1进入ANSYS 【开始】→【程序】→ANSYS 10.0→ANSYS Product Launcher →change the working directory →Job Name: shiti1→Run 1.2设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK 1.3选择单元类型 单元是三节点常应变单元，可以用4节点退化表示。 ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4 node 42 →OK (back to Element Types window)→Options… →select K3: Plane Strain→OK→Close (the Element Type window) 1.4定义材料参数
北航2010-2011年研究生数值分析期末模拟试卷1-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,)(2 323x 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 ________，=)(2 x 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）此迭代的收敛阶是多少？证明你的结论. 四、(16分) 试确定常数A ，B ，C 和，使得数值积分公式 有尽可能高的代数精度. 试问所得的数值积分公式代数精度是多少？它是否为Gauss 型的？
2020年度北航材料力学试题
2003北京航天大学材料力学试题 【一】、选择题,从所给答案中选择一个正确答案(本题共10分,每小题5分) 1、在下列四种工程材料中,_________不可应用各向同性假设。 A.铸铁；B.玻璃；C.松木；D.铸铜。 2、设图示任意平面图形对该平面内的1Z 、2Z 、3Z 轴的惯性矩分别为1I 、2I 、3I ,对点的极惯性矩为p I ,在列关系式中,_________是正确的。(1Z 轴垂直于3Z 轴) A.2I =1I ＋3I ； B.p I =1I ＋2I ； C.p I =1I ＋3I ； D.p I =2I ＋3I 。 题一.2图 【二】、填空题(本题共20分,每小题5分) 1、如题二(1)图a 所示圆轴承受扭距T,在沿轴线成45°处贴有电阻应变片1R 和2R 。将1R 和2R 接到题二(1)图b 所示电桥上,电桥中3R 和4R 是阻值相同的固定电阻。设电阻应变仪上的读数应变为ε,则应变片1R 的应变1ε=________。 题一.1图 2、杆1、2和3的横截面积及长度均相等,其材料的应力应变去向如题二(2)图所示。则______强度最高。 ______刚度最高。______塑性最好。
3、已知各向同性线弹性材料常数为E ,泊松比为μ,材料内某点主应变 1σ>2σ>3σ=0。则此点第三主应变3ε(1ε>2ε>3ε)的大小为______，此点最大切应变m ax γ的大小为______。 4、影响构件疲劳的主要因素包括___________________，___________________，和___________________， 【三】、(15分)画出题三所示梁的剪力图和弯矩图。 题三图 【四】、(20分)图示桁架两杆材料相同,拉压许用应力相等,为[]σ。两杆夹角为α,2杆长为l 。节点B 作用向下的载荷F 。不考虑稳定条件。1.设计两杆的横截面积1A 和2A ；2.将1A 、2A 、α作为可设计量,保持结构其余参数不变,求结构重量最轻时的α值。
有限元分析试题(同济)
同济大学本科课程期终考试统一命题纸A卷 2007—2008学年第二学期 一．是非题（认为该题正确，在括号中打√；该题错误，在括号中打×。）(每小题2分) （1）用加权余量法求解微分方程，其权函数V和场函数u的选择没有任何限制。（）（2）四结点四边形等参单元的位移插值函数是坐标x、y的一次函数。（）（3）在三角形单元中，其面积坐标的值与三结点三角形单元的结点形函数值相等。（）（4）二维弹性力学问题的有限元法求解，其收敛准则要求试探位移函数C1连续。（）（5）有限元位移法求得的应力结果通常比应变结果精度低。（）（6）等参单元中Jacobi行列式的值不能等于零。（）（7）在位移型有限元中，单元交界面上的应力是严格满足平衡条件的。（）（8）四边形单元的Jacobi行列式是常数。（）（9）利用高斯点的应力进行应力精度的改善时，可以采用与位移插值函数不同结点的形函数进行应力插值。（）（10）一维变带宽存储通常比二维等带宽存储更节省存储量。（）二．单项选择题（共20分，每小题2分） 1 在加权余量法中，若简单地利用近似解的试探函数序列作为权函数，这类方法称为 ________________。 （A）配点法（B）子域法（C）伽辽金法 2 等参变换是指单元坐标变换和函数插值采用______的结点和______的插值函数。 （A）不相同，不相同（B）相同，相同（C）相同，不相同（D）不相同，相同 3 有限元位移模式中，广义坐标的个数应与___________相等。 （A）单元结点个数（B）单元结点自由度数（C）场变量个数 4 采用位移元计算得到应力近似解与精确解相比较，一般___________。 （A）近似解总小于精确解（B）近似解总大于精确解（C）近似解在精确解上下震荡（D）没有规律 5 如果出现在泛函中场函数的最高阶导数是m阶，单元的完备性是指试探函数必须至少 是______完全多项式。 （A）m-1次（B）m次（C）2m-1次 6 与高斯消去法相比，高斯约当消去法将系数矩阵化成了_________形式，因此，不用进 行回代计算。 （A）上三角矩阵（B）下三角矩阵（C）对角矩阵 7 对称荷载在对称面上引起的________________分量为零。 （A）对称应力（B）反对称应力（C）对称位移（D）反对称位移 8 对分析物体划分好单元后，__________会对刚度矩阵的半带宽产生影响。 （A）单元编号（B）单元组集次序（C）结点编号 9 n个积分点的高斯积分的精度可达到______阶。 （A）n-1 （B）n（C）2n-1 （D）2n 10 引入位移边界条件是为了消除有限元整体刚度矩阵K的__________。 （A）对称性（B）稀疏性（C）奇异性 三．简答题（共20分，每题5分）
北航材料力学试题往年试卷-01-答案
班号 学号 姓名 成绩 《 材 料 力 学 A 》期末试卷 一、选择题 (单选，共12分，每小题3分) 1、下列说法正确的是： D 。 A 、各向同性材料是指材料内部任意位置的力学性能都是相同的； B 、材料的刚度越大越不容易发生屈服变形； C 、塑性材料的强度极限高于脆性材料的强度极限； D 、脆性材料的单向压缩许用应力一般高于单向拉伸许用应力。 2、下列说法正确的是 C 。 A 、薄壁圆管的扭转切应力公式只适用于线弹性、各向同性材料； B 、任意截面形状的闭口薄壁杆横截面上的扭转切应力为常数； C 、开口薄壁杆的抗扭性能很差，对于受扭构件，一般不要采用开口薄壁杆； D 、圆轴扭转的刚度条件是圆轴的抗扭模量不能超过许用值。 3、下图所示两均质矩形截面等厚梁，材料相同、总长度相同、AB 段的长度与横截面形状相同。下面结论中正确的是： B 。 A 、两梁 B 截面的挠度和转角不相同； B 、两梁B 截面的挠度和转角相同； C 、两梁C 截面的挠度和转角相同； D 、图（a ）所示梁C 截面挠度和转角均小于图 （b ）所示梁C 截面的挠度和转角。 F 图（b ）
4、下图所示正方形截面杆件，横截面ABCD 上拉应力最大的点是 B 。 二、填空题（8分，每空1分） 1、由于截面急剧变化引起的应力局部增大现象，称为 应力集中 。 2、 几个载荷同时作用时产生的效果，等于各个载荷单独作用时产生的效果的总和， 这一原理 称为载荷叠加原理。 当构件的受力和变形满足 线弹性（物理线性） 条件与 小变形（几何线性）条件时， 载荷叠加原理方能适用。 3、在进行纯弯梁横截面上正应力分析时，除了运用静力学方程外，还运用了 几何 方程与 物理 方程，并根据梁的变形特点作了 平面 假设与 单向受力 假设。 F
北航数值分析报告第三次大作业
数值分析第三次大作业 一、算法的设计方案： （一）、总体方案设计： x y当作已知量代入题目给定的非线性方程组，求（1）解非线性方程组。将给定的(,) i i
得与(,)i i x y 相对应的数组t[i][j],u[i][j]。 （2）分片二次代数插值。通过分片二次代数插值运算，得到与数组t[11][21],u[11][21]]对应的数组z[11][21]，得到二元函数z=(,)i i f x y 。 （3）曲面拟合。利用x[i],y[j],z[11][21]建立二维函数表，再根据精度的要求选择适当k 值，并得到曲面拟合的系数矩阵C[r][s]。 （4）观察和(,)i i p x y 的逼近效果。观察逼近效果只需要重复上面（1）和（2）的过程，得到与新的插值节点(,)i i x y 对应的(,)i i f x y ，再与对应的(,)i i p x y 比较即可，这里求解 (,)i i p x y 可以直接使用（3）中的C[r][s]和k 。 （二）具体算法设计： （1）解非线性方程组 牛顿法解方程组()0F x =的解* x ，可采用如下算法： 1）在* x 附近选取(0) x D ∈，给定精度水平0ε>和最大迭代次数M 。 2）对于0,1, k M =执行 ① 计算() ()k F x 和()()k F x '。 ② 求解关于() k x ?的线性方程组 () ()()()()k k k F x x F x '?=- ③ 若() () k k x x ε∞∞ ?≤，则取*()k x x ≈，并停止计算；否则转④。 ④ 计算(1) ()()k k k x x x +=+?。 ⑤ 若k M <，则继续，否则，输出M 次迭代不成功的信息，并停止计算。 （2）分片双二次插值 给定已知数表以及需要插值的节点，进行分片二次插值的算法： 设已知数表中的点为： 00(0,1,,) (0,1,,)i j x x ih i n y y j j m τ=+=???=+=?? ，需要插值的节点为(,)x y 。 1) 根据(,)x y 选择插值节点(,)i j x y ： 若12h x x ≤+ 或12 n h x x ->-，插值节点对应取1i =或1i n =-，
北京科技大学有限元试题及答案
一 判断题（20分） （×）1. 节点的位置依赖于形态，而并不依赖于载荷的位置 （√）2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元 （×）3. 不能把梁单元、壳单元和实体单元混合在一起作成模型 （√）4. 四边形的平面单元尽可能作成接近正方形形状的单元 （×）5. 平面应变单元也好，平面应力单元也好，如果以单位厚来作模型化 处理的话会得到一样的答案 （×）6. 用有限元法不可以对运动的物体的结构进行静力分析 （√）7. 一般应力变化大的地方单元尺寸要划的小才好 （×）8. 所谓全约束只要将位移自由度约束住，而不必约束转动自由度 （√）9. 同一载荷作用下的结构，所给材料的弹性模量越大则变形值越小 （√）10一维变带宽存储通常比二维等带宽存储更节省存储量。 二、填空（20分） 1．平面应力问题与薄板弯曲问题的弹性体几何形状都是 薄板 ，但前者受力特点是： 平行于板面且沿厚度均布载荷作用 ，变形发生在板面内； 后者受力特点是： 垂直于板面 的力的作用，板将变成有弯有扭的曲面。 2．平面应力问题与平面应变问题都具有三个独立的应力分量： σx ，σy ，τxy ，三个独立的应变分量：εx ，εy ，γxy ，但对应的弹性体几何形状前者为 薄板 ，后者为 长柱体 。3．位移模式需反映 刚体位移 ，反映 常变形 ，满足 单元边界上位移连续 。 4．单元刚度矩阵的特点有：对称性 ， 奇异性 ，还可按节点分块。 5．轴对称问题单元形状为：三角形或四边形截面的空间环形单元 ，由于轴对称的特性，任意一点变形只发生在子午面上，因此可以作为 二 维问题处理。 6．等参数单元指的是：描述位移和描述坐标采用相同的形函数形式。等参数单元优点是：可以采用高阶次位移模式，能够模拟复杂几何边界，方便单元刚度矩阵和等效节点载荷的积分运算。 7．有限单元法首先求出的解是 节点位移 ，单元应力可由它求得，其计算公式为 {}{} [][]e D B σδ=。（用符号表示即可） 8．一个空间块体单元的节点有 3 个节点位移： u ，v ，w 9．变形体基本变量有位移应变应力 基本方程 平衡方程 物理方程 几何方程 10.实现有限元分析标准化和规范化的载体就是单元
有限元分析及应用大作业
有限元分析及应用大作业 作业要求： 1）个人按上机指南步骤至少选择习题中3个习题独立完成，并将计算结果上交； 也可根据自己科研工作给出计算实例。 2）以小组为单位完成有限元分析计算； 3）以小组为单位编写计算分析报告； 4）计算分析报告应包括以下部分： A、问题描述及数学建模； B、有限元建模（单元选择、结点布置及规模、网格划分方案、载荷及边界 条件处理、求解控制） C、计算结果及结果分析（位移分析、应力分析、正确性分析评判） D、多方案计算比较（结点规模增减对精度的影响分析、单元改变对精度的 影响分析、不同网格划分方案对结果的影响分析等） 题一：图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；（注意ANSYS中用四边形单元退化为三节点三角形单元） 2）分别采用不同数量的三节点常应变单元计算； 3）当选常应变三角单元时，分别采用不同划分方案计算。 解：1.建模： 由于大坝长度>>横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看作无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作
用于大坝上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力，满足平面应变问题的载荷条件。因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况P=98000-9800*Y；建立几何模型，进行求解；假设大坝的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3； 2：有限元建模过程： 2.1 进入ANSYS : 程序→ANSYS APDL 15.0 2.2设置计算类型: ANSYS Main Menu: Preferences →select Structural →OK 2.3选择单元类型: ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182，六节点三角形单元选择Solid Quad 8node 183)→OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window) 2.4定义材料参数： ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK 2.5生成几何模型： 生成特征点： ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS →依次输入四个点的坐标：input:1(0,0),2(10,0),3(1,5),4(0.45,5) →OK 生成坝体截面： ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点，1(0,0),2(6,0),3(0,10) →OK 2.6 网格划分： ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window) 2.7 模型施加约束： 给底边施加x和y方向的约束： ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK 给竖直边施加y方向的分布载荷： ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量；2) 在Result窗口中出现{X},写入所施加的载荷函数： 98000-9800*{Y}；3) File>Save(文件扩展名：func) →返回：Parameters →Functions →Read from file：将需要的.func文件打开，参数名取meng，它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边；OK →在下拉列表框中，选择：Existing table →OK →选择需要的载荷为meng参数名→OK 2.8 分析计算： ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load
北航数值分析大作业第二题
数值分析第二次大作业 史立峰 SY1505327
一、 方案 （1）利用循环结构将sin(0.50.2)() 1.5cos( 1.2)() {i j i j ij i j i j a +≠+==(i,j=1,2,……,10)进行赋值，得到需要变换的 矩阵A ； （2）然后，对矩阵A 利用Householder 矩阵进行相似变换，把A 化为上三角矩阵A (n-1)。 对A 拟上三角化，得到拟上三角矩阵A (n-1)，具体算法如下： 记A(1)=A ，并记A(r)的第r 列至第n 列的元素为()n r r j n i a r ij ,,1,;,,2,1) ( +==。 对于2,,2,1-=n r 执行 1. 若 ()n r r i a r ir ,,3,2) ( ++=全为零，则令A(r+1) =A(r),转5；否则转2。 2. 计算 () ∑+== n r i r ir r a d 1 2 )( ()( )r r r r r r r r r r d c a d a c ==-=++则取,0sgn ) (,1)(,1若 )(,12r r r r r r a c c h +-= 3. 令 () n T r nr r r r r r r r r R a a c a u ∈-=++) ()(,2)(,1,,,,0,,0 。 4. 计算 r r T r r h u A p /)(= r r r r h u A q /)(= r r T r r h u p t /= r r r r u t q -=ω T r r T r r r r p u u A A --=+ω)()1( 5. 继续。 （3）使用带双步位移的QR 方法计算矩阵A (n-1)的全部特征值，也是A 的全部特征值，具体算法如下： 1. 给定精度水平0>ε和迭代最大次数L 。 2. 记n n ij n a A A ?-==][) 1()1()1(，令n m k ==,1。
lhz有限元分析及应用-课程试卷
有限元分析与应用试题 1.有限元求解问题的主要思路是什么？并做简要介绍。 ● 将连续系统分割成有限个分区或单元（离散化） 离散化 将直杆划分成n 个有限段，有限段之间通过一个铰接点连接。两段之间的连接点称为节点，每个有限段称为单元。第i 个单元的长度为L i ，包含第i ，i+1个节点。 ● 用标准方法对每个单元提出一个近似解（单元分析） 单元分析 用单元节点位移表示单元内部位移－第i 个单元中的位移用所包含的结点位移来表示。 线性插值所得到的 第i 结点的位移 第i 结点的坐标 第i 个单元的 应变 ) ()(1i i i i i x x L u u u x u --+=+i u i x i i i i L u u dx du -== +1εi i i i i L u u E E )(1-= =+εσ
应力 内力 将所有单元按标准方法组合成一个与原有系统近似的系统（整体分析） 首先把外载荷集中到节点上： 把第i 单元和第i+1单元重量的一半，集中到第i+1结点上 建立结点的力平衡方程：对于第i+1结点，由力的平衡方程可得 （i=1，n-1） i i i i i L u u EA A N ) (1-= =+σ2 ) (11+++= -i i i i L L q N N ) (2 )()(11121++++++=---i i i i i i i i L L q L u u EA L u u EA
令 对于第n+1个结点，第n 个单元的内力与第n+1个结点上的外载荷平衡， 再加上约束条件 因此可以得到n+1个方程构成的方程组，可解出n+1个结点的位移。 有限元方法的基本思想和原理是“简单”而“朴素”的，在发展初期，许多学术权威对该方法的学术价值有所鄙视，国际著名刊物Journal of Applied Mechanics 许多年来拒绝刊登有关有限元方法的文章，其理由是没有新的科学实质。 现在完全不同了，由于有限元方法在科学研究和工程分析中的地位， 1 += i i i L L λ22 1)11(2)1(i i i i i i i L EA q u u u λλλ+=-++-++1()2 n n n n n n EA u u qL N A L σ+-== = EA qL u u n n n 221= +-+0 1=u
北航数值分析报告大作业第八题
北京航空航天大学 数值分析大作业八 学院名称自动化 专业方向控制工程 学号 学生姓名许阳 教师孙玉泉 日期2014 年11月26 日
一．题目 关于x , y , t , u , v , w 的方程组(A.3) ???? ?? ?=-+++=-+++=-+++=-+++79 .0sin 5.074.3cos 5.007.1cos sin 5.067.2cos 5.0y w v u t x w v u t y w v u t x w v u t (A.3) 以及关于z , t , u 的二维数表（见表A-1）确定了一个二元函数z =f (x , y )。 表A-1 二维数表 t z u 0 0.4 0.8 1.2 1.6 2 0 -0.5 -0.34 0.14 0.94 2.06 3.5 0.2 -0.42 -0.5 -0.26 0.3 1.18 2.38 0.4 -0.18 -0.5 -0.5 -0.18 0.46 1.42 0.6 0.22 -0.34 -0.58 -0.5 -0.1 0.62 0.8 0.78 -0.02 -0.5 -0.66 -0.5 -0.02 1.0 1.5 0.46 -0.26 -0.66 -0.74 -0.5 1. 试用数值方法求出f (x , y ) 在区域}5.15.0,8.00|), {≤≤≤≤=y x y x D （上的近似表达式 ∑∑===k i k j s r rs y x c y x p 00 ),( 要求p (x , y )以最小的k 值达到以下的精度 ∑∑==-≤-=10020 7210)],(),([i j i i i i y x p y x f σ 其中j y i x i i 05.05.0,08.0+==。 2. 计算),(),,(* ***j i j i y x p y x f (i =1,2,…,8 ; j =1,2,…,5) 的值，以观察p (x , y ) 逼 近f (x , y )的效果，其中j y i x j i 2.05.0,1.0**+==。
材料力学试题及答案)汇总
2010—2011材料力学试题及答案A 一、单选题（每小题2分，共10小题，20分） 1、 工程构件要正常安全的工作，必须满足一定的条件。下列除（ ）项，其他各项是必须满足的条件。 A 、强度条件 B 、刚度条件 C 、稳定性条件 D 、硬度条件 2、内力和应力的关系是（ ） A 、内力大于应力 B 、内力等于应力的代数和 C 、内力是矢量，应力是标量 D 、应力是分布内力的集度 3、根据圆轴扭转时的平面假设，可以认为圆轴扭转时横截面（ ）。 A 、形状尺寸不变，直径线仍为直线。 B 、形状尺寸改变，直径线仍为直线。 C 、形状尺寸不变，直径线不保持直线。 D 、形状尺寸改变，直径线不保持直线。 4、建立平面弯曲正应力公式z I My =σ,需要考虑的关系有（ ）。 A 、平衡关系,物理关系，变形几何关系； B 、变形几何关系，物理关系，静力关系； C 、变形几何关系，平衡关系，静力关系； D 、平衡关系, 物理关系，静力关系； 5、利用积分法求梁的变形，不需要用到下面那类条件（ ）来确定积分常数。 A 、平衡条件。 B 、边界条件。 C 、连续性条件。 D 、光滑性条件。 6、图示交变应力的循环特征r 、平均应力m σ、应力幅度a σ分别为（ ）。 A -10、20、10； B 30、10、20； C 31- 、20、10； D 31-、10、20 。
7、一点的应力状态如下图所示，则其主应力1σ、2σ、3σ分别为（）。 A 30MPa、100 MPa、50 MPa B 50 MPa、30MPa、-50MPa C 50 MPa、0、-50Mpa、 D -50 MPa、30MPa、50MPa 8、对于突加载的情形，系统的动荷系数为（）。 A、2 B、3 C、4 D、5 9、压杆临界力的大小，（）。 A 与压杆所承受的轴向压力大小有关； B 与压杆的柔度大小有关； C 与压杆材料无关； D 与压杆的柔度大小无关。 10、利用图乘法计算弹性梁或者刚架的位移，要求结构满足三个条件。以下那个条件不是必须的（） A、EI为常量 B、结构轴线必须为直线。 C、M图必须是直线。 D、M和M至少有一个是直线。 二、按要求作图(共12分) 1、做梁的剪力、弯矩图(10分)
中南大学有限元分析及应用(2008研)考试试卷
中南大学考试试卷(试卷共2页) 2007 -- 2008学年下学期时间110分钟有限元分析及应用课程 36 学时 2 学分 专业年级： 07级研究生总分100分考试形式：大型作业注：此页不作答题纸，请将答案写在答题纸上A4 一．简答题（共40分,每题10分） 1．简述非节点载荷移置的缘由及遵循的原则(p49)？写出集中力移置的普遍公式(p50)。 2．任何一个有限元分析问题都是空间问题，什么情况下可以简化为平面问题(p84)、轴对称问题(p128)？ 3．简述有限元几何模型(关键点、线、面、体)、网格模型(节点、单元)、有限元模型层次关系(p204)及所包含的主要内容，并说明约束的意义(p265)？ 4．写出3节点平面三角形单元(p90)、8节点六面体单元(p168)的广义位移函数，并说明单元位移模式选取的原则(p32)？ 二．综合题（20分） 1.对于如图所示结构模型，若以角速度W绕中心轴线旋转，求其在惯性力的作用下的 最大应力与变形。（1）概述该分析模型的简化方法及理由(p363)？（2）说明约束的施加方法及理由(p369)？(3)用图表示简化后的物理分析模型 三计算与软件操作题（本题共40分）
据提取等关键命令 点应力、形变、位移数束与载荷、求解以及节成、单元生成、施加约要求：简明概述节点生及位移。 号节点）的应力、应变（即软件，求形心，用对于图出刚度矩阵 细节不要描述，直接写注意：单元刚度的求解；）的应变、应力、位移，（手工计算求板形心点，按有限元的解题步骤对于图求解下列问题： 构离散，节点三角形单元进行结。用泊松比的拉力，弹性模量承受。左端固定，右端角点，板厚，宽已知矩形薄板，长5P ANSYS 22)2(50100P 12)1(33.0,102E N 5001010020011--=?====u Pa cm t cm w cm l 图2-1 矩形薄板2单元结构离散方式 图2-2 矩形薄板4单元结构离散方式 X X Y Y P
数值分析
习 题 1. 指出有效数49×102，0.0490，490.00的绝对误差限、相对误差限和有效数字位数． 2. 将 3.142作为π的近似值，它有几位有效数字，相对误差限和绝对误差限各为多少？ 3. 要使101的近似值x * 的相对误差限不超过4102 1?×，问查开方表时x * 需要保留几位有效数字？ 4. 已知近似数x * 有两位有效数字，试估计其相对误差限． 5. 设x * 为x 的近似数, 证明n x * 的相对误差大约为x * 相对误差的n 1倍． 6. 某矩形的长和宽大约为100cm 和50cm, 应该选用最小刻度为多少cm 的测量工具, 才能保证计算出的面积误差（绝对值）不超过0.15cm 2． 7. 已知三角形面积c ab S sin 2 1=，测量a , b , c 时产生的相对误差为)(*a e r ，)(*b e r ，)(*c e r ，其中2 ,0*π<
北航材料力学在线作业三 附答案
北航《材料力学》在线作业三 一、单选题（共 20 道试题，共 80 分。） 1. 对于不同柔度的塑性材料压杆，其最大临界应力将不超过材料的（） A. 比列极限 B. 弹性极限 C. 屈服极限 D. 强度极限 满分：4 分 2. 在平面图形的几何性质中，（）的值可正、可负、也可为零。 A. 静矩和惯性矩 B. 极惯性矩和惯性矩 C. 惯性矩和惯性积 D. 静矩和惯性积 满分：4 分 3. 图示简支梁，截面C的挠度与（）成反比例关系。 A. 跨度L B. 荷载集度q C. q的作用区域a D. 材料的弹性模量E 满分：4 分 4. 某机器的圆轴用45号钢制成，在使用中发现弯曲刚度不够，改善抗弯刚度的有效措施是（） A. 对轴进行调质热处理 B. 改用优质合金钢
C. 加粗轴径 D. 增加表面光洁度 满分：4 分 5. 如图所示，在平板和受啦螺栓之间垫上一个垫圈，可以提高（）强度。 A. 螺栓的拉伸 B. 螺栓的剪切 C. 螺栓的挤压 D. 平板的挤压 满分：4 分 6. 表示扭转变形程度的量（）。 A. 是扭转角，不是单位扭转角 B. 是单位扭转角，不是扭转角 C. 是扭转角和单位扭转角 D. 不是扭转，也不是单位扭转角 满分：4 分 7. 挠曲线近似微分方程不能用于计算（）的位移。 A. 变截面直梁 B. 等截面曲梁 C. 静不定直梁 D. 薄壁截面等直梁 满分：4 分
8. 在水平压缩冲击问题中，曾得到这样一个结论，杆件体积越大，相应的冲击应力越小，该结 论（） A. 只适用于等截面直杆，不适用于变截面直杆 B. 只适用于变截面直杆，不适用于等截面直杆 C. 既适用于等直杆，也适用于变截面直杆 D. 既不适用于等直杆，也不适用于变截面直杆 满分：4 分 9. 图示刚性槽内嵌入一个铝质立方块，设铝块与钢槽间既无间隙，也无摩擦，则在均布压力p 作用下铝块处于（） A. 单向应力状态，单向应变状态 B. 平面应力状态，平面应变状态 C. 单向应力状态，平面应变状态 D. 平面应力状态，单向应变状态 满分：4 分 10. 非对称薄壁截面梁只发生平面弯曲、不发生扭转的横向力作用条件是（）。 A. 作用面平行于形心主惯性平面 B. 作用面重合于形心主惯性平面 C. 作用面过弯曲中心 D. 作用面过弯曲中心且平行于形心主惯性平面 满分：4 分 11. 在下列关于轴向拉压杆轴力的说法中，（）是错误的。 A. 拉压杆的内力只有轴力
有限元考试试题及答案
一、 简答题（共40分，每题10分） 1. 论述单元划分应遵循的原则。 2. 说明形函数应满足的条件。 3. 说明四边形等参数单元中“等参数”的含义，即为什么要引入等参数单元。 4. 阐述边界元法的主要优缺点。 二、 计算题（共60分，每题20分） 1. 一杆件如图3所示，杆件上方固定后，在下方受垂直向下的集中力作用，已 知：杆件材料的杨氏模量2 721/100.3in lbf E E ?==，截面积2125.5in A =， 2275.3in A =，长度in L L 1221==，集中力lbf P 100=，用有限元方法求解B 点 和C 点位移。备注：（1）1 lbf （磅力，libra force ） = N 。（2）杨氏模量、弹性 模量、Young 氏弹性模量具有相同含义（10分） 2. 如图2 所示，有一正方形薄板，沿对角承受压力作用，厚度t=1m ，载荷 F=20KN/m ，设泊松比μ=0，材料的弹性模量为E ，试求它的应力分布。（15分） 学院 专业 学号 姓名 y 图1
图2 3. 图示结点三角形单元的124边作用有均布侧压力q，单元厚度为t，求单元的等效结点荷载。 图3
一、简答题 1. 答： 1）合理安排单元网格的疏密分布 2）为突出重要部位的单元二次划分 3）划分单元的个数 4）单元形状的合理性 5）不同材料界面处及荷载突变点、支承点的单元划分 6）曲线边界的处理，应尽可能减小几何误差 7）充分利用结构及载荷的对称性，以减少计算量 2. 答： 形函数应满足的三个条件： a.必须能反映单元的刚体位移，就是位移模式应反映与本单元形变无关的由 其它单元形变所引起的位移。 b.能反映单元的常量应变，所谓常量应变，就是与坐标位置无关，单元内所 有点都具有相同的应变。当单元尺寸取小时，则单元中各点的应变趋于相 等，也就是单元的形变趋于均匀，因而常量应变就成为应变的主要部分。 c.尽可能反映位移连续性；尽可能反映单元之间位移的连续性，即相邻单元 位移协调。 3. 答： 含义：所谓的等参数单元，就是在确定单元形状的插值函数和确定单元位移场的插值函数中采用了完全相同的形函数。 意义：构造出一些曲边地高精度单元，以便在给定地精度下，用数目较少地单元，解决工程实际地具体问题。 4. 答： 有限单元法是基于变分原理的里兹（Ritz）法的另一种形式，从而使里兹法分析的所有理论基础都适用子有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法．利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，面且事先不要求满足任何边界条件，因此它可以用来处理很复杂的连续介质问题。有限单元法中所利用的主要是伽辽金（Galerkin）法。它可以用于已经知道问题的微分方程和
北京航空航天大学2015春《材料力学》在线作业一满分答案
北京航空航天大学2015春《材料力学》在线作业一满分答案
15春北航《材料力学》在线作业一满分答案 一、单选题（共20道试题，共80分。） 1. 图示平面刚架AB段的内力分量（）为零 M 和Q A. B. M 和N c.N和Q d.N 正确答案：A 2.在轴向拉压杆和受扭圆轴的横截面上分别产生（）。 A.线位移、线位移 B.角位移、角位移
C.线位移、角位移
D.角位移、线位移 正确答案：C 3.梁的挠度是（）。 A.横截面上任一点沿梁轴垂直方向的线位移 B.横截面形心沿梁轴垂直方向的线位移 C.横截面形心沿梁轴方向的线位移 D.横截面形心的位移 正确答案：B 4.中性轴是梁的（）的交线。 A.纵向对称面与横截面 B.纵向对称面与中性层 C.横截面与中性层 D.横截面与顶面或底面 正确答案：C 5. 图示单元体（）无线应变。 A.仅沿X方向
仅沿y方向 B. c.沿x,y两个方向 d.沿任意方向 正确答案：C 5.在横截面面积相等的条件下，（）截面杆的抗扭强度最高。 A.正方形 B.矩形 C.实心圆形 D.空心圆形 正确答案：D 6.在冲击应力和变形实用计算的能量法中，因 为不计被冲击物的质量，所以计算结果与实 际情况相比，（） A.冲击应力偏大，冲击变形偏小 B.冲击应力偏小，冲击变形偏大 C.冲击应力和变形均偏大 D.冲击应力和变形均偏小
正确答案：C 7.长度和受载形式均相同的两根悬臂梁，若其抗弯截面刚度EI相同，而截面形状不同， 则两梁的（） A.最大正应力相等，最大挠度不等 B.最大正应力不等，最大挠度相等 C.最大正应力和最大挠度都不等 D.最大正应力和最大挠度都相等 正确答案：B 8.在三向压应力接近相等的情况下，脆性材料和塑性材料的破坏方式（）。 A.分别为脆性断裂、塑性流动 B.分别为塑性流动、脆性断裂 C.都为脆性断裂 D.都为塑性流动 正确答案：D 10. 下列四根圆轴，横截面面积相同，单位长度扭转角
北航数值分析大作业第二次
《数值分析》计算实习作业 （第二题）
算法设计方案： 1、对矩阵A 赋值，取计算精度ε=1×10-12； 2、对矩阵A 进行拟上三角化，得到A (n-1)，并输出A (n-1)； 对矩阵A 的拟上三角化，通过直接调用子函数inftrianglize(A)来实现；拟上三角化得到的矩阵A (n-1)输出至文件solution.txt 中。 3、对A (n-1)进行QR 分解并输出Q 、R 及RQ 矩阵； QR 分解通过直接调用子函数QRdescom(A,Q,R, n)实现。 4、运用QR 方法求所有的特征值，并输出； (1)初始时令m=n ，在m>2的条件下执行； (2)判断如果|A mm-1|<ε，则得到一个特征值，m=m-1，转(4)；否则转(3)； (3)判断如果|A m-1m-2|<ε，则得到两个特征值，m=m-2，转(4)； (4)判断如果m ≤2，转(6)；否则转(5)； (5)执行相似迭代，转(2)； k k T k k k k k k k k k k Q A Q A R Q M I D A D tr A M ==+-=+1)2)det(( (6)求出最后的一个或两个特征值； (7)输出全部的特征值至文件solution.txt 中。 5、输出QR 分解法迭代结束之后的A (n-1)至文件solution.txt 中； 6、通过反幂法求出所有实特征值的特征向量并输出。 首先令B=(A-λi I)，其中λi 是实特征值；反幂法通过调用子函数Bpowmethod(B,x1)实现，最终λi 对应的特征向量就是x1；最后将所有的实特征值的特征向量输出。