钢管受到冲击屈曲分析

一个钢管受到外力的撞击而发生变形，假定每次的撞击力大小相等，变形已知，但是撞击次数未知，现在反求撞击力，不知道是否可行？
1、在次数已知的时候（可以假定为10次），如何操作
2、如果次数未知，是否力的大小可以为任意值了？
下面是基本命令流，不过没有实现我的想法，请大家帮忙指正一下


!*********************************************************************************
!                        薄壁钢管的冲击屈曲分析
!******************************************************************************
!分析环境设置
/filname,pipe-buckling
/title,analysis of pipe-buckling
!定义单元类型及算法选项，实参数
/prep7
et,1,shell163
keyopt,1,1,7
r,1,5/6,5,8e-3,,,,
!定义材料模型
mp,dens,1,7800
mp,ex,1,2.07e11
mp,nuxy,1,0.27
tb,biso,1,,,,
tbdat,1,3.0e8
tbdat,2,1.0e10
edmp,hgls,1,5,0.1,1.5,0.06
!建立几何模型
k,1,0,0,0
k,2,0.5864,0,1.3808
lstr,1,2
wprot,0,0,23
FLST,2,2,8  
FITEM,2,0,0,0   
FITEM,2,0.305,0,0  
CIRCLE,P51X, , , ,360, ,
adrag,2,3,4,5,,,1
CYL5,0.305,0,-0.305,0  
wpoff,0,0,1.5
CYL5,0.305,0,-0.305,0
wpoff,0,0,-1.5
NUMMRG,ALL, , , ,LOW
wprot,0,0,-23
!网格的划分
aatt,1,1,1,0,
esize,0.05
amesh,all
                                   
               
!施加位移约束
LSEL,S, , ,       7
nsll,s,1
d,all,ux
csys,2
nsel,s,loc,x,0.30499,0.30501
csys,0
nsel,u,loc,x,-0.306,0.0001
d,all,all
csys,2
nsel,s,loc,x,1.531,1.5305
csys,0
nsel,u,loc,x,0.5866,0.2
d,all,all
allsel
cskp,100,0,1,3,4                     
csys,100   
!定义冲击载荷数组或曲线
nsel,s,loc,z,0.01,0.3
nsel,r,loc,y,-0.1,0.1
nsel,r,loc,x,-0.16,-0.4
cm,load,nodes,
*dim,time,array,2,1,1,,,
time(2,1,1)=0.01
*dim,compress,array,2,1,1,,,
compress(2,1,1)=-50000
!施加冲击载荷
!以载荷时间数组形式施加显式分析荷载
edload,add,fx,0,load,time,compress,0,,,,,
allsel,all
edmp,hgls,1.5,0.1,1.5,0.06,,,
!定义表面接触信息
!定义钢管表面的单面自动接触
EDCGEN,ASSC, , ,0.1,0.1,0,0,0, , , , ,0,10000000
fini
!进入求解器
/solu
!设置求解时间
time,0.005,
!结果文件输出类型
edopt,add,blank,both
!结果文件输出步数的设置
edrst,100
edhtime,100
!设置壳单元厚度方向积分点输出个数
edint,5,0,
!设置能量控制选项
edopt,add,blank,both
!递交求解
solve
fini

                               
登录/注册后可看大图