一致缺陷模态法.docx

《一致缺陷模态法.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、!一致缺陷模态法!finish$/clear$js=1$*use,g:static_stabiiity_analysis(wq).mac!!长度单位为m，力单位为N，质量单位为kg!!!!定义参数/config,nres,99999!!结果集或者子步的最大数目，AssignsvaluestoANSYSconfigurationparameters.k='lh'!!凯威特几型网壳q='qk'a='ds's=101!!跨度f=3!!矢高比倒数js_stop=2d=300!!初始几何缺陷取跨度的1/300、1/400、1/500……elme_div1=1!!杆划分的单元个数elme_div2

2、=1!!杆划分的单元个数nb=30!!双重非线性时，施加的nb倍荷载ibmax=3!!双重非线性时，弧长半径因子上限值ibmin=0.00001!!双重非线性时，弧长半径因子下限值nbstepsub=100!!双重非线性时，荷载子步数uzb=0.5!!Z向位移绝对值到达uzb时，停止计算!!改变工作目录、文件名、图标名/cwd,g:N_Eigenvalue_Results!!改变工作目录/filename,%k%-%s%-%f%!!工作名称%k%-%s%-%f%/title,%k%-%s%-%f%/prep7*use,g:N_Eigenvaluemacrocolor.mac!==

3、============================================================================!!建立节点*use,g:N_EigenvaluemacroWQ%k%-%s%-%f%orinodecoodin.mac*USE,g:N_EigenvaluemacroWQ%k%-%s%-%f%real.mac*get,kp_num,kp,0,count!返回所有关键点总数的到kp_num这个变量。NKPT,,all!将所有的关键点定义为节点*USE,g:N_EigenvaluemacroWQ%k%-%s%-%

4、f%Creat_Lines.mac！创建线并且，直接给先定义属性，则线就具备了属性。allsel！！选择所有的单元以及节点lsel,s,line,,all!!将所有的单元选出来作为一个新的选择集，此行原本是用real来选择的lesize,all,,,elme_div1!!将一根线划分为elme_div1个单元，线的单元尺寸定义lmesh,all!!网格划分，在几何线上生成线单元nummrg,node!!重合节点合并allsel!边界条件施加*use,g:N_EigenvaluemacroWQ%k%-%s%-%f%boundaryconditions.mac!06节点荷载及重力

5、加速度施加FCUM,ADD!!以累加的形式施加荷载*use,g:N_EigenvaluemacroWQ%k%-%s%-%f%equiforce_%q%.mac,1.0!!施加工况acel,,,9.8!!重力加速度eplotallsel!定义静力求解选项/solu!!进入求解层pstres,on!!打开预应力开关，必须在第一个荷载步内指定nsubst,1!!每一个荷载步分为一个子步solvefinish/post1pldisp,1!显示静力求解之后的结构变形的形状finish!定义特征值屈曲求解选项/solu!!进入求解层antype,buckle!!定义求解类型为特征值屈曲分析

6、bucopt,lanb,20!!定义特征值提取方法为lanb，提取特征值数为20outres,all,all!!定义输出全部荷载步结果这一步是定义向数据库及结果文件中写入的数据，solvefinish!施加缺陷/post1set,listset,1,js!!循环读取每阶模态，寻找各阶模态的最大位移pldispALLSELNSORT,U,SUM,0,0,ALL!!将节点结果按照合位移顺序排列*GET,MAX_U,SORT,0,IMAX!!取MAX_U=节点数据中最大值对应的节点编号*GET,USUM,NODE,MAX_U,U,SUM!!USUM_%b%=编号为MAX_U的节点的合位移(寻找

7、各阶模态最大和位移，为施加缺陷做准备)，强制置换《参数化》P26/PREP7allsellsel,s,real,,1,3ksll,snslk,supgeom,s/d/USUM,1,1,%k%-%s%-%f%,rst!!模态阶数与大小。施加缺陷(《参数化》P174)ALLSEL!获取缺陷荷载坐标（此时生成关键点）（统一）*use,g:N_Eigenvaluemacro8_GetDefectCoodinates.mac!将变量保