首页
发现
课程
培训
文章
案例
问答
需求
服务
行家
赛事
热门搜索
发布
消息
注册
|
登录
首页
/
文章
/
详情
弱形式到底“弱”在哪?一次讲清它与强形式的关系
兵心依旧
5天前
浏览2
关注
假设某连续体结构在给定荷载和边界条件下,其真实位移场为
u
。在连续介质力学中,
u
需要满足几何关系和材料本构关系,并使动量平衡方程在空间域内
逐点成立
。这组控制方程构成问题的
强形式
,并要求位移场具有足够的光滑性。
有限元计算的基础是“离散化”。然而,面对复杂的几何、非线性材料以及多变的边界条件,直接对要求如此苛刻的强形式进行离散几乎无法实现。因此,我们需要引入一种等价的全局积分表达——
“
弱形式
”
。
那么,弱形式究竟“弱”在哪里?它是否降低了物理要求?有限元误差又从哪一步开始产生?本文将围绕这些问题展开讨论。
PART 01
强形式:逐点平衡
在完全拉格朗日(
)格式下,强形式是一个完整的边界值问题,它要求物理守恒方程在初始构型
的每一个材料点上绝对成立。
数学表达
动量平衡(控制方程)
:
位移边界条件:
受力边界条件:
物理含义:域内任意一点的应力散度与体力之和,等于该点的惯性力。
“强”的代价:有限元的导数困境
“强”在哪里:
(1)由于应力
本身包含位移的一阶导数,计算散度
必然要求位移场具备二阶连续导数。(2)物体内部的
每一个材料点
,都必须满足动量平衡,差一丝一毫都不行。
有限元的痛点:
常规位移型有限元使用低阶形函数,仅能保证单元
间的基本连续(连续),根本无法算出二阶导数。因此,直接对强形式进行常规网格离散在数学上是行不通的。
物理逻辑的闭环
强形式绝非几条孤立的方程,而是一条单向传递的物理逻辑链:
即:未知位移场主导变形,变形通过材料特性激发出应力,应力最终代回动量方程以检验力学平衡。
图1 强形式的逐点平衡
PART 02
从强形式到弱形式:数学降阶与虚功原理
强形式要求动量平衡方程逐点成立,这对位移场的光滑性提出了较高要求。为降低这一要求,可以通过加权残值和分部积分,将强形式改写为弱形式。
引入测试函数:建立加权残值
将动量平衡方程写成残差形式:
强形式要求
在域内逐点成立。为了将其转化为积分形式,我们引入任意的测试函数
。在结构力学中,它具有明确的物理意义——
虚位移
。
作为允许的虚位移,它必须满足位移边界条件上的约束,即:
将残差与测试函数相乘并在全域积分:
核心逻辑:
这并非一次简单的“总体平均”。由于该积分等式必须对所有任意允许
的成立,它实际上等价于要求残差
在任何局部区域和任何方向上都必须为零。
分部积分:降低导数要求
直接积分并未解决应力散度(
)带来的高阶求导难题。关键的一步是对散度项应用分部积分(散度定理):
将边界积分拆分:在位移边界
上,因
使得积分为零;在受力边界
上,代入牵引力边界条件
。代回原式,整理得到动力学弱形式:
数学降阶:
分部积分后,体积分中不再出现
,从而降低了对位移场高阶导数的要求。这使得
未知位移场只需满足一阶导数存在即可,契合常规有限元的低阶形函数。
物理意义与边界条件
推导出的弱形式方程,在物理上正是经典的
虚功原理
:
等式左侧
:
内部虚功 + 惯性力虚功(静力问题无此项);
等式右侧
:
体力虚功 + 表面力虚功(外部虚功)。
此外,基于角动量平衡带来的柯西应力对称性,若定义虚应变
,则有
。这深刻表明:内部虚功仅由实际变形产生,刚体虚转动不作功。
经过转换,两类边界条件在弱形式中承担不同角色:
自然边界条件(受力边界):
受力边界条件以边界积分形式自动进入弱方程。
本质边界条件(位移边界):
位移条件无法通过积分化解,必须预先施加到允许的位移场和测试函数空间中,其中测试函数在该边界上必须为零。
因此,弱形式并没有改变平衡规律,而是通过积分表达和分部积分降低了对解的光滑性要求。
图2 强形式到弱形式的推导链
PART 03
弱形式“弱”在哪里?
“弱”仅仅指代对数学光滑性要求的降低,而绝非对物理定律的妥协。
弱化的是光滑性要求
强形式:
要求应力场能够被继续求导。位移场具有二阶或更高阶连续导数。
弱形式:
体积分只包含
和
。
位移和测试函数只需具备一阶弱导数(属于
空间),单元之间只保持位移连续的
形函数便可以用于离散。
图3 强形式到弱形式光滑性要求
弱化的不是物理规律
若强形式的解足够光滑,将它乘以任意测试函数并积分,必然得到弱形式。因此,强解一定也是弱解。
反之,根据变分学基本引理,若弱解具有足够的正则性,且弱形式对
所有
允许的测试函数均成立,则可以精确恢复域内的强形式平衡方程:
总结:对于足够光滑的解,强形式与弱形式均可以通过有限元法离散求解,且两者等价;但解不够光滑(无二阶连续导数)时,采用强形式无法用有限元离散求解,一般通过弱形式离散求解。
图4 强形式到弱形式关系及有限元离散
弱形式本身不是近似
强形式推导弱形式是没有损失任何信息的,两者反映物理规律一致。那为什么我们算出来的有限元结果总是存在误差呢?
答案:误差根本不产生于“强形式➡️弱形式”这一步,而是产生于“弱形式➡️离散化求解”这一步(伽辽金投影)。
在连续弱形式中,真实的位移
与测试函数均定义在无限维函数空间中。然而,计算机无法直接处理无限维问题。因此,有限元离散引入了两个关键的妥协:
解空间截断(降维)
:
强制将解空间和测试空间限制在由有限个网格节点与形函数
张成的
有限维子空间
中。离散解
和测试函数
只能表达为:
正交性弱化(误差诞生)
:
由于测试函数被降维至
,丧失了在无限维空间中的绝对“任意性”。此时,变分学基本引理不再全局成立——方程的残差仅能保证对有限组离散形函数保持正交,而无法再保证域内处处为零。
总结
:
有限元的本质是用有限维空间去逼近无限维空间。误差的诞生,正是由于我们将自由的无限维测试函数,局限在了有限维的形函数基底之中。
图5 有限元误差不来自于强形式与弱形式转换
来源:兵哥讲力学
登录后免费查看全文
立即登录
非线性
材料
控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-08-21
最近编辑:5天前
兵心依旧
博士
兵哥出品,必是精品
关注
获赞 81
粉丝 521
文章 51
课程 3
点赞
收藏
作者推荐
可试听
跟着兵哥学ABAQUS子程序二次开发(0基础入门)
¥100
5.0
立即查看
0/200
清空
提交
还没有评论
课程
培训
服务
行家
2024赛前研讨会(五):清华大学康进武老师《深度学习AI与数值模拟》报告会
线性强化弹塑性umat子程序教程
感应加热炉熔池技巧(考虑电磁搅拌力、浮升力、液冷等)
matlab矩形薄板有限元编程/小挠度弯曲/板壳力学/板单元/壳单元/四节点/刚度矩阵/
相关推荐
ABAQUS求解器从入门到精通
【2024年最新】Mechanical三维橡胶压缩变形非线性收敛调试仿真分析
Ansys workbench结构仿真分析
Abaqus材料模型-随动硬化kinematic hardening
最新文章
Qt嵌入式开发优选:QGraphicsView打造高性能、可定制的工业级图形界面
一年发射50多次、升至2万米高空:美军在波兰测试“微型高空气球”有何用意?
【历史上的今天】1836年8月21日:你天天解的 N–S 方程,名字里那位纳维,今天忌辰
【CAE】特斯拉Model Y整车碰撞仿真K文件碰撞模型及报告下载
GMW14116悬置设计及仿真标准
热门文章
通用/泛亚/蔚来/弗迪/上海电驱动等百位新能源专家确认发言,易贸三电活动6月苏州见
机械知识库:轴承市场趋势与未来发展
Compensators for Control of Switching Supplies
华为顶牛硬件女工程师:等我老了可以跟孩子说:“你娘当年,干过一票大的!”
锂电设备技术发展趋势!
其他人都在看
美太空军首次组建专职征兵中队:人才竞争开始向“太空作战准备”延伸
水下、太空、AI排在最前:白宫新科技战略释放美军技术投资新信号
从“有多少卫星”到“多快交付能力”:美太空军开始重新定义“太空战备”
中考高考落榜之后,普通家庭孩子该如何破局?
罗永浩x朱江明 丨 零跑凭什么成了销冠?
VIP会员
学习计划
福利任务
下载APP
联系我们
微信客服
联系客服
人工服务时间为周一至周五的9:30-19:30
非工作时间请在微信客服留言
客服热线:
4000-969-010
邮箱:
service@fangzhenxiu.com
地址:
北京市朝阳区莱锦创意园CN08座
帮助与反馈
返回顶部