作者ahuie (阿辉)
站内Mechanical
标题[问题] 有限元素法 任意四边形板元素
时间Thu Jun 11 22:52:26 2015
各位前辈好
小弟我研究所做的研究是有关FEM及振动学
我作的分析是要自行写程式(MATLAB)
求出刚性矩阵与质量矩阵
再用F=Kq 解出位移q
目前的Plate bending element板元素
有导出矩形元素及三角形元素
与ANSYS的振动频率结果比较差异不大
因为我题目是做圆板及椭圆板
必须用三角形或任意四边形切割
三角形元素做出来的结果与ANSYS差异不大
但我推导出的任意四边形做出来的结果
在与ANSYS同样Mesh下 答案比较起来有明显的差距
任意四边形是用自然座标系导的
Jocobian矩阵与一般四节点元素相同
高斯积分也取很多点了 看起来应该是收敛
找了很多书板的任意四边形元素都很少提到
期刊跟论文找过很多 正在努力看懂@@
有没有前辈做FEM 希望可以指导一下> <
PS. 用我推导的任意四边形元素做矩形板的振动模态
MESH都切割成正方形 答案会与ANSYS接近
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 1.34.51.102
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Mechanical/M.1434034349.A.AD2.html
※ ahuie:转录至看板 Cad_Cae 06/11 22:53
1F：推 choral: 你是把元素切到多细? 06/11 23:00
2F：→ ahuie: 1000多个元素 再增加看起来答案也收敛了@@ 06/11 23:11
3F：→ ahuie: 1200 06/11 23:13
4F：推 choral: 喔喔！元素量不多~ 那就跟初始设定有关了(大问题惹) 06/11 23:23
5F：→ ahuie: QQ 我在想应该是形状函数的问题 06/11 23:29
※ 编辑: ahuie (1.34.51.102), 06/11/2015 23:31:42
6F：推 choral: 任意四边形没啥涉猎的说 抱歉啦～ 06/11 23:31
7F：→ ahuie: 没关系QQ 还有很多时间可以做... 06/11 23:32
8F：推 lolitass: 有专版Cad_Cae 06/12 05:37
9F：推 gbcowandy: 你的元素编号是哪个？几个节点？ 06/12 07:57
10F：→ gbcowandy: 喔喔看到了4节点，你可以改成五节点看看 06/12 07:58
11F：推 jnth5566: cadcae很多高手 06/12 10:21
12F：推 ATand: 其实简单想你用四个点去舖圆形，你的边界不会好看到哪 06/12 21:39
13F：→ ATand: 三角点去舖面，自然与原来形状近似度最高 06/12 21:40
14F：→ ahuie: 我一开始想法也跟楼上一样 三角形本来就适合斜边与曲面 06/12 22:45
15F：→ ahuie: 边界分割 但是我们教授叫我试看看用四边形割 06/12 22:46
16F：→ ahuie: ANSYS用四边形做答案也满准 只是不知道他的方法... 06/12 22:46
17F：→ gamer: 我不懂任意四边形元素，不过标准元素三角形的收敛性远差於 06/14 23:32
18F：→ gamer: 四边形，至於为什麽，修过有限元素法你就会知道了。和数值 06/14 23:32
19F：→ gamer: 积分时转换到高斯积分域的形状函数有关，四边形元素在转换 06/14 23:33
20F：→ gamer: 时的误差会比三角形小。 06/14 23:33