我想看一级黄色大片_久久亚洲国产精品一区二区_久久精品免视看国产明星_91久久青青青国产免费

您的位置:網站首頁 > Ansys教程

ANSYS根據后處理結果來修改數據庫中的某些參數

時間:2010-11-14 11:34:51 來源:未知

重啟動在非線性和動力學分析中經常用到,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#p#分頁標題#e#