线膨胀系数求救

dongfengye

本人做了一简单的线膨胀方面的例题，所求结果值特别小，请高手指点原因，是否因为命令流设置有问题还是对线膨胀理解有问题
/clear
/FILENAME,example,1
l=5
b=0.5
/prep7
k,1,0,0,0
k,2,l,0,0
k,3,0,b,0
k,4,l,b,0
a,1,2,4,3
et,1,plane13
mp,ex,1,2.1e10
mp,alpx,1,1.2e-5
allsel
lesize,1,,,10
lesize,3,,,10
lesize,2,,,1
lesize,4,,,1
type,1
mat,1
allsel
amesh,all
/solu
antype,static
allsel
nsel,s,loc,x,0
d,all,ux,,,,,uy
allsel
nsel,s,loc,x,l
d,all,ux,,,,,uy
allsel
ic,all,temp,20
allsel
nsel,s,loc,y,0
bf,all,temp,15
allsel
nsel,s,loc,y,b
bf,all,temp,25
allsel
solve

ghost1903

你用的是什么长度单位？如果是毫米的话，弹性模量是2.1e5MPa，如果是米的话，弹性模量是2.1e11Pa，而你给的弹性模量是2.1e10，所以求出的变形回很小的

dongfengye

谢谢楼上的朋友回答,不过还是有问题,混凝土的弹模是E10量级,我用的单位是米,关键是求出来的应力量级是E-8次方,用结构力学法可以算出是不可能的,所以我想可能是我在计算中某个命令理解出错了!

nwpuyl

运行一下，先找到一个问题，约束没加上！
发现是单元定义不到位，plane13默认状态只有AZ自由度！
所以可以把单元定义那一句：“et,1,plane13”改为
et,1,plane13,4或者在其后面添加一句命令：keyopt,1,1,4
楼主先看看能解决问题不，其它问题再仔细检查！

xhcugb2005

版主一针见血 佩服

dongfengye

先谢谢楼上的朋友,你的回答确实是我的一个错误,不过按照你的建议,算出来的依然答案很小,所以问题可能还有其他的,希望各路朋友帮一把,谢谢!

dongfengye

本人最后用SOLID65试算了一个,结果X方向的应力是1.41MPA,与材料力学法算出来的1.26MPA相差不多,我想应该是对的吧,把命令流贴上来给大家参考一下:
/clear
/FILENAME,example,1
l=5
b=0.5
/prep7
k,1,0,0,0
k,2,l,0,0
k,3,0,b,0
k,4,l,b,0
a,1,2,4,3
vext,all,,,,,0.5
et,1,solid65
mp,ex,1,2.1e10
mp,prxy,1,0.17
mp,alpx,1,1.2e-5
allsel
lsel,s,,,1,7,2
lesize,all,,,10
allsel
lsel,s,,,2,8,2
lesize,all,,,2
allsel
lsel,s,,,9,12
lesize,all,,,2
mat,1
allsel
vmesh,all
/solu
antype,static
allsel
nsel,s,loc,x,0
d,all,ux,,,,,uy,uz,
allsel
nsel,s,loc,x,l
d,all,ux,,,,,uy,uz,
!allsel
!bf,all,temp,20
allsel
nsel,s,loc,y,0
bf,all,temp,-5
allsel
nsel,s,loc,y,b
bf,all,temp,5
allsel
solve

csloy

alpx=1.2E-5 好大的值哦，是什么材料哦