一致缺陷模态法.docxVIP

!一致缺陷模态法 ! finish$/clear $js=1 $*use,g:\static_stabiiity_analysis(wq).mac !!长度 单位为 m ，力单位为 N，质量单位为 kg !! !!定义参数 /config,nres,99999 !!结果集或者子步的最大数目，Assigns values to ANSYS configuration parameters. k=lh !!凯威特几型网壳 q=qk a=ds s=101 !!跨度 f=3 !!矢高比倒数 js_stop=2 d=300 !!初始几何缺陷取跨度的1/300、1/400、1/500…… elme_div1=1 !!杆划分的单元个数 elme_div2=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_Eigenvalue\macro\color.mac !==============================================================================! !建立节点 *use,g:\N_Eigenvalue\macro\WQ\%k%-%s%-%f%\orinodecoodin.mac *USE,g:\N_Eigenvalue\macro\WQ\%k%-%s%-%f%\real.mac *get,kp_num,kp,0,count !返回所有关键点总数的到kp_num这个变量。 NKPT,,all !将所有的关键点定义为节点 *USE,g:\N_Eigenvalue\macro\WQ\%k%-%s%-%f%\Creat_Lines.mac ！创建线并且，直接给先定义属性，则线就具备了属性。 allsel ！！选择所有的单元以及节点 lsel,s,line,,all !!将所有的单元选出来作为一个新的选择集，此行原本是用real来选择的 lesize,all,,,elme_div1 !!将一根线划分为elme_div1个单元，线的单元尺寸定义 lmesh,all !!网格划分，在几何线上生成线单元 nummrg,node !!重合节点合并 allsel !边界条件施加 *use,g:\N_Eigenvalue\macro\WQ\%k%-%s%-%f%\boundaryconditions.mac ! 06 节点荷载及重力加速度施加 FCUM,ADD !!以累加的形式施加荷载 *use,g:\N_Eigenvalue\macro\WQ\%k%-%s%-%f%\equiforce_%q%.mac,1.0 !!施加工况 acel,,,9.8 !!重力加速度 eplot allsel !定义静力求解选项 /solu !!进入求解层 pstres,on !!打开预应力开关，必须在第一个荷载步内指定 nsubst,1 !!每一个荷载步分为一个子步 solve finish /post1 pldisp,1 !显示静力求解之后的结构变形的形状 finish !定义特征值屈