用ANSYS进行滞回分析

用ANSYS进行滞回分析/PREP7!定义单元类型，实常数，材料特性ET,1,SHELL143R,1,12,,,,,MP,EX,1,196784MP,NUXY,1,0.3!双线性随动强化模型TB,BKIN,1,1,2,1TBDATA,,310,600,,,,!定义关键点、线、面K,1,54,0,0K,2,-54,0,0K,3,54,0,1000K,4,-54,0,1000A,1,2,4,3!定义边界荷强迫位移,划分网格AESIZE,ALL,27,MSHAPE,0,2DMSHKEY,0CM,_Y,AREAASEL,,,,1CM,_Y1,AREACMSEL,S,_YAMESH,_Y1*do,i,1,5D,i,ALL,0*enddoOUTPR,BASIC,ALL,OUTRES,ALL,ALL,D,46,ux,30TIME,1AUTOTS,0NSUBST,10,,,1KBC,0LSWRITE,01,!第2荷载步D,46,ux,-30TIME,3AUTOTS,0NSUBST,20,,,1KBC,0LSWRITE,02,!第3荷载步D,46,ux,30TIME,5AUTOTS,0NSUBST,20,,,1KBC,0LSWRITE,03,!第4荷载步D,46,ux,-30TIME,7AUTOTS,0NSUBST,20,,,1KBC,0LSWRITE,04,!第1荷载步D,46,ux,40TIME,1AUTOTS,0NSUBST,10,,,1KBC,0LSWRITE,05,!第2荷载步D,46,ux,-40TIME,3AUTOTS,0NSUBST,20,,,1KBC,0LSWRITE,06,!第3荷载步D,46,ux,40TIME,5AUTOTS,0NSUBST,20,,,1KBC,0LSWRITE,07,!第4荷载步D,46,ux,-40TIME,7AUTOTS,0NSUBST,20,,,1KBC,0LSWRITE,08,!求解FINISH/SOLULSSOLVE,1,8,1,!画出荷载位移曲线FINISH/POST26NSOL,2,46,U,X,RFORCE,3,46,F,X,XVAR,2PLVAR,3,,,,,,,,,,回复：【分享】用ANSYS进行滞回分析请教这位高手，本人也要做一个滞回分析，是个软钢圆柱，而我采用的是实体建模，采用SOLID45单元，双线性随动强化，可结果在大位移的情况就是出现蝶形曲线，而实验的情况则是出现一个梭形，跟你所画的图一样，我换了其他可用的SOLID单元，但结果还是一样。希望你能给我提些建议。能给个联系方式吗？以下是本人的一个命令流/prep7et,1,solid45mp,ex,1,2.01e5mp,prxy,1,0.26TB,BKIN,1,1,2TBTEMP,0TBDATA,,185,0,!建立单元及划分网格block,,15,,15,,200esize,5vmesh,all!施加底部约束nsel,s,loc,z,0d,all,all,0ALLSEL,ALL!定义加载的位移数组*dim,disp,,6wlw1=80disp(1)=0disp(2)=wlw1disp(3)=0disp(4)=-wlw1disp(5)=0disp(6)=wlw1!进入solution阶段/solunlgeom,onsstif,offautots,onoutres,all,alloutpr,all,all!施加位移荷载time,0.0001nsel,s,loc,z,200nsubst,1,0,0d,all,ux,disp(1)ALLSEL,ALLsolve*do,i,2,6time,insel,s,loc,z,200d,all,ux,disp(i)nsubst,40,0,0allsel,allsolve*enddofinish那么怎么提取滞回数据呢/PREP7ET,1,BEAM3ET,2,COMBIN14KEYOPT,2,1,0KEYOPT,2,2,0KEYOPT,2,3,2R,1,0.16,0.00213333,0.4,,,,R,2,0.18,0.0054,0.6,0,0,600,R,3,,5000000,,!阻尼器线性系数C1MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,3e10MPDATA,PRXY,1,,0.2MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,2500K,1,,,,K,2,6,,,K,3,,6,,K,4,6,6,,KPLOTLSTR,1,3LSTR,3,4LSTR,2,4LSEL,s,LINE,,1,3,2,LATT,1,1,1,,,,LSEL,S,LINE,,2LATT,1,2,1,,,,LSEL,all,LESIZE,ALL,1,,,,,,,1LMESH,ALL,/SHRINK,0/ESHAPE,1.0/EFACET,1/RATIO,1,1,1/REPLOTTYPE,2MAT,1REAL,3ESYS,0E,2,14D,1,all,,,14,13,FINISH*SET,NT,1001*SET,DT,0.02*DIM,AC,,NT*VREAD,AC(1),RECORD,TXT(F8.3)/SOLU!模态分析ANTYPE,2MODOPT,SUBSP,8MXPAND,8,,,1SOLVEFINI!得到自振频率1*GET,FREQ1,MODE,1,FREQ/CONFIG,NRES,20000/SOLUANTYPE,TRANSTRNOPT,FULLALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)*DO,I,1,500ACEL,AC(I),0,0TIME,I*0.02OUTRES,ALL,ALLSOLVE*ENDDOFINISH单柱滞回曲线问题（命令流＋图）建立一个单柱模型，进行位移加载。分别采用了随动强化和等向强化两种强化准则。材料的应力应变曲线（见命令流中）为三折线，均有明显的下降段，但是计算后的柱顶位移－柱底剪力滞回曲线上没有发现结构有明显的刚度退化现象，反而呈现一种理想弹塑性的滞回曲线样式，不得其解！命令流如下：fini/clear/prep7n,1,n,16,1.5n,17,0,1000fill,1,16et,1,beam188mp,ex,1,3E10mp,nuxy,1,0.167mp,dens,1,0!随动强化TB,KINH,1,1,3,PLASTICTBTEMP,0TBPT,,6.6e-3,37.23e6TBPT,,0.034,28.034e6TBPT,,0.051,0!等向强化!TB,MISO,1,1,3!TBTEMP,0.0!TBPT,DEFI,0.001,3E7!TBPT,DEFI,6.6e-3,37.23e6!TBPT,DEFI,0.034,28.034e6r,1,SECTYPE,1,BEAM,RECT,pier,0SECOFFSET,CENTSECDATA,0.25,0.25,10,10,0,0,0,0,0,0type,1mat,1real,1*do,ii,1,15e,ii,ii+1,17*enddod,1,all/soluantype,staticnropt,fulloutpr,all,alloutres,all,all*do,tt,1,20,2time,ttnsubst,10,,d,16,,tt*0.01,,,,uy!lswrite,ttsolvetime,tt+1nsubst,10,,d,16,,-1*tt*0.01,,,,uy!lswrite,tt+1solve*enddo!lssolve,1,20,1savefini采用随动强化时的滞回曲线如下图：先把命令流贴一下：/PREP7K,,0,0,0,K,,0,10,0,K,,60,0,0,K,,60,10,0,FLST,2,4,3FITEM,2,2FITEM,2,1FITEM,2,3FITEM,2,4A,P51XFLST,2,1,5,ORDE,1FITEM,2,1VEXT,P51X,,,0,0,3,,,,/VIEW,1,1,1,1/ANG,1/REP,FASTSAVEET,1,SOLID45MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,206000MPDATA,PRXY,1,,0.29TB,BISO,1,1,2,TBTEMP,0TBDATA,,300,12000,,,,/prep7MSHAPE,0,3DMSHKEY,1VMESH,all/SOLUDA,3,ALL,*DIM,dis,TABLE,9,1,,TIME,,DIS(1,0)=0,1,2,3,4,5,6,7,8DIS(1,1)=0,3,0,-3,0,4,0,-4,0D,22,,%DIS%,,,,UZ,,,,,NSUBST,40,0,0OUTRES,BASIC,-40TIME,9/STATUS,SOLUSOLVEFINISH/post26NSOL,2,22,U,z,UzRFORCE,3,22,F,z,FzPROD,3,3,,,,,,0.001,1,1,VARNAM,3,LOADPLTIME,0,0XVAR,2SPREAD,0PLCPLX,0PLVAR,3,,,,,,,,,,/AXLAB,X,displacement(cm)/AXLAB,Y,load(N)其中定义施加往复位移的命令：*DIM,dis,TABLE,9,1,,TIME,,DIS(1,0)=0,1,2,3,4,5,6,7,8DIS(1,1)=0,3,0,-3,0,4,0,-4,0D,22,,%DIS%,,,,UZ,,,,,各位朋友：我在分析一悬臂板的滞回曲线时，底边固定，顶部节点X方向的自由度耦合(编号为1)，顶部节点采用位移加载，命令流如下。请问我如何得到顶部节点的力(即顶部所有节点的合力)-水平位移(顶部所有节点耦合后节点位移相同)关系，请指教。fini/clear/PREP7ET,1,SHELL63R,1,3,,,,,,MP,EX,1,2.06E+005MP,PRXY,1,0.3TB,BKIN,1,1,2,1TBDATA,,235,3000,,,,k,1k,2,50k,3,50,200k,4,0,200a,1,2,3,4TYPE,1MAT,1REAL,1ESIZE,5,0,MSHAPE,0,2DMSHKEY,1AMESH,1NSEL,S,LOC,Y,200CP,1,UX,ALLSAVEFINISH施加的位移为：0,5,0,-5,0,10,0,-10,0,15,0,-15,0提供一个本人做支撑位移控制低周往复荷载下命令流，希望能对你有帮助/solu！求解选项设置ANTYPE,STATICnlgeom,onpred,offnropt,full,,onsstif,offNSUBST,100,10000,10cnvtol,u,,0.03,0autots,onOUTRES,ALL,1kbc,0！斜坡荷载*dim,disp,,number*VREAD,disp(1),disp2,txt,,IJK,number(E3.0)！读荷载数组*DO,I,1,NUMBERtime,iDK,s_num+1,,-235*s_length*disp(i)/6/2.06e5,,0,UX,,,,,,！位移加载ACEL,0,9.8,0,SOLVE*ENDDOsave