!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!参考同济大学郭小农教授的mooc!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
0 基本教程
0.1 界面介绍
0.1.1 preprocessor
前处理模块
负责建模?约束?加力?
0.1.2 solution
求解器模块
0.1.3 general poseproc
常规的后处理模块
0.1.4 timehist postpro
时间处理相关模块
0.2 关于mac文件的一些写法注意
0.2.1 基本介绍
像txt文件一样,里面写着ansys的一些命令。
可以直接在命令窗口里面输入mac文件的名字,比如是ex01a.mac
ex01a
或者
/input,ex01a,mac
- 同时保存之后还会有个.db文件,是数据库。
- 还有个.rst文件,里面保存了计算结果。
- file.log文件可以用来查询运行过程中的错误。
0.2.2 文件规范
(1)注释
在mac文件中,!表示注释
一般开始的时候都会有如下命令
finish ! 结束当前命令
/clear ! 清除当前的所有操作
/filename,ex01a ! 命名当前文件名为ex01a
0.3 plotctrls
0.3.1 numbering
可以打开nodes or elements的编号
0.3.2 pan zoom(放大镜?)
1 求解简支梁的挠度
要注意在ansys里面没有单位,所以要换算成国际单位制,也就是m、N、Pa这种
1.1 简单注释和最开始的三句命令
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!第01讲:求解简支梁的挠度!!!!!!
!!!!!!算例ex01a !!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
finish !! 结束当前命令
/clear !! 清除当前的所有操作
/filename,ex01a !! 命名当前文件名为ex01a
1.2 参数输入
像excel一样去定义试件的一些信息
/prep7
表示进入preprocessor模块
以下就要定义一些参数,比如说,基本几何性质(截面性质等)、物理本构、力的大小。
!!!!参数输入!!!!
/prep7 !! 表示进入preprocessor模块
L=6 !!跨度
NN=20 !!有限元的分段
Ex=2.00e11 !!弹性模量
h=0.4 !!方管截面宽度
b=0.2 !!方管截面高度
t=0.01 !!方管壁厚
ep=0.002 !!边缘纤维应变
Area=b*h-(b-2*t)*(h-2*t) !!截面面积
Iz=(b*h**3-(b-2*t)*(h-2*t)**3)/12 !!截面惯性矩
Wz=Iz/(h/2) !!截面模量
Mz=ep*Ex*Wz !!根据边缘纤维应变计算得到弯矩
输入之后可以打开Parameter-scalar parameter查看之前定义的参数。
1.3 前处理模块(preprocessor)
在这个模块中,我们去定义了单元类型、材料特性、单元、节点、弯矩等等
这里的1都表示第一种。
- 单元类型(element type et):(beam3似乎在19.2中已经被淘汰,但是还是不影响,所以可以用)
19.2中好像只有beam4了,这是一个三维梁,而beam3是一个平面梁,显然更简单一些。
- 材料特性(material proporty mp):ex表示线弹性材料?mp里面还有很多材料特性,可以去看看。
- 实常数(realconstant r):注意这里实常数的输入的顺序要和单元类型中规定好的顺序要一样。
定义好之后,在preprocessor里面的相关模块可以看到定义好的东西
!!!!前处理!!!!
/prep7
et,1,beam3 !!定义单元类型(element type)
mp,ex,1,Ex !!定义材料特性
r,1,Area,Iz,h !!定义单元的实常数
- 采用循环来定义节点(nodes n): 注意节点数量
循环的形式:
*do,i,1,n !! 就相当于for i in range(1,n+1) 这里的最后一个数是取的到的
STATEMENT
*enddo
dx=L/NN !!每段单元长度
*do,i,1,NN+1 !!定义节点
n,i,(i-1)*dx,0,0 !!对第i个节点,设定它的x y z
*enddo
- 采用循环来定义单元(elements e)
*do,i,1,NN !!定义单元
e,i,i+1 !!每个单元都是从第i个节点到第i+1个节点,所以这里的i其实指的就是节点的编号
*enddo
如果element消失,可以在plot中点选elements让他显示出来
同时,在plotctrls里面的numbering可以打开编号
- 施加约束(其实就是让节点按我们想要的去变形,所以是deformation d)
d,1,ux,0 !!对1号节点施加约束(x方向位移为0)
d,1,uy,0
d,NN+1,uy,0
- 施加荷载(force f)
注意这里的弯矩是右手定则,向外为正。
f,1,mz,-Mz !!施加弯矩,对1号节点施加z方向弯矩,大小(包括正负)为-Mz
f,NN+1,mz,Mz
1.4 求解器模块(solution)
就是有限元迭代求解
!!!!求解!!!!
/solu !!进入求解器模块
solve !!开始求解
在大变形的前提下,需要打开大变形开关,启用时间子步。
time,1 !!总时间步1
nsubst,500 !!子步为500步 n substep
nlgeom,on !!打开几何大变形开关 nonlinear geometry?
1.5 后处理模块(postprocessor)
想要查看变形图?弯矩图?就需要进入到这个模块来
为了查看形状,需要打开形状开关,/eshape,1
但是显示的话,还需要右键,replot一下
plnsol 可以画变形图?plot n solution
!!!!后处理!!!!
/post1 !!进入到后处理模块
/eshape,1 !!打开形状开关(可以看到截面)
plnsol,u,y,0,1.0 !!画y方向的变形u,比例是1.0