資源描述：
《一致缺陷模態(tài)法.docx》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、!一致缺陷模態(tài)法!finish$/clear$js=1$*use,g:static_stabiiity_analysis(wq).mac!!長(zhǎng)度單位為m，力單位為N，質(zhì)量單位為kg!!!!定義參數(shù)/config,nres,99999!!結(jié)果集或者子步的最大數(shù)目，AssignsvaluestoANSYSconfigurationparameters.k='lh'!!凱威特幾型網(wǎng)殼q='qk'a='ds's=101!!跨度f(wàn)=3!!矢高比倒數(shù)js_stop=2d=300!!初始幾何缺陷取跨度的1/300、1/400、1/500……elme_div1=1!!桿劃分的單元個(gè)數(shù)elme_div2
2、=1!!桿劃分的單元個(gè)數(shù)nb=30!!雙重非線性時(shí)，施加的nb倍荷載ibmax=3!!雙重非線性時(shí)，弧長(zhǎng)半徑因子上限值ibmin=0.00001!!雙重非線性時(shí)，弧長(zhǎng)半徑因子下限值nbstepsub=100!!雙重非線性時(shí)，荷載子步數(shù)uzb=0.5!!Z向位移絕對(duì)值到達(dá)uzb時(shí)，停止計(jì)算!!改變工作目錄、文件名、圖標(biāo)名/cwd,g:N_Eigenvalue_Results!!改變工作目錄/filename,%k%-%s%-%f%!!工作名稱(chēng)%k%-%s%-%f%/title,%k%-%s%-%f%/prep7*use,g:N_Eigenvaluemacrocolor.mac!==
3、============================================================================!!建立節(jié)點(diǎn)*use,g:N_EigenvaluemacroWQ%k%-%s%-%f%orinodecoodin.mac*USE,g:N_EigenvaluemacroWQ%k%-%s%-%f%real.mac*get,kp_num,kp,0,count!返回所有關(guān)鍵點(diǎn)總數(shù)的到kp_num這個(gè)變量。NKPT,,all!將所有的關(guān)鍵點(diǎn)定義為節(jié)點(diǎn)*USE,g:N_EigenvaluemacroWQ%k%-%s%-%
4、f%Creat_Lines.mac！創(chuàng)建線并且，直接給先定義屬性，則線就具備了屬性。allsel！！選擇所有的單元以及節(jié)點(diǎn)lsel,s,line,,all!!將所有的單元選出來(lái)作為一個(gè)新的選擇集，此行原本是用real來(lái)選擇的lesize,all,,,elme_div1!!將一根線劃分為elme_div1個(gè)單元，線的單元尺寸定義lmesh,all!!網(wǎng)格劃分，在幾何線上生成線單元nummrg,node!!重合節(jié)點(diǎn)合并allsel!邊界條件施加*use,g:N_EigenvaluemacroWQ%k%-%s%-%f%boundaryconditions.mac!06節(jié)點(diǎn)荷載及重力
5、加速度施加FCUM,ADD!!以累加的形式施加荷載*use,g:N_EigenvaluemacroWQ%k%-%s%-%f%equiforce_%q%.mac,1.0!!施加工況acel,,,9.8!!重力加速度eplotallsel!定義靜力求解選項(xiàng)/solu!!進(jìn)入求解層pstres,on!!打開(kāi)預(yù)應(yīng)力開(kāi)關(guān)，必須在第一個(gè)荷載步內(nèi)指定nsubst,1!!每一個(gè)荷載步分為一個(gè)子步solvefinish/post1pldisp,1!顯示靜力求解之后的結(jié)構(gòu)變形的形狀finish!定義特征值屈曲求解選項(xiàng)/solu!!進(jìn)入求解層antype,buckle!!定義求解類(lèi)型為特征值屈曲分析
6、bucopt,lanb,20!!定義特征值提取方法為lanb，提取特征值數(shù)為20outres,all,all!!定義輸出全部荷載步結(jié)果這一步是定義向數(shù)據(jù)庫(kù)及結(jié)果文件中寫(xiě)入的數(shù)據(jù)，solvefinish!施加缺陷/post1set,listset,1,js!!循環(huán)讀取每階模態(tài)，尋找各階模態(tài)的最大位移pldispALLSELNSORT,U,SUM,0,0,ALL!!將節(jié)點(diǎn)結(jié)果按照合位移順序排列*GET,MAX_U,SORT,0,IMAX!!取MAX_U=節(jié)點(diǎn)數(shù)據(jù)中最大值對(duì)應(yīng)的節(jié)點(diǎn)編號(hào)*GET,USUM,NODE,MAX_U,U,SUM!!USUM_%b%=編號(hào)為MAX_U的節(jié)點(diǎn)的合位移(尋找
7、各階模態(tài)最大和位移，為施加缺陷做準(zhǔn)備)，強(qiáng)制置換《參數(shù)化》P26/PREP7allsellsel,s,real,,1,3ksll,snslk,supgeom,s/d/USUM,1,1,%k%-%s%-%f%,rst!!模態(tài)階數(shù)與大小。施加缺陷(《參數(shù)化》P174)ALLSEL!獲取缺陷荷載坐標(biāo)（此時(shí)生成關(guān)鍵點(diǎn)）（統(tǒng)一）*use,g:N_Eigenvaluemacro
