在ADINA中如何实现节点坐标的update?

rongtian

在分析中想将变形后的位置作为结构的初始位置重新进行分析，即新的节点坐标＝原节点坐标＋节点位移，在ADINA中应如何实现，此功能在膜分析等方面很重要。
下面是一个简单的命令流
  
*
* Command file created from session file information stored within AUI database
*
*--- Database created 16 February 2004, 00:00:00 ---*
*--- by ADINA: AUI version 8.0.6 ---*
*
DATABASE NEW SAVE=NO PROMPT=NO
FEPROGRAM ADINA
CONTROL FILEVERSION=V80
*
MATERIAL ELASTIC NAME=1 E=2.00000000000000E+09 NU=0.300000000000000,
     DENSITY=0.00000000000000 ALPHA=0.00000000000000 MDESCRIP='NONE'
*
COORDINATES POINT SYSTEM=0
@CLEAR
1 0.00000000000000 0.00000000000000 0.00000000000000 0
2 0.00000000000000 10.0000000000000 0.00000000000000 0
3 0.00000000000000 10.0000000000000 10.0000000000000 0
4 0.00000000000000 0.00000000000000 10.0000000000000 0
@
*
SURFACE VERTEX NAME=1 P1=1 P2=2 P3=3 P4=4
*
FIXBOUNDARY LINES FIXITY=ALL
@CLEAR
4  'ALL'
@
*
FIXBOUNDARY POINTS FIXITY=ALL
@CLEAR
4  'ALL'
1  'ALL'
2  'ALL'
3  'ALL'
@
*
SUBDIVIDE MODEL MODE=LENGTH SIZE=0.500000000000000 NDIV=1,
    &nbspROGRESS=GEOMETRIC MINCUR=1
*
EGROUP SHELL NAME=1 DISPLACE=DEFAULT MATERIAL=1 RINT=DEFAULT,
     SINT=DEFAULT TINT=2 RESULTS=STRESSES STRESSRE=GLOBAL PRINTVEC=0,
     NLAYERS=1 INITIALS=NONE FAILURE=0 SECTION=0 CMASS=DEFAULT,
     STRAINS=DEFAULT RUPTURE=ADINA TIME-OFF=0.00000000000000,
     OPTION=NONE
*
GSURFACE NODES=9 PATTERN=AUTOMATIC NCOINCID=BOUNDARIES NCEDGE=1234,
     NCVERTEX=1234 NCTOLERA=1.00000000000000E-05 SUBSTRUC=0 GROUP=1,
     PREFSHAP=AUTOMATIC MESHING=MAPPED SMOOTHIN=NO DEGENERA=NO,
     COLLAPSE=NO MIDNODES=CURVED METHOD=ADVFRONT FLIP=NO
@CLEAR
1
@
*
LOAD PRESSURE NAME=1 MAGNITUD=2000000.00000000 BETA=0.00000000000000,
     LINE=0
*
APPLY-LOAD BODY=0
@CLEAR
@
*
APPLY-LOAD BODY=0
@CLEAR
1  'PRESSURE' 1  'SURFACE' 1 0 1 0.00000000000000 0 -1 0 0 0  'NO',
     0.00000000000000 0.00000000000000 1 0
@
*--- Database saved 16 February 2004, 00:00:00 ---*

rongtian

原始构形

rongtian

变形后的构形如下，想将变形后的构形作为初始构形

jinjinhao

restart好像可以,我没研究过这个.楼下的详细说说,呵呵

feizhang

版主，我感觉restart好像不行，因为restart不仅继承了第一步的变形，而且继承了第一步的内力，楼主的意思应该是将变形后的位置当作初始位置，而不考虑由此引起的应力，在ADINA中可以实现，在第一次计算后输出每一个节点的position(变形后的位置坐标），作一下格式转换后作为节点坐标输入即可。

rongtian

正如feizhang兄所理解的，能否说的更详细点，此问题已经困扰我好久了，不解决它我的工作就无法进展，急！！

feizhang

过程如下：
1. 对计算结果进行后处理，输出节点变形后的位置（＝原始坐标＋节点位移），并输出为文件保存（此处假设文件名为g1)。如下图

feizhang

2、对文件进行修改。
    原始文件如下：
  
  ADINA: AUI version 8.0.6, 17 February 2004:  *** NO HEADING DEFINED ***
  Licensed from ADINA R&D, Inc.
  Finite element program ADINA, response range type load-step:
  Listing for zone WHOLE_MODEL:
     &nbspOINT         X-POSITION      Y-POSITION      Z-POSITION
  
  Time 0.00000E+00
  
  Node 1            0.00000E+00     0.00000E+00     0.00000E+00
  Node 2            0.00000E+00     5.00000E-01     0.00000E+00
  Node 3            0.00000E+00     1.00000E+00     0.00000E+00
  Node 4            0.00000E+00     1.50000E+00     0.00000E+00
  Node 5            0.00000E+00     2.00000E+00     0.00000E+00
  Node 6            0.00000E+00     2.50000E+00     0.00000E+00
  Node 7            0.00000E+00     3.00000E+00     0.00000E+00
  Node 8            0.00000E+00     3.50000E+00     0.00000E+00
  Node 9            0.00000E+00     4.00000E+00     0.00000E+00
......
Node 1670        -5.60560E-01     4.25000E+00     9.75000E+00
  Node 1671        -6.13770E-01     4.75000E+00     9.75000E+00
  Node 1672        -6.50054E-01     5.25000E+00     9.75000E+00
  Node 1673        -6.67455E-01     5.75000E+00     9.75000E+00
  Node 1674        -6.64802E-01     6.25000E+00     9.75000E+00
  Node 1675        -6.41729E-01     6.75000E+00     9.75000E+00
  Node 1676        -5.98703E-01     7.25000E+00     9.75000E+00
  Node 1677        -5.37063E-01     7.75000E+00     9.75000E+00
  Node 1678        -4.59092E-01     8.25000E+00     9.75000E+00
  Node 1679        -3.68077E-01     8.75000E+00     9.75000E+00
  Node 1680        -2.68132E-01     9.25000E+00     9.75000E+00
  Node 1681        -1.61746E-01     9.75000E+00     9.75000E+00
  
  *** End of list.
  
将其修改为如下（掐头去尾）并保存。
Node 1            0.00000E+00     0.00000E+00     0.00000E+00
  Node 2            0.00000E+00     5.00000E-01     0.00000E+00
  Node 3            0.00000E+00     1.00000E+00     0.00000E+00
  Node 4            0.00000E+00     1.50000E+00     0.00000E+00
  Node 5            0.00000E+00     2.00000E+00     0.00000E+00
  Node 6            0.00000E+00     2.50000E+00     0.00000E+00
  Node 7            0.00000E+00     3.00000E+00     0.00000E+00
  Node 8            0.00000E+00     3.50000E+00     0.00000E+00
  Node 9            0.00000E+00     4.00000E+00     0.00000E+00
......
Node 1670        -5.60560E-01     4.25000E+00     9.75000E+00
  Node 1671        -6.13770E-01     4.75000E+00     9.75000E+00
  Node 1672        -6.50054E-01     5.25000E+00     9.75000E+00
  Node 1673        -6.67455E-01     5.75000E+00     9.75000E+00
  Node 1674        -6.64802E-01     6.25000E+00     9.75000E+00
  Node 1675        -6.41729E-01     6.75000E+00     9.75000E+00
  Node 1676        -5.98703E-01     7.25000E+00     9.75000E+00
  Node 1677        -5.37063E-01     7.75000E+00     9.75000E+00
  Node 1678        -4.59092E-01     8.25000E+00     9.75000E+00
  Node 1679        -3.68077E-01     8.75000E+00     9.75000E+00
  Node 1680        -2.68132E-01     9.25000E+00     9.75000E+00
  Node 1681        -1.61746E-01     9.75000E+00     9.75000E+00

feizhang

使用excel读入该文件，并删除第一列，如下图。

feizhang

4、已文本文件（制表符分隔）的格式江文件保存为g3. 用记事本打开该文件，并将文件的第一行修改为空行，保存该文件。

feizhang

回到ADINA前处理，通过meshing>elements>element nodes打开节点列表，如下图。

feizhang

将上表clear,并用import输入文件g3。
大功告成，说起来复杂，其实作起来不难。

feizhang

前面的步骤1中有点小错误，在输出变形后的节点位置时最好选择single response.如下图

BosFEA

看来feizhang兄用ADINA是非常深入了，感谢feizhang的贡献。