板壳单元计算示例
图8-12圆弧旋转角度设置对话框
在ARC选项栏中填写旋转角度180°，生成圆弧。 建立关键点，点号分别为4，5，6 （ 0.05,0,-0.1 ），（ -0.05,0,-0.1 ），（ 0,0.05,-0.1 ）。 建立完关键点后沿关键点生成线，将图形连接如图 8-13所示图形。运行Modeling>Operate>Booleans>Add> Lines，利用布尔运算将线相加：将CD，DE相加为一条 线，另外将圆弧的两段线加为一条线。 生成面：运行 Preprocessor>Modeling>Create>Areas>Arbitrary>By lines ，拾取相交后的线可得到如图 8-14 所示的图形。
ANSYS板壳单元计算示例
1 问题描述
开口的圆弧为 180 °，半径为 0.05m ， CA 长为 0.1m ，壳体厚 0.001m ， CDE =90 °，弹性模量 E=210GPa，μ =0.3。约束：边CD和边DE全约束。承 受载荷：B点作用集中载荷F＝10N，方向水平向右。
图8-5 壳体示意图
图8-10 设置实常数对话框
（4）建立模型 生成圆环，选择Preprocessor> Modeling>Create>Lines>Arcs>By Cent&Radius ，弹出如图 8-11 所示 对话框，填写数据0,0,0，单击 Apply 按钮，接着对话框提示再选 取一点，填入0.05作为半径，单击 OK 按钮，弹出如图 8-12 所示对话 框。
图8-13 壳体线框模型
图8-14 壳体模型
（5）划分网格 设置网格大小：运行MeshTool，弹出如图8-15所示 对话框，在Smart Size前面打勾，将下面的条形框拉到 最右边。点击Mesh，拾取壳体。划分网格完 毕，如图8-16所示。
图8-16 壳体有限元模型
（6）施加约束 选择菜单Solution>Define Loads>Apply > Structure>Displacement>On Lines，拾取线 段CD和DE，选择All DOF在这两条线段上 施加全约束。
2 ANSYS求解操作过程
（1）选择单元类型 运行Preprocessor>Element Type>Add/Edit/Delete，弹 出Element Types对话框如图8-6所示，然后单击Add，弹 出 Library of Element Types 窗口，如图 8-7 所示，选择 Shell Elastic 4node 63，选择完毕单击OK按钮。图8-18 云图显示对话框
3 结论
从图8-19壳体的位移云图可知，最大位移发生在B点 处，且最大位移值为0.000691m。 从图8-20壳体的应力云图可知，最大应力发生在B点 处，且最大应力值为40.1MPa。
图8-19 位移变形云图
图8-20 等效应力云图
图8-17 壳体约束与载荷
（9）后处理 运行 General Postproc>Plot Results>Contour Plot>Nodal Solu，弹出如图8-18所示对话框。运行DOF Solution>Displacement vector sum和Stress>von Mises stress ，分别显示壳体位移和应力云图。结果如图 8-19 和图8-20所示。
图8-15 设置网格尺寸对话框
（7）施加载荷 执行命令Solution>Define Loads>Apply >Structure> Force/Moment>On Keypoints，拾取B点, 施加载荷FX＝ 10N，如图8-17所示。
（8）求解 选择Solution>Solve> Current LS ，开始计算， 计算结束会弹出计算完毕 对话框，单击 Close 关闭 对话框，计算完毕。
图8-5 壳体示意图
图8-6 单元类型对话框
（2）设置材料属性 运行 Preprocessor> Material Props>Material Models，弹出图8-8所示 对话框，双击Isotropic， 弹出如图8-9所示对话框， 在EX 选项栏中设置数值 2.1e11，在PRXY选项栏 中设置数值0.3。设 置完毕单击OK按钮。
图8-7 单元类型库对话框
图8-8 选择材料属性对话框
图8-9 设置材料属性对话框
（3）设置壳厚参数 运行Preprocessor> Real Constants>Add/Edit/ Delete，弹出如图8-10所示 对话框，设置 TI(k) 选项栏 为0.001。设置完毕单击OK 按钮完成设置。
ANSYS有限元分析——ANSYS板壳问题实例
2、分析类型静力分析3、问题描述板壳问题4、ANSYS单位m N Pa5、单元SHELL63：4节点（每个节点6个自由度）6、材料弹性模量和泊松比7、实常数厚度8、建模取1/4模型41进入ANSYS程序→ANSYSED 9.0→ANSYS Product Launcher→change the working directory into yours →input Initial jobname:board→Run2 设置计算类型ANSYS Main Menu: Preferences →select Structural →OK3 定义材料参数ANSYS Main Menu: Preprocessor →Material Props→Material Models →Structural →Linear →Elastic→Isotropic →input EX:2e11, PRXY:0.3 →OK54 选择单元类型ANSYS Main Menu:Preprocessor →Element Type→Add/Edit/Delete→Add →select Shell Elastic 4node 63→OK (back to Element Types window) →Close (the Element Type window)65、定义实常数ANSYS Main Menu: Preprocessor →Real Constants…→Add…→select Type 1→OK→input Shell thickness at node I :0.01→OK →Close (the Real Constants Window)76、创建矩形ANSYS Main Menu:Preprocessor →Modeling →Create →Areas→Rectangle→By Dimensions →依次按下图输入：→OK87、划分单元ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool →OK(Size Controls) Global: Set →input NDIV:20→Mesh Tool →Mesh :select Areas→Shape:Quad→Free→Mesh →Pick All →Close( the Mesh Tool window)98、施加固定边界ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Displacement→On Lines→拾取上边和右边→OK→select Lab2:ALL DOF →OK109、施加对称约束边界ANSYS Main Menu: Solution→Define Loads →Apply→Structural →Displacement →Symmetry B.C →On Lines→拾取下边和左边→OK1111 分析计算ANSYS Main Menu: Solution →Solve →Current LS→OK(to closethe solve Current Load Step window) →OK12、挠度ANSYS Main Menu: General Postproc→Plot Results→Contour Plot→Nodal Solu→select: DOF solution, Z-Component of displacement→OK13按右侧的等侧视图141617板壳单元的压力载荷面1 (I-J-K-L)，（底部，+Z方向）面2 (I-J-K-L)，（顶部，-Z方向）面3 (J-I)，面 4 (K-J)，面 5 (L-K)，面 6 (I-L)19显示壳单元的厚度ANSYS命令菜单栏：PlotCtrls>Style >Size and Style→弹出对话框中（Display of element项设置为On)控制壳单元的结果输出ANSYS Main Menu: General Postproc→Options for OutpTop layer: 顶面Middle layer:中面Bottom layer：底面2021作业三要求1、属于力学的那类问题？2、单位制；单元类型；单元描述；实常数；材料参数3、划分网格的方案，施加载荷和边界条件4、计算结果的体现：MISIS 应力和位移，需要知道最大值以及位置。
板壳分析
∂w ∂w ∂ w [ χ ] = − 2 ; − 2 ; −2 [ N] de = [ B] de ∂x ∂y ∂xy [ B] = [ B]1 [ B]2 [ B]3 [ B]4
弹性薄板矩形（R12） 弹性薄板矩形（R12）单元
薄板的单元列式 在前面已导出
∂w ∂w ∂ w [ χ ] = − 2 ; − 2 ; −2 ∂y ∂xy ∂x
2 2 2 T
Q1 1 My1
Mx1
4 x w3
θx3 θy3
2 y z
3
将位移模式代入可得形变矩阵为
2 2 2 T
2
∂x
;v = −z
∂y
2
; w = w( x, y)
y向曲率 向曲率
扭率
2
他们完全确定板的变形， 他们完全确定板的变形，因此称他们组成的矩 阵为形变矩阵，记作[χ ， 阵为形变矩阵，记作 χ]，也即
∂w ∂w ∂w ε x = −z 2 ;ε y = −z 2 ;γ xy = −2z ∂x ∂y ∂xy
2) 形函数的确定 龙驭球院士指出利用对 薄板的形函数可以用广 称性可减少确定广义坐标 义坐标法， 义坐标法，也可以用试凑 的工作量， 的工作量 法得到。 ，但是仍然不很 法得到。由于单元自由度 方便。为此 介绍试凑法， 介绍试凑法， 方便，因此可有12个广义 。为此,介绍试凑法 为12，因此可有 个广义 2 首先引入自然坐标ξ=x/a, 坐标， 坐标，位移模式可设为如 η=y/b。 。 下不完全四次多项式
∂w ∂w ∂ w [χ] = − 2 ;− 2 ;−2 ∂y ∂xy ∂x
2 2 2
T
弹性薄板基本知识
由此可得薄板应变矩阵为[ 由此可得薄板应变矩阵为 ε]=z[χ]。 χ。 薄板内力和总势能 1) 设平面应力弹性矩阵为 设平面应力弹性矩阵为[D]’，则薄板应力矩阵 ， 为 [σ]=-z[D]’ [χ]。 2) 薄板内力 χ。 微元体如图所示。 微元体如图所示。 h/2 x 由图可得 τ xy zdydz ∫-h/2 h/2 ∫-yh/2 σ x zdydz τ xydydz h/2 σ x dydz τ yx zdxdz ∫
ansys板壳计算
4 板壳问题当比较薄的构件承受横向力作用时，它将发生弯曲变形，此时需要利用板壳理论来计算。
Ａｎｓｙｓ中板壳单元是同一种单元，最常用的是弹性壳体单元Ｓｈｅｌｌ６３，它可以模拟板，也可以模拟薄壳，或者厚壳。
壳体单元上的荷载是通过ＳＦ（Ｓｕｒｆａｃｅ Ｆｏｒｃｅ）来定义压力的。
ＳＦＡ和ＳＦＥ分别是对面积或者单元定义表面压力的。
壳体单元计算结果可以直接给出内外表面的应力值。
下面几个实例的代码均在Ａｎｓｙｓ５．６的ＥＤ版中调试通过。
４．１ ＤｒａｗＰｉｐｅ．ｔｘｔ 利用拉伸操作建立膨胀弯管模型 化工管道的弯管，主要用来释放温度和内压引起的轴向变形的，这里我们通过拉伸（Ｅｘｔｒｕｄｅ）操作建立一个膨胀弯管模型的一半。
由于ＥＤ版的限制，这里的直线管段长度较短。
为了方便起见，在直线段的端部采用了固定端约束（图中用彩色箭头标识），弯管上部的对称位置采用了对称面约束（图中用Ｓ表示）。
同样的原因，这里单元划分的尺寸取的较大，管道的应力计算结果不够理想． 下图给出的是管道在内压作用下的Ｍｉｓｅｓ应力。
全部代码如下： ＦＩＮＩＳＨ ！退出以前模块 ／ＣＬＥＡＲ，ＮＯＳＴＡＲＴ ！清除内存中的所有数据，不读入初始化文件 ／ＦＩＬＥＮＡＭＥ，ＤＲＡＷＰＩＰＥ ／ＴＩＴＬＥ，ＤｒａｗＰｉｐｅ．ｔｘｔ， Ａ ｓａｍｐｌｅ ｔｏ Ｇｅｎｅｒａｔｅ Ｅｌｂｏｗ ｂｙ Ｅｘｔｒｕｄｅ Ｏｐｅｒａｔｉｏｎ． ／ＰＲＥＰ７ ！进入前处理模块ＰＲＥＰ７ ｋ，１ ！定义位于管道轴线上的控制点１至６ ｋ，２，，，１ ｋ，３，，０．８，１ ｋ，４，，０．８，１．５ ｌ，１，２ ！定义通过这些控制点的直线 ｌ，２，３ ｌ，３，４ ｌｆｉｌｌｔ，１，２，．２８ ！对这些相交的直线定义倒角 ｌｆｉｌｌｔ，２，３，．２８ ｋ，２１！定义第２１号关键点 ｃｉｒｃｌｅ，２１，．１ ！以２１号关键点为圆心，０．１为半径作圆，圆弧线编号依次是６，７，８，９ ／ＶＩＥＷ，１，１，１，１ ＡＤＲＡＧ，６，７，８，９，，，１，４，２，５，３ ！开始沿前面定义的路径用圆扫描形成管道，管道轴线是１，４，２，５，３号线 ｅｔ，１，ｓｈｅｌｌ６３ ！定义第一类单元为壳单元ＳＨＥＬＬ６３ ｒ，１，．０１ ！定义单元的第一类实常数：Ｔｈｉｎｋｎｅｓｓ＝０．０１ ｍｐ，ｅｘ，１，３ｅ１１ ！定义第一类材料的弹性模量ＥＸ＝３ｅ１１ ＤＬ，６，ＡＬＬ，ＡＬＬ，０ ！定义６，７，８，９号线固定端条件 ＤＬ，７，ＡＬＬ，ＡＬＬ，０ ＤＬ，８，ＡＬＬ，ＡＬＬ，０ ＤＬ，９，ＡＬＬ，ＡＬＬ，０ ＤＬ，４２，，ＳＹＭＭ ！定义４２，４５，４７，４９号线对称条件 ＤＬ，４５，，ＳＹＭＭ ＤＬ，４７，，ＳＹＭＭ ＤＬ，４９，，ＳＹＭＭ ＳＦＡ，ＡＬＬ，１，ＰＲＥＳＳ，１ ！对所有面定义压力荷载，大小为１ ＥＳＩＺＥ，０．０８，０， ！单元划分尺寸０．０８ ＡＭＥＳＨ，ＡＬＬ ！单元划分 ／ＳＯＬＵ ！进入求解模块ＳＯＬＵＴＩＯＮ ＳＯＬＶＥ ！开始求解 ／ＰＯＳＴ１ ！进入后处理模块ＰＯＳＴ１ ＳＥＴ，１ ！读入第一个荷载步的计算结果 ＰＬＮＳＯＬ，Ｓ，１ ！在图形窗口显示结点上的主应力 ＦＩＮＩＳＨ ！退出后处理模块 ４．２ Ｊｕｅｚｕ．ｔｘｔ 管道噘嘴模型受内压作用 噘嘴是管道加工过程中的一类主要缺陷，是指在卷边时边缘部分的曲率没有达到管道的曲率要求，是圆截面管道在焊缝位置处出现一个凸起。
2014-计算力学-9-板壳问题
的不同，它除了弯曲变形外还存在着中面变形，所以壳体中 的内力包括有弯曲内力和中面内力。
壳体弯曲问题
在壳体理论中，有以下几个计算假定： ① 垂直于中面方向的正应变极其微小，可以不计。
② 中面的法线总保持为直线，且中面法线及其垂直线段之间 的直角也保持不变，即这两方向的剪应变为零。
③ 与中面平行的截面上的正应力（即挤压应力），远小于其 垂直面上的正应力，因而它对变形的影响可以不计。 ④ 体力及面力均可化为作用在中面的载荷。
板壳问题
平板弯曲问题
矩形单元
壳体弯曲问题
平板弯曲问题
在弹性力学里，把两个平行面和垂直于这两个平行面的柱面或
棱柱面所围成的物体称为平板，简称为板，如图9-1所示。两个 板面之间的距离t称为板的厚度，而平分厚度t的平面称为板的中 间平面，简称中面。如果板的厚度 t 远小于中面的最小尺寸 b （如小于b/8～b/5），该板就称为薄板，否则就为厚板。
N yi a i (1 0 )(1 0 )(1 2 ) / 8
式中 0 i ，0 i。
矩形单元
矩形单元的刚度矩阵 矩形单元的刚度矩阵可以写成如下形式：
k11 k k 21 k 31 k 41
h/2
其中子矩阵为：
T
k12 k 22 k 32 k 42
0
,
y
z 0
0
,
xy
z 0
0
这就是说，中面的任意一部分，虽然弯曲成为弹性曲面的一 部分，但它在 xy 面上的投影形状却保持不变。
矩形单元
按薄板弯曲的基本假定，板内各点的位移为：
u z w x
w v z y
w w( x, y )
7_板壳问题有限元分析
1 1 2 h 1 1 2
h
BiT DB j abd d dz
(6.17)
21 /44
薄板问题的有限元法
代入 D 、 Bi 和 B j 于是有
D 1 1 b2 T kij N i , N j , uN iT, N T, uN iT, N T, j j 1 1 a 2 ab +2(1- )N
2
24 /44
薄板问题的有限元法
k23 15H ab(i j )(i j ) b2 b2 k31 3Ha (2 3 5 2 ) j0 15 2 j 5i0 a a k32 15H ab(i j )(i j )
23 /44
薄板问题的有限元法
其中
b2 a2 a2 b2 k11 3H 0 15( 2 0 2 0 ) (14 4 5 2 5 2 ) 00 b b a a a2 a2 k12 3Hb (2 3 5 2 ) 0i 15 2 i 5 0i b b b2 b2 k13 3Ha (2 3 5 2 )i0 15 2 i 50 j a a a2 a2 k21 3Hb (2 3 5 2 ) 0 j 15 2 j 5 0i b b a2 k22 Hb 2(1 ) 0 (3 50 ) 5 2 (3 0 )(3 0 ) b
1 E D 2 1 0
薄板问题的有限元法
图 6.2 平板内力
10 /44
薄板问题的有限元法
设 M x 、 M y 和 M xy 表示单位宽度上的内力矩，于是有
2w 2 x Mx h h3 2 w h3 M M y h2 z dz D DC D 'C (6.5) 2 12 y 12 2 M xy 2w 2 xy
4 .板壳问题的有限元法(4学时)
第五章 板壳问题的有限元法
章节内容： 5.1 薄板弯曲的基本理论 5.2 薄板单元：矩形单元和三角形单元 5.3 薄壳有限元分析的简介
车辆工程教研室
机电工程学院
5.1 薄板弯曲的基本理论
5.1.1 薄板（thin plate）
工程实际中，存在大量的板壳构件（plate and shell） 几何特点：厚度远远小于其它两个方向的尺寸。 薄板：t/b < 1/15 中面：平分板厚度的平面 坐标系oxyz ：xy轴在中面上，z轴垂直于中面 z 载荷 作用于中面内的载荷：平面应力问题 垂直于中面的载荷：板弯曲
其中
车辆工程教研室
机电工程学院
5.5 薄壳有限元分析
局部坐标系
局部坐标系对整体 坐标系的方向余弦 矩阵（从整体坐标 到局部坐标）
局部坐标系与整体坐标系的关系
车辆工程教研室
机电工程学院
5.5 薄壳有限元分析
坐标变换矩阵
车辆工程教研室
机电工程学院
5.5 薄壳有限元分析
单元刚度矩阵
转换矩阵：
3.
应力
引起的形变很小，在计算变形时可以忽略。
车辆工程教研室
机电工程学院
5.1 薄板弯曲的基本理论
5.1.2 位移
位移分量：薄板中面的挠度 w 根据挠度，可以计算：在x和y轴方向上的位移分量和绕x和y轴方 向的转角。
y
z
b
o
车辆工程教研室
t
x
机电工程学院
5.1 薄板弯曲的基本理论
5.1.3 应变及几何方程
机电工程学院
5.1.5 平衡方程
板壳问题的有限元
N i Ni
N xi N yi
8
(i 1,2,3,4)
8
形函数
N i (1 0 )(1 0 )(2 0 0 2 2 ) N xi bi (1 0 )(1 0 )(1 2 ) N yi a i (1 0 )(1 0 )(1 )
D D 12 1 o o 1 o 1 o 2
Et 3 12 1 2
第五章 板壳问题有限单元法
一、薄板弯曲基本假定和基本方程
zy y xy z y x
考虑薄板上、下板面的边界条件
zy z t
0
2
zx z t 0
2
解得横向剪应力，为
E 2 t2 2 z zx w 2 4 x 21
E 2 t2 2 z zy w 2 4 y 21
2
上式中，利用板下面的边界条件 z z t q ，得：
Et 3 4 wq 2 12(1 )
D w q
4
Et 3 D＝ 12(1 v 2 )
D是板的弯曲刚度，板厚的三次方成正比，与弹模成正比，与梁的弯曲刚度类似
第五章 板壳问题有限单元法
一、薄板弯曲基本假定和基本方程
w a1 a2 a3 a4 2 a5 a6 2 a7 3 a8 2 a9 2 a10 3 a11 3 a12 3
x y 其中： , a b
第五章 板壳问题有限单元法
二、矩形薄板单元
2. 单元位移插值函数：
2
8
i 1,2,3,4
2014-计算力学-9-板壳问题
（9-7）
矩形单元
其中 a和b分别是单元的长和宽。将单元的四个节点坐标 分别代入（9-6）和（9-7）式，即可求得位移模式中的12个 参数，再代入（9-6）式，得
w
N w N
i i i 1
4xi xi 来自 yi yi N i i N e
由此得到：
x
（9-6）
w w 1 3 5 2 6 8 2 2 9 310 2 11 3 312 2 y b b w w 1 y 2 2 4 5 3 7 2 2 8 9 2 311 2 12 3 x a a
v w z y
故有
,
由于z =0， zx 0 , zy 0 ，所以中面的法线在薄板弯曲时保持不伸缩，成 为弹性曲面的法线。此外，由于不计z 所引起的应变，故其物理方程为
1 x y E 1 y y x E 2(1 ) xy xy E
（9-9）
式中 D
矩形单元
矩形单元的等效节点力
当平板单元受有分布横向载荷q时，其相应的等效节点力为
Qi
e
Wi M xi M yi
1
1
1 1
qN i T abdd （i = 1,2,3,4）
（9-10）
若q = q0 为常量时，有
1 1
k13 k 23 k 33 k 43
k14 k 24 k 34 k 44
kij Bi DB j dxdydz h / 2 1 1Bi T DB j abdd
D ab b2 a
2
4.板壳单元-Tang
e
3b ( 1 ) 0 b ( 1 3 )( 1 ) i i i i a i z 3a i (1 i ) ai (1 i )(1 3i ) 0 4ab b 2 2 2 2 ai (3 2 i 1) ii (3 3 4) b i (3 2i 1)
y w x ( 2 2 4 x 5 y 3 7 x 2 28 xy 9 y 2 311 x 2 y 12 y 3 )
17
利用12个节点位移值可待定12个系数，整理w(x,y)为插值
函数形式：
w( x, y ) Ni wi N xi xi N yi yi N e Nl wl N xl xl N yl yl
[S’]分块矩阵形式
' ' [ S ' ] [ S i' S 'j S m Sn ]
t 1 1 ~ 薄膜 t 厚度 b 80 100 1 1 t 1 1 ~ ~ 薄板 80 100 b 5 8 t 1 1 ~ 厚板 b 板长宽最小值 b 5 8
2
3
薄板基础理论知识
如图所示平板，取其中性面为坐标面，z轴垂直于 中性面。其中 t 为板厚。当板受有垂直于板中性面的外 力时，板的中性面将发生弯扭变形，从而变成一个曲面。 板变形的同时，在板的横截面上将存在内力——弯矩和 扭矩。
其中，形函数：
1 x y x x y y N r (1 )(1 )[2 (1 ) (1 )] 8 xr yr xr xr yr yr 1 x y 2 y N xr yr (1 )(1 ) (1 ) 8 xr yr yr 1 x y x N yr xr (1 ) 2 (1 )(1 ) 8 xr yr xr r i, j , m, l
【2019年整理】ANSYS壳体计算示例
图8 壳体线框模型
9. 运行Modeling>Operate>Booleans>Add>Lines 10. 利用布尔运算将线相加： 将CD，DE相加为一条线；将圆弧的两段线加为一条线 11. 生成面：运行Preprocessor>Modeling>Create>Areas>Arbitrary>By lines 12. 拾取相交后的线可得到如图9所示的图形
图5 设置实常数对话框
（4）建立模型
1. 选择Preprocessor>Modeling>Create>Lines>Arcs>By Cent&Radius，生成圆环 2. 弹出如图6所示对话框，填写数据0,0,0，单击Apply按钮
3. 接着对话框提示再选取一点，填入0.05作为半径 4. 单击OK按钮
图9 壳体模型
（5）划分网格
1. 设置网格大小：运行Meshing>Size Cntrls> Manual Size>Areas>All Areas 2.弹出如图10所示对话框，设置Element edge length为0.005
图10 设置网格尺寸对话框 3. 设置完毕点击OK
4. 运行Meshing>Mesh>Areas>Free，拾取壳体 5. 划分网格完毕，如图11所示
2.弹出Element Types对话框如图1所示
3.单击Add
图1
4.弹出Library of Element Types窗口，如图2所示
5.选择Shell Elastic 4node 63 6.单击OK按钮. 单元设置完毕
图2
（2）设置材料属性 1. 运行Preprocessor>Material Props>Material Models 2. 弹出图3所示对话框