布西涅斯克问题（半无限弹性体的点荷载） — 故障排除
更丰富的内容请在 boussinesq-halfspace.html 查看。
模型边界的影响
如果模型的外形太小会怎样？
固定边界的反射会导致应力被过大评估。如果底面固定太浅，荷载点正下方的$\sigma_{zz}$可能会比理论值大20%以上。症状表现为模型尺寸变化时结果也随之变化，此时应怀疑这种影响。
对策：
- 逐步增大模型尺寸，采用结果不再变化的尺寸
- 用无限单元近似远场条件
- 经验值为 $R_{model} \geq 20 z_{eval}$
无限单元效果不好怎么办？
无限单元的衰减函数可能与问题特性不匹配。Abaqus的无限单元假设指数衰减，但布西涅斯克解是$1/r$衰减，在严格意义上并不完全适配。尽管如此，实际应用中通常能获得足够的精度。如果对残差有顾虑，更安全的做法是扩大有限单元区域（远离无限单元边界）。
荷载施加的问题
在单个节点施加集中荷载时，该节点的应力会变得无意义地很大，对吗？
这是正常行为。荷载点的应力随着网格细化而无限增大。这反映了物理上正确的奇异性，不是有限元分析的缺陷。
实际对策：
1. 不计算荷载点处的应力。在 $z > 5h_{elem}$ 处计算
2. 如果将荷载分散为面压分布，半径$a$如果小于计算距离的1/10，则与理论值的差异小于1%
3. 对于Hertz接触等物理上确定面压分布的情况，直接给定该分布，与布西涅斯克解的积分结果比较
分散荷载的具体方法是什么？
通过验证反力总和等于施加荷载来确认等价性。
数值伪影
结果中会出现奇怪的波纹吗？
用低阶单元（TET4、HEX8）绘制应力等值线时，会在单元之间出现不连续性，形成锯齿状外观。这是因为C0连续的单元在单元边界处的应力（应变的导数）会出现间断。
对策：
- 改用二阶单元（TET10、HEX20）可大幅改善
- 用节点平均化方法（SPR法：Zienkiewicz-Zhu超收敛补丁恢复）进行平滑化
- 平滑化前后差异较大的区域表示网格不足
SPR法是什么？第一次听说。
这是Zienkiewicz和Zhu（1992）提出的事后误差估计方法。在每个节点的补丁内，用最小二乘法从积分点应力恢复高精度节点应力。恢复前后的差值是局部离散化误差的估计值。许多商用求解器都将其作为选项实现。Abaqus中可通过*ERROR INDICATOR使用。
验证数据的可视化
理论值与计算值的定量对比。合格标准为误差5%以内。
| 评估项目 | 理论值/参考值 | 计算值 | 相对误差 [%] | 判定 |
|---|---|---|---|---|
| 最大位移 | 1.000 | 0.998 | 0.20 | 合格 |
| 最大应力 | 1.000 | 1.015 | 1.50 | 合格 |
| 固有振动数（1阶） | 1.000 | 0.997 | 0.30 | 合格 |
| 反力合计 | 1.000 | 1.001 | 0.10 | 合格 |
| 能量守恒 | 1.000 | 0.999 | 0.10 | 合格 |
判定标准：相对误差 < 1%：■ 优秀，1～5%：■ 可接受，> 5%：■ 需检查
布西涅斯克问题（半无限弹性体的点荷载） — 故障排除的CAE实务质量检查
布西涅斯克问题（半无限弹性体的点荷载） — 故障排除不是单独的公式，而应作为验证和有效性确认中的工程模型来处理。要获得可信的结果，需要将支配物理、材料值、边界条件、离散化、求解器设置和后处理标准作为一个整体来描述。在用于设计判断前，请明确哪些量是输入、哪些是计算结果、哪些是诊断指标。
建模检查清单
- 用途明确化：确定布西涅斯克问题（半无限弹性体的点荷载） — 故障排除是用于概算、详细设计、故障调查还是其他分析的验证。
- 单位统一：内部计算采用SI单位，记录荷载、形状、材料常数、时间和频率尺度的换算。
- 假设明文化：确认线性性、定常/非定常、小变形、连续体假设、对称条件、理想边界条件的适用范围。
- 与基准解比对：手算、极限情况、网格收敛或独立求解器结果对比后再采用。
验证中应观察的信号
| 确认项目 | 应查看的内容 | 应警惕的征兆 |
|---|---|---|
| 输入条件 | 形状、材料、荷载、约束是否与验证和有效性确认问题相符。 | 图形看起来合理，但数量级或单位不符。 |
| 数值设置 | 网格、时间步、收敛容差、求解器设置是否对布西涅斯克半空间问题足够。 | 轻微改变设置就导致结果大幅变化。 |
| 物理适用范围 | 使用的理论在应力、温度、速度、频率范围内是否有效。 | 将结果外推到超出模型假设的条件。 |
实务中，应将输入表、模型文件、结果图、审查意见用同一单位保存。这样可追踪布西涅斯克问题（半无限弹性体的点荷载） — 故障排除的计算依据，避免将此页面当作黑盒答案的风险。
相关主题
帮助
更详细
错误