马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
就做一悬臂梁的动载荷加载,看相关结点的加速度响应时程,现在看来总是出现抖动的现象(见图),请教如和消除、、、、
finish
/clear
/filnam,try2
/title,try2
/prep7
block,0,0.2,0,0.4,0,10
et,1,solid185 !设置分析单元
mp,ex,1,2.06e5 !设置材料属性
mp,prxy,1,0.3 !设置泊松比
mp,dens,1,7.8e-6 !设置密度
!划分网格
mshape,0,3D !选择3D六面体网格
mshkey,1 !选择映射划分
esize,0.2,0.2,0.2
vmesh,all
!nlist
!低端四边加约束
dl,1,,all
dl,2,,all
dl,3,,all
dl,4,,all
!加载重力加速度
timent,off
nlgeom,on
time,1E-8
nsubst,2
kbc,1
acel,0,0,9.8 !加载重力加速度
solve
timent,on
nlgeom,off
nt=150
*DIM,acx_1,,nt,1 !定义数组
*DO,T,1,50,1
acx_1(T,1)=T
*ENDDO
*DO,T,51,100,1
acx_1(T,1)=50
*ENDDO
*DO,T,101,150,1
acx_1(T,1)=150-T
*ENDDO
/INPUT,ansuitmp
!****************求解****************
!设置求解选项
FINI
/SOLU
ANTYPE,4 !指定分析类型为瞬态动力学分析
TRNOPT,FULL !瞬态动力学分析采用FULL法
OUTRES,ALL,All !不输出任何结果、、只有(all,all)才能输出加速度
tim=0
*DO,T,1,nt,1 !循环读入载荷数据并求解
TIME,T*0.5 !设置时间步
KBC,0 !指定载荷为递增载荷
NSUB,15 !设定子步数为
ALPHAD,0.05 !设定质量阻尼ALPHA为0.05
BETAD,0.01 !设定刚度阻尼BETA为0.01
f,12,FX,*acx_1(T,1)
ALLS !选中所有元素
OUTRES,NSOL,LAST
SOLVE !求解
*ENDDO
SAVE !保存
!!!!!!!后处理观察位移响应曲线
/post26
nsol,2,282,u,x,ux_12 !定义节点的输出,显示位移响应为U,速度响应改为V,加速度响应改为A
nsol,3,282,a,x,ax_12
prvar,2,3
/axlab,X,time !设置图形坐标轴标志
/axlab,y,ux_12
plvar,2,3 !显示变化曲线
finish !退出时间历程后处理器
[ 本帖最后由 fkts2007 于 2010-4-18 16:43 编辑 ] |