(2025年)高等教育有限元分析复习题(附答案)

(2025年)高等教育有限元分析复习题(附答案)一、选择题（每题3分，共18分）1.有限元法的核心思想是将连续体离散为有限个单元的集合，其数学本质是（）。A.加权余量法的点匹配形式B.变分原理与分片插值的结合C.差分法的高阶近似D.边界元法的区域扩展答案：B2.平面应力问题中，三角形三节点单元的应变矩阵B是（）阶矩阵。A.3×3B.3×6C.6×3D.6×6答案：B（每个节点2个位移自由度，3个节点共6个自由度；平面应力问题应变有3个分量，故B为3×6）3.等参元的“等参”指（）。A.单元几何形状与位移插值使用相同阶数的形函数B.单元节点数与自由度参数相等C.单元材料参数在积分中保持恒定D.单元应变与应力呈线性关系答案：A4.采用高斯积分计算单元刚度矩阵时，若形函数为二次多项式，通常需要（）积分点。A.1点B.2点C.3点D.4点答案：B（二次多项式在1维需2点高斯积分，2维需2×2=4点，但实际计算中常取2点以平衡精度与效率）5.对于位移边界条件u=0的处理，以下方法中误差最小的是（）。A.直接代入法（主元法）B.罚函数法C.拉格朗日乘子法D.等效节点力法答案：A（直接修改刚度矩阵主对角线元素，严格满足边界条件）6.有限元解的收敛性要求形函数必须满足（）。A.仅完备性B.仅协调性C.完备性与协调性D.连续性与周期性答案：C（完备性保证常应变/常应力状态，协调性保证单元间位移连续）二、填空题（每空2分，共24分）1.加权余量法中，伽辽金法的权函数取为__________。答案：试探函数（或形函数）2.平面应变问题中，应变分量包括εₓ、εᵧ和__________。答案：γₓᵧ（剪切应变）3.三维实体单元每个节点的自由度为__________。答案：3（x、y、z方向位移）4.单元刚度矩阵的物理意义是：当单元某节点产生单位位移时，在__________节点上需要施加的力。答案：所有5.等参元雅可比矩阵J的行列式|J|需满足__________，否则单元会发生畸变。答案：|J|>0（保证坐标变换可逆）6.应力恢复的超收敛Patch法通过__________拟合单元边界应力，提高整体精度。答案：最小二乘法7.集中质量矩阵与一致质量矩阵的主要区别在于__________。答案：集中质量矩阵假设质量集中于节点，一致质量矩阵通过形函数积分得到8.热应力分析中，温度载荷转化为等效节点力时需考虑__________和热膨胀系数。答案：弹性模量与温度变化量三、简答题（每题8分，共40分）1.简述形函数的基本性质。答：形函数Nᵢ需满足：（1）在节点i处Nᵢ=1，其他节点j≠i时Nᵢ=0；（2）所有形函数之和为1（∑Nᵢ=1），保证刚体位移模式；（3）形函数是坐标的连续函数，且具有足够阶数的可微性（至少一阶），以保证应变矩阵B的连续性；（4）对于高阶单元，形函数需满足协调性条件，即相邻单元公共边上的位移插值函数一致。2.说明单元刚度矩阵的组装过程及物理意义。答：组装过程：（1）确定整体节点与单元节点的对应关系（节点编号映射）；（2）将单元刚度矩阵kᵉ的元素kᵉ_ij按整体节点自由度编号i’、j’，叠加到整体刚度矩阵K的K_i’j’位置（即“对号入座，同号相加”）。物理意义：整体刚度矩阵K描述了结构所有节点位移与节点力的线性关系（F=Kδ），其元素K_ij表示节点j产生单位位移时，节点i所需施加的力，体现了结构的整体刚度分布。3.数值积分在有限元计算中的作用是什么？高斯积分为何被广泛使用？答：作用：有限元中单元刚度矩阵、质量矩阵等需通过积分计算（如kᵉ=∫BᵀDBdV），直接解析积分复杂或不可行，数值积分通过离散积分点近似求解。高斯积分优势：（1）精度高，n点高斯积分可精确积分2n-1次多项式；（2）计算效率高，相同精度下所需积分点少于牛顿-科特斯公式；（3）适用于等参元，通过自然坐标积分避免复杂的几何变换积分。4.对比位移边界条件（本质边界条件）与力边界条件（自然边界条件）的区别及处理方法。答：区别：位移边界条件直接约束节点位移（如u=u₀），是强制满足的几何条件；力边界条件描述边界上的外力（如∫TᵀNdS=F），通过虚功原理自然包含在方程中。处理方法：位移边界条件常用直接代入法（修改刚度矩阵主元）、罚函数法（引入大刚度项）或拉格朗日乘子法（增广方程）；力边界条件通过等效节点力公式（Fᵉ=∫TᵀNdS）转化为节点力向量，直接叠加到整体载荷向量中。5.有限元解的收敛性需要满足哪些条件？工程中如何验证收敛性？答：必要条件：（1）完备性：形函数包含常数项（刚体位移）和线性项（常应变）；（2）协调性（C⁰连续性）：相邻单元公共边界上的位移连续。充分条件：当单元尺寸h→0时，解趋近于精确解（h收敛）或单元阶数p→∞时解收敛（p收敛）。工程验证方法：（1）h收敛性检验：逐步细化网格，观察关键节点位移或应力的变化，当变化量小于阈值时认为收敛；（2）p收敛性检验：提高单元阶数（如从线性到二次），对比结果；（3）能量范数误差估计：计算有限元解与精确解的能量差，判断是否满足收敛准则。四、计算题（共58分）1.（15分）推导平面杆单元（长度L，截面积A，弹性模量E）在局部坐标系下的刚度矩阵。解：平面杆单元仅轴向（x’轴）有位移，节点1位移u₁，节点2位移u₂。形函数N₁=1-ξ，N₂=ξ（ξ=x’/L，0≤ξ≤1）。应变ε=du/dx’=d(N₁u₁+N₂u₂)/dx’=(-u₁+u₂)/L应力σ=Eε=E(u₂-u₁)/L单元内力F=σA=EA(u₂-u₁)/L根据虚功原理，δW=δu₁F₁+δu₂F₂=∫δεσAdx’δε=d(δu)/dx’=(-δu₁+δu₂)/L代入得：δu₁F₁+δu₂F₂=∫₀ᴸ[(-δu₁+δu₂)/L][E(u₂-u₁)/L]Adx’=(EA/L²)(-δu₁+δu₂)(-u₁+u₂)L=(EA/L)(δu₁u₁δu₁u₂δu₂u₁+δu₂u₂)比较系数得单元刚度方程：[F₁;F₂]=(EA/L)[1-1;-11][u₁;u₂]故局部坐标系下刚度矩阵kᵉ=(EA/L)[[1,-1],[-1,1]]2.（18分）已知三角形三节点单元节点坐标为(0,0)、(2,0)、(0,1)，求其形函数表达式及应变矩阵B。解：三角形单元节点i(0,0)、j(2,0)、k(0,1)，面积A=(1/2)|(xj-xi)(yk-yi)-(xk-xi)(yj-yi)|=(1/2)(2×1-0×0)=1。形函数公式：Nᵢ=(aᵢ+bᵢx+cᵢy)/(2A)，其中：aᵢ=xjykxkyj=2×10×0=2bᵢ=yjyk=01=-1cᵢ=xkxj=02=-2同理，aⱼ=xkyixiyi=0×00×0=0，bⱼ=ykyi=10=1，cⱼ=xixk=00=0aₖ=xiyjxjyi=0×02×0=0，bₖ=yiyj=00=0，cₖ=xjxi=20=2因此：Nᵢ=(2x2y)/(2×1)=10.5xyNⱼ=(0+x+0y)/2=0.5xNₖ=(0+0x+2y)/2=y应变矩阵B=[BᵢBⱼBₖ]，其中Bᵣ=[∂Nᵣ/∂x0;0∂Nᵣ/∂y;∂Nᵣ/∂y∂Nᵣ/∂x]（平面应力问题）计算偏导：∂Nᵢ/∂x=-0.5，∂Nᵢ/∂y=-1∂Nⱼ/∂x=0.5，∂Nⱼ/∂y=0∂Nₖ/∂x=0，∂Nₖ/∂y=1故Bᵢ=[-0.50;0-1;-1-0.5]，Bⱼ=[0.50;00;00.5]，Bₖ=[00;01;10]整体B矩阵为3×6矩阵：B=[[-0.5,0,0.5,0,0,0],[0,-1,0,0,0,1],[-1,-0.5,0,0.5,1,0]]3.（25分）四节点矩形单元（长2a，宽2b，弹性模量E，泊松比ν，平面应力状态），采用2×2高斯积分计算其刚度矩阵。解：单元在自然坐标(ξ,η)下，ξ∈[-1,1]，η∈[-1,1]，节点坐标：1(-a,-b)、2(a,-b)、3(a,b)、4(-a,b)。形函数Nᵢ=(1/4)(1+ξᵢξ)(1+ηᵢη)（i=1,2,3,4，ξ₁=-1,ξ₂=1,ξ₃=1,ξ₄=-1；η₁=-1,η₂=-1,η₃=1,η₄=1）。雅可比矩阵J=∂(x,y)/∂(ξ,η)=[∑Nᵢ,ξxᵢ∑Nᵢ,ξyᵢ;∑Nᵢ,ηxᵢ∑Nᵢ,ηyᵢ]计算Nᵢ,ξ=(1/4)ξᵢ(1+ηᵢη)，Nᵢ,η=(1/4)ηᵢ(1+ξᵢξ)代入节点坐标xᵢ=(-a,a,a,-a)，yᵢ=(-b,-b,b,b)，得：J=[(1/4)(-1)(1+η)(-a)+(1/4)(1)(1+η)(a)+(1/4)(1)(1-η)(a)+(1/4)(-1)(1-η)(-a),(1/4)(-1)(1+η)(-b)+(1/4)(1)(1+η)(-b)+(1/4)(1)(1-η)(b)+(1/4)(-1)(1-η)(b);(1/4)(-1)(1+ξ)(-b)+(1/4)(1)(1+ξ)(-b)+(1/4)(1)(1-ξ)(b)+(1/4)(-1)(1-ξ)(b),(1/4)(-1)(1+ξ)(-b)+(1/4)(1)(1+ξ)(-b)+(1/4)(1)(1-ξ)(b)+(1/4)(-1)(1-ξ)(b)]化简后J=[[a,0],[0,b]]，故|J|=ab，J⁻¹=[[1/a,0],[0,1/b]]应变矩阵B=LN，其中L为微分算子矩阵（平面应力：L=[∂/∂x,0;0,∂/∂y;∂/∂y,∂/∂x]）。N为形函数矩阵（4节点×2自由度，N=[N₁0N₂0N₃0N₄0;0N₁0N₂0N₃0N₄]）B=LN=[∂N₁/∂x0∂N₂/∂x0∂N₃/∂x0∂N₄/∂x0;0∂N₁/∂y0∂N₂/∂y0∂N₃/∂y0∂N₄/∂y;∂N₁/∂y∂N₁/∂x∂N₂/∂y∂N₂/∂x∂N₃/∂y∂N₃/∂x∂N₄/∂y∂N₄/∂x]由∂Nᵢ/∂x=∂Nᵢ/∂ξ·∂ξ/∂x=(Nᵢ