Question

我的ANSYS命令流检查了无数遍了，就是没法求解，求各位大神帮我看看是哪里出了问题，跪求啊。。。。

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!第06讲：认识材料非线性 !!!
!!!!算例ex06b !!!
!!!!泰勒实验 !!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

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和导向点，然后再梁上一点施加扭矩，但是模型建出来，位移约束和荷载都加了，就是运行求解命令后没反应，没有结果。。。。。哎，弄了两天了 ，求大神帮忙啊。。急急急。

Answer 1

运行了一下你的程序,程序有警告!(问题不是很大)
我去掉了你接触设置,扭转载荷施加,运行了你的程序,能运行的,也就是说在单纯拉伸载荷作用下是有结果的。那就说明问题出在接触单元设置及扭转施加这两个方面。
这里推荐你看一下这个链接:http://www.jxcad.com.cn/topic_1166053.html
这里详细地讲了施加扭转载荷的具体过程……
希望能帮到你……

Answer 2

beam189的定位点错误，定位点与梁轴线在一条直线上。如看改为：
k,1002,LL,1.02*LL

追问

谢谢了，不过没分了。。

追答

哈哈，无所谓分不分的，重在交流。

Answer 3

哥们，你这个命令流有警告，你得先把警告解决了。

追问

那个警告不是问题 知识提示，我查过以前也遇到过相似问题，所以主要问题不是那个警告

Answer 4

看看节点合并了没，也许结构没固定住