2026年有限元期末试题库及答案

2026年有限元期末试题库及答案一、选择题（每题2分，共20分）1.有限元法的核心思想是将连续体离散为有限个单元的组合，其数学基础是（）。A.差分法B.变分原理与加权残值法C.边界元法D.有限体积法答案：B2.线性三角形单元的形函数N_i(x,y)在节点i处的值为（）。A.0B.1C.x_iD.y_i答案：B3.计算二维四边形单元的刚度矩阵时，若采用高斯积分，通常选择2×2积分点可精确积分（）次多项式。A.1B.2C.3D.4答案：C4.位移模式的构造需满足“完备性条件”，即必须包含（）。A.常数项和一次项B.二次项C.三角函数D.指数函数答案：A5.总刚度矩阵K的主要特性不包括（）。A.对称性B.奇异性C.稀疏性D.非对角占优答案：D6.对于本质边界条件（位移边界条件），常用的处理方法是（）。A.直接代入法B.罚函数法C.拉格朗日乘子法D.以上均是答案：D7.平面应力问题与平面应变问题的弹性矩阵D的区别在于（）。A.平面应力的D中μ替换为μ/(1-μ)B.平面应变的D中E替换为E/(1-μ²)C.平面应力的D中ν替换为ν/(1-ν)D.平面应变的D中E替换为E/(1-ν²)，ν替换为ν/(1-ν)答案：D8.空间梁单元的每个节点通常具有（）个自由度。A.3B.4C.5D.6答案：D9.有限元计算中的误差主要来源于（）。A.离散误差与数值积分误差B.网格畸变误差C.材料参数误差D.以上均是答案：D10.非协调单元是指（）。A.形函数不满足连续性要求B.单元刚度矩阵不对称C.无法通过分片试验D.节点自由度不匹配答案：A二、填空题（每题2分，共20分）1.形函数N_i的基本性质之一是：在所有节点上，ΣN_i=______。答案：12.单元刚度矩阵k的表达式为k=∫_ΩB^TDBdΩ，其中B是______矩阵，D是______矩阵。答案：应变-位移；弹性3.高斯积分的精度与积分点数量有关，1点高斯积分可精确积分______次多项式。答案：14.处理自然边界条件（力边界条件）时，通常通过______将面力或集中力等效到节点荷载向量中。答案：虚功原理5.平面三角形单元（3节点）每个节点有______个自由度，单元总自由度为______。答案：2；66.总刚度矩阵的稀疏性是由于______。答案：节点仅与相邻单元的节点相关，非相邻节点的刚度系数为07.协调单元要求形函数在单元交界面上满足______连续性。答案：C0（位移连续）8.稳态热传导问题的控制方程为∇·(k∇T)+Q=0，其中k是______，Q是______。答案：热导率；内热源密度9.有限元误差估计的主要方法有______和______（任写两种）。答案：残差法；能量范数法10.杆单元（2节点）的应变ε与节点位移u1、u2的关系为ε=______（设杆长为l）。答案：(u2-u1)/l三、简答题（每题6分，共30分）1.简述有限元法的主要分析步骤。答案：（1）连续体离散化：将求解区域划分为有限个单元，确定节点坐标和单元连接关系；（2）选择形函数：假设单元内场变量（如位移）的插值函数，通常为多项式；（3）建立单元方程：通过变分原理或加权残值法，推导单元刚度矩阵和节点荷载向量；（4）组装总刚度矩阵：将所有单元的刚度矩阵和荷载向量按节点编号组装为整体方程组KU=F；（5）处理边界条件：施加位移边界条件，消除总刚度矩阵的奇异性；（6）求解线性方程组：得到节点未知量（如位移）；（7）后处理：计算单元应力、应变等派生量，分析结果合理性。2.形函数需要满足哪些基本性质？答案：（1）节点值特性：N_i(x_j,y_j)=δ_ij（δ为克罗内克函数，i=j时为1，否则为0）；（2）求和为1：ΣN_i=1（保证刚体位移模式的正确性）；（3）完备性：形函数应包含足够低次的多项式（至少常数项和一次项），以反映单元的刚体位移和常应变状态；（4）协调性（对协调单元）：相邻单元的形函数在公共边界上连续，避免位移间断。3.比较平面应力问题与平面应变问题的异同。答案：相同点：均为二维问题，应力或应变分量仅与x、y坐标有关，且σ_xy=τ_xy，σ_xz=σ_yz=0（平面应力）或ε_xz=ε_yz=0（平面应变）。不同点：（1）平面应力：厚度方向应力σ_z≈0，应变ε_z≠0；适用于薄板（厚度远小于平面尺寸）。（2）平面应变：厚度方向应变为0（ε_z=0），应力σ_z≠0；适用于长柱体（长度远大于截面尺寸）。（3）弹性矩阵D不同：平面应力的D中E→E，ν→ν；平面应变的D中E→E/(1-ν²)，ν→ν/(1-ν)。4.总刚度矩阵为何具有稀疏性和带状性？答案：稀疏性：每个节点仅与相邻单元的节点产生刚度耦合，非相邻节点的刚度系数为0，因此总刚度矩阵中大部分元素为0，呈现稀疏性。带状性：若节点编号按“先局部后整体”的顺序（如带宽最小化编号），则非零元素集中在主对角线附近的带状区域内，带宽由单元连接的最大节点编号差决定。5.说明高斯积分在单元刚度矩阵计算中的作用。答案：单元刚度矩阵的计算涉及对B^TDB的体积分，直接积分复杂。高斯积分通过选择特定积分点（高斯点）和权重系数，将积分转化为离散点的函数值加权和，可高效计算多项式积分。例如，2×2高斯积分可精确计算3次多项式的积分，减少计算量；高阶高斯积分可提高精度，平衡计算效率与结果准确性。四、计算题（每题10分，共30分）1.推导两节点杆单元的刚度矩阵。已知杆长l，弹性模量E，截面积A，假设位移模式为u(x)=N1u1+N2u2，其中N1=1-x/l，N2=x/l。答案：（1）应变ε=du/dx=(-1/l)u1+(1/l)u2=Bu^e，其中B=[-1/l,1/l]；（2）单元刚度矩阵k=∫_0^lB^TEABdx；（3）代入B得：k=EA/l∫_0^l[1-1;-11]dx=(EA/l)[1-1;-11]。物理意义：k_11表示节点1产生单位位移（u1=1，u2=0）时，节点1所需的力；k_12表示节点2产生单位位移（u2=1，u1=0）时，节点1所需的力（负号表示反向）。2.构造三节点三角形单元（节点1(0,0)、节点2(2,0)、节点3(0,2)）的形函数N1、N2、N3，并验证其性质。答案：（1）单元面积A=1/2|(x2-x1)(y3-y1)-(x3-x1)(y2-y1)|=1/2|2×2-0×0|=2；（2）形函数表达式：N1=(a1+b1x+c1y)/(2A)，其中a1=x2y3-x3y2=2×2-0×0=4，b1=y2-y3=0-2=-2，c1=x3-x2=0-2=-2，故N1=(4-2x-2y)/4=1-x/2-y/2；同理，N2=(a2+b2x+c2y)/(2A)，a2=x3y1-x1y3=0×0-0×2=0，b2=y3-y1=2-0=2，c2=x1-x3=0-0=0，故N2=(0+2x+0y)/4=x/2；N3=(a3+b3x+c3y)/(2A)，a3=x1y2-x2y1=0×0-2×0=0，b3=y1-y2=0-0=0，c3=x2-x1=2-0=2，故N3=(0+0x+2y)/4=y/2；（3）验证性质：节点值：N1(0,0)=1，N1(2,0)=1-1-0=0，N1(0,2)=1-0-1=0（符合δ_ij）；求和：N1+N2+N3=1-x/2-y/2+x/2+y/2=1（满足ΣN_i=1）；完备性：包含常数项（1）和一次项（x/2、y/2），可反映刚体位移和常应变。3.已知某结构总刚度矩阵K和荷载向量F如下，节点1的位移u1=0，节点3的位移u3=u3（未知），节点2的位移u2未知。试处理边界条件并求解u2和u3（假设K=[[4,-2,0],[-2,5,-3],[0,-3,6]]，F=[0,5,10]^T）。答案：（1）原始方程：KU=F→4u1-2u2+0u3=02u1+5u2-3u3=50u1-3u2+6u3=10（2）代入u1=0，方程简化为：2u2=0→矛盾（说明需采用直接代入法修改K和F）；正确处理：将已知位移u1=0对应的行和列修改为K(1,1)=1，其余K(1,j)=K(j,1)=0，F(1)=0；修改后K=[[1,0,0],[-2,5,-3],[0,-3,6]]，F=[0,5,10]^T；（3）方程变为：1×u1+0×u2+0×u3=0→u1=0（满足）2u1+5u2-3u3=5→5u2-3u3=50u1-3u2+6u3=10→-3u2+6u3=10（4）联立后两式：由第二式：5u2=5+3u3→u2=1+0.6u3；代入第三式：-3(1+0.6u3)+6u3=10→-3-1.8u3+6u3=10→4.2u3=13→u3=13/4.2≈3.095；则u2=1+0.6×3.095≈2.857；（5）支反力计算（取原始方程第一式）：F1=4u1-2u2+0u3=0-2×2.857≈-5.714（负号表示方向与假设相反）。五、综合题（20分）某简支梁跨度L=4m，承受均布荷载q=10kN/m，材料弹性模量E=2×10^4MPa，截面惯性矩I=1×10^-4m^4。要求：（1）将梁离散为2个等长单元（单元长度l=2m），选择欧拉-伯努利梁单元（每个节点2个自由度：w（挠度）和θ（转角））；（2）推导单个梁单元的刚度矩阵（已知形函数对应的B矩阵为[-6/l²4/l6/l²2/l]，其中l为单元长度）；（3）组装总刚度矩阵K，施加边界条件（左端点w1=0，θ1=0；右端点w3=0）；（4）求解节点2的挠度w2和转角θ2；（5）讨论网格加密（如4个单元）对结果精度的影响。答案：（1）离散化：节点1(0)、节点2(2m)、节点3(4m)，单元1（1-2）和单元2（2-3），l=2m。（2）梁单元刚度矩阵推导：欧拉-伯努利梁的应变能U=1/2∫EI(w'')²dx，位移模式w(x)=N1w1+N2θ1l+N3w2+N4θ2l（N为三次形函数）；曲率κ=w''=Bu^e，其中B=[-6/l²4/l6/l²2/l]（对应节点自由度u^e=[w1,θ1,w2,θ2]^T）；单元刚度矩阵k=∫_0^lEIB^TBdx=EI∫_0^l[36/l^4-24/l^3-36/l^412/l^3;24/l^316/l^224/l^3-8/l^2;36/l^424/l^336/l^4-12/l^3;12/l^3-8/l^2-12/l^34/l^2]dx=(EI/l^3)[12-6l-126l;6l4l²6l-2l²;126l12-6l;6l-2l²-6l2l²]。（3）组装总刚度矩阵：单元1的节点1-2对应总自由度[1(w1),2(θ1),3(w2),4(θ2)]；单元2的节点2-3对应总自由度[3(w2),4(θ2),5(w3),6(θ3)]；总刚度矩阵K为6×6矩阵，非零块为：K(1,1)=12EI/l³（单元1的k11），K(1,2)=-6lEI/l³（k12），K(1,3)=-12EI/l³（k13），K(1,4)=6lEI/l³（k14）；K(2,2)=4l²EI/l³（k22），K(2,3)=6lEI/l³（k23），K(2,4)=-2l²EI/l³（k24）；K(3,3)=12EI/l³（单元1的k33）+12EI/l³（单元2的k11）=24EI/l³；K(3,4)=-6lEI/l³（单元1的k34）+6lEI/l³（单元2的k12）=0；K(3,5)=-12EI/l³（单元2的k13）；K(3,6)=6lEI/l³（单元2的k14）；同理填充其他项（过程略）。（4）边界条件：w1=0（自由度1），θ1=0（自由度2），w3=0（自由度5）。将对应行/列置零，主对角线置1，荷载向量F考虑均布荷载等效节点力（q=10kN/m，单元等效荷载为[ql/2,ql²/12,ql/2,-ql²/12]^T）。总荷载向量F=[0,0,ql/2+ql/2,ql²/12-