本人刚接触索寞结构及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

无人问津啊，期待有人指点下啊
论坛毕竟不是百科全书，耐心等待，也许会有答案

EM2=8e8/100000     !膜的虚拟弹性模量
-----------------------------------------------

膜的虚拟弹性模量降低太多了，降低100-1000倍比较合适，索的弹性摸量也要降低的，楼主可以试试

索的模量降低似乎没有影响，因为前面应变给的很大（0.99），事实上，我改过索的模量，影响不是很大。
另外，膜的模量我改到过100-1000（书上都这么说），最后也没有成功。我之所以改到100000，是因为有些论文上说降温幅度DETT在-0.1°~-1°比较合适。还是没有成功，希望多些人帮忙改改......

楼主的模型有点问题，第一好像忘记了生成边界索单元，第二，有些数据可能搞错了单位，比如，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

1、如何生成索？
2、楼上改完后，初次找形的结果是不错了，为何我在多次迭代求解之后，反而不及第一次找形后的应力均匀？
3、楼上仁兄如果方便，请贴出后面迭代部分....

在索膜结构找形中，索膜的弹性模量设置，边界设置（刚性边界或者柔性边界），对于最终的找形结果影响很大。出现问题的时候，多对上述参数进行调整，同时最好用另外的软件进行验证对比。
索膜找形是个比较复杂的过程，慢慢试试就好了