作业帮我的ANSYS命令流检查了无数遍了,就是没法求解,求各位大神帮我看看是哪里出了问题,
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/25 02:16:28
我的ANSYS命令流检查了无数遍了,就是没法求解,求各位大神帮我看看是哪里出了问题,
finish
/clear
/filname,ex06b,1 定义文件名
/prep7
参数输入
PI=3.1415926535
DD=100e-3 管径
tt=4.0e-4 壁厚
RR=(DD-tt)/2 壁中面半径
LL=0.2 长度
Ex=2.06e11 弹性模量
Nu=0.3 泊松比
Fy=235e6 屈服强度
qy=Fy*tt/2 对轴向产生应力Fy的轴向拉力,N/m
MY=Fy*tt*2*PI*RR*RR/2 对环向产生应力Fy/2的切应力的扭矩,N*m
alfy=2.0 轴向荷载系数
alfr=2.0 径向荷载系数
定义单元类型,实常数,材料特性
et,1,shell181
et,2,beam189
et,3,conta175
et,4,targe170
keyopt,3,2,2
keyopt,3,4,2
keyopt,3,12,5
mp,ex,1,Ex
mp,prxy,1,Nu
TB,BKIN,1,1,2,1
TBDATA,Fy,0,
r,1,tt
r,2
sectype,1,beam,ctube
secdata,DD/2-tt,DD/2
建立几何模型
k,1,0,0,0
k,2,0,LL,0
k,3,DD/2-tt,0,0
k,4,DD/2-tt,LL,0
l,3,4
arotat,1,1,2,360,24
allsel,all
单元划分
lsel,s,length,LL
lesize,all,20*tt
lsel,all
aatt,1,1,1,0
aesize,all,PI*RR/24 设置单元划分尺寸
amesh,all
建立梁单元
k,1000,LL,0
k,1001,1.01*LL
k,1002,1.02*LL
l,1000,1001
l1=_return
lsel,s,l1
latt,1,1,2,1002,1
esize,tt
lmesh,all
allsel,all
定义接触单元和目标单元的导向节点,以方便施加施加扭矩
nsla,s,1 选择面上所有节点
nsel,r,loc,Y,LL
type,3
real,2
esurf
allsel,all
nbeam=node(0,LL,0)
type,4
tshap,pilo
e,nbeam
施加约束
allsel,all
nsel,s,loc,y,0,0
d,all,all
施加荷载
allsel,all
lsla,s
lsel,R,loc,y,LL,LL
SFL,all,PRES,-qy*alfy
allsel,all
f,node(0,1.01*LL,0),my,MY*almy
allsel,all
首先进行静力解
/solu
antype,0
nsubst,10
outres,all,all
arclen,1 弧长法
solve
finish
后处理
/post1
PLNSOL,s,eqv,0,1.0
这是一个薄壁筒,高为LL,在它的上面建一个长为0.01的梁,然后在梁和筒之间设置接触单元conta175和导向点,然后再梁上一点施加扭矩,但是模型建出来,位移约束和荷载都加了,就是运行求解命令后没反应,没有结果.哎,弄了两天了 ,急.
finish
/clear
/filname,ex06b,1 定义文件名
/prep7
参数输入
PI=3.1415926535
DD=100e-3 管径
tt=4.0e-4 壁厚
RR=(DD-tt)/2 壁中面半径
LL=0.2 长度
Ex=2.06e11 弹性模量
Nu=0.3 泊松比
Fy=235e6 屈服强度
qy=Fy*tt/2 对轴向产生应力Fy的轴向拉力,N/m
MY=Fy*tt*2*PI*RR*RR/2 对环向产生应力Fy/2的切应力的扭矩,N*m
alfy=2.0 轴向荷载系数
alfr=2.0 径向荷载系数
定义单元类型,实常数,材料特性
et,1,shell181
et,2,beam189
et,3,conta175
et,4,targe170
keyopt,3,2,2
keyopt,3,4,2
keyopt,3,12,5
mp,ex,1,Ex
mp,prxy,1,Nu
TB,BKIN,1,1,2,1
TBDATA,Fy,0,
r,1,tt
r,2
sectype,1,beam,ctube
secdata,DD/2-tt,DD/2
建立几何模型
k,1,0,0,0
k,2,0,LL,0
k,3,DD/2-tt,0,0
k,4,DD/2-tt,LL,0
l,3,4
arotat,1,1,2,360,24
allsel,all
单元划分
lsel,s,length,LL
lesize,all,20*tt
lsel,all
aatt,1,1,1,0
aesize,all,PI*RR/24 设置单元划分尺寸
amesh,all
建立梁单元
k,1000,LL,0
k,1001,1.01*LL
k,1002,1.02*LL
l,1000,1001
l1=_return
lsel,s,l1
latt,1,1,2,1002,1
esize,tt
lmesh,all
allsel,all
定义接触单元和目标单元的导向节点,以方便施加施加扭矩
nsla,s,1 选择面上所有节点
nsel,r,loc,Y,LL
type,3
real,2
esurf
allsel,all
nbeam=node(0,LL,0)
type,4
tshap,pilo
e,nbeam
施加约束
allsel,all
nsel,s,loc,y,0,0
d,all,all
施加荷载
allsel,all
lsla,s
lsel,R,loc,y,LL,LL
SFL,all,PRES,-qy*alfy
allsel,all
f,node(0,1.01*LL,0),my,MY*almy
allsel,all
首先进行静力解
/solu
antype,0
nsubst,10
outres,all,all
arclen,1 弧长法
solve
finish
后处理
/post1
PLNSOL,s,eqv,0,1.0
这是一个薄壁筒,高为LL,在它的上面建一个长为0.01的梁,然后在梁和筒之间设置接触单元conta175和导向点,然后再梁上一点施加扭矩,但是模型建出来,位移约束和荷载都加了,就是运行求解命令后没反应,没有结果.哎,弄了两天了 ,急.
运行了一下你的程序,程序有警告!(问题不是很大)
我去掉了你接触设置,扭转载荷施加,运行了你的程序,能运行的,也就是说在单纯拉伸载荷作用下是有结果的.那就说明问题出在接触单元设置及扭转施加这两个方面.
这里推荐你看一下这个链接:http://www.jxcad.com.cn/topic_1166053.html
这里详细地讲了施加扭转载荷的具体过程……
希望能帮到你……
我去掉了你接触设置,扭转载荷施加,运行了你的程序,能运行的,也就是说在单纯拉伸载荷作用下是有结果的.那就说明问题出在接触单元设置及扭转施加这两个方面.
这里推荐你看一下这个链接:http://www.jxcad.com.cn/topic_1166053.html
这里详细地讲了施加扭转载荷的具体过程……
希望能帮到你……