本文介紹了ANSYS后處理結果修改數據庫中的某些參數的若干技巧。

　　重啟動在非線性和動力學分析中經常用到，ANSYS的重啟動分析能力非常強大。在非線性分析中，經常會遇到如下情況：

　　在一個非線性或動力分析過程中，隨著載荷逐漸增大，應力或者變形也將增大，而材料的特性(如剛度)以及實常數(如材料厚度等)也會隨著應力或者變形的變化而改變，材料特性又不能夠通過定義應力應變關系來確定，需要有人工根據計算結果來調整。

　　對于這樣的情況，可以采用逐步重啟動的方法來模擬。具體過程是，首先作一個初始分析，然后進入后處理器，對計算結果進行處理，對應力(或位移)超過一個規定值的單元或者單元組，按照一定的規則修改材料常數(相當于調整局部區域的剛度)以及實常數，再在前一步分析的基礎上，采用修改了的材料常數、實常數、載荷等進行重啟動分析。在這里，需要修改的數據可以采用參數的方式定義，根據計算結果修改這些參數，然后作為一個外部文件保存在磁盤上，在重啟動時，將這些數據文件以參數的形式讀入，用這些參數修改數據庫的相關數據，然后求解。為此，可以用APDL語言來編制程序，實現這個過程。下面就是實現這個過程的流程：

　　/Title,根據前一步計算結果，修改參數進行重啟動分析

　　*set,neg,, !屬性(材料、實常數)組數

　　/prep7

　　*dim,tk,array,neg

　　*dim,em,array,neg

　　*dim,pr,array,neg

　　*set,tk(1),,,, ! 設置實常數的初值(Thickness)

　　*set,em(1),,,, ! 設置材料常數的初值

　　*set,pr(1),,,,

　　et,1,181

　　*do,i,1,neg

　　mp,ex,i,em(i)

　　mp,prxy,i,pr(neg)

　　r,i,tk(i)

　　*enddo

　　建立幾何模型

　　指定屬性

　　分網

　　建立單元

　　finish

　　/solution

　　antype,,

　　邊界條件

　　載荷

　　求解選項(非線性或瞬態)

　　求解控制

　　輸出控制

　　solve ! 初始求解

　　*do,i,1,ns

　　/POST1

　　*get,emax,elem,,num,max

　　*dim,smax,array,emax

　　*dim,dmax,array,emax

　　etable,seqv,s,eqv

　　*do,ie,1,emax

　　*get,smax(i1),etable,ie,seqv

　　if,smax(i1),gt,slim,then

　　*get,je,elem,ie,attr,real

　　tk(je)=tk(je) tk(je)*0.10 !修改實常數 (Thickness)

　　em(je)=em(je)-ex(je)*0.10 ! 修改材料常數 (Ex)

　　pr(je)=pr(je) ! 修改材料常數(Prxy)

　　*endif

　　*enddo

　　*CREATE,outp_data

　　*CFOPEN,'inpu_ex',' ',' '

　　*VWRITE,em(1)

　　(f20.5)

　　*CFCLOS

　　*END

　　/INPUT,outp_data

　　*CREATE,outp_data

　　*CFOPEN,'inpu_pr',' ',' '

　　*VWRITE,pr(1),

　　(f20.5)

　　*CFCLOS

　　*END

　　/INPUT,outp_data

　　*CREATE,outp_data

　　*CFOPEN,'inpu_tk',’’

　　*VWRITE,tk(1)

　　(f20.5)

　　*CFCLOS

　　*END

　　/INPUT,outp_data

　　FINISH

　　/SOLU

　　ANTYPE,,REST,,,

　　*dim,em,array,neg

　　*dim,pr,array,neg

　　*dim,tk,array,neg

　　*CREATE,read_data

　　*VREAD,em(1),inpu_ex

　　(f20.5)

　　*end

　　/INPUT,read_data

　　*CREATE,read_data

　　*VREAD,pr(1),inpu_pr

　　(f20.5)

　　*END

　　/INPUT,read_data

　　*CREATE,read_data

　　*VREAD,tk(1),inpu_tk

　　(f20.5)

　　*END

　　/INPUT,read_data

　　*do,j,1,neg

　　mp,ex,j,em(j)

　　mp,prxy,j,pr(j)

　　r,j,tk(j)

　　*enddo

　　!求解選項

　　!求解控制

　　!載荷

　　......

　　solve

　　finish

　　*enddo