本文介紹了ANSYS后處理結(jié)果修改數(shù)據(jù)庫(kù)中的某些參數(shù)的若干技巧。

　　重啟動(dòng)在非線性和動(dòng)力學(xué)分析中經(jīng)常用到，ANSYS的重啟動(dòng)分析能力非常強(qiáng)大。在非線性分析中，經(jīng)常會(huì)遇到如下情況：

　　在一個(gè)非線性或動(dòng)力分析過(guò)程中，隨著載荷逐漸增大，應(yīng)力或者變形也將增大，而材料的特性(如剛度)以及實(shí)常數(shù)(如材料厚度等)也會(huì)隨著應(yīng)力或者變形的變化而改變，材料特性又不能夠通過(guò)定義應(yīng)力應(yīng)變關(guān)系來(lái)確定，需要有人工根據(jù)計(jì)算結(jié)果來(lái)調(diào)整。

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

　　/Title,根據(jù)前一步計(jì)算結(jié)果，修改參數(shù)進(jìn)行重啟動(dòng)分析

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

　　/prep7

　　*dim,tk,array,neg

　　*dim,em,array,neg

　　*dim,pr,array,neg

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

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

　　*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

　　建立幾何模型

　　指定屬性

　　分網(wǎng)

　　建立單元

　　finish

　　/solution

　　antype,,

　　邊界條件

　　載荷

　　求解選項(xiàng)(非線性或瞬態(tài))

　　求解控制

　　輸出控制

　　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 !修改實(shí)常數(shù) (Thickness)

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

　　pr(je)=pr(je) ! 修改材料常數(shù)(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

　　!求解選項(xiàng)

　　!求解控制

　　!載荷

　　......

　　solve

　　finish

　　*enddo