有限单元法部分课后题答案
- 格式:docx
- 大小:24.84 KB
- 文档页数:4
-----好资料学习有限单元法中“离散”的含义是什么?有限单元法是如何将具有无限自由度的连续介1.1质问题转变成有限自由度问题的?位移有限元法的标准化程式是怎样的?)离散的含义即将结构离散化,即用假想的线或面将连续体分割成数目有限的单元,并1(数的节在其上设定有限个节点;用这些单元组成的单元集合体代替原来的连续体,而场函点值将成为问题的基本未知量。
)给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律,即2(无限自通过插值以单元节点位移表示单元内任意点的位移。
因节点位移个数是有限的,故由度问题被转变成了有限自由度问题。
)有限元法的标准化程式:结构或区域离散,单元分析,整体分析,数值求解。
(3 ?单元刚度矩阵和整体刚度矩阵各有哪些性质?各自的物理意义是什么?两者有何区别1.3整体刚度矩阵的性单元刚度矩阵的性质:对称性、奇异性(单元刚度矩阵的行列式为零)。
个自j Kij 即单元节点位移向量中第稀疏性。
单元 Kij 物理意义质:对称性、奇异性、整体刚度 j 个自由度方向引起的节点力。
由度发生单位位移而其他位移分量为零时,在第中每一列元素的物理意义是:要迫使结构的某节点位移自由度发生单位位移,而其 K 矩阵他节点位移都保持为零的变形状态,在所有个节点上需要施加的节点荷载。
什么叫应变能?什么叫外力势能?试叙述势能变分原理和最小势能原理,并回答下述2.2问题:势能变分原理代表什么控制方程和边界条件?其中附加了哪些条件?,外力所做的功将以变形能的形式储存εσ和应变(1)在外力作用下,物体内部将产生应力起来,这种能量称为应变能。
(2)外力势能就是外力功的负值。
势能变分原理可叙述如下:在所有满足边界条件的协调位移中,那些满足静力平衡条件(3) 的位移使物体势能泛函取驻值,即势能的变分为零V=0 +δp=δ Uεδ∏此即变分方程。
对于线性弹性体,势能取最小值,即02V≥ε+δδ2∏P=δ2U 此时的势能变分原理就是著名的最小势能原理。
2 弹性力学问题的有限单元法思考题2.1 有限元法离散结构时为什么要在应力变化复杂的地方采用较密网格,而在其他地方采用较稀疏网格?答:在应力变化复杂的地方每一结点与相邻结点的应力都变化较大,若网格划分较稀疏,则在应力突变处没有设置结点,而使得所求解的误差很大,若网格划分较密时,则应力变化复杂的地方可以设置更多的结点,从而使得所求解的精度更高一些。
2.2 因为应力边界条件就是边界上的平衡方程,所以引用虚功原理必然满足应力边界条件,对吗?答:对。
2.3 为什么有限元只能求解位移边值问题和混合边值问题?弹性力学中受内压和外压作用的圆环能用有限元方法求解吗?为什么?答:有限元法是一种位移解法,故只能求解位移边值问题和混合边值问题。
而应力边值问题没有确定的位移约束,不能用位移法求解,所以也不能用有限元法求解。
2.4 矩形单元旋转一个角度后还能够保持在单元边界上的位移协调吗?答:能。
矩形单元的插值函数满足单元内部和单元边界上的连续性要求,是一个协调元。
矩形的插值函数只与坐标差有关,旋转一个角度后各个结点的坐标差保持不变,所以插值函数保持不变。
因此矩形单元旋转一个角度后还能够保持在单元边界上的位移协调。
2.5 总体刚度矩阵呈带状分布,与哪些因素有关?如何计算半带宽? 答:因素:总体刚度矩阵呈带状分布与单元内最大结点号与最小结点号的差有关。
计算:设半带宽为B ,每个结点的自由度为n ,各单元中结点整体码的最大差值为D ,则B=n(D+1),在平面问题中n=2。
2.6 为什么单元尺寸不要相差太大,如果这样,会导致什么结果? 答:由于实际工程是一个二维或三维的连续体,将其分为具有简单而规则的几何单元,这样便于网格计算,还可以通过增加结点数提高单元精度。
在几何形状上等于或近似与原来形状,减小由于形状差异过大带来的误差。
若形状相差过大,使结构应力分析困难加大,误差同时也加大。
2.7 剖分网格时,在边界出现突变和有集中力作用的地方要设置结点或单元边界,试说明理由。
有限单元法智慧树知到课后章节答案2023年下山东科技大学山东科技大学绪论单元测试1.有限元法的核心思想是“数值近似”和“离散化”。
( )A:错 B:对答案:对第一章测试1.下列属于平面应力问题的是()。
A:挡土墙 B:受内水压力作用的圆管 C:平板坝的平板支墩 D:重力水坝答案:平板坝的平板支墩2.平衡方程研究的是()之间关系的方程式。
A:应力和位移 B:应力和应变 C:应变和位移 D:应力和体力答案:应力和体力3.弹性力学的边界条件有()。
A:应力边界条件 B:位移边界条件 C:混合边界条件 D:应变边界条件答案:应力边界条件;位移边界条件;混合边界条件4.弹性力学的基本假定有()。
A:假设物体是连续的 B:假设物体的变形是很小的 C:假设物体是完全弹性的 D:假设物体内无初应力 E:假设物体是均匀的和各向同性的答案:假设物体是连续的;假设物体的变形是很小的;假设物体是完全弹性的;假设物体内无初应力;假设物体是均匀的和各向同性的5.在体力为常量时,平衡方程、相容方程及应力边界条件中均不含弹性常数E和μ,故我们可以由一种材料替代另一种材料,用平面应力问题替代平面应变问题作实验,得到的应力是完全一样的。
()A:对 B:错答案:对第二章测试1.一维变带宽存储的存储量()。
A:与结点编号有关 B:与结点编号和单元编号有关 C:与单元编号有关 D:与存储上三角或者下三角有关答案:与结点编号有关2.应变矩阵与()。
A:材料参数有关 B:单元几何尺寸和材料参数都有关 C:单元几何尺寸和材料参数都无关 D:单元几何尺寸有关答案:单元几何尺寸有关3.单元刚度矩阵建立了单元的与之间的关系。
()A:应力,结点位移 B:应力,应变 C:结点力,结点位移 D:应变,结点位移答案:结点力,结点位移4.为了保证有限元解的收敛性,位移函数要满足()条件。
A:位移函数应能反映单元的常应变状态 B:位移函数应包含刚体位移 C:位移函数在单元内要连续,在单元边界上要协调。
1、有限元是近似求解一般连续场问题的数值方法2、有限元法将连续的求解域离散为若干个子域,得到有限个单元,单元和单元之间用节点连接3、直梁在外力的作用下,横截面的内力有剪力和弯矩两个.4、平面刚架结构在外力的作用下,横截面上的内力有轴力、剪力、弯矩 .5、进行直梁有限元分析,平面刚架单元上每个节点的节点位移为挠度和转角6、平面刚架有限元分析,节点位移有轴向位移、横向位移、转角。
7、在弹性和小变形下,节点力和节点位移关系是线性关系。
8、弹性力学问题的方程个数有15个,未知量个数有15个。
9、弹性力学平面问题方程个数有8,未知数8个。
10、几何方程是研究应变和位移之间关系的方程11、物理方程是描述应力和应变关系的方程12、平衡方程反映了应力和体力之间关系的13、把经过物体内任意一点各个截面上的应力状况叫做一点的应力状态14、9形函数在单元上节点上的值,具有本点为_1_.它点为零的性质,并且在三角形单元的任一节点上,三个行函数之和为_1_15、形函数是_三角形_单元内部坐标的_线性_函数,他反映了单元的_位移_状态16、在进行节点编号时,同一单元的相邻节点的号码差尽量小.17、三角形单元的位移模式为_线性位移模式_-18、矩形单元的位移模式为__双线性位移模式_19、在选择多项式位移模式的阶次时,要求_所选的位移模式应该与局部坐标系的方位无关的性质为几何_各向同性20、单元刚度矩阵描述了_节点力_和_节点位移之间的关系21、矩形单元边界上位移是连续变化的1. 诉述有限元法的定义答:有限元法是近似求解一般连续场问题的数值方法2. 有限元法的基本思想是什么答:首先,将表示结构的连续离散为若干个子域,单元之间通过其边界上的节点连接成组合体。
其次,用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。
3. 有限元法的分类和基本步骤有哪些答:分类:位移法、力法、混合法;步骤:结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移。
有限元法理论及应用大作业1、试简要阐述有限元理论分析的基本步骤主要有哪些?答:有限元分析的主要步骤主要有:(1)结构的离散化,即单元的划分;(2)单元分析,包括选择位移模式、根据几何方程建立应变与位移的关系、根据虚功原理建立节点力与节点位移的关系,最后得到单元刚度方程;(3)等效节点载荷计算;(4)整体分析,建立整体刚度方程;(5)引入约束,求解整体平衡方程。
2、有限元网格划分的基本原则是什么?指出图示网格划分中不合理的地方。
题2图答:一般选用三角形或四边形单元,在满足一定精度情况,尽可能少一些单元。
有限元划分网格的基本原则:1.拓扑正确性原则。
即单元间是靠单元顶点、或单元边、或单元面连接2.几何保持原则。
即网络划分后,单元的集合为原结构近似3.特性一致原则。
即材料相同,厚度相同4.单元形状优良原则。
单元边、角相差尽可能小5.密度可控原则。
即在保证一定精度的前提下,网格尽可能的稀疏一些。
(a)(b)中节点没有有效的连接,且(b)中单元边差相差很大。
(c)中没有考虑对称性,单元边差很大。
3、分别指出图示平面结构划分为什么单元?有多少个节点?多少个自由度?题3图答:(a )划分为杆单元, 8个节点,12个自由度。
(b )划分为平面梁单元,8个节点,15个自由度。
(c )平面四节点四边形单元,8个节点,13个自由度。
(d )平面三角形单元,29个节点,38个自由度。
4、什么是等参数单元?。
答:如果坐标变换和位移插值采用相同的节点,并且单元的形状变换函数与位移插值的形函数一样,则称这种变换为等参变换,这样的单元称为等参单元。
5、在平面三节点三角形单元中,能否选取如下的位移模式,为什么?(1).⎪⎩⎪⎨⎧++=++=26543221),(),(y x y x v yx y x u αααααα (2). ⎪⎩⎪⎨⎧++=++=2652423221),(),(yxy x y x v yxy x y x u αααααα 答:(1)不能,因为位移函数要满足几何各向同性,即单元的位移分布不应与人为选取的 坐标方位有关,即位移函数中的坐标x,y 应该是能够互换的。
有限元法基础习题答案有限元法是一种常用的工程分析方法,广泛应用于结构力学、热传导、流体力学等领域。
它通过将复杂的物理问题离散化为一系列简单的子问题,并利用数值方法求解这些子问题,从而得到整体问题的近似解。
在学习有限元法的过程中,习题是必不可少的一环。
本文将给出一些有限元法基础习题的答案,希望能够帮助读者更好地理解和掌握这一方法。
习题一:一维线性弹性力学问题考虑一根长度为L的弹性杆,杆的截面积为A,杨氏模量为E。
在杆的一端施加一个沿杆轴向的拉力F,另一端固定。
假设杆轴向变形u(x)满足以下方程:EAu''(x) = -F,0 < x < Lu(0) = 0, u(L) = 0其中,u''(x)表示u(x)对x的二阶导数。
解答:根据上述方程,我们可以得到杆的位移函数u(x)的表达式。
首先,对方程两边进行积分,得到:EAu'(x) = -Fx + C1其中,C1为积分常数。
再次对方程两边进行积分,得到:EAu(x) = -F/2*x^2 + C1*x + C2其中,C2为积分常数。
根据边界条件u(0) = 0,可得C2 = 0。
代入边界条件u(L) = 0,可得:EAu(L) = -F/2*L^2 + C1*L = 0由此可得C1 = F/2*L。
将C1代入上式,可得:EAu(x) = -F/2*x^2 + F/2*L*x最终得到杆的位移函数u(x)的表达式为:u(x) = (-F/2*E)*(x^2 - L*x),0 < x < L习题二:二维平面弹性力学问题考虑一个正方形薄板,边长为L,板的厚度为h。
假设薄板的杨氏模量为E,泊松比为ν。
在薄板的一侧施加一个沿法向的均匀表面压力P,另一侧固定。
求薄板的位移和应力分布。
解答:根据平面弹性力学理论,我们可以得到薄板的位移和应力分布。
首先,根据杨氏模量E、泊松比ν和薄板的厚度h,可以计算出薄板的弹性模量D:D = E*h^3 / (12*(1-ν^2))接下来,根据薄板的边界条件和平衡方程,可以得到薄板的位移和应力分布。
《有限元法及其应用》课后习题目录第1章绪论 (3)第2章有限单元法理论基础 (4)第3章杆系结构单元 (5)第4章平面三角形单元 (7)第5章平面四边形等参数单元 (9)第6章常用有限元软件及其在岩土工程中的应用 (10)第1章绪论1-1试说明有限元法解题的基本思路。
1-2试说明用有限元法解题的主要步骤。
1-3有限元法主要有哪些优点?第2章有限单元法理论基础2-1 何为虚功,虚功原理的具体思路是什么?2-2 虚功原理的适用条件有哪些?2-3 位移模式的概念是什么?2-4 如何构造位移模式?2-5 弹性力学问题的求解需要满足哪些条件?第3章 杆系结构单元3-1 推导横截面积为A 的一维桁架结构的单元刚度矩阵。
3-2 图示(见题图3-1)为一平面超静定桁架结构,在载荷P 作用下,求各杆件的轴力。
此结构可看成由14、24、34三个杆单元组成,每个杆单元的两端为杆单元的结点,各结点的水平、铅直位移分别用u 、v 表示。
题图3-1 平面超静定桁架结构a —平面结构;b —单元组成;c —各结点位移3-3 图示(见题图3-2)刚架中,两杆为尺寸相同的等截面杆件,横截面面积为20.5m A =,截面惯性矩为41m 24I =,弹性模量7310kPa E =⨯,求解此结构。
题图3-2 等截面刚架结构第4章平面三角形单元4-1 按位移求解的有限单元法中:(1)应用了哪些弹性力学的基本方程?(2)应力边界条件及位移边界条件是如何反映的?(3)力的平衡条件是如何满足的?(4)变形协调条件是如何满足的?4-2 在有限单元法中,如何应用虚功原理导出单元内的应力和结点力的关系式,并将外荷载静力等效地变换为结点荷载?4-3 为了保证有限单元法解答的收敛性,平面三角形单元位移模式应满足哪些条件?μ=,记杨氏弹性模4-4 题图4-1所示等腰直角三角形单元,设14量为E,厚度为t,求形函数矩阵[]N、应变矩阵[]B、应力矩阵[]S与单元刚度矩阵[]eK。
1.1有限单元法中“离散”的含义是什么?有限单元法是如何将具有无限自由度的连续介质问题转变成有限自由度问题的?位移有限元法的标准化程式是怎样的?(1)离散的含义即将结构离散化,即用假想的线或面将连续体分割成数H有限的单元,并在其上设定有限个节点;用这些单元组成的单元集合体代替原来的连续体, 而场函数的节点值将成为问题的基本未知量。
(2)给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律,即通过插值以单元节点位移表示单元内任意点的位移。
因节点位移个数是有限的,故无限自由度问题被转变成了有限自由度问题。
(3)有限元法的标准化程式:结构或区域离散,单元分析,整体分析,数值求解。
1.3单元刚度矩阵和整体刚度矩阵各有哪些性质?各自的物理意义是什么?两者有何区别?单元刚度矩阵的性质:对称性、奇异性(单元刚度矩阵的行列式为零)。
整体刚度矩阵的性质:对称性、奇异性、稀疏性。
单元Kij物理意义Kij即单元节点位移向量中第j个自山度发生单位位移而其他位移分量为零时,在第j个自111度方向引起的节点力。
整体刚度矩阵K中每一列元素的物理意义是:要迫使结构的某节点位移自由度发生单位位移,而其他节点位移都保持为零的变形状态,在所有个节点上需要施加的节点荷载。
2.2什么叫应变能?什么叫外力势能?试叙述势能变分原理和最小势能原理,并回答下述问题:势能变分原理代表什么控制方程和边界条件?其中附加了哪些条件?(1)在外力作用下,物体内部将产生应力C和应变£,外力所做的功将以变形能的形式储存起来,这种能量称为应变能。
⑵外力势能就是外力功的负值。
(3)势能变分原理可叙述如下:在所有满足边界条件的协调位移中,那些满足静力平衡条件的位移使物体势能泛函取驻值,即势能的变分为零8 Up二8 U e + 5 V二0此即变分方程。
对于线性弹性体,势能取最小值,即6 2UP二§2U £ + 6 2VN0此时的势能变分原理就是著名的最小势能原理。
实用文档有限元法基础与程序设计李亚智课后答案1.诉述有限元法的定义?答:有限元法是近似求解一般连续场问题的数值方法2.有限元法的基本思想是什么?答:首先,将表示结构的连续离散为若干个子域,单元之间通过其边界上的节点连接成组合体。其次,用.每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。3.有限元法的分类和基本步骤有哪些?答:分类:位移法、力法、混合法;步骤:结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移。4.有限元法有哪些优缺点?答:优点:有限元法可以模拟各种几何形状复杂的结构,得出其近似解;通过计算机程序,可以广泛地应用于各种场合;可以从其他CAD软件中导入建好的模型;数学处理比较方便,对复杂形状的结构也能适用;有限元法和优化设计方法相结合,以便发挥各自的优点。缺点:有限元计算,尤其是复杂问题的分析计算,所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。对无限求解域问题没有较好的处理办法。尽管现有的有限元软件多数使用了网络自适应技术,但在具体应用时,采用什么类型的单元、多大的网络密度等都要完全依赖适用者的经验。5.梁单元和平面钢架结构单元的自由度由什么确定?答:由每个节点位移分量的总和确定6.简述单元刚度矩阵的性质和矩阵元素的物理意义?答:单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵7.有限元法基本方程中的每一项的意义是什么P14?答:Q整个结构的节点载荷列阵(外载荷、约束力) ; 整个结构的节点位移列阵;结构的整体刚度矩阵,又称总刚度矩阵。8.位移边界条件和载荷边界条件的意义是什么?答:由于刚度矩阵的线性相关性不能得到解,引入边界条件,使整体刚度矩阵求的唯一解。9.简述整体刚度矩阵的性质和特点P14?答:对称性;奇异性;稀疏性;对角线上的元素恒为正。10简述整体坐标的概念P25?答:在整体结构上建立的坐标系叫做整体坐标,又叫做统一坐标系。。
传热学_燕山大学中国大学mooc课后章节答案期末考试题库2023年1.求解导热问题时,常用的数值计算法有:()参考答案:有限单元法_有限差分法2.导热问题数值求解的基本步骤包括:()参考答案:代数方程组的求解_问题的数学描述_区域离散_节点离散方程的建立3.求解传热问题的基本方法有哪些?()参考答案:理论分析_实验测量_数值计算4.解析法的特点:()参考答案:能够得到精确解_解有普遍性5.采用有限差分法计算导热问题,可以得到节点处温度值的精确解。
参考答案:错误6.有限差分法较有限元法更适合求解不规则物体边界的导热问题。
参考答案:错误7.高斯-赛德尔迭代法比简单迭代法计算精度高。
参考答案:错误8.下面()种说法是不正确的。
参考答案:同类物理现象就是彼此相似的现象9.流体横掠单管时,相似准则中的定型尺寸取为()。
参考答案:管外径10.流体在圆管内流动时,相似准则中的定型尺寸取为()。
参考答案:管内径11.对流换热系数不是流体的物性参数,是一个反应对流换热复杂过程的多变量的函数。
参考答案:正确12.自然对流换热与强制对流换热都是对流换热,属于是同类物理现象。
参考答案:错误13.不符合肋片描述的是()参考答案:减小换热量作用14.对温度梯度描述错误的是()参考答案:沿温度增加方向为负15.黑体是指能够吸收投射到其表面上的全部辐射能的物体。
参考答案:正确16.导热是固体内部或固体之间的传热过程。
参考答案:错误17.工业上,多数工程材料作为灰体处理不会引起较大的误差。
参考答案:正确18.导热问题数值求解的基本思想是:()参考答案:用求解区域上离散点的温度集合代替连续的温度场,计算难度小_用代数方程代替微分方程19.一般而言,与非金属材料相比,金属材料的导热系数()参考答案:较大20.五种具有实际意义的换热过程为:导热、对流换热、复合换热、传热过程和()参考答案:辐射换热21.热量传递的三种基本方式是()参考答案:热对流、导热、热辐射22.由表面 1 和表面 2 组成的封闭系统中,【图片】()【图片】参考答案:可能大于、等于或小于23.流体纵掠平板时,对流动边界层的描述,下列第()种说法是不准确的。
1.1 有限单元法中“离散”的含义是什么?有限单元法是如何将具有无限自由度的连续介质问题转变成有限自由度问题的?位移有限元法的标准化程式是怎样的?(1)离散的含义即将结构离散化,即用假想的线或面将连续体分割成数目有限的单元,并在其上设定有限个节点;用这些单元组成的单元集合体代替原来的连续体,而场函数的节点值将成为问题的基本未知量。
(2)给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律,即通过插值以单元节点位移表示单元内任意点的位移。
因节点位移个数是有限的,故无限自由度问题被转变成了有限自由度问题。
(3)有限元法的标准化程式:结构或区域离散,单元分析,整体分析,数值求解。
1.3 单元刚度矩阵和整体刚度矩阵各有哪些性质?各自的物理意义是什么?两者有何区别?单元刚度矩阵的性质:对称性、奇异性(单元刚度矩阵的行列式为零)。
整体刚度矩阵的性质:对称性、奇异性、稀疏性。
单元 Kij 物理意义 Kij 即单元节点位移向量中第 j 个自由度发生单位位移而其他位移分量为零时,在第 j 个自由度方向引起的节点力。
整体刚度矩阵 K 中每一列元素的物理意义是:要迫使结构的某节点位移自由度发生单位位移,而其他节点位移都保持为零的变形状态,在所有个节点上需要施加的节点荷载。
2.2 什么叫应变能?什么叫外力势能?试叙述势能变分原理和最小势能原理,并回答下述问题:势能变分原理代表什么控制方程和边界条件?其中附加了哪些条件?(1)在外力作用下,物体内部将产生应力σ和应变ε,外力所做的功将以变形能的形式储存起来,这种能量称为应变能。
(2)外力势能就是外力功的负值。
(3)势能变分原理可叙述如下:在所有满足边界条件的协调位移中,那些满足静力平衡条件的位移使物体势能泛函取驻值,即势能的变分为零δ∏p=δ Uε+δV=0此即变分方程。
对于线性弹性体,势能取最小值,即δ2∏P=δ2Uε+δ2V≥0此时的势能变分原理就是著名的最小势能原理。
势能变分原理代表平衡方程、本构方程和应力边界条件,其中附加了几何方程和位移边界条件。
2.3 什么是强形式?什么是弱形式?两者有何区别?建立弱形式的关键步骤是什么?等效积分形式通过分部积分,称式∫ΩCT(v)D(u)dΩ+∫ΓET(v)F(u)dΓ为微分方程的弱形式,相对而言,定解问题的微分方程称为强形式。
区别:弱形式得不到解析解。
建立弱形式的关键步骤:对场函数要求较低阶的连续性。
2.4 为了使计算结果能够收敛于精确解,位移函数需要满足哪些条件?为什么?只要位移函数满足两个基本要求,即完备性和协调性,计算结果便收敛于精确解。
2.6 为什么采用变分法求解通常只能得到近似解?变分法的应用常遇到什么困难?Ritz 法收敛的条件是什么?(1)在 Ritz 法中,N 决定了试探函数的基本形态,待定参数使得场函数具有一定的任意性。
如果真实场函数包含在试探函数之内,则变分法得到的解答是精确的;如果试探函数取自完全的函数序列,则当项数不断增加时,近似解将趋近于精确解。
然而,通常情况下试探函数不会将真实场函数完全包含在内,实际计算时也不可能取无穷多项。
因此,试探函数只能是真实场函数的近似。
可见,变分法就是在某个假定的范围内找出最佳解答,近似性就源于此。
(2)采用变分法近似求解,要求在整个求解区域内预先给出满足边界条件的场函数。
通常情况下这是不可能的,因而变分法的应用受到了限制。
(3)Ritz 法的收敛条件是要求试探函数具有完备性和连续性,也就是说,如果试探函数满足完备性和连续性的要求,当试探函数的项数趋近于无穷时,则 Ritz 法的近似解将趋近于数学微分方程的精确解。
3.1 构造单元形函数有哪些基本原则?形函数是定义于单元内坐标的连续函数。
单元位移函数通常采用多项式,其中的待定常数应该与单元节点自由度数相等。
为满足完备性要求,位移函数中必须包括常函数和一次式,即完全一次多项式。
多项式的选取应由低阶到高阶,尽量选择完全多项式以提高单元的精度。
若由于项数限制而不能选取完全多项式时,也应使完全多项式具有坐标的对称性,并且一个坐标方向的次数不应超过完全多项式的次数。
有时为了使位移函数保持一定阶次的完全多项式,可在单元内部配置节点。
然而,这种节点的存在将增加有限元格式和计算上的复杂性,除非不得已才加以采用。
形函数应保证用它定义的位移函数满足收敛要求,即满足完备性要求和协调性条件。
3.1 构造单元形函数有哪些基本原则?试采用构造单元的几何方法,构造 T10 单元的形函数,并对其收敛性进行讨论。
通常单元位移函数采用多项式,其中的待定常数由节点位移参数确定,因此其个数应与单元节点自由度数相等。
根据实体结构的几何方程,单元的应变是位移的一次导数。
为了反映单元刚体位移和常应变即满足完备性要求,位移函数中必须包含常数项和一次项,即完全一次多项式。
3.3 何谓面积坐标?其特点是什么?为什么称其为自然坐标或局部坐标?(1)三角形单元中,任一点 P(x,y)与其 3 个角点相连形成 3 个子三角形,其位置可以用下述称为面积坐标的三个比值来确定:L1=A1/A L2=A2/A L3=A3/A其中 A1,A2,A3 分别为 P23,P31,P12 的面积。
(2)面积坐标的特点:a T3 单元的形函数 Ni 就是面积坐标 Lib 面积坐标与三角形在整体坐标系中的位置无关。
c 三个节点的面积坐标分别为节点 1(1, 0, 0)、节点 2(0, 1, 0)、节点 3(0, 0, 1),形心的面积坐标为(1/3, 1/3, 1/3)。
d 单元边界方程为 Li=0(i=1,2,3)e 在平行于 23 边的一条直线上,所有点都有相同的面积坐标 L1(L1 对应的三角形具有相同的高和底边),而且 L1 就等于此直线至 23 边的距离与节点 1 至 23 边的距离之比值。
f 面积坐标与直角坐标互为线性关系。
(3)面积坐标与三角形在整体坐标系中的位置无关,故称为局部坐标或自然坐标。
4.1 与平面问题相比,轴对称问题有何特点?在有限元表达格式上有何区别?轴对称问题是空间问题的一种特殊情况,结构的几何形状、约束条件及荷载分布都对称于某个轴,其位移、应变、应力等也对称于此轴,而与环向坐标无关。
4.2 试用体积坐标构造 10 节点四面体单元的形函数并讨论收敛性。
5.1 何谓等参单元?等参单元具有哪些优越性?在等参单元计算中,数值积分的阶次是否越高越好?为什么?等参单元(简称等参元)就是坐标变换和单元内的等变量(通常是位移函数)采用相同的节点参数和相同的插值函数进行变换而设计出的一种单元。
优越性:一,有些工程结构的形状比较复杂,如果用直边单元离散这些结构将需要大量的单元才能得到较好的近似,而曲边的等参单元可非常方便的离散复杂结构。
二,如果在单元内多取些节点,单元便具有较多的位移自由度,从而就能够插值表示较复杂的单元内部位移场,这样也就提高了单元本身的精度。
三,等参单元刚度矩阵、荷载矩阵的计算是在规则单元域内进行的,因此不管被积函数多么复杂都可方便的采用标准化数值分析。
在等参单元计算中,数值积分的阶次并不是越高越好,5.6 何谓位移的零能模式?在什么条件下会发生零能模式?对应于某种非刚体位移模式,减缩积分时高斯点上的应变正好等于零,此时的应变能当然也为零,这种非刚体位移模式称为零能模式。
采用减缩积分时会发生零能模式。
6.1 对于杆系结构单元,为什么要在局部坐标系内建立单元刚度矩阵?为什么还要坐标变换?(1)在局部坐标系内可以更方便的建立单元刚度矩阵。
(2)在整体分析中,对所有单元都应采用同一个坐标系即整体坐标系 X Y,否则围绕同一节点的不同单元对节点施加的节点力不能直接相加。
因此,在进行整体分析之前,还需要进行坐标转换工作,把局部坐标系中得出的单元刚度方程转换成整体坐标系中的单元刚度方程,从而得出整体坐标系中的单元刚度矩阵。
6.2 有哪几种梁弯曲理论?如何用中性轴位移确定梁内任一点的位移?工程梁理论、剪切梁理论、通用梁理论、空间梁理论。
梁弯曲理论(包括工程梁理论和剪切梁理论)在弹性力学基本假定的基础上引入了某些附加假定,将问题归结为求解中性轴位移,而梁内任一点的位移都可以通过中性轴位移来表示。
7.1 在薄板弯曲理论中做了哪些假设?如何用中面位移确定板内任一点的位移?假设:(1)板厚度方向的挤压变形可忽略不计,即εZ=0。
(2)在板弯曲变形中,中面法线保持为直线,且仍为弹性曲面(挠度曲面)的法线,即直法线假设。
(3)薄板中面只发生弯曲变形,没有面内的伸缩变形,即中面水平位移。
(u)z=0=(v)z=0 =0 薄板的全部位移、应力和应变分量都可以用板的挠度ω来表示,而薄板小挠度弯曲被简化为中面的弯曲问题,只要中面挠度ω确定,任何点的位移都可确定。
薄板内不等于零的应变分量有如下三个:εx=бu/бx=-z б2ω/бx2εy=бv/бy=-б2ω见P116,式(7.3a) r xy=бu/бy+бv/бx=-2z б2ω/бxбy7.2 薄板单元和厚板单元的基本假设有什么不同?各自是怎样选择节点位移参数的?不同点:薄板单元假设横向纤维无挤压,板的中面法线变形后仍保持为直线,该直线垂直于变形后的中面,但是厚板单元的假设考虑横向变形的影响,板的中面法线变形后仍基本保持为直线,但该直线不再垂直于变形后的中面,法线绕坐标轴的转角不再是挠度的导数,而是独立的变量。
7.3 在薄板单元中,节点力矩与薄板内力有何区别?节点力矩 Mxi, Myi 是集中力矩,而板内力矩 Mx, My 是分布力矩,此外,两者的正负号规定也不相同,因为 Mx, My 与应力正负号的规定相应。
8.1 薄壳理论有哪些假设?与薄板理论的假设有何异同?厚壳分析中引入了何种假设?与厚板理论的假定有何异同?薄壳理论的假设:薄壳发生微笑变形时,忽略沿壳体厚度方向的挤压变形;且认为直法线假设成立,即变形后中面法线保持为直线且仍为中面的法线;壳体变形时中面不但发生弯曲,而且面内也将产生面内伸缩变形;折板假设;非耦合假设。
与薄板理论的假设的相同点:直法线假设和法向(板厚度方向)的纤维无挤压假设均成立。
不同点:薄板中面只发生弯曲变形,没有面内的伸缩变形,即中面水平位移为零,而壳体变形时中面不但发生弯曲,而且也将产生面内伸缩变形。
厚壳分析的假设:变形前后的中曲面法线变形后仍基本保持为直线,但因横向剪切变形的缘故,该直线不再垂直于变形后的中曲面,此外,壳体厚度方向的挤压变形可以忽略。
与厚板理论的假设的相同点:中面法线变形后仍基本保持为直线,但因横向剪切变形的缘故,该直线不在垂直于变形后的中面。
厚度方向的挤压变形忽略不计。
不同点:厚板理论的假设中,中面内的线位移可以忽略,而厚壳理论的假设中,中面内的位移不可忽略,并且厚壳的位移场可用中面位移表示。