问题如下：将实测的加速度数据输入到模型中，输入的加速度时程曲线数据不大，为10-2m/s2左右的数，可是计算完的输出加速度结果却异常的大，不知是什么方面出了问题，得 从几个方面考虑，对模型进行调整呢？下面附上命令流，和计算图，求高人指点迷津！谢谢啦！
finish/clear/PREP7*SET,dxyz,0.3          !网格尺寸*SET,KbN,22517000*SET,CbN,2430000et,1,45                  !基础mp,ex,1,3.3e10  mp,prxy,1,0.2mp,dens,1,2500k, ,dxyz,10*dxyz,k, ,63*dxyz,10*dxyz,k, ,63*dxyz,20*dxyz,k, ,dxyz,20*dxyz,a,1,2,3,4,vext,1,,,,,4*dxyz    !中间实体1k, ,dxyz,12*dxyz,4*dxyzk, ,63*dxyz,12*dxyz,4*dxyzk, ,7*dxyz,13*dxyz,4*dxyzk, ,57*dxyz,13*dxyz,4*dxyzk, ,57*dxyz,14*dxyz,4*dxyzk, ,7*dxyz,14*dxyz,4*dxyzk, ,7*dxyz,17*dxyz,4*dxyzk, ,57*dxyz,17*dxyz,4*dxyzk, ,57*dxyz,18*dxyz,4*dxyzk, ,7*dxyz,18*dxyz,4*dxyza,7,8,9,10a,11,12,13,14a,15,16,17,18asba,7,8asba,10,9vext,7,,,,,dxyz       !开槽后的落在实体1上的实体2LESIZE,ALL,dxyz, , , ,1, , ,1,  VMESH,1 VSWEEP,2  allsel,all            !实体1和实体2划分网格k, ,0,9*dxyz,k, ,dxyz,9*dxyz,k, ,dxyz,22*dxyz,k, ,0,22*dxyz,a,31,32,33,34AGEN,2,21, , ,63*dxyz, , , ,0   a,4,3,38,33a,1,32,35,2vext,21,24,1,,,7*dxyz      !周围实体k, ,10*dxyz,,k, ,53*dxyz,,k, ,53*dxyz,9*dxyz,k, ,10*dxyz,9*dxyz,a,51,52,53,54vext,45,,,,,7*dxyz       !下方的实体LESIZE,ALL,dxyz, , , ,1, , ,1,  MSHAPE,0,3d MSHKEY,1VMESH,3,7,1allsel,all                  !划分单元nummrg,allnumcmp,all                  !采用合并节点方法连接各体 !*******************************************************************************!以下建立底边界法向和切向弹簧阻尼单元et,2, combin14,          !法向r,2, KbN, CbN*do,iq,1,9nsel,s,loc,z,0  !选中基础底面所有节点nsel,r,loc,y,(iq-1)*dxyz  !y=0所有节点*get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node((ip+9)*dxyz,(iq-1)*dxyz,0)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x,y,z-dxyz/10    !定义底边界法向结点以便与边界点形成法向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************下底面粘弹性边界建成*do,iq,10,23nsel,s,loc,z,0  !选中基础底面所有节点nsel,r,loc,y,(iq-1)*dxyz  !y=0所有节点*get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node((ip-1)*dxyz+0,(iq-1)*dxyz,0)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x,y,z-dxyz/10    !定义底边界法向结点以便与边界点形成法向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************上底面粘弹性边界建成*do,iq,1,8nsel,s,loc,x,0  !选中基础底面所有节点nsel,r,loc,z,(iq-1)*dxyz  !y=0所有节点*get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node(0,(ip-1)*dxyz+2.7,(iq-1)*dxyz)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x-dxyz/10,y,z    !定义底边界切向结点以便与边界点形成切向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************左1侧面粘弹性边界建成*do,iq,1,8,nsel,s,loc,x,64*dxyz  nsel,r,loc,z,(iq-1)*dxyz  *get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node(64*dxyz,(ip-1)*dxyz+2.7,(iq-1)*dxyz)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x+dxyz/10,y,z    !定义底边界切向结点以便与边界点形成切向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************左4侧面粘弹性边界建成*do,iq,1,8nsel,s,loc,x,10*dxyz  !选中基础底面所有节点nsel,r,loc,y,0,9*dxyznsel,r,loc,z,(iq-1)*dxyz  !y=0所有节点*get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node(10*dxyz,(ip-1)*dxyz+0,(iq-1)*dxyz)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x-dxyz/10,y,z    !定义底边界切向结点以便与边界点形成切向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************左2侧面粘弹性边界建成*do,iq,1,8nsel,s,loc,x,53*dxyz  !选中基础底面所有节点nsel,r,loc,y,0,9*dxyznsel,r,loc,z,(iq-1)*dxyz  !y=0所有节点*get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node(53*dxyz,(ip-1)*dxyz+0,(iq-1)*dxyz)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x+dxyz/10,y,z    !定义底边界切向结点以便与边界点形成切向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************左3侧面粘弹性边界建成*do,iq,1,8nsel,s,loc,y,0  !选中基础底面所有节点nsel,r,loc,z,(iq-1)*dxyz  !y=0所有节点*get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node((ip+9)*dxyz,0,(iq-1)*dxyz)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x,y-dxyz/10,z    !定义底边界切向结点以便与边界点形成切向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************前1面粘弹性边界建成*do,iq,1,8nsel,s,loc,y,9*dxyz  !选中基础底面所有节点nsel,r,loc,x,0,10*dxyznsel,r,loc,z,(iq-1)*dxyz  !y=0所有节点*get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node((ip-1)*dxyz,9*dxyz,(iq-1)*dxyz)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x,y-dxyz/10,z    !定义底边界切向结点以便与边界点形成切向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************前2面左粘弹性边界建成*do,iq,1,8nsel,s,loc,y,9*dxyz  !选中基础底面所有节点nsel,r,loc,x,53*dxyz,63*dxyznsel,r,loc,z,(iq-1)*dxyz  !y=0所有节点*get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node((ip-1)*dxyz+16.2,2.7,(iq-1)*dxyz)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x,y-dxyz/10,z    !定义底边界切向结点以便与边界点形成切向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************前2面右粘弹性边界建成*do,iq,1,8nsel,s,loc,y,22*dxyz  !选中基础底面所有节点nsel,r,loc,z,(iq-1)*dxyz  !y=0所有节点*get,np,node,,count      !得到选中的结点数，存入np*get,npmax,node,,num,maxd    !得到已经定义的最大结点数，存入npmax*do,ip,1,npnpnum=node((ip-1)*dxyz,22*dxyz,(iq-1)*dxyz)x=nx(npnum)y=ny(npnum)z=nz(npnum)npmax=npmax+1n,npmax,x,y+dxyz/10,z    !定义底边界切向结点以便与边界点形成切向单元type,2real,2e,npnum,npmaxd,npmax,all,0       !约束新生成的点*enddo*enddoallsel,all!***********************************************************************前3面粘弹性边界建成/SOLANTYPE,2MODOPT,LANB,15  EQSLV,SPAR  MXPAND,15, , ,0 LUMPM,0 PSTRES,0MODOPT,LANB,15,0,0, ,OFF  solve/SOLU    ANTYPE,2 !设置分析类型为模态分析MXPAND,15,0,0,1,0,  !设置扩展模态阶数为15阶SOLVEFINISH/SOLU    ANTYPE,8 !设置分析类型为谱分析SPOPT,PSD,15,1  !设置幅的激励为基础激励!定义激励方向/SOL!*  ALPHAD,0,  BETAD,0, DMPRAT,0.07,   MDAMP,1, , , , , , ,nsel,s,loc,y,9*dxyz-dxyz/10nsel,a,loc,y,-dxyz/10D,all,UY,1.0 allsel,all!设置频率与谱值关系PSDFRQ,1, ,1,3,3.5,4,5,5.1,5.5,5.9,6.033333,6.666667PSDFRQ,1, ,8,10.1,10.7,11,13,15,18,20,20.3,20.5PSDFRQ,1, ,20.9,22,23.433333,23.733333,24.2,24.7,26,27.3,28,29PSDFRQ,1, ,30,31,31.966667,33,34,36,38,41,43,45,46,48,50!*  PSDVAL,1,1.3549E-05,4.61063E-05,0.000401264,0.000128004,0.000363115,0.000564224,0.000357226,0.000308852,0.001051814PSDVAL,1,0.002002561,0.000636296,0.01042194,0.001984788,0.000551114,0.000131586,0.000243743,0.000148657,0.000279433,0.00278082PSDVAL,1,0.001451292,0.000462117,0.000153696,5.50298E-05,0.003670355,0.000202857,0.000234703,0.000138016,0.000978275,0.000240522PSDVAL,1,0.000354334,8.37025E-05,0.000120601,0.00046909,0.000184796,0.001429391,0.000147629,0.00004309,1.96368E-05,4.58687E-05,1.4166E-05,2.13832E-05,1.2542E-05,5.94052E-06OUTRES,ALL,ALL,  !设置输出所有子步PSDUNIT,1,acel,386.4, !定义谱类型为位移谱PFACT,1,BASE, PSDRES,DISP,ABS!定义输出控制PSDRES,VELO,ABSPSDRES,ACEL,ABS  solvefinish/SOLU   ANTYPE,8 !设置分析类型为谱分析PSDCOM,0.002,15, !设置合并模态的阀值和阶数SOLVE finish*DIM,acex,ARRAY,6000,1,1, , , *CREATE,ansuitmp*VREAD,acex(1,1),'ACEX','txt',' ',JIK,1,6000,1,0,   (pe13.7)*END/INPUT,ansuitmpfinish/config,nres,20000/soluantype,transtrnopt,fullACEL,0,0,-9.81, alphad,0.099betad,0.025*do,i,1,6000p=12251.7*acex(i)   !质量乘以加速度得到作用于模型迎波面的面力nsel,s,loc,y,9*dxyz-dxyz/10nsel,a,loc,y,-dxyz/10F,all,FY,pallsel,alltime,i*0.005outres,all,allsolve*enddofinish

 返回小木虫查看更多