如题
求命令流

命令流暂时没有，不过我用ansys做过很多算例，效果还是不错的。具体过程大致如下：
1。建立几何平面投影。
2。选择单元，膜用shell41,索用link10，设定材料常数。通常将两种单元的弹摸将低100倍。将索的应变放大100倍，设定膜的热膨张系数，参考温度。
3。化分单元，膜用三角形单元。
4。固定支撑点，提升支坐到设定高度。
5。打开大变形和应力刚化，求解。
6。更新座标，改材料常数为真实值，重新固定支撑点。
7。关掉大变形，求解。
8。循换若干次，应力分布均匀度可达95%以上。
还不清楚，联系emai:yzxzrq@yahoo.com.cn

非常感谢楼上这位大哥的回答！我试着写了一下，却发现如果我的模型是一个矩形的话，我用mshape,1,2D划分是，总是有一个警告
   Structural elements without mid nodes usually produce much more accurate results in quda or brick shape.
   这会影响计算结果吗？
   盼相告！！！

      发一个膜结构找形的算例吧，LOG文件做的修改，供参考。复杂些的也可以照此比划！！！！我做了很多。
!!!定义单元与材料特性
/PREP7   
ET,1,SHELL41                  
KEYOPT,1,1,2
ET,2,LINK10   
R,1,0.001, , , , , ,
RMORE, ,  
R,2,0.000314,0.9999999,   
MP,EX,1,8e6
MP,PRXY,1,0.3
MP,alpx,1,1
UIMP,1,REFT,,,0  
MP,EX,2,2e8
MP,PRXY,2,0.3   
K,1,,,,  
K,2,15,,,   
K,3,7.5,6,,  
K,4,7.5,-6,,
  
!!!定义结构几何模型
LSTR,       1,       3        
LSTR,       3,       2   
LSTR,       2,       4   
LSTR,       4,       1   
FLST,2,4,4   
FITEM,2,4   
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
AL,P51X  
AATT,       1,       1,   1,       0,     
LATT,2,2,2, , , ,      
LESIZE,ALL, , ,15, ,1, , ,1,       !!!定义单元边界长度
MSHAPE,1,2D  
MSHKEY,1  
CM,_Y,AREA   
ASEL, , , ,       1  
CM,_Y1,AREA  
CHKMSH,'AREA'   
CMSEL,S,_Y   
AMESH,_Y1                          !!!三角划分面单元
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2
  
FLST,2,4,4,ORDE,2                   !!!划分索单元
FITEM,2,1   
FITEM,2,-4   
LMESH,P51X      
FINISH
  
!!!设定边界并且第一次求解
/SOLU   
FLST,2,4,3,ORDE,2   
FITEM,2,1   
FITEM,2,-4   
/GO  
DK,P51X, , , ,0,ALL, , , , , ,   
NLGEOM,1
NROPT,AUTO, ,   
LUMPM,0  
EQSLV, , ,0,
PRECISION,0  
MSAVE,0  
PIVCHECK,1   
SSTIF,ON  
FLST,2,2,3,ORDE,2   
FITEM,2,3   
FITEM,2,-4     
/GO  
DK,P51X, ,5, ,0,UZ, , , , , ,   
/STATUS,SOLU
SOLVE   
EPLOT   
UPCOORD,1,OFF   
FLST,2,2,3,ORDE,2   
FITEM,2,3   
FITEM,2,-4   
/GO  
DK,P51X, ,0, ,0,ALL, , , , , ,   
FINISH   
/PREP7   
MP,EX,2,2E+11   
MP,PRXY,2,0.3     
R,2,0.000314,0.9999999e-3,   
FINISH   
/SOL
FLST,2,1,5,ORDE,1   
FITEM,2,1     
BFA,P51X,TEMP,-0.1   
/STATUS,SOLU
SOLVE   
EPLOT   
FINISH   
/POST1   
PLDISP,0
/VIEW, 1 ,,1
/ANG, 1  
/REP,FAST   
FINISH  
   
!!!第一次循环更新坐标，释放不平衡应力
n=15
*do,i,1,n
/solu
antype,static
sstif,on
/solu
solve
upcoord,0.5,off,lx,rst
*enddo
  
!!!施加真实应力
/SOL
FLST,2,1,5,ORDE,1   
FITEM,2,1   
!*   
BFA,P51X,TEMP,-2
/STATUS,SOLU
SOLVE   
FINISH   
/POST1   
AVPRIN,0, ,  
!*   
PLESOL,S,EQV,0,1
FINISH
  
!!!第二次循环更新坐标，释放不平衡应力
n=15
*do,i,1,n
/solu
antype,static
sstif,on
/solu
solve
upcoord,0.5,off,lx,rst
*enddo  
  
!!!还原膜的真实弹性模量
/PREP7     
MP,EX,1,8E+8   
MP,PRXY,1,0.3   
FINISH   
/SOL
FLST,2,1,5,ORDE,1   
FITEM,2,1      
BFA,P51X,TEMP,-2e-2  
/STATUS,SOLU
SOLVE   
FINISH   
/POST1   
AVPRIN,0, ,   
PLESOL,S,EQV,0,1
  
!!!第三次循环更新坐标，释放不平衡应力,找形结束
n=3
*do,i,1,n
/solu
antype,static
sstif,on
/solu
solve
upcoord,0.5,off,lx,rst
*enddo  
   
  
  

请问civil fans兄，在更新几何坐标时，到底用upcoord还是用upgeom,两者有什么区别！

不错，CIVIL FANS的用ANSYS找形的方法对我用动力松弛法和非线形有限元找形也很有启发。我们知道对伞锥的膜结构，很多是不存在等应力分布的。
但是我们在找形的初始赋值的时候一般是赋给等应力，然后通过小杨士模量来使应变的变化反映到应力的变化上去，从而找出非等应力的膜面。
但是有个问题是如果这种伞锥的控制点不是变化给出的话，预应力在达到符合要求的非等应力分布之前，伞锥就颈缩了。
CIVIL FANS这种让一些控制点逐步提升的方法实际上就象牛顿-拉费孙法中的荷载步，让荷载是分步加上去的从而给预应力的变化有一个充分的缓冲！
对我很有启发，在此表示谢意！

回答fatgao
  两个命令实际上是有差别的，其差别可以参看help文件，在此我就不多说了。不过这两个命令用于索膜结构找形时，其功能是一致的。
   关于jinsongxia所论，确如其言！有两点补充，实际索膜结构找形，理论上是可以得到等应力曲面的，在ansys中通过增加坐标更新的循环次数可以实现，索网结构一般循环5--10次，膜结构则可能要循环20---50次，根据结构复杂程度而定。分步提升和循环更新坐标形式不一样，但有异曲同工之秒。实际工程结构的应力分布差别控制在20---30%内一般就可以了。

请问civil fans兄：
上面列出的命令流中，你所说的关掉大变形nlgeom,off怎么没有体现呢？
是不是不关上也可以？若关上的话要好收敛些？
我在你的命令流中试着关上，找形结果有些变化。
望能告知 ，多谢了！

请各位大侠看下贴并请恳请提供一个命令流。(关于ＡＮＳＹＳ找形的难题)
  
http://okok.org/forum/viewthread.php?tid=47943&pg=2&bpg=1

多谢这位大侠了！真是救我于水火

对于ansys和索膜我都是初学者,有几个比较弱的问题向各位请教
1:一开始为什么要把弹性模量缩小一百倍
2:第一次更新坐标之后为什么要加一个温度力
3:第一个循环中为什么更新坐标的系数是0.5,就是为了多算几次!为什么多算几次后能实现等应力曲面
先谢谢能给我指点的大侠

我作的 感觉挺麻烦的  还没有恢复E  大家提意见

我们需要该伞型结构设计的兼职人员，但要求在北京，有感兴趣的吗？
请联系我：13070153010