两端简支 ！求解索力与第一阶自振频率关系命令流 ！单位：kg-m-s FINISH /CLEAR *DIM,FREQL,ARRAY,16 *DIM,SUOLI,ARRAY,16 /PREP7 PI=ACOS(-1) LL=150 EA=0.0026 EE=1.95E11 EI=EA*EA/4/PI FORCE=600000 ISTRAIN=FORCE/EA/EE ET,1,BEAM3 R,1,EA,EI,0.1, ,ISTRAIN MP,EX,1,EE MP,PRXY,1,0 MP,DENS,1,8250 *AFUM,DEG K,1 K,2,LL*COS(45),LL*SIN(45) L,1,2 LATT,1,1,1 LESIZE,1,1,5 LMESH,ALL FINISH /SOL D,1,UX,0, , , ,UY D,2,UY FINISH *D0,I,1,16 /PREP7 T0=400000+I*200000 NMODIF,2, , , ,45 F,2,FX,T0 ALLSEL,ALL ACEL,0,10,0 ISTRAIN=T0/EA/EE !定义自振频率的数组 !定义索力的数组 !进入前处理模块 !定义圆周率大小 !定义索长 !定义截面积 !定义弹性模量 !定义截面惯性矩 !定义初始索力 !定义初应变 !定义单元类型 !定义单元实常数 !定义材料常数 !定义关键点 !连接关键点 !赋予单元特征值 !网络划分尺寸设置 !进行划分网格 !前处理结束 !进入求解模块 !定义边界条件 !求解模块结束 !进入循环求解 !索力增量设置 !旋转端点的局部坐标体系 !在节点 2 施加索力 !选择所有项目 !定义重力加速度 !定义初应变 

!修改所有单元的实常数为 I !前处理结束 !进行静力分析 !计入初应力 !求解 !求解模块结束 !进行模态分析 !定义模态分析方法、模态书目 !计入初应力 !求解 !求解模块结束 !进入通用后处理器 !将各索力输入索数组 !得到自振频率 !循环结束 !绘制一阶振型图 !画出力和频率的曲线 !列表显示频率值 R,I,EA,EI,0. 1,0,ISTRAIN,0 EMODIF,ALL,REAL,I FINISH /SOL ANTYPE,0 PSTRES,ON SOLVE FINISH /SOL ANTYPE,2 MODOPT,SUBSP,1 PSTRES,ON SOLVE FINISH /POST1 SUOLI(I)=400000+200000*I *GET,FREQL(I),MODE,1,FREQL *ENDDO PLDISP,1 *VPLOT,SUOLI(1),FREQL(1) *STATUS,FREQL FINISH 两端固支 finish /clear /PREP7 D=0.15 L=5 N=10000 PI=ACOS(-1) A0=PI*D*D/4 I0=PI*D**4/64 TM=1.95E11 BSB=0.3 MD=7920 PZ=1.2E-4 YINGLI=N*1000/A0 YINGBIAN=YINGLI/TM WENCHA=-YINGBIAN/PZ 

!泊松比 !密度 !热膨胀系数 !面积、IZZ、IYY JDSL=101 ET,1,BEAM3 MP,EX,1,TM !弹性模量 MP,NUXY,1,BSB MP,DENS,1,MD MP,ALPX,1,PZ R,1,A0,I0,I0 N,1,0,0,0 N,JDSL,L,0,0 FILL,1,JDSL,JDSL-2, TYPE,1 MAT,1 REAL,1 *DO,I,1,JDSL-1 E,I,I+1 *ENDDO TREF,0 BFE,ALL,TEMP,,WENCHA,WENCHA,WENCHA,WENCHA finish /SOL ANTYPE,STATIC D,1,ALL D,JDSL,all ! D,11,,,,,,uy ! D,89,,,,,,uy PSTRES,ON SOLVE FINISH /SOL ANTYPE,MODAL MODOPT,LANB,13 MXPAND,13 PSTRES,ON SOLVE FINISH /post1 pldisp etable,mforce,smisc,1 etable,mstress,ls,1 plls,mforce,mforce,1 plls,mstress,mstress,1 prrsol finish /POST1 

SET, LIST SET, FIRST PLDI ! SET, NEXT ! PLDI FINISH