北航《数值分析》总复习题(2011年)

北航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 、α作为可设计量,保持结构其余参数不变,求结构重量最轻时的α值。

北航材料力学试题往年试卷-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

北航数值分析大作业第一题幂法与反幂法

《数值分析》计算实习题目 第一题： 1. 算法设计方案 （1）1λ，501λ和s λ的值。 1)首先通过幂法求出按模最大的特征值λt1，然后根据λt1进行原点平移求出另一特征值λt2，比较两值大小，数值小的为所求最小特征值λ1，数值大的为是所求最大特征值λ501。 2)使用反幂法求λs ，其中需要解线性方程组。因为A 为带状线性方程组，此处采用LU 分解法解带状方程组。 （2）与140k λλμλ-5011=+k 最接近的特征值λik 。 通过带有原点平移的反幂法求出与数k μ最接近的特征值 λik 。 （3）2cond(A)和det A 。 1）1=n λλ2cond(A)，其中1λ和n λ分别是按模最大和最小特征值。 2）利用步骤（1）中分解矩阵A 得出的LU 矩阵，L 为单位下三角阵,U 为上三角阵，其中U 矩阵的主对角线元素之积即为det A 。 由于A 的元素零元素较多，为节省储存量，将A 的元素存为6×501的数组中，程序中采用get_an_element()函数来从小数组中取出A 中的元素。 2.全部源程序 #include #include void init_a();//初始化A double get_an_element(int,int);//取A 中的元素函数 double powermethod(double);//原点平移的幂法 double inversepowermethod(double);//原点平移的反幂法 int presolve(double);//三角LU 分解 int solve(double [],double []);//解方程组 int max(int,int); int min(int,int); double (*u)[502]=new double[502][502];//上三角U 数组 double (*l)[502]=new double[502][502];//单位下三角L 数组 double a[6][502];//矩阵A int main() { int i,k; double lambdat1,lambdat2,lambda1,lambda501,lambdas,mu[40],det;

北航数值分析大作业第二题

数值分析第二次大作业 史立峰 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。

北航数值分析报告第三次大作业

数值分析第三次大作业 一、算法的设计方案： （一）、总体方案设计： 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 =-，

北航数值分析报告大作业第八题

北京航空航天大学 数值分析大作业八 学院名称自动化 专业方向控制工程 学号 学生姓名许阳 教师孙玉泉 日期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分)

数值分析

习 题 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*π<>2时的情形．用所设计的算法以及二次方程求根公式计算05.240=p ，00.1=q 时方程根的近似值(计算过程保留2位小数)，并给出它们的相对误差限（根的准确值为L 0916683.4801?=x , L 002082935.02?=x ）．

北航材料力学在线作业三 附答案

北航《材料力学》在线作业三 一、单选题（共 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. 拉压杆的内力只有轴力

北航数值分析课程第一次大作业讲解

《数值分析A》计算实习题目第一题 一．算法设计方案： 1.矩阵A的存储与检索 将带状线性矩阵A[501][501]转存为一个矩阵MatrixC[5][501] . 由于C语言中数组角标都是从0开始的，所以在数组MatrixC[5][501]中检索A的带内元素a ij的方法是： A的带内元素a ij=C中的元素c i-j+2,j 2.求解λ1，λ501，λs ①首先分别使用幂法和反幂法迭代求出矩阵按摸最大和最小的特征值λmax和λmin。λmin即为λs； 如果λmax>0,则λ501=λmax；如果λmax<0,则λ1=λmax。 ②使用带原点平移的幂法（mifa（）函数），令平移量p=λmax，求出对应的按摸最大的特征值λ，max， 如果λmax>0,则λ1=λ，max+p；如果λmax<0,则λ501=λ，max+p。 3.求解A的与数μk=λ1+k（λ501-λ1）/40的最接近的特征值λik （k=1,2，…，39）。 使用带原点平移的反幂法，令平移量p=μk，即可求出与μk最接近的特征值λik。 4.求解A的（谱范数）条件数cond（A）2和行列式d etA。 ①cond（A）2=|λ1/λn|，其中λ1和λn分别是矩阵A的模最大和最小特征值。 ②矩阵A的行列式可先对矩阵A进行LU分解后，detA等于U所有

对角线上元素的乘积。 二．源程序(VS2010环境下，C++语言) #include #include #include #include #include #include #include #define E 1.0e-12 /*定义全局变量相对误差限*/ int max2(int a,int b) /*求两个整型数最大值的子程序*/ { if(a>b) return a; else return b; } int min2(int a,int b) /*求两个整型数最小值的子程序*/ { if(a>b) return b; else return a; } int max3(int a,int b,int c) /*求三整型数最大值的子程序*/ { int t; if(a>b) t=a; else t=b; if(t

北京航空航天大学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；最后将所有的实特征值的特征向量输出。

北航材料力学各章重点(选择题)

材料力学各章重点 一、绪论 1．各向同性假设认为，材料沿各个方向具有相同的 A 。 (A)力学性质； (B)外力； (C)变形； (D)位移。 2．均匀性假设认为，材料内部各点的 C 是相同的。 (A)应力； (B)应变； (C)位移； (C)力学性质。 3．构件在外力作用下 B 的能力称为稳定性。 （A）不发生断裂；（B）保持原有平衡状态； （C）不产生变形；（D）保持静止。 4．杆件的刚度是指 D 。 （A）杆件的软硬程度；（B）件的承载能力； （C）杆件对弯曲变形的抵抗能力；（D）杆件对弹性变形的抵抗能力。 二、拉压 1．低碳钢材料在拉伸实验过程中，不发生明显的塑性变形时，承受的最大应力应当小于 D 的数值， （A）比例极限；（B）许用应力；（C）强度极限；（D）屈服极限。 2．对于低碳钢，当单向拉伸应力不大于 C 时，虎克定律ζ=Eε成立。 (A) 屈服极限ζs；(B)弹性极限ζe；(C)比例极限ζp；(D)强度极限ζb。 3．没有明显屈服平台的塑性材料，其破坏应力取材料的 B 。 （A）比例极限ζp；（B）名义屈服极限ζ0.2；

（C ）强度极限ζb ；（D ）根据需要确定。 4．低碳钢的应力~应变曲线如图所示，其上 C 点的纵坐标值为该钢的强度极限σb 。 (A)e ； (B)f ； (C)g ； (D)h 。 5 、三种材料的应力—应变曲线分别如图所示。其中强度最高、刚度最大、塑性最好的材料分别是 A 。 (A)a 、b 、c ； (B)b 、c 、a ； (C)b 、a 、c ； (D)c 、b 、a 。 5．材料的塑性指标有 C 。 (A)ζs 和δ； (B)ζs 和ψ； (C)δ和ψ； (D)ζs ,δ和ψ。 6．确定安全系数时不应考虑 D 。 (A)材料的素质；(B)工作应力的计算精度；(C)构件的工作条件；(D)载荷的大小。 8．系统的温度升高时，下列结构中的____A______不会产生温度应力。 9、图示两端固定阶梯形钢杆，当温度升高时 D 。 (A)AC 段应力较大，C 截面向左移； A B C D 3题图

北航数值分析作业第一题题解

北航数值分析作业第一题： 一、算法设计方案 1.要求计算矩阵的最大最小特征值，通过幂法求得模最大的特征值，进行一定 判断即得所求结果； 2.求解与给定数值接近的特征值，可以该数做漂移量，新数组特征值倒数的绝 对值满足反幂法的要求，故通过反幂法即可求得； 3.反幂法计算时需要方程求解中间过渡向量，需设计Doolite分解求解； 4.|A|=|B||C|，故要求解矩阵的秩，只需将Doolite分解后的U矩阵的对角线相 乘即为矩阵的Det。 算法编译环境：vlsual c++6.0 需要编译函数：幂法，反幂法，Doolite分解及方程的求解 二、源程序如下： #include #include #include #include int Max(int value1,int value2); int Min(int value1,int value2); void Transform(double A[5][501]); double mifa(double A[5][501]); void daizhuangdoolite(double A[5][501],double x[501],double b[501]); double fanmifa(double A[5][501]); double Det(double A[5][501]); /***定义2个判断大小的函数，便于以后调用***/ int Max(int value1,int value2) { return((value1>value2)?value1:value2); } int Min(int value1,int value2) { return ((value1

北航2010-2015年研究生数值分析期末模拟试卷与真题

数值分析模拟卷A 一、填空（共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 型的？ 五、（15分） 设有常微分方程的初值问题???=='00 )() ,(y x y y x f y ，试用Taylor 展开原理构造形如 )()(11011--++++=n n n n n f f h y y y ββα的方法，使其具有二阶精度，并推导其局部截断误 差主项.

2015北航《材料力学》在线作业二

北航《材料力学》在线作业二 试卷总分：100 测试时间：-- 试卷得分：100 1. T 形截面铸铁梁，设各个截面的弯矩均为正值，则将其截面按图（ ）所示的方式布置，梁的强度最高。 A. A B. B C. C D. D 满分：4 分 得分：4 2. 三根杆的横截面面积及长度均相等，其材料的应力—应变曲线分别如图所示，其中强度最高、刚度最大、塑性最好的杆分别是（ ）。 A. a,b,c B. b,c,a C. b,a,c D. c,b,a 满分：4 分 得分：4 3. 一拉伸钢杆，弹性模量E=200GPa ，比例极限σp=200MPa ，今测得其轴向应变ε=0.0015，则横截面上的正应力（）。 A. σ=Eε=300MPa B. σ>300MPa C. 200MPa<σ<300MPa D. σ<200MPa 单选题 判断题 一、单选题（共 20 道试题，共 80 分。） 得分：80V n m l k j i n m l k j n m l k j n m l k j n m l k j i n m l k j n m l k j n m l k j n m l k j n m l k j n m l k j i n m l k j

满分：4 分 得分：4 4. 如图所示，在平板和受啦螺栓之间垫上一个垫圈，可以提高（ ）强度。 A. 螺栓的拉伸 B. 螺栓的剪切 C. 螺栓的挤压 D. 平板的挤压 满分：4 分 得分：4 5. 某机轴的材料为45号钢，工作时发生弯曲和扭转组合变形。对其进行强度计算时，宜采用（ ）强度理论。 A. 第一或第二 B. 第二或第三 C. 第三或第四 D. 第一或第四 满分：4 分 得分：4 6. 在下列关于轴向拉压杆轴力的说法中，（）是错误的。 A. 拉压杆的内力只有轴力 B. 轴力的作用线与杆轴重合 C. 轴力是沿杆轴作用的外力 D. 轴力和杆的横截面和材料无关 满分：4 分 得分：4 7. 若某低碳钢构件危险点的应力状态为近乎三向等值拉伸，进行强度校核时宜采用（）强度理论。 A. 第一 B. 第二 C. 第三 D. 第四 满分：4 分 得分：4 8. 图示简支梁，截面C 的挠度与（ ）成反比例关系。 A. 跨度L n m l k j n m l k j n m l k j n m l k j i n m l k j n m l k j n m l k j i n m l k j n m l k j n m l k j n m l k j i n m l k j n m l k j i n m l k j n m l k j n m l k j n m l k j

北航数值分析第三次大作业

数值分析第三次大作业 一、算法的设计方案 1、求解非线性方程组 将题目中给出的(,)i i x y 当作已知量代入题目给定的非线性方程组，求出与 (,)i i x y 相对应的数组te[i][j]，ue[i][j]，此处采用的是牛顿法解非线性方程组，其 算法如书上91页所示。 2、分片二次代数插值 对所求出的数组te[i][j]，ue[i][j]，通过分片二次代数插值运算，得到与数组te[11][21]，ue[11][21]对应的数组ze[11][21]，从而得到二元函数z=(,)i i f x y ，此处采用如书上101页例2中所示的分片二次代数插值。 3、曲面插值 利用x[11]，y[21]，ze[11][21]建立二维函数表，进行曲面插值计算，逐步提高k 值，计算其精度，看其是否满足要求，以此来确定循环结束的时刻，并得到曲面拟合的系数矩阵C[r][s]，此处的算法如书142页所示，只需将所需矩阵给出，然后按公式进行计算即可。 4、比较 观察和),(j i y x p 逼近(,)i i f x y 的效果。观察逼近效果只需要利用新给的点列 (,)i i 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 。 5、几点说明 分片二次插值的结果x[i]，y[j]，ze[i][j]输出到一个文件shubiao.txt 中，方便结果的复制与粘贴。 曲面插值的结果输出到一个文件xishu.txt 中，包括循环中每一次的k 值以及误差平方和sigma 的值，还有最后满足误差要求时曲面插值的系数C[r][s]。 观察逼近效果的结果输出到一个文件shubiao1.txt 中，方便结果的复制与粘贴。

07-08年度北航材料力学试题第一学期期末试卷

2007~2008第1学期期末考试A 卷 一、选择题（每小题4分，共20分） 1. 为计算图示拉压杆m-m 截面的内力， 建立平衡方程。 A 只能对m-m 截面左段； B 只能对m-m 截面右段； C 既能对m-m 截面左段也能对m-m 截面右段； D 需取整体为研究对象。 2．低碳钢拉伸、压缩屈服极限分别为+ s σ，- s σ，灰口铸铁拉伸、压缩强度极限分别为+ b σ，- b σ，则有 。 A -+-+==b b s s σσσσ,； B - +-+=

4．设图示简支梁的弯矩为M （x ），则挠曲轴方程可以写为??++=D Cx dxdx x M EIw )(，其 中EI 为弯曲刚度，w 为挠度，对于积分常数有 。 A 0,0==D C ； B 0,0≠=D C ； C 0,0=≠ D C ； D 0,0≠≠D C 。 5．中性轴是梁的 的交线。 A 纵向对称面与横截面； B 横截面与中性层； C 纵向对称面与中性层； D 横截面与顶面。 二、作图题（5分） 试绘制图示梁的挠曲轴的大致形状，并标明凹凸性和拐点。 q 29 qa M 图

材料力学试题及答案汇总

1、图示刚性梁AB 由杆1和杆2支承，已知两杆的材料相同，长度不等，横截面积分别为A 1和A 2，若载荷P 使刚梁平行下移，则其横截面面积（）。 A 、A 1〈A 2B 、A 1 〉A 2 C 、A 1=A 2D 、A 1、A 2为任意 2、建立圆周的扭转应力公式τρ =M ρρ/I ρ时需考虑下列因素中的哪几个？答：（ ） （1） 扭矩M T 与剪应力τ ρ 的关系M T =∫A τρ ρdA （2） 变形的几何关系（即变形协调条件） （3） 剪切虎克定律 （4） 极惯性矩的关系式I T =∫A ρ2dA A 、（1） B 、（1）（2） C 、（1）（2）（3） D 、全部 3、二向应力状态如图所示，其最大主应力σ1=（） A 、σ B 、2σ C 、3σ D 、4σ 4、高度等于宽度两倍(h=2b)的矩形截面梁，承受垂直方向的载荷，若仅将竖放截面改为平放截面，其它条件都不变，则梁的强度 A 、提高到原来的2倍 B 、提高到原来的4倍 C 、降低到原来的1/2倍 D 、降低到原来的1/4倍 5. 已知图示二梁的抗弯截面刚度EI 相同，若二者自由端的挠度相等，则P 1/P 2=（） A 、2 B 、4C 、8 D 、16 6、下列结论中正确的是 （ ） A 、材料力学主要研究各种材料的力学问题 B 、材料力学主要研究各种材料的力学性质 C 、材料力学主要研究杆件受力后变形与破坏的规律 D 、材料力学主要研究各种材料中力与材料的关系 7、有两根圆轴，一根为实心轴，直径为D 1，另一根为空心轴，内外径比为d 2/D 2=0.8。若两轴的长度、材料、轴内扭矩和产生的扭转角均相同，则它们的重量之比W 2/W 1为（ ） A 、0.74 B 、0.62 C 、0.55 D 、0.47 8、材料的失效模式 B 。 A 只与材料本身有关，而与应力状态无关； B 与材料本身、应力状态均有关； C 只与应力状态有关，而与材料本身无关； 题一、3图 题一、5图 题一、 4 题一、1