标题:
APDL索膜找形-程序纠错
[打印本页]
作者:
LEEQING
时间:
2010-3-23 11:08
标题:
APDL索膜找形-程序纠错
本人刚接触索寞结构及APDL找形,最近尝试对在找形论文中经常出现的一例子进行找形,进过反复调试,还是有很多错误,现求助于索膜结构中的高手们,忘不吝赐教,你的每一句话,也许对我都是莫大的帮助.....
!马鞍形膜面找形分析
!鞍形膜面找形2010/3/26
finish
/clear
/prep7
EM1=1.6E11 !索的弹性模量
A1=49/10000 !索的面积
ISTRAN=0.999 !索的应变
EM2=8e8/100000 !膜的虚拟弹性模量
H=2 !节点提升高度
T=1/1000 !膜厚
BOS=0.3 !泊松比
SIG0=3000 !膜的预应力
F=30000 !索的预应力
REF=1.0 !热膨胀系数
DETT=-(1-BOS)*SIG0/(EM2*REF) !降低温度
ET,1,LINK10
ET,2,SHELL41
MP,EX,1,EM1
MP,PRXY,1,BOS
MP,EX,2,EM2
MP,PRXY,2,BOS
MP,ALPX,2,REF
R,1,A1,ISTRAN
R,2,T
K,1,-8,0,0$K,2,8,0,0
K,3,0,8,0$K,4,0,-8,0
L,1,4$L,4,2
L,2,3$L,3,1
A,1,4,2,3 !
LATT,1,1,1
AATT,2,2,2
DK,1,UX,,,,,UY,UZ !在关键点上定义约束条件。
DK,2,UX,,,,,UY,UZ
DK,3,UX,,,,,UY,UZ
DK,4,UX,,,,,UY,UZ
DK,1,UZ,H
DK,2,UZ,H
DK,3,UZ,-H
DK,4,UZ,-H
lesize,all,,,20
lsel,all
ASEL,ALL
mshape,1 !指定单元为三角形单元
mshkey,1 !指定映射网格划分
amesh,all
TREF,0
ASEL,ALL
BFA,ALL,TEMP,DETT !温度荷载
!第一次找形分析
/solu !保存模型及求解数据
antype,0 !申明分析类型。0表示静力分析
nlgeom,on !打开大变形
sstif,on !打开应力刚度选项
nsubst,15 !子步数
outres,all,all
KBC,0
time,1
solve
finish
!有形.............
!更新坐标
/prep7
upcoord,1 !更新所有节点坐标
finish
!多次迭代
*do,i,1,5
/prep7
EM1=1.6E11 !索的真实弹性模量
EM2=8E8 !膜的真实弹性模量
MP,EX,1,EM1
MP,EX,2,EM2
F=2E4
A1=49/10000
ISTRAN=F/A1*EM1 !索的应变量
R,1,A1,ISTRAN
R,2,T
DETT=-(1-BOS)*SIG0/(EM2*REF)
/solu
DK,1,UX,,,,,UY,UZ !在关键点上定义约束条件。
DK,2,UX,,,,,UY,UZ
DK,3,UX,,,,,UY,UZ
DK,4,UX,,,,,UY,UZ
TREF,0 !参考温度0°
BFA,ALL,TEMP,DETT
ASEL,ALL
solve
finish
/prep7
upcoord,1
finish
*enddo
附件:
101file002.jpg
(2010-3-23 11:08, 46.29 K) / 下载次数 2
http://bbs.okok.org/attachment.php?aid=121154
作者:
LEEQING
时间:
2010-3-26 10:44
标题:
APDL索膜找形-程序纠错
无人问津啊,期待有人指点下啊
论坛毕竟不是百科全书,耐心等待,也许会有答案
作者:
xxqxxq2002
时间:
2010-3-27 15:15
EM2=8e8/100000 !膜的虚拟弹性模量
-----------------------------------------------
膜的虚拟弹性模量降低太多了,降低100-1000倍比较合适,索的弹性摸量也要降低的,楼主可以试试
作者:
LEEQING
时间:
2010-3-28 10:10
索的模量降低似乎没有影响,因为前面应变给的很大(0.99),事实上,我改过索的模量,影响不是很大。
另外,膜的模量我改到过100-1000(书上都这么说),最后也没有成功。我之所以改到100000,是因为有些论文上说降温幅度DETT在-0.1°~-1°比较合适。还是没有成功,希望多些人帮忙改改......
作者:
xxqxxq2002
时间:
2010-3-30 20:50
楼主的模型有点问题,第一好像忘记了生成边界索单元,第二,有些数据可能搞错了单位,比如,SIG0=3000 !膜的预应力,其实应该为3000000才符合实际,我对楼主的命令稍微修改了 ,效果还不错,就是不知道怎么发送图片。晕死。后半部分楼主自己完善了。
finish
/clear
/prep7
EM1=1.6E8 !索的弹性模量
A1=0.0002 !索的面积
EM2=8e8/1000 !膜的虚拟弹性模量
H=2 !节点提升高度
T=1/1000 !膜厚
BOS=0.0 !泊松比
SIG0=3e6 !膜的预应力 !SIG0=3000 (施加的力太小了,注意单位)
F=30000 !索的预应力
ISTRAN=F/(A1*EM1) !索的应变
REF=1.0 !热膨胀系数
DETT=-SIG0/(EM2*REF) !降低温度
ET,1,LINK10
ET,2,SHELL41
MP,EX,1,EM1
MP,PRXY,1,BOS
MP,EX,2,EM2
MP,PRXY,2,BOS
MP,ALPX,2,REF
R,1,A1,ISTRAN
R,2,T
K,1,-8,0,0$K,2,8,0,0
K,3,0,8,0$K,4,0,-8,0
L,1,4$L,4,2
L,2,3$L,3,1
A,1,4,2,3 !
LATT,1,1,1
AATT,2,2,2
DK,1,UX,,,,,UY,UZ !在关键点上定义约束条件。
DK,2,UX,,,,,UY,UZ
DK,3,UX,,,,,UY,UZ
DK,4,UX,,,,,UY,UZ
DK,1,UZ,H
DK,2,UZ,H
DK,3,UZ,-H
DK,4,UZ,-H
lesize,all,,,20
lsel,all
ASEL,ALL
mshape,1 !指定单元为三角形单元
mshkey,1 !指定映射网格划分
amesh,all
lmesh,all
TREF,0
ASEL,ALL
BFA,ALL,TEMP,DETT !温度荷载
!第一次找形分析
/solu !保存模型及求解数据
antype,0 !申明分析类型。0表示静力分析
nlgeom,on !打开大变形
sstif,on !打开应力刚度选项
nsubst,15 !子步数
outres,all,all
KBC,0
time,1
solve
finish
作者:
LEEQING
时间:
2010-4-2 10:21
标题:
APDL索膜找形-程序纠错
1、如何生成索?
2、楼上改完后,初次找形的结果是不错了,为何我在多次迭代求解之后,反而不及第一次找形后的应力均匀?
3、楼上仁兄如果方便,请贴出后面迭代部分....
作者:
lianrenweiman
时间:
2010-5-18 21:14
在索膜结构找形中,索膜的弹性模量设置,边界设置(刚性边界或者柔性边界),对于最终的找形结果影响很大。出现问题的时候,多对上述参数进行调整,同时最好用另外的软件进行验证对比。
索膜找形是个比较复杂的过程,慢慢试试就好了
欢迎光临 中华钢结构论坛 China Structure Forum (http://bbs.okok.org/)
Powered by Discuz! 5.0.0