力学有限元的基石:虚功原理的推导

推导虚功方程的过程

弹性力学的平衡方程

在张量形式中,平衡方程为:

∇ ⋅ σ + b = 0 \nabla \cdot \sigma + b = 0 σ+b=0

用下标表示为:

∂ σ i j ∂ x j + b i = 0 \frac{\partial \sigma_{ij}}{\partial x_j} + b_i = 0 xjσij+bi=0

其中, σ i j \sigma_{ij} σij 是应力张量的分量, b i b_i bi 是体积力的分量。

平衡方程弱形式

我们乘以一个任意的虚位移场 v i v_i vi 并在整个体积 Ω \Omega Ω 上积分:

∫ Ω ( ∂ σ i j ∂ x j + b i ) v i d Ω = 0 \int_\Omega \left( \frac{\partial \sigma_{ij}}{\partial x_j} + b_i \right) v_i \, d\Omega = 0 Ω(xjσij+bi)vidΩ=0

展开得到:

∫ Ω ∂ σ i j ∂ x j v i d Ω + ∫ Ω b i v i d Ω = 0 \int_\Omega \frac{\partial \sigma_{ij}}{\partial x_j} v_i \, d\Omega + \int_\Omega b_i v_i \, d\Omega = 0 ΩxjσijvidΩ+ΩbividΩ=0

对第一个积分项使用分部积分

对第一个积分项 ∫ Ω ∂ σ i j ∂ x j v i d Ω \int_\Omega \frac{\partial \sigma_{ij}}{\partial x_j} v_i \, d\Omega ΩxjσijvidΩ 使用分部积分和微分乘法法则,将导数从应力张量 σ i j \sigma_{ij} σij 转移到虚位移场 v i v_i vi 上。

∫ Ω ( v i ∂ σ i j ∂ x j + σ i j ∂ v i ∂ x j ) d Ω = ∫ Ω ∂ ( σ i j v i ) ∂ x j d Ω \int_\Omega \left( v_i \frac{\partial \sigma_{ij}}{\partial x_j} + \sigma_{ij} \frac{\partial v_i}{\partial x_j} \right) d\Omega = \int_\Omega \frac{\partial (\sigma_{ij} v_i)}{\partial x_j} \, d\Omega Ω(vixjσij+σijxjvi)dΩ=Ωxj(σijvi)dΩ

根据高斯散度定理有:

∫ Ω ∂ ( σ i j v i ) ∂ x j d Ω = ∫ ∂ Ω σ i j v i n j d S \int_\Omega \frac{\partial (\sigma_{ij} v_i)}{\partial x_j} \, d\Omega = \int_{\partial\Omega} \sigma_{ij} v_i n_j \, dS Ωxj(σijvi)dΩ=ΩσijvinjdS

这里, ∂ Ω \partial\Omega Ω 是体积 Ω \Omega Ω 的边界, n j n_j nj 是边界上的单位外法向量。因此,我们可以将第一个积分项写成:

∫ Ω v i ∂ σ i j ∂ x j d Ω = ∫ ∂ Ω σ i j v i n j d S − ∫ Ω σ i j ∂ v i ∂ x j d Ω \int_\Omega v_i \frac{\partial \sigma_{ij}}{\partial x_j} \, d\Omega = \int_{\partial\Omega} \sigma_{ij} v_i n_j \, dS - \int_\Omega \sigma_{ij} \frac{\partial v_i}{\partial x_j} \, d\Omega ΩvixjσijdΩ=ΩσijvinjdSΩσijxjvidΩ

代入平衡方程弱形式

将上面分部积分得到的公式带入到平衡方程弱形式可得:

∫ Ω σ i j ∂ v i ∂ x j d Ω = ∫ ∂ Ω σ i j v i n j d S + ∫ Ω b i v i d Ω \int_\Omega \sigma_{ij} \frac{\partial v_i}{\partial x_j} \, d\Omega = \int_{\partial\Omega} \sigma_{ij} v_i n_j \, dS + \int_\Omega b_i v_i \, d\Omega ΩσijxjvidΩ=ΩσijvinjdS+ΩbividΩ

根据边界条件:

  • 在边界 ∂ Ω t \partial\Omega_t Ωt 上有应力边界条件 σ i j n j = t i \sigma_{ij} n_j = t_i σijnj=ti
  • 在边界 ∂ Ω u \partial\Omega_u Ωu 上有位移边界条件 u i = u 0 i u_i = u_{0i} ui=u0i

最终虚功原理(虚位移原理)

结合体积力项和边界条件,得到弱形式:

∫ Ω σ i j ∂ v i ∂ x j d Ω = ∫ ∂ Ω t t i v i d S + ∫ Ω b i v i d Ω \int_\Omega \sigma_{ij} \frac{\partial v_i}{\partial x_j} \, d\Omega = \int_{\partial\Omega_t} t_i v_i \, dS + \int_\Omega b_i v_i \, d\Omega ΩσijxjvidΩ=ΩttividS+ΩbividΩ

由应力张量的对称性

σ i j ∂ v i ∂ x j = 1 2 ( σ i j + σ j i ) ∂ v i ∂ x j = σ i j 1 2 ( ∂ v i ∂ x j + ∂ v j ∂ x i ) \sigma_{ij} \frac{\partial v_i}{\partial x_j} = \frac{1}{2}(\sigma_{ij} + \sigma_{ji})\frac{\partial v_i}{\partial x_j} = \sigma_{ij} \frac{1}{2}(\frac{\partial v_i}{\partial x_j} + \frac{\partial v_j}{\partial x_i}) σijxjvi=21(σij+σji)xjvi=σij21(xjvi+xivj)

最终,平衡方程的弱形式为:

∫ Ω σ i j ϵ i j d Ω − ∫ Ω b i v i d Ω − ∫ ∂ Ω t t i v i d S = 0 \int_\Omega \sigma_{ij} \epsilon_{ij} \, d\Omega - \int_\Omega b_i v_i \, d\Omega - \int_{\partial\Omega_t} t_i v_i \, dS = 0 ΩσijϵijdΩΩbividΩΩttividS=0

如果将上面的 v i v_i vi写成位移的变分形式可以得到如下的虚功方程

∫ Ω σ : δ ϵ d Ω = ∫ Ω b ⋅ δ u d Ω + ∫ ∂ Ω t t ⋅ δ u d S \int_\Omega \sigma : \delta \epsilon \, d\Omega = \int_\Omega b \cdot \delta u \, d\Omega +\int_{\partial\Omega_t} t \cdot \delta u \, dS Ωσ:δϵdΩ=ΩbδudΩ+ΩttδudS

总结

虚位移原理是平衡方程和力的边界条件的等效积分“弱”形式。虚位移原理的力学意义是:如果力系(包括内力和外力)是平衡的,则它们在虚位移和虚应变上所作之功的总和为零。反之,如果力系在虚位移及虚应变上所作之功的和等于零,则它们一定是满足平衡的,所以虚位移原理表述了力系平衡的必要而充分的条件。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/diannao/43538.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

知识图谱入门笔记

自学参考: 视频:斯坦福CS520 | 知识图谱 最全知识图谱综述 详解知识图谱的构建全流程 知识图谱构建(概念,工具,实例调研) 一、基本概念 知识图谱(Knowledge graph):由结…

Redis管理禁用命令

在redis数据量比较大时,执行 keys * ,fluashdb 这些命令,会导致redis长时间阻塞,大量请求被阻塞,cpu飙升,严重可能导致redis宕机,数据库雪崩。所以一些命令在生产环境禁止使用。 Redis 禁用命令…

【C语言】指针(3):探索-不同类型指针变量

目录 一、字符指针变量 二、数组指针变量 三、二维数组传参的本质 四、函数指针变量 4.1 函数指针变量 4.2 函数指针变量的使用 4.3 函数指针变量的拓展 五、函数指针数组 六、转移表的应用 通过深入理解指针(1)和深入理解指针(2&am…

67.SAP FICO-凭证类型学习

目录 SAP凭证类型 凭证类型的作用 - OBA7 SAP默认的凭证类型更改 FI相应事务代码默认凭证类型 - OBU1 对FB50、60、70默认凭证类型的更改 - OBZO 后勤货物移动默认凭证类型 - OMBA 发货凭证类型 收货凭证类型 自动移动凭证类型 存货盘点凭证类型 发票默认的凭证类…

深度学习Day-24:ResNeXt-50算法思考

🍨 本文为:[🔗365天深度学习训练营] 中的学习记录博客 🍖 原作者:[K同学啊 | 接辅导、项目定制] 要求: 阅读给出代码,判断是否存在错误,正确与否都请给出你的思考;查找…

如何减少开发过程中的bug-数据库篇

1.1慢查询 1.1.1 是否命中索引 提起慢查询,我们马上就会想到加索引。如果一条SQL没加索引,或者没有命中索引的话,就会产生慢查询。 索引哪些情况会失效? 查询条件包含or,可能导致索引失效 如果字段类型是字符串&am…

LeetCode 0724.寻找数组的中心下标:前缀和(时空复杂度O(n)+O(1))

title: 724.寻找数组的中心下标 date: 2024-07-08 13:22:58 tags: [题解, LeetCode, 简单, 数组, 前缀和] 【LetMeFly】724.寻找数组的中心下标:前缀和(时空复杂度O(n)O(1)) 力扣题目链接:https://leetcode.cn/problems/find-pi…

数据结构--二叉树相关习题5(判断二叉树是否是完全二叉树 )

1.判断二叉树是否是完全二叉树 辨别: 不能使用递归或者算节点个数和高度来判断。 满二叉树可以用高度和节点来判断,因为是完整的。 但是完全二叉树前面是满的,但是最后一层是从左到右连续这种 如果仍然用这种方法的话,如下图…

暑期备考2024小学生古诗文大会:吃透真题和知识点(持续)

2024年上海市小学生古诗文大会的自由报名初赛将于10月19日(星期六)正式开始,还有3个多月的时间。 为帮助孩子们备考,我持续分享往年上海小学生古诗文大会真题,这些题目来自我去重、合并后的1700在线题库,每…

加密与安全_密钥体系的三个核心目标之完整性解决方案

文章目录 Pre机密性完整性1. 哈希函数(Hash Function)定义特征常见算法应用散列函数常用场景散列函数无法解决的问题 2. 消息认证码(MAC)概述定义常见算法工作原理如何使用 MACMAC 的问题 不可否认性数字签名(Digital …

SketchUp Pro 2024:现代科技之诗意体验

在那遥远的唐朝,李白曾以诗酒为伴,游历山川,挥洒才情。而今,若李白穿越时空,手握现代科技之利器——SketchUp Pro 2024,定会以诗意之笔,描绘这款软件的神奇与魅力。 初识SketchUp Pro 2024 初…

Vue Router:History 模式 vs. Hash 模式

在开发 SPA(单页应用程序)时,路由管理是不可或缺的一部分。Vue.js 框架中的 Vue Router 提供了两种主要的路由模式:History 模式和 Hash 模式。理解这两种模式的区别及其实现方式,对于开发和部署 Vue 应用至关重要。 …

k8s record 20240708

一、PaaS 云平台 web界面 资源利用查看 Rancher 5台 CPU 4核 Mem 4g 100g的机器 映射的目录是指docker重启后,数据还在 Rancher可以创建集群也可以托管已有集群 先docker 部署 Rancher,然后通过 Rancher 部署 k8s 想使用 kubectl 还要yum install 安…

如何分析前后端bug?

如何分析前后端bug,我来支你1️⃣招 📝一般通过查看接口的方式分析前后端bug。 . 【方法】: web项目,用浏览器自带的F12抓包看接口请求。 app客户端,一般用fiddler等工具进行抓包接口。 . ✅用经典的电商项目举例&…

应用软件受到网络攻击怎么办?

大家都知道在目前的互联网社会中,大型的网络游戏与电商网站企业是网络攻击的重要对象,同时软件应用也无法避免地会受到各种网络攻击,那么当我们的软件应用被攻击时,该怎么办呢? 首先我们可以使用高防CDN,安…

2. 年龄问题

年龄问题 题目描述 本题为填空题,只需要算出结果后,在代码中使用输出语句将所填结果输出即可。 S 夫人一向很神秘。这会儿有人问起她的年龄,她想了想说: "2020 年前,我丈夫的年龄刚好是我的 22 倍,…

ATA-8035射频功率放大器在声动力疗法中的应用

声动力疗法是一种基于声波能量的治疗方法,广泛应用于医疗和美容领域。它利用高强度聚焦的声波来实现切割、破碎或加热组织,以治疗各种疾病和美容问题。在声动力疗法中,射频功率放大器起着至关重要的作用,它负责提供足够的能量来激…

达梦数据库的系统视图v$auditrecords

达梦数据库的系统视图v$auditrecords 在达梦数据库(DM Database)中,V$AUDITRECORDS 是专门用来存储和查询数据库审计记录的重要系统视图。这个视图提供了对所有审计事件的访问权限,包括操作类型、操作用户、时间戳、目标对象等信…

详解 | 什么是GeoTrust

GeoTrust是一家全球知名的数字证书颁发机构(Certificate Authority,简称CA),专注于提供SSL/TLS证书和其他相关的网络安全产品。 1、历史背景: GeoTrust成立于2001年,最初作为一个独立的公司运营。2006年&a…

js+spring boot实现简单前后端文件下载功能

jsboot项目实现自定义下载 一、前端页面 1、先导入axios的js包 2、注意axios响应的格式:result.data.真实的数据内容 3、这里请求的url就是你boot项目的getMapping的url,保持一致即可 4、如果想在后端设置文件名,那么后端生成后&#xf…