abaqus黑盒猜测和理论学习： S4R 单元的实现方式 板壳、梁理论
对称的壳单元（位于theory介绍壳体的较前部分）
对称的壳单元的6个假设
abaqus关于AXISYMMETRIC shell element的theory6个假设，可以看到实际上对称的壳单元满足Kir直法线假设
但以实验为准，见博主这篇博文，S4R单元（并不属于thin ,按abaqus属于general的范畴）的3轴应力应变都是0。
事实上，应力为0是因为S4R单元是结合了平面应力假设+Kir直法线假设的壳元.
abaqus 有限应变壳单元理论学习 重点学习S4R
S4R理论
看了一下abaqus 《theory》 3.6.5。
里面实际是连续介质力学的内容，讲了各个量——如法向量、曲率随时间的更新公式等。
此外还有上图开始讲了能预防剪切自锁的strain field formulation，并推导出了S4R这一类finite strain单元的
stiffness matrix，这是我真正需要的。但是Ks并不是刚度，它的单位是应力。Cs is the transverse shear stiffness
finite strain单元 3.6最后：
Element S4 is a fully integrated finite-membrane-strain shell element. Since the element’s stiffness is fully integrated, no spurious membrane or bending zero energy modes exist and no membrane or bending mode hourglass stabilization is used. Drill rotation control, however, is required. Element S4 uses the same drill stiffness formulation as used for element S4R. Similarly, element S4 assumes that the transverse shear strain (and force, since the transverse shear treatment is elastic based on the initial elastic modulus of the material) is constant over the element. Therefore, all four stiffness integration locations will have the same transverse shear strain, transverse shear section force, and transverse shear stress distribution. The transverse shear treatment for S4 is identical to that for S4R.
It is well known that a standard displacement formulation will exhibit shear locking for applications dominated by in-plane bending deformation. However, a standard displacement formulation for the out-of-plane bending stiffness is not subject to similar locking response.Hence, S4 uses a standard displacement formulation for the element’s bending stiffness, and the theory presented above for the rotation kinematics and bending strain measures applies to S4. The primary difference between the element formulations for S4 and S4R is the treatment of the membrane strain field. This formulation is the topic of the following discussion.
The membrane formulation used for S4 does not rely on the fact that S4 is a shell element. Hence, the discussion below details the formulation from the point of view that the membrane response is governed by the equilibrium for a three-dimensional body in a state of plane stress.
1.standard displacement formulation?
2.why “a standard displacement formulation for the out-of-plane bending stiffness is not subject to similar locking response”?
3.What is a membrane formulation, what is finite-membrane-strain ?
薄膜应变就是面内的那三个应变。剩下三个叫transverse shear strain。
abaqus theory前面40页也很重要，很基本的.
initial stress:预应力
我的手稿地址：
https://github.com/BraveDrXuTF/S4R_theoryPDF
通过这个锻炼，找到一个abaqus乃至其他CAE软件复杂文档的组成逻辑，帮助以后的阅读查找工作。
想要参考abaqus壳更多的内容，可以看analysis 29章，里面对于各类单元的讲解很详细。
经典的板壳力学教材和上课讲的东西能否成为突破口？结合看一下
图片来自于徐芝纶的《板壳力学》
我或许还需要N1,N2,M1,M2等数据
找到了RF,RM系列数据
什么叫薄膜内力？
轴对称壳元
轴对称壳元（s,θ\thetaθ两个坐标即可确定壳体内一点位置）的薄膜状态和弯曲状态，壳元的内力是有分类的，分为薄膜内力和弯曲内力：
轴对称斜锥壳元：
轴对称截锥壳元：以直代曲
几何方程—B—K
考虑横向剪切的轴对称壳体理论基本原理：
比不考虑横向剪切多了个γ\gammaγ.
考虑剪切变形的梁称为铁木辛柯梁。——有限单元法9.2.4
挠度和截面转角独立插值的称为铁木辛柯梁单元。——有限单元法9.2.4
保证K2奇异，K1非奇异，这样才能使极限情况下方程仍然有非零解。怎么理解？——有限单元法8.5节
我的理解（不一定正确）：极限情况下K2占据主导地位，方程(K1+K2)a=F退化为K2a=0，K2奇异才能有非零解。
位移和转动独立插值的截锥壳元
原理和铁木辛柯梁单元一样，也需要保证K2奇异，K1非奇异。
本科学的梁单元是不考虑剪切变形的梁单元，即《有限单元法》9.2.3节
考虑剪切变形有两种方式
1.叠加
wb部分方程和本科学的梁单元一样,ws代表剪切部分引起的附加法向位移，它的方程比较简单
2.铁木辛柯梁单元
分别插值。
考虑层间应力无法直接应用S4R，但可以加弹簧
http://www.doc88.com/p-9921475056045.html
声明：
近日作者发现，有部分网站：
https://wenku.baidu.com/view/70282ae40f22590102020740be1e650e52eacf01.html
https://www.pudn.com/news/623f187f2cfc38172120b885.html
在未通知作者的情况下以不恰当的方式转载了本文，将本文变为下载需要收费的文章，没有注明原作者和原始来源。并且第一个网址抹除了图片的水印，在此作者给予强烈的谴责，并保留追责的权利。
根据CC BY-SA 4.0协议，转载本文，必须附上原始作品的来源及本许可证的链接。
20220906
这一块的话，没想到S4R如此神通广大，竟然有如此的formulation。
3D单元复合材料建模，在文档abaqus analysis 4。
20220908
没想到啊没想到，S4R的层中是有TSHR13和23的。在上下边界这两个值则为0.
在下边界，我试了一阶剪切的strain formulation是符合abaqus的结果的，在里面-0.6667的位置，我的小算例printabaqus.py是符合的，而大规模test_data.py则暂时不符合：
20221031：S4R的厚度方向默认是3个积分点（亲测），而不是5个。
以4层为例，将section_num_tuple改为(1,2，…13)倒也不会报错，但是，odb中只会有12个积分点的结果。
20221104
成功验证abaqus S4R 单元在多层变刚度平板下使用的是一般壳理论。
20221124
成功校对abaqus S4R平板大部分位置的应变能。我还计算了特定工况下的外力做功，应变能-外力功=总势能。
但是有一点很奇怪，边界的平衡方程量严重的对不上，我想abaqus在其中定有计较，有可能涉及到沙漏控制等等。
上图是abaqus平板的平衡方程量，不是应变能哈。如图，大部分都对上了（平衡方程量=0就是对上了），边界的严重对不上。
abaqus黑盒猜测和理论学习： S4R 单元的实现方式 板壳、梁理论相关推荐
- 应用Abaqus有限元软件中的cohesive单元模拟压头侵入地层随机断裂过程
应用 Abaqus有限元软件中的 cohesive单元模拟压头侵入地层随机断裂过程. 建立的压头侵入地层的有限元模型如图2所示.其中球形压头直径为1mm,高为3mm,厚度为0.5mm.将其设置为刚体, ...
- ABAQUS中Cohesive粘聚力模型的2种定义方式（附案例操作步骤）
附赠仿真学习包,包含结构.流体.电磁.热仿真等多学科视频教程,点击领取: 仿真秀粉丝专属礼包 导读:大家好,我是仿真秀专栏作者-烧仙草.消费电子行业仿真,擅长胶材等材料的本构模型研究和构建 ...
- php excel 右对齐,excel中单元格对齐方式在哪里设置？
excel单元格对齐方式可以在"设置单元格格式"框的"文本对齐方式"中设置.设置方法:1.选中要设置的单元格或单元格区域:2.鼠标点击右键,点击"设置 ...
- python3-xlwt-Excel设置表格基础(字体颜色 背景颜色 设置边框 边框颜色 冻结窗口 字体大小 字体高度 字体加粗 下划线 斜体字 单元格对齐方式 自动换行 删除线 超链接 插入公式)
文章目录 引入xlwt和创建workboos对象 初始化样式和创建设置字体,赋值给style 保存文件 自定义样式 字体颜色 背景颜色 设置边框 边框颜色 冻结窗口 字体大小 字体高度 字体加粗 下划 ...
- Abaqus中C3D8R单元和C3D8I单元的区别
C3D8R单元: 8节点六面体线性减缩积分单元.与完全积分单元相比,线性减缩积分单元仅在单元中心包含一个积分点,而二次减缩积分单元的积分点数量与线性完全积分单元的相同. 线性减缩积分单元有以下优点: ...
- Abaqus怎么关闭和显示节点单元编号
Mesh----View-assembly display option–mesh(show node labels) View--Part Display Options--Mesh--Show n ...
- asp.net gridview删除 获取到第一行第一列的单元格内容_VBA中的常用单元格引用方式...
VBA编程经常和"对象"打交道,其中最频繁的对象大概就是"单元格"了.(听说您还没有对象?那--我想你大概需要一份Excel,包邮988--) 今天我们就来聊一 ...
- 纹理对象纹理单元纹理目标_网页设计理论：纹理
纹理对象纹理单元纹理目标 Texture has become an indispensable element in web design. It is not only a trend but a ...
- 设置单元格填充方式_单元格的选择及设置单元格格式
数据输入完毕,接下来可以设置字体.对齐方式.添加边框和底纹等方式设置单元格格式,从而美化工作表.要对单元格进行设置,首先要选中单元格. 选择单元格 选择单元格是指在工作表中确定活动单元格以便在单元格中 ...
最新文章
- Soul网关发布2.2.0：让高性能网关变得如此简单！
- mysql分页查询limit_MySQL查询语句(where,group by,having,order by,limit)
- ANTLR和Jetbrains MPS：解析文件并以树符号显示AST
- Python语法糖——遍历列表时删除元素
- python集合运算_从零开始学Python - 第014课：常用数据结构之集合
- return 输出为空php,thinkphp5 返回json数据的方法---以及返回json为空的原因
- springboot项目实例_Springboot项目的接口防刷的实例
- UnitTest测试框架全栈详解
- 通俗易懂的随机森林模型讲解
- 使用python解析pdf文档中的单词导入到有道云词典中
- ios storyboard 传参 返回传参（segue）
- 3Ds Max 精简材质编辑器 一片黑问题
- 开发者应知的 如何“推销自己”？
- 2020年8月中国编程语言排行榜
- 一些常见的处理器如arm，arduino，stm32，51，树莓派的联系和区别，还有各自的长短板？
- String字符串截取方式
- 极客日报第 73 期：Twitter 正谈判收购印度初创公司 ShareChat：欲打造 TikTok 竞品；华为 Mate X2 正式发布，17999 元起；爱奇艺员工10万“买”了个北京
- 微信公众号怎么集赞服务器,【微信点赞】_微信怎么点赞如何集赞微信点赞活动新手必备点赞教程...
- 当在浏览器中输入一个域名后，会发生什么
- Python 爬取财务报表