【BiliBili】 【知乎】 【QQ空间】 【微博】


博学慎思,坚持创新,求知若飢,虛心若愚                       关键字     类型    
  • 博主
  • 全部
  • 日志
  • 研究
  • 论文
  • 自编程序
  • 工程
  • Dino笔记
  • DinoBox
  • OpenSEES
  • Dino公开课
  • 每周挑战
  • Dino爱编程
  • 留言板

DinoBox

ETE_Frame3D
结构分析程序
ETE_Frame3D
2022-01-19
ETE-section截面分析程序
ETE-section
2021-08-19
ETO_OpenSEES前后处理程序
程序通过邮件申请。
★注意,软件不能下载,是通过邮箱发送软件的。

申请ETO时,注意EMAIL 的写法:(格式一定要对,否则电脑收不到)
邮件题目: ETO (版本号) 下载
注册人: (真实姓名)
学校/单位: (大学或工作单位)
邮件回复,大约需要1~4天的时间。
ETO_程序
2014-03-26
上一篇 | 下一篇
Dino爱编程 订阅所有Dino爱编程的日志

[Dino爱编程]有限元概念形函数ShapeFunction

作者:dinochen1983 日期:2021-06-17

字体大小: 小 中 大

【Dino爱编程】有限元理论中用到的形函数

关键词:数学,有限元,生物学,形函数
【程序下载】点击下载:程序的下载地址:ShapeFunctionProgram.rar
形函数(Shape Function),在知乎有一篇关于形函数的介绍,写得不错,“有限元简单科普之形函数”,
https://zhuanlan.zhihu.com/p/92684111
但笔者想通过编程的方法及展示一本书的内容来介绍我对形函数的理解。

1. 以四结点四边形单元为例,由于有4个结点,所以形函数写成如下公式。
当然,公式的推导过程可以看有限元的相关教科书。


基于上述公式,我们了解到这个N_i就叫做形函数,其实每个形函数N_i,其实是一个关于ξ与η的函数,完整的形态应该写成N_i (ξ,η),所以在参考坐标系ξη中,随着坐标(ξ,η)点的变化,这个参数是发生变化的。比如我们要求求解高斯积分点(0.57735,0.57735)的形函数,我们就把ξ=0.57735,η=0.57735 代入公式中计算得到全部的形函数。

2. 以四结点四边形单元为例,可以通过以下求和公式,把参考坐标系的点变成真实二维坐标系的点。这相当于映射变换,其中u_i与v_i为真实坐标的四个边角点(已确定),参考坐标点(ξ,η)可以确定全部的形函数N_i值,通过以下公式就可以把(ξ,η)找到对应真实坐标点(u,v),这个点的对应关系叫做映射,随着4点角点的变化,映射后的形状也会发生变化。参考坐标内的图形,由真实坐标的四个角点控制形状,所以这个函数叫做形函数(Shape Function)。


通过四个角点来控制形状,如果大家使用过一些图形软件如Rhino, 3dsmax, photoshop就一定很熟悉这个功能,这个就是几何拉伸(控制点变形),如下图所示。
 
图,3DSMAX的三维图形拉伸功能FFD2X2X2

图 PhotoShop的二维图形拉伸功能


形函数也体现了自然界的生物形状相关性,这涉及到动物遗传理论及进化论的延伸,虽然很多现象是很明显的,但现代科学还不能完全给出解析,这里列举了一本书的内容,是一个数学家达西·汤普森的著作《On Growth and Form》,达西·汤普森是一位英国苏格兰数学家、生物学家和古典学家,是数理生物学的早期先驱之一,主要作品有《生长和形态》(On Growth and Form),我拿到了书的简读版,如下图所示。好书是从孙一非同学那里借来的。本书就是让我大开眼界,没错,我又来推荐大家读数学书了。
   

神书《On Growth and Form》, 作者:达西·汤普森

书上展示了动物相似性的规律,最典型的例子就是二维形态的生物:鱼。不同种类的鱼通过一次(4结点)或二次(8结点)的形函数变换就可以得到对应的形状。还有人的头骨、猩猩的头骨,狒狒的头骨通过二次形函数即可变换得到。可见遗传基因控制着某种发育机制最终通过形函数的表达出来。即为什么形函数叫做形函数,控制形态的函数。
 
图 不同种类的鱼的几何变换关系

  

图 不同人类,猩猩与狒狒头骨的几何变换关系
图 不同昆虫外形几何变换关系


笔者通过编程的方法实现4结点与8结点的几何映射计算,首先列出8点的形函数如下图所示。

注意,8个结点的顺序每个教科书都有所不同,顺度会影响最终形函数的表达式。

那么八结点四边形单元的形函数,边界与内部点之前的关系是二次曲线关系。形函数公式如下所示。



笔者通过编程实现图形的形函数4结点与8结点的变换,特别是8结点变换在图形软件中比较少见。
关键程序的源代码可见附件。

 

图,8结点形函数对正方形的变换



示例, Argyropelecus olfersi与Sternoptyx diaphana 两种鱼的变换关系(斜切)

示例, Polyprion 与 Pseudopriacanthus altus 的几何变换关系。


示例, Diodon 与 orthagoriscus 鱼的几何变换关系。

有限元中形函数的应用就是把4结点8结点的任意四边单元,规则化成正方形单元,那么两都的应力场,位移场需要通过形函数建立起对应关系,那么求解过程中的积分点,只需要通过正方形的坐标点表达,对于四个高斯积分点的情况就是高斯积分点(0.57735,0.57735)等,这个积分点的坐标就是参考坐标了。那么应变(对位移求导)的关系,需要通过一个叫做雅克比矩阵建立联系,这个以后再聊。

本篇讨论”形函数“,讲述一个横跨数学几何,生物学(物种进化),编程与有限元都相关的题目,这也是跨学科学习的一种乐趣。
【程序下载】点击下载:程序的下载地址:ShapeFunctionProgram.rar




[本日志由 dinochen1983 于 2021-11-08 04:08 PM 编辑]
文章来自: 本站原创
引用通告: 查看所有引用 | 我要引用此文章
Tags: 有限元有限元 数学数学
评论: 0 | 引用: 0 | 查看次数: 4248
发表评论
昵 称:
密 码: 游客发言不需要密码.
验证码: 验证码
Dino:如果看不到验证码,再刷几下就可以了
内 容:
选 项:
虽然发表评论不用注册,但是为了保护您的发言权,建议您注册帐号.
字数限制 200 字 | UBB代码 关闭 | [img]标签 关闭

信息(Information)

DinoChen 视频频道

My GitHub Page



扫描以下二维码,加入公众号!


dinochen.com需要你的支持


日历(Calendar)

2022年7月
  • 日
  • 一
  • 二
  • 三
  • 四
  • 五
  • 六
  • 26
  • 27
  • 28
  • 29
  • 30
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20
  • 21
  • 22
  • 23
  • 24
  • 25
  • 26
  • 27
  • 28
  • 29
  • 30
  • 31
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6

查找(Search)

关键字
 
类 型

文章(Recent Article)

[OpenSEES]OpenSEES实例教程[汇总...[Dino笔记]粒子群算法在桁架结构中的应用[获奖]WSP亚洲2021年卓越奖获得一等奖与创新...[Dino笔记]基于遗传算法对剪力墙结构优化设计方...[Dino笔记]基于OpenSEES虚功原理的钢结...[OpenSEES] 新版ETO更新及相关下载方法...[新闻]祝韩老师荣获广东省科技进步一等奖[程序]ETE-section功能更新:可修改划分...[App开发]采用Swfit语言编制经典的推箱子游...[程序]ETE-section更新_考虑混凝土保护...

评论(Comments)

陈博士您好,这个小工具好像有点问题,有效翼板宽度b...打开的时候显示“你没有许可证”,怎么弄啊陈博士您好,这个是需要购买的软件还是可以试用的怎么查看文档计算过程呢,陈博士一直显示出现一个错误字符无法计算是什么原因呢?为什么一直出错计算不了?陈博士,能共享下文件嘛?第二版没有光盘陈博士您好,我换了好几台电脑ETE还是需要许可,陈博士您好,现在程序还是打不开,需要许可请问如何获得会员许可证?

留言(GuestBook)

Thanks regarding furnishi...陈博士您好,想请问一下零长度单元用在桥梁的支座上,...陈博,新版ETO下载以后界面打不开是什么原因呢,旧...陈博士您好,请问可以出一期画易损性曲线的教程吗陈博士您好,我在学习您4.33侧向多自由度简化模型...陈博士您好,計算wind2019的program ...陈博士,您好,我下载使用了您编写的压型钢板组合楼板...陈博士您好,我是opensee初学者,想向您请教两...陈博士你好,请问开洞混凝土截面,ETE-secti...陈博士您好,感谢您的回复,我的刚性隔板设置了转动惯...

TagBar

Loading...

用户面板(User Panel)

登录用户注册

统计(Statistics)

日志: 408 篇
评论: 893 个
留言: 413 个
访问: 14024155 次
在线: 10 人
建站时间: 2008-02-08

Powered By PJBlog3 V2.7 Build 05 CopyRight 2005 - 2011, 陈学伟的博客

Processed in 0.281250 second(s) , 6 queries 粤ICP备15071789号-1