有限元分析试题(同济)
- 格式:doc
- 大小:307.40 KB
- 文档页数:2
有限元考试试题有限元考试试题在工程学领域中,有限元分析是一种常用的数值计算方法,用于解决结构力学、热传导、流体力学等问题。
有限元方法的应用广泛,因此在相关领域中的考试中,有限元试题是非常重要的一部分。
本文将探讨一些有限元考试试题,以帮助读者更好地理解和应用这一方法。
1. 问题描述:一根长度为L的杆件,两端固定,如何确定杆件上各个位置的位移?解答:这是一个典型的弹性力学问题,可以通过有限元方法进行求解。
首先,将杆件分割成若干个小单元,每个小单元内部的位移近似为线性。
然后,根据杆件的边界条件,建立相应的刚度矩阵和载荷向量。
最后,通过求解线性方程组,得到杆件上各个位置的位移。
2. 问题描述:如何确定一个结构的应力分布情况?解答:有限元分析可以用来计算结构的应力分布情况。
首先,将结构分割成若干个小单元,每个小单元内部的应力近似为线性。
然后,根据结构的边界条件和加载情况,建立相应的刚度矩阵和载荷向量。
最后,通过求解线性方程组,得到结构上各个位置的应力分布情况。
3. 问题描述:如何确定一个结构的固有频率?解答:固有频率是指结构在没有外界激励下自由振动的频率。
有限元分析可以用来计算结构的固有频率。
首先,将结构分割成若干个小单元,每个小单元内部的位移近似为线性。
然后,根据结构的边界条件,建立相应的刚度矩阵和质量矩阵。
最后,通过求解特征值问题,得到结构的固有频率和相应的振型。
4. 问题描述:如何考虑非线性材料的影响?解答:有限元分析可以考虑非线性材料的影响。
在材料的应力-应变关系中,通常存在非线性现象,如材料的屈服、硬化、蠕变等。
为了考虑这些非线性现象,可以采用增量形式的有限元分析方法。
在每个增量步骤中,根据当前应力状态和材料的非线性特性,更新刚度矩阵和载荷向量。
通过迭代求解,可以得到结构的非线性响应。
5. 问题描述:如何考虑流体结构耦合问题?解答:有限元分析可以考虑流体结构耦合问题。
在流体结构耦合问题中,结构的变形会影响流体的流动,而流体的流动又会对结构施加载荷。
同济大学有限元试卷-CAL-FENGHAI.-(YICAI)-Company One1
1
2
3
4
5
07有限元试卷
(每题10分)
1、为什么说变形体虚位移原理与势能原理彼此等价?
2、
3、试述有限单元法单元分析步骤和核心问题?
4、
5、试求左端铰接单元的单元刚度矩阵。
6、试述有限单元法整体分析步骤。
7、试求四边形12节点单元形函数N3,N8。
(图画的不好,各位参见书上图片)
8、试述协调元收敛准则。
证明矩形双线性单元为协调元。
9、什么是等参元?试述母元和子元的关系。
10、
11、常应变三角元和三角形环单元的区别为什么
12、
13、试述非完全协调元收敛准则。
Wilson非协调元一定收敛么?
14、试述影响有限元分析精度的因素及解决方法。
2。
有限元期末考试题及答案一、选择题1. 有限元方法是一种数值分析方法,主要用于求解什么类型的数学问题?A. 线性代数方程B. 微分方程C. 积分方程D. 代数方程答案:B2. 在有限元分析中,单元的划分是基于什么原则?A. 单元数量B. 单元形状C. 问题域的几何特性D. 计算资源答案:C3. 下列哪项不是有限元分析中常用的单元类型?A. 三角形单元B. 四边形单元C. 六面体单元D. 圆形单元答案:D二、填空题4. 有限元方法中,______是指将连续的物理域离散成有限数量的小区域,这些小区域称为单元。
答案:离散化5. 在进行有限元分析时,通常需要定义材料属性,包括______、密度和弹性模量等。
答案:泊松比三、简答题6. 简述有限元方法的基本步骤。
答案:有限元方法的基本步骤包括:定义问题域、离散化问题域、选择单元类型、定义材料属性、构建全局刚度矩阵、施加边界条件、求解线性代数方程、提取结果。
7. 解释什么是有限元分析中的收敛性,并说明影响收敛性的因素。
答案:收敛性是指随着单元数量的增加,有限元分析结果逐渐接近真实解的性质。
影响收敛性的因素包括单元的类型、形状、大小以及网格的布局等。
四、计算题8. 假设有一个长度为2米的杆,两端固定,中间施加了一个向下的力F=1000N。
如果杆的材料是钢,其弹性模量E=210 GPa,泊松比ν=0.3,请计算杆的弯曲位移。
答案:首先,根据Euler-Bernoulli梁理论,可以写出弯曲位移的方程为:\[ w(x) = \frac{F}{384EI} L^3 \]其中,\( w(x) \) 是位移,\( F \) 是施加的力,\( L \) 是杆的长度,\( E \) 是弹性模量,\( I \) 是截面惯性矩。
对于一个矩形截面,\( I \) 可以表示为:\[ I = \frac{bh^3}{12} \]假设杆的截面宽度为b,高度为h,代入上述公式,可以计算出位移。
有限元分析及应用作业报告试题10一、问题描述确定图示扳手中的应力, E=210Gpa,μ=0.3, 假设厚度为10mm;并讨论采用何种处理可降低最大应力或改善应力分布。
图1为扳手的基本形状和基本尺寸图二、数学建模与分析由图1及问题描述可知,板手的长宽尺寸远远大于厚度,研究结构为一很薄的等厚度薄板,满足平面应力的几何条件;作用于薄板上的载荷平行于板平面且作用在沿厚度方向均匀分布在办手柄的左边缘线,而在两板面上无外力作用,满足平面应力的载荷条件。
故该问题属于平面应力问题,薄板所受的载荷为面载荷,分布情况及方向如图1所示,建立几何模型,并进行求解。
薄板的材料为钢,则其材料参数:弹性模量E=2.1e11,泊松比σ=0.3三、有限元建模1、单元选择:选取三节点常应变单元来计算分析薄板扳手的位移和应力。
由于此问题为平面应力问题,:三节点常应变单元选择的类型是PLANE42(Quad 4node42),该单元属于是四节点单元类型,在网格划分时可以对节点数目控制使其蜕化为三节点单元。
2、定义材料参数:ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK3、生成几何模型:a.创建关键点点:ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入16个点的坐标→OKb、将这16个关键点有直线依次连起来,成为线性模型4、生成实体模型:ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →连接特征点→生成两个area→Operate→Subtract→拾取整个扳手区域→OK→生成扳手模型5、结点布置及规模6、网格划分方案ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →Mesh: Areas, Shape: Tri,Free →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window)7、载荷及边界条件处理8、求解控制A、模型施加约束给模型施加x方向约束ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On Lines →拾取模型左部的竖直边→OKB、给模型施加载荷ANSYS Main Menu: Solution →Define Loads →Apply →Structural →force→on keypoints→拾取上面左端关键点→700N/mm→okC、分析计算:ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load Step window) →OK6)结果显示:ANSYS Main Menu: General Postproc →Plot Results →Deformed Shape… →select Def + Undeformed →OK (back to Plot Results window) →Contour Plot →Nodal Solu →select: DOF solution →displacement vector sum,von mises stress→OK四、计算结果及结果分析1、三节点常应变单元1)三节点单元的网格划分图2 常应变三节点单元的网格划分平面图图3 常应变三节点单元的网格划分立体图2)三节点单元的约束受载情况图4 常应变三节点单元的约束受载图3)三节点单元的位移分析图5 常应变三节点单元的位移分布图4)三节点单元的应力分析图6 常应变三节点单元的应力分布图2、六节点三角形单元1)六节点三角形单元网格划分图7 六节点三角形单元网格划分图2)六节点三角形单元约束和受载情况分析图8 六节点三角形单元约束受载图3)六节点三角形单元位移分析图9 六节点三角形单元的变形分布图4) 六节点三角形单元的应力分析图9 六节点三角形单元的应力分布图图10 六节点三角形单元的局部应力分布图根据以上位移和应力图,可以得出常应变三节点单元和六节点三角形单元的最小最大位移应力如表1-1所示。
一、如图所示的1D 杆结构,试用取微单元体的方法建立起全部基本方程和边界条件,并求出它的所有解答。
注意它的弹性模量为E 、横截面积A解:如图1.1所示的1D 杆结构,其基本变量为 位移 x u 应变 x ε 应力 x σ取微单元体Adx ,其应力状态如图1.2,由泰勒展开式知()⋅⋅⋅⋅⋅+∂∂+⋅∂∂+=+22221dx x dx x dx x x x x σσσσ略去2阶以上的商阶微量知()dx xdx x xx ⋅∂∂+=+σσσ 由力的平衡知0=∑i x :0=-⎪⎭⎫ ⎝⎛∂∂+A A dx x x x x σσσ即力的平衡方程为:⋅⋅⋅⋅=0dxd xσ① 位移由图1.3知(泰勒展开,略去商阶微量)()dx xu u dx x u xx ⋅∂∂+=+ dxu dxdxdx u dx x uu ABABB A xx x x x ∂=-+-∂∂+=-=∴)(''ε应变 即几何方程为:⋅⋅⋅⋅=dxdu xx ε② 根据虎克定律知⋅⋅⋅⋅⋅⋅⋅=⋅=dxdu E E xx x εσ③ 由①、②、③知该1D 杆的基本方程为⎪⎪⎪⎩⎪⎪⎪⎨⎧====dx du E E dx du dx d x x xx xxεσεσ0 在节点1时位移:00==x x u 在节点2时应力:APlx x==σ即其边界条件为00==x x u on u SAPlx x==σ on P S 由①式知⋅⋅⋅⋅⋅=0c x σ ④ ④代入③解得:dxdu Ec x=0 ⋅⋅⋅⋅⋅⋅⋅+=10c x Ec u x ⑤ 0c 、1c 为待定系数结合边界条件知⎪⎪⎩⎪⎪⎨⎧==+A P c c x Ec 010解知得APc =0,01=c ∴⎪⎪⎪⎩⎪⎪⎪⎨⎧==⋅==EA P E x EA P u A P x xx x σεσ二、设平面问题中的应力问题y a x a a x 321++=σy a x a a y 654++=σ y a x a a xy 987++=τ其中i a (1、2、………9)为常数,令所有体积力为零,对下面特殊情况说明平衡是否满足?为什么?或者i a 之间有什么关系才满足平衡。
一、填空(共10个空,每空2分,共20分)11、有限元法是近似求解连续场问题的数值方法。
2、有限元法将连续的求解域离散,得到有限个单元,单元和单元之间用节点相连。
3、直梁在外力作用下,横截面上的内力有剪力和弯矩两个。
4、平面刚架结构在外力作用下,横截面上的内力有剪力、弯矩和轴力。
5、进行直梁的有限元分析,梁单元上每个节点的节点位移为挠度和转角。
、平面刚架结构中,已知单元e的坐标变换矩阵[T e]及局部坐标系x´O´y ´下的单元刚度矩阵[K´]e,则单元在整体坐标系xOy下的单元刚度矩阵为 P31 。
7、平面刚架结构中,已知单元e的坐标变换矩阵[T e]及整体坐标系xOy下的单元节点力矩阵{p}e,则单元在局部坐标系x´O´y´下的单元节点力矩阵为 P30 。
8、在弹性范围和小变形的前提下,节点力和节点位移之间是线性系。
9、弹性力学问题的方程个数有 15个,未知量个数有 15 个。
10、弹性力学平面问题的方程个数有个,未知量个数有个。
11、把经过物体内任意一点各个截面的应力状况叫做一点的应力状态。
12、形函数在单元节点上的值,具有本点为 1 、它点为零的性质,并且在三角形单元的任一节点上,三个形函数之和为 1 。
13、形函数是定义于元内部坐标连续函数。
14、在进行节点编号时,要尽量使同一单元的相邻节点的号码差尽可能小,以便最大限度地缩小刚度矩阵带宽,节省存储、提高计算效率。
15、三角形单元的位移模式为。
16、矩形单元的位移模式为。
17、在选择多项式位移模式的阶次时,要求所选的位移模式应该与局部坐标系的方位无关,这一性质称为几何各向同性。
18、单元刚度矩阵描述了节点力和节点位移之间的关系。
19、在选择多项式作为单元的位移模式时,多项式阶次的确定,要考虑解答的收敛性,即要满足单元的完备性和协调性的要求。
20、三节点三角形单元内的应力和应变是常数,四节点矩形单元内的应力和应变是线性变化的。
同济大学本科课程期终考试统一命题纸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分)
1、简述有限单元法结构刚度矩阵的特点。
2、简述有限元法中选取单元位移函数(多项式)的一般原则。
3、简述有限单元法的收敛性准则。
4、考虑下列三种改善应力结果的方法(1)总体应力磨平、(2)单元应力磨平和(3)分片应力磨平,请分别将它们按计算精度(高>低)和计算速度(快>慢)进行排序。
四.计算题(共40分,每题20分)
1、如图1所示等腰直角三角形单元,其厚度为t ,弹性模量为E ,泊松比0ν=;单元的边
长及结点编号见图中所示。
求 (1) 形函数矩阵N
(2) 应变矩阵B 和应力矩阵S
(3) 单元刚度矩阵e K
2、图2(a )所示为正方形薄板,其板厚度为t ,
四边受到均匀荷载的作用,荷载集度为21/N m ,同时在y 方向相应的两顶点处分别承受大小为2/N m 且沿板厚度方向均匀分布的荷载作用。
设薄板材料的弹性模量为E ,泊松比0ν=。
试求
(1) 利用对称性,取图(b )所示1/4结构作为研究对象,并将其划分为4个面积大小
相等、形状相同的直角三角形单元。
给出可供有限元分析的计算模型(即根据对称性条件,在图(b )中添加适当的约束和荷载,并进行单元编号和结点编号)。
(2) 设单元结点的局部编号分别为i 、j 、m ,为使每个单元刚度矩阵e K 相同,试在
图(b )中正确标出每个单元的合理局部编号;并求单元刚度矩阵e K 。
(3) 计算等效结点荷载。
(4) 应用适当的位移约束之后,给出可供求解的整体平衡方程(不需要求解)。
123
图1
图
2 (a ) (b )。