《有限单元法》复习参考题

1、有限单元法复习参考题一、简答题：1、简述应用有限单元法解决具体问题的要点。（1）将一个表示结构或者连续体的求解域离散为若干个子域（单元），并通过他们边界上的结点相互结合为组合体。（2）用每个单元内所假设的近似函数来分片地表示全求解域内待求的未知场变量。而每个单元内的近似函数由未知场函数（或及其导数，为了叙述方便，后面略去此加注）在单元各个节点上的数值与其对应的插值函数来表达。（3）通过和原问题数学模型（基本方程、边界条件）等效的变分原理或者加权余量法，建立求解基本未知量（场函数的结点值）的代数方程或者常微分方程组。2、等效积分形式和等效积分“弱”形式的区别何在？为什么等效积分“弱”形式在数值分

2、析中得到更多的应用？在很多情况下对微分方程的等效积分形式进行分部积分可以得到等效积分的弱形式，如下式CYD）D（u）dQ+jET（6）F（u）dT=0,其中C、D、E、F是微分算子。像这种逋过适当提高对任意函数u和6的连续性要求，以降低对微分方程场函数u的连续性要求所建立的等效积分形式称为微分方程的等效积分“弱”形式。值得指出的是，从形式上看“弱”形式对函数u的连续性要求降低了，但对于实际的物理问题却常常较原始的微分方程更逼近真正的解，因为原始微分方程往往对解提出了过分的要求。所以等效积分“弱”形式在数值分析中得到更多的应用。3、什么是Ritz（里兹）方法？其优缺点是什么？收敛的条件是什么？基

3、于变分原理的近似解法称为Ritz（里兹），解法如下：1.求解步麻：d（“）4 Ktat = Ft=0da. t3）求解线性代数方程组勺 一A u的近似解1）假设近似解：“=且=£一一3为待定参数，满足强制边界条件.2）将百代人中。W4I泛函/（«）的极值问题（求函数U）,转化为求多元（"，0,）函数的极值问题。优缺点：一般来说，使用里兹方法求解，当试探函数族的范围扩大以及待定参数的数目增多时，近似解的精度将会提高。局限性：（1）在求解域比较复杂的情况下，选取满足边界条件的试探函数，往往会产生难以克服的困难。（2）为了提高近似解的精度，需要增加待定参数，即增加试探函

4、数的项数，这就增加了求解的复杂性，而且由于试探函数定义于全域，因此不可能根据问题的要求在求解域的不同部位对试探函数提出不同精度的要求，往往由于局部精度的要求使整个问题求解增加许多困难。收敛的条件：试探函数NrN”.，Nn应取自完备函数系列，试探函数NN2，.,Nq应满足Cm_1连续性要求。4、什么是最小住能原理？该原理在有限单元分析中的作用是什么？对场函数的试探函数有什么要求？如此公式所示Blip=0,Blip是系统的总位能，它是弹性体变形位能和外力位能之和。该式表明，在所有区域内连续可导的并在边界上满足给定位移条件的可能位移中，真实位移使系统的总位能取驻值。在所有的可能位移中，真实位移使系统

5、总位能取得最小值，因此Blip=0所表达的称为最小位能原理。利用最小位能原理求得位移近似解的弹性变形能是精确解变性能的下界，即近似的位移场在总体上偏小，也就是说结构的计算模型显得偏于刚硬。要求：最小位能原理的试探函数-位移，应事先满足几何方程和给定的位移的边界条件。5、有限单元法中单元的位移模式为什么通常采用多项式作为近似函数？选择广义坐标有限元位移模式的一般原则是什么？因为多项式运算简便，并且随着项数的增多，可以逼近任何一段光滑的函数曲线，多项式的选择应由低次到高次。1）大多数函数可用泰勒级数展开，根据需要取前几项逼近真实的位移函数解；2）多项式函数可保持各向同性，不偏惠某一坐标方向；3）多

6、项式函数便于积分和微分，使有限元公式简单、直观。4）多项式函数很容易满足收敛准则。一般原则：（1）广义坐标是由结点场变量确定的，因此它的个数应与结点自由度数相等。（2）选取多项式时，常数项和坐标的一次项必须完备。（3）多项式的选取应由低阶到高阶，尽量选取完全多项式以提高单元的精度。6、在有限单元法中，保证有限元解收敛有哪些准则？六节点三角形单元是收致的单元吗？为什么？完备性要求。如果出现在泛函中的场函数的最高阶导数是m阶，则有限元解收敛的条件之一是单元内场函数的试探函数至少是m次完全多项式。协调性要求。如果出现在泛函中的最高阶导数是m阶，则试探函数在单元交界面上必须具有Cm连续性，即在相邻单元

7、的交界面上函数应有直至阶的连续导数。不是收敛单元，因为不满足完备性要求和协调性要求。7、何谓位移元？为什么位移元解具有下限性？请给出力学上的解释。位移元：以位移为基本未知量，并基于最小位能原理建立的有限元称之为位移元。位移元解具有下限性可以解释如下：单元原是连续体的一部分，具有无限多个自由度。在假定了单元的位移函数后，自由度限制为只有以结点位移表示的有限自由度，即位移函数对单元的形变进行了约束和限制，使单元的刚度较实际连续体加强了，因此连续体的整体刚度随之增加，离散后的K较实际的K为大,因此求得的位移近似解总体上将小于精确解。8、什么是拉格朗日单元和Seiendipity单元？比较这两种单元的

8、各自特点。拉格朗日单元特点：（1）插值函数构造方便；（2）内部结点较多，单元的次数越高相应自由度越高；（3）单元阶次增高，非完全高次项增加。Seienchpity单元作用是：不改变精度的条件下，减少内部结点，即对Lagrange单元简化。9、什么是阶谱单元？如何在有限单元法中采用阶谱单元？相对于通用的标准单元有何好处？阶谱单元：特点：（1）插值函数（阶谱函数）不再具有“0-1特性工（2）高阶单元的单元特性矩阵可承袭低阶单元的单元特性矩阵。在用于自适应分析中可以节省编程的工作量。10、什么是等参变换？在有限单元法中，等参数单元的主要优点是什么？等参变换是指单元的几何形状和单元内的场函数采用相同数

9、目的结点参数及相同的插值函数进行变换。优点是借助于等参元可以对于一般的任意几何形状的工程问题方便地进行有限元离散。等参元的插值函数是用自然坐标给出的，等参元的一切计算（如单元刚度矩阵、单元载荷列阵等）都是在自然坐标系中规格化的母单元内进行，相关运算大大简化。不管各个积分形式的矩阵的被积函数如何复杂，都可以采用标准化的数值积分方法计算，从而使工程问题的有限元分析纳入了统一的通用化程序。11、等参元计算中数值积分阶次的选择应遵循哪些原则？如何检查所采用的积分方案是否满足所述的原则？1.保证积分的精度。2.保证结构总刚度矩阵k是非奇异的。对于一个给定形式的单元，如果采用精确积分，则插值函数中所有项次

10、在1人常数的条件下能被精确积分，并能保证刚度矩阵的非奇异性。如果采用减缩积分，因为插值函数中只有完全多项式的项次能被精确积分，因此需要进行刚度矩阵非奇异必要条件的检查。若能通过检查，则可以考虑采用减缩积分方案，以减少计算工作量，并可能对计算结果有所改进。12、简述有限元网格划分的基本原则。网格疏密的布置，不连续处的网格自然划分，不同密度划分网格过渡13、什么是自适应分析方法？用什么方法进行自适应的重分析？自适应有限元技术是一种根据中间计算结果自动控制计算过程的求解偏微分方程的方法。它主要利用中间计算结果自动计算所需的网格，选取最佳离散方式，从而逐步对误差做适当调节以达到所需精度。h型改进，p型

11、改进14、为什么双线性四边形单元用于弯曲应力分析时表现出较差的性能？不能有重节点不能出现内角大于180。的情况内角最好介于30。-150。之间（有限变形的情况）15、什么是罚函数法？罚函数法求解近不可压缩弹性力学问题时的有限元方程系数矩阵应具有什么性质？如何保证它具有这样的性质？kl非奇异，k2奇异，kl+ak2非奇异二、计算分析题：1、试写出下述定解问题的等效积分形式和等效积分弱形式，并说明构造“弱”形式的意义。(需+器)="")加Qu(x,y) = u(x,y)on r =（提示：利用Green公式：.vdxdy=，dS-口,n为so的单位外法向1、解答：原定解问题的等

12、效积分形式为JJ（-）】必丫办+（“（X,了）一u（x,y）yc/S=0Q6其中,v=v（x,>9r。=WXy）为任意函数,且分别在区域C及边界为上可机由Green公式可得;厂力一】煞+L(羽刃 u(x9y)fvdS = 0这里，要求口和在区域Q上可积。上式即为等效积分“弱”形式。构造“弱”形式的意义：看课件PPT2、已知：A(。)=普+Q(x)=0(0<X<L)?其中(aJxQ(M=1(0<x<172),zAxn期。(L/2<X4L),边界条件沏阿尸°；瓦=10x=L假设近似函数为。5)=。0+。逮+42/+%/,试用配点法，子域法和伽辽金法求解

13、。2、解答：设吠幻二&十十十口3,余量R(x)=-+Q(x)dr由边界条件0(0)=0,可得劭=o；由"ax=10可得A1+2a2L+3aL2-10=0(a)配点法：取x=L/3和2L/3为配点,要求:/?(£/?>=0(2)/?(2Z73)=0(3)町+2啊L+3町L=1。329+2与匚+1=0,29+4叫1=0,«1=10+1.1,=JL«2=4，22L解方程组(1>(3),可得为=10+L/2M2=-1必=:NL(b)子域法：取04W2和为子域，则山2|R(x)cix=0；M=°(5)2Rai02旬+2与L+3由L=1

14、0,与L+一包L+_L=0,«2+-43=0424.=10+!么6=_与=三23L4解方程组,(4),(5),可得卬=10+£/2,%=-3/4,%=；（b）伽辽金法：取权函数%=x，匕=/,匕=/,则-10)=0x-£-10)=0=0Y-f.(6)(8)）也2d3L3+L(2d2+l)L2+1a2L2-Le1+2a2L+3d3-10徵+2 与 L + 3 的 1? - 10) = 0,6丫5,1c4上15t4t35sL+瑟(2与+1)L+_«2-L齐广+孑虫+1）八"己1"+2微L+3由-10）=015.2317Txing=?=一L

15、+10P百T为37132175解方程组(6)(8),可得=10十一匕出=-23/32,«3=3、某问题的微分方程是至+至+=+。=0mQ,边界条件为/-。on=1=qon匚，。-其中，C和。仅是坐标的函数，证明此方程的微分算子是自伴随的，并建立相应的自然变分原理。3、解答；微分算子为2=冒十紫十项，取任意函数小，公十不t+ cv dxdydu dv du dv十刖，卜力山"打.十哈人d2u d2ndx2 dy2vnt + vn it (is 4- &，dsex cy J J dn=JJvL（u）dxdyvs 十?da故算子是自伴随的。原问题等价于：（假设。己满足上的

16、边界条件）(av国7+而十 c 十。dxdy 一ft 的G（谢）dA 3（邮）d迎中誉望4“邮+。叼"3_* ds山嚏樽卜黑+ c轮dxdy + J军师杰-/-c(jf2 + Q( dxdy + f 前陋2%=6哈（沙聘+c© 4- 0, dxdy + J cj(j)cis i2v1=0等价的自然变分原理为:二51但© 一4翌/.十0y办+/押$2Jr2砥氏舒圉- c4 - 2QG dxdy - £(j(/fds4、弹性薄板的控制方程为：装+ 2 及4d4w d4wdx2dy2 + dy幺，建立周边固支时的自然D变分原理。必富+2嘉+券匏/，=°

17、;,其中“，满足强迫边界条件.JL署苏s=TL笠"年）加”+f翳兄振"JJnexJJndvexJex三JL零6詈）心仍一答?况当）“几票输心力TL等陪必曲+f需小小=乩察次鲁地心-一城（勖八JJ°dydydy=f£82w小|加”"-JL卷蹊）2+于舂科门MWp-rb（即）dTJexey,ex工上加小小=f也母次直乙以包）“dxdy22改办dxdy，小力勿相加得：山富C券+翳-舒”坊dxexdyrdyD17rd"'"d2卬济卬e'v0*1,fc2wu/We=JJ/f”）十十2T6T7r卜-杰凶9Fe（=）/rJ

18、）。exexoyaydxdxeyfDJexdxT吆冏把）"7丘"（旦）女",s伫）打Jdy2ydyJdxeyydxJdxdyxdy"一Mj山岩）工（会）、（西牙一生四加7g曳沆电）打卢丝5（包）/J"2fix"2-dxyD丁Gxdndxdydndy=/乩:（彳彳尸十：（不）.十（二）.一点时小'=oJ也2/2uyoxeyD故原问题的自然变分原理为:n<w，>=山5（获+2（r+W3中阳,或者口的J片喘才+万（羽）”（e%dxd）-qwdxdy5、如有一问题的泛函为+kw22+ qwdx,其中，民/,k是常数，q是给

19、定函数，W是未知函数，试导出原问题的微分方程和边界条件。5、解答：泛函的一阶变分为矶=演1：（与+殍小）=j仪+当+g“去=J（Eiw+上”曲仙+q（iv）dx两次利用分部枳分公式，我们有fEhi8vdx=Elw8>vElir+Elw'SwdxEQIoJu所以，加I=£7/6w-Elh6闻十(£Av<41+hv+wdx由汨=0可得Ei iv t> w =0. E! vv3 二0及Io£（£/h>C4'+h，+wdx=0上述各式对任给的变分及3”，成立，由变分法预备定理可得:微分方程为：£Av<4）+

20、夕=0边界条件为：E/M*l=0，£7M1/iJ=0（II）6、考虑如图1所示悬臂梁，设其跨长为/,抗弯刚度为七/,在梁的中部x及2端点x=/处受集中荷载尸作用。（1）若用Ritz（里兹）法计算粱的挠度曲线方程，试问：是否可取如下表达式？卬=a（lcos?F）,其中，。为待定常数。1可取挠度曲线表达式为w=a（l-cos宁），仅需要验证该挠度函数满足位移边界条件即可，即显然满足。（2）若是可以，试利用最小位能原理求出相应的挠度曲线方程。图1（2）用最小势能原理确定待定常数0。弯矩为：，M（x）=-E/会二cZrII梁的应变能为：u=，一口，3去二空二2Efh4"外力势能为：

21、V=一。（用/一P（M，）I二一3X-系统的总势能为：nw”二需由<sn=s（u+v）=o,得：型9=0.即de6"户由此解得；=也7£7/6PI'八7TX：.W=（1-COS）E/I7、证明三节点三角形单元的形状函数满足乂(七,力)=可=:",及M+N,+N/1乂 =（q+3+纣）（工时一%力）+（力-几）+（% 一）为 2A1NM2-4-1y-(.vAVw-a;jv)+(x9ry.-avw)+-x.y.)1<N+N+N-3+"/+4)+3+”+S+T=r.，EF8、设有一弹性平面问题，厚度为f,弹性模量为E,泊松比4=0,对于如图

22、2所示的三节点三角形单元，试计算其单元刚度矩阵。图29、对如图3所示四边形单元，试计算其单元刚度矩阵，写出其基本思路即可。图39、解答:基本思路：（利用等参元计算）步h将此四边形单元变换到边长为2的正方形单元（标准单元），写出相应的儿何变换及其雅可比矩阵；步2:计算在标准单元上的应变矩阵，进而计算原单元上的应变矩阵；步3：将单元刚度矩阵的计算化归为标准单元上进行，即灯=Jj卬螃10、试用“试凑法”构造图4各节点形函数，要求写出详细的计算过程。10、解答，详见课件PPT11、解答，详见课件PPT图411、利用构造变结点数单元插值函数的方法，构造如图5所示8结点单元的插值函数。4（口1）73（L1

23、）l（h-l）523-1）图5(因为单元的边是直线，可用4个结点定义单12、(1)图6所示为二次四边形单元，试计算。乂/&和。N?/在自然坐标为(1/2,1/2)的点Q的数值元的几何形状)。12、Cae1:几何变换采用4个点;在单元内有x = / Nr ,二| 11Jacobi矩阵为J=J（f，7）=此=，1 +矮）（1+ ,川）"123,4一办正加丸ax-ax-即1 . I一工 0 + 07,）% £ £ （1+，7,7r=ii=iX乙a+M况E后）x 1=1J=1一 一 o O 2 6 5 5 62.12. 1 - 4 1X7-91 - 2-11751

24、 一14 12一175-1一170-a CN =N；N；N&,1215缗aNx i Z1 八 i/i u，、i /I - =-(1)-(1-<)+-7(1)当=:,飒飒 9比训 16cN2湾二10+# +抑+ )+；("/)，当占=_1, 如=工- 1 2 结 161的 16购 网 j av,一-=a-L + ft-Lox觉 OIdN- 4力呜当弓冈二9三二双'识二 69-7 = 4760Case2:几何变换采用8个点:叫=0 0257 劄乂=(1+/(1+矶尹川4M=(17)(1+X3+-1”4%3=（1-4）（1-）（一岁一一1）/4N"(l+&l

25、t;)XW-/7-W4M.("3Xl+”2-功/2N?=(1-X1-”2M=(l-/f)(l+G/2Q (1/2.1“沙,4kQ(121/2)I.1=15.625尸54HINXN3N3.125V,NtMV15Z*1Ul =15.625513.12515j15.6253.1255T1fO.562511O.O2571150.5625/"1o.O32Ij用由优H噂5T： Jo. 18751J0.01I4115J O.O625J to.(X)l8j(2)图7所示为二次三角形单元，试计算6NJ永和。NJ。，在点P(L5,2.0)的数值。 。1(1.520,(2)%=20.2458 1

26、11520)49”之-0.347 49耳=«24-1)&=L；(2LD乂=4(24-1)十4L-(2A,b)+2AlA)2A2ACAf-12塾=以冬,义区二4L&4乙工6 第力，也力 -2A 气 2A (2A):- '1 -号=p(2A?ba 2A也)=A|dx(-2)5x 4j= 0.245t(2A,c. + 2Alr,力(2Ar1 2)=yr|4K(-3)4.5x(-1)1=-0.U713、对如图8所示的四边形单元（见左图）（1）写出将该单元变换到一边长为2的正方形单元（见右图）的坐标变换o13、解答，（1）坐标变换为:>=M8/7)=M'+

27、N/z+N/3+T-加2十加打十加加十/十沁4川十加3y=9(*)=Mm+"2+n3y3+'必小8）吐力2+抄加一加3+3+纨1+机9+33）（1+办8（2）计算该单元的雅可比（Jacobi）矩阵。要求写出详细计算过程。<2>雅可比矩阵为:立即办-即r I-> 9也可以写为班-"ar丽14、有一个三角形单元，受有如图9所示的分布载荷，试计算该单元的等效节点载荷列阵。要求写出详细的计算过程。14、解答：设单元等效节点载荷列阵为等"它由12边的等效节点载荷列阵为鸿万及23边的等效节点载荷列阵为浦f叠加而成，用'RJ十坛。首先，计算4/

28、、依据公式，有寓 J0=JA/M o M0 N1 0对三节点三角形单元，恒有N?=JN.=L.这样，K,=lOOzjM杰=100/JL.ds利用面积坐标的耳函数在三角形某一边上的积分公式可求得;5OJ类似地,有'=50，。因此，H°，50f，0,50/,0,0r，其中r为厚度。其次,计算（“JZ完全类似于上述方法，可得R13r=|0，0，100/,0,200gor，其中，为厚度。最后，计算即="4+小，=（0,50/,100r,50r,200/,0r15、图10为一给定的六节点三角形单元，在力边上作用有线性分布的面载荷（x方向），假设单元厚度为1，试用两种不同的方法

29、求单元等效节点载荷列阵。要求写出详细的计算过程。（给出均布测压的等效节点载荷计算的例子：图示结点三角形单元的124边作用有均布侧压力q,单元厚度为3求单元的等效结点荷载。3(xjJ解答;Ni=Li(2Lrl)； N2=L2(2L2-I)； N3=Lj(2I-vI)；N4=4UL2： N<=4L2L5； N产IL山;在三知形142边上建立局部坐:标示如图:-1010ee142L 41 2节点在局部坐标系下的形函数为，的把；包=（1打1 +0）；%=；（1 +六所以土/»="/第：r Ij- 2.fi ="J_|W< = -5/BJ.人=14忸灌=，;其中

30、，I为三龟形142边的边k7cos1/sin 06rZr4 小0叫 2 pcosl 2 J% -M 4J /. sin <9j = 3/sin= 3 "k . x,f2 cos a s sin。/cos/sin &4其中,U为q与水平方向的夹角°）将上述应用于该题可得到计算结果oocooo16、图11为一边长为2的八节点正方形单元，它的边界平行于整体坐标轴，在边152上受有均布表面荷载，P、=P。,假设单元厚度为试求单元等效节点载荷。要求写出详细的计算过程。图1116.解答：寻效节点我荷为:=（P）=（PJ=P=P9=，。其中,取a=1即为本腮答案c17、考虑

31、如图12所示受均布载荷作用的悬臂梁，将其剖分成两个单元，单元和节点编号如图所示，设节点位移和节点力分别为（匕，2）t和,j=l,2,3,已知平面梁单元单元刚度矩阵为1261-1261614/2-6/2广-12-6112-61612/2-614/2EI灯=皆其中，/为长度，E/为抗弯刚度，试求节点2和节点3的位移值，要求写出详细计算过程。图1217、解答:依题意可知，单元的位移列阵为/产=（匕节点力列阵为用=（居,氏,%）了,它们满足平衡方程：勺如包一二仍m单元（2）的位移列阵为M产;（孙内,与必产，节点力列阵为尸二仁,但，玛,外了,它们满足平衡方程：K-/产=尸产整个结构的位移列阵为，=（%©死鸟,将网了，力列阵为尸二（aM/22，FM）T，它们满足平衡方程:勺=（尸经计算可得，两个单元的刚度矩阵为126-126-126-12Kn)=£764-62,Kfy=El64-6-12-612-6-12-61262-6462662一64单元等效节点载