（一） 初始應力場施加實例實例1要在分層的SHELL181單元中的一指定層上施加初始應力,可以使用指定層的材料ID號。 本文介紹了ANSYS施加初始應力場的相關方法在用ANSYS做結構分析時,可以把初始應力指定為一項載荷，只在靜態分析和全瞬態分析中被承認(分析可以是線性,也可以是非線性), 初始應力載荷只能在分析的第一載荷步中施加.恒應力可以用ISTRESS命令指定.運用ISFILE命令從輸入文件中獲取初始應力值,為了列表或是刪除初始應力,命令變元允許初始應力被限制為具體的材料類型.要在分層的SHELL181單元中的一指定層上施加初始應力,可以使用指定層的材料ID號。在初始應力命令中使用MATx 字段可以根據層的材料ID號施加需要的應力，或者可以在SHELL181單元的所有層上施加初始應力，然后刪除那些不希望有初始應力的層上的初始應力。用delete, elementID, 層號命令，還可以用初始應力命令列表或刪除分層SHELL181單元上指定層上的初始應力值。初始應力功能只支持以下單元類型:PLANE2, PLANE42, SOLID45, PLANE82, SOLID92, SOLID95, LINK180, SHELL181, PLANE182, PLANE183, SOLID185, SOLID186, SOLID187, BEAM188,還有BEAM189.ISFILE命令可以將合成初始應力寫入文件。這些命令可以在/SOLUTION處理器中運用，在ANSYS Commands Reference(ANSYS命令參考手冊)中可以看到所有初始應力命令的描述。【1】首先將初始應力數據用一個外部的ASII文件給出，假設文明名為istress.ist，內容如下/CSYS,0! ELEM ID ELEM INTG LAY/CELL SECT INTG SX SY SZ SXY SYZ SX1,all,all,all,100,0,0,0,0,02,all,all,all,100,0,0,40,0,03,all,all,all,100,0,0,0,0,04,all,all,all,100,0,0,0,0,405,all,all,all,100,20,0,0,0,0【2】讀取應力場文件/prep7/title, Example of Initial stress import into ANSYSet,1,182! Plane stress PLANE182 elementmp,ex,1,1.0e9mp,nuxy,1,0.3! Define the nodes!n,1n,2,2.0n,3,4.0n,4,6.0n,5,8.0n,6,10.0n,7,1.0n,8,2.0,1.0n,9,4.0,1.0n,10,6.0,1.0n,11,8.0,1.0n,12,10.0,1.0! Define the 5 elements!e,1,2,8,7e,2,3,9,8e,3,4,10,9e,4,5,11,10e,5,6,12,11! Constrain all dofs on all nodes at x=0 to be zeronsel,s,loc,x,d,all,allnallfinish!/solu!讀應力文件inis,read,istress,ist! List the initial stressesinis,listoutres,all,allsolvefinish!/post1set,lastPLNSOL,S,EQV,0,1.0finish實例2 一個巖土工程中常用的巖層簡化平面模型同時受重力和地表均壓作用，比較以下兩種方法所獲的結果。 1.常規方法，同時施加重力加速度和均壓。 2.初應力載荷法，先計算模型在重力作用下的結果，同時輸出初應力文件，然后在全新的模型中調入此初應力文件，計算等效節點力反力；施加等效節點力，進行求解，得到與只有重力場作用下的模型同樣的受力和變形;最后施加均壓，求解。finish/filename,try01,1/prep7et,1,42mp,ex,1,2e11mp,nuxy,1,0.3mp,dens,1,7800blc4,1,2esize,0.1amesh,all /solu nsel,s,loc,y,0d,all,uy,0nsel,s,loc,x,0nsel,a,loc,x,1d,all,ux,0allsel acel,10,inistate,write,1,ssolve!得到只有重力作用下的初應力文件finish/clear,nostart /filename,try02,1/prep7et,1,42mp,ex,1,2e11mp,nuxy,1,0.3mp,dens,1,7800blc4,1,2esize,0.1amesh,all /solu d,all,all isfile,read,try01,ist,2isfile,list outres,all,all allsel solve!計算等效節點力反力finish/filename,try03,1/solu lsclear,all nsel,s,loc,y,0d,all,uy,0nsel,s,loc,x,0nsel,a,loc,x,1d,all,ux,0allsel ldread,reac,1,1,try02,rstsolvesfl,3,pres,1e5allsel solve/post1plnsol,u,y 實例3ANSYS中關于初應力和荷載步設置的算例算例為一個地基+地基上面的一塊方塊墻吧，先通過一次計算僅加邊界條件和自重，計算得到自重應力場，并輸出初應力文件用來模擬初應力場。然后施加應力和自重并進行計算，此時的位移基本為0，即消除了初位移，所得應力場即為自重應力場。在些基礎上進一步施加墻上法向面荷載，并進行第二步計算，得到的位移應該是僅有面荷載引起的位移。關于初應力想說明的幾點：1、初應力只能加在第一個荷載步，用命令流和GUI方式均可。但在求解前不能退出求解器，如果加了初應力后退出求解器到前處理器或者后處理器后再回到前處理器，剛才施加的初應力就沒有了，必須再次施加。2、也可以不用初應力而直接分兩個荷載步進行計算，如第一步僅計算自重，第二步再加面荷載后進行計算，在后處理中用工況組合來得到“凈位移”，但工況組合中的應力結果似乎是不正確的。3、如果用LSSOLVE從荷載文件進行求解，在寫荷載文件時初應力的設置并不會寫入荷載文件，所以，在命令流或者GUI方式下在求解前必須顯示指定加載初應力。關于荷載步設置的幾點建議：1、在GUI方式下，每次進入求解器進行求解似乎都是開始一個新的分析（這一點偶也不是很明白）。如果不退出求解器，即便不改變約束和荷載，只要求解一次，就會多一個荷載步結果，但所有結果是一樣的；如果退出求解器后再進來，求解就重新開始（根據時間值）。2、對點、線、面、體荷載都有替換和疊加兩種方式，在替換方式下，在同一位置重復加荷載，只有最后一次加的荷載有效；在疊加方式下，在同一位置重復加荷載，所有荷載會疊加后共同作用在結構上。/BATCHKEYW,PR_SET,1 KEYW,PR_STRUC,1KEYW,PR_THERM,0KEYW,PR_FLUID,0KEYW,PR_ELMAG,0KEYW,MAGNOD,0 KEYW,MAGEDG,0 KEYW,MAGHFE,0 KEYW,MAGELC,0 KEYW,PR_MULTI,0KEYW,PR_CFD,0 /GO!*/prep7!* define material protertiesmp,dens,1,2490mp,ex,1,20e9mp,prxy,1,.167mp,dens,2,2660mp,ex,2,10e9mp,prxy,2,.28!* define element typeet,1,solid45!* define geometryblock,-5,-0.5,-5,5,-5,0block,-0.5,0.5,-5,5,-5,0block,0.5,5,-5,5,-5,0block,-0.5,0.5,-5,5,0,5vglue,allvsel,s,loc,z,-6,0vatt,1,1,0vsel,invevatt,2,1,0!* mesh geometryESIZE,0,10lsel,s,loc,x,-0.4,0.4lesize,all,4allsel,allvmesh,allASEL,S,LOC,Z,-6,-0.1ASEL,U,LOC,X,-0.6,-0.4ASEL,U,LOC,X,0.4,0.6cm,uxyz,areaasel,s,loc,z,0,5asel,r,loc,x,0.4,0.6cm,sfload,areaALLSEL,ALLFINISH/SOLU !*ANTYPE,0!* the zero load step only for get initial stress filelsclear,allTIME,1AUTOTS,-1 NSUBST,1, , ,1KBC,0!* define constrainscmsel,s,uxyzda,all,all,0.!* deifine gravityacel,0.,0.,10.allsel,alliswrite,on !* open option for write initial stress filesolveFINISH/SOLU !*ANTYPE,0!* the first load steplsclear,allTIME,1AUTOTS,-1 NSUBST,1, , ,1KBC,0!* define constrainscmsel,s,uxyzda,all,all,0.!* deifine gravityacel,0.,0.,10.allsel,allISFILE,READ,AINIStress,ist, ,1!*初應力僅在第一個荷載步施加iswrite,offsolve!* the second load steplsclear,allTIME,2AUTOTS,-1 NSUBST,1, , ,1KBC,0!* define constrainscmsel,s,uxyzda,all,all,0.!* deifine gravityacel,0.,0.,10.!* define surface loadcmsel,s,sfloadsfa,all,1,pres,10000allsel,alliswrite,offsolve/prep7set,firstplnsol,u,z,0,1.0 !* 查看初位移應該基本為0plnsol,s,z,0,1.0 !* 可以查到初應力是存在的set,lastplnsol,u,z,0,1.0 !* 查看僅面荷載引起的位移plnsol,s,z,0,1.0 !* 查看應力結果(二)多工況任意組合實例對單層或二層框架進行彈性分析，需要考慮四種荷載 恒荷載，活荷載，風荷載和吊車荷載 1,幾何模型(beam3和beam54)建立后，定義所需的element table，主要包括桿端力和最大應力，最小應力等。 然后保存數據庫。分別施加四種荷載的標準值（不乘分項系數），并分別存成四個load step file。 2,使用solutionfrom ls files，求解四種荷載 3,荷載組合，命令流如下： /post1 lcdef,1,1 lcdef,2,2 lcdef,3,3 lcdef,4,4 !定義四種工況，分別為四種荷載下的計算結果 lcfact,1,1.2 lcfact,2,1.4 lcfact,3,1.19 lcfact,4,1.4 !指定各工況的組合系數 lcase,1 !讀入工況1，database1 sumtype,prin !指定加操作的對象 lcoper,add,2 !荷載組合，databasedatabase2 lcoper,add,4 !荷載組合，databasedatabase4 lcoper,lprin !計算線性主應力 lcwrite,11 !把database結果寫到工況11,即恒荷載活荷載吊車荷載的結果 lcase,1 lcfact,2,1.19 lcfact,4,1.19 !改變組合系數 sumtype,prin lcoper,add,2 lcoper,add,3 lcoper,add,4 lcoper,lprin lcwrite,12 !把database結果寫到工況12,即恒荷載活荷載吊車荷載風荷載的結果 !. .其他荷載組合 !之后使用lcase,n 就可調入工況n，并查看它的變形和內力 !可使用如下命令流得到工況11和12，13的較大者99,進而查看最大應力 lcase,11 lcase,min,12 lcase,min,13 lcwrite,98 lcase 98 !查看工況98的應力分布. . lcase,11 lcase,max,12 lcase,max,13 lcwrite,99 lcase 99 !查看工況99