【新人】做了一个拱结构

wang780738

导师让做一个拱结构的实例，让我分析平面内失稳和平面外失稳的临界荷载有多大，现在的基本参数有高11米，跨度25米，让我找了一个frp的本构模型，我找的是弹性模量是38GPA，抗拉强度是800MPA，我自己算了一下，单位荷载作用下（竖向）屈曲系数是10N，单位线荷载作用下是280n/m，我觉得这个也太小了吧，还是说我做的不对，还有，在ansys中怎么实现平面外稳定和平面内稳定，下面把命令流发上来

wang780738

/BATCH  
/COM,ANSYS RELEASE 15.0    UP20131014       15:08:57    03/09/2019
/TITLE,asd  
/REPLOT,RESIZE  
!*  
/NOPR   
KEYW,PR_SET,1   
KEYW,PR_STRUC,1
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,0
/GO
!*  
/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Structural   
!*  
/PREP7  
!*  
ET,1,BEAM188
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,38e9   
MPDATA,PRXY,1,,0.3  
tb,biso,1,1,2,  
tbtemp,0
tbdata,,8e8,,,,
SECTYPE,   1, BEAM, RECT, asd, 0
SECOFFSET, CENT
SECDATA,0.01,0.03,0,0,0,0,0,0,0,0,0,0   
csys,1  
k,1,12.5
k,2,11,90   
k,3,-12.5   
l,1,2   
l,2,3   
csys,0  
!*  
LATT,1, ,1, , , ,1  
ESIZE,0,10,
FLST,2,2,4,ORDE,2   
FITEM,2,1   
FITEM,2,-2  
LMESH,P51X  
FINISH  
/SOL
!*  
ANTYPE,0
PSTRES,1
FLST,2,2,3,ORDE,2   
FITEM,2,1   
FITEM,2,3   
!*  
/GO
DK,P51X, , , ,0,ALL, , , , , ,  
*SET,qh,1   
*get,ne,elem,,count
*do,i,1,ne  
*SET,ni,nelem(i,1)  
*SET,nj,nelem(i,2)  
*SET,dy,ny(nj)-ny(ni)   
*SET,dx,nx(nj)-nx(ni)   
*SET,dyx,dy/dx  
*SET,dyx2,1+dyx*dyx
*SET,qsv,qh/dyx2
*SET,qsp,qh*dyx/dyx2
sfbeam,i,1,pres,-qsv
sfbeam,i,2,pres,qsp
*enddo  
/STATUS,SOLU
SOLVE   
FINISH  
/SOLU   
!*  
ANTYPE,1
!*  
BUCOPT,LANB,6,0,0,CENTER
MXPAND,6,0,0,1,0.001,   
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  
SET,LIST
!*  
*GET,load,ACTIVE, ,SET,FREQ
FINISH  
/PREP7  
UPGEOM,0.005,1,1,'asd','rst',' '
!*  
FINISH  
/SOL
!*  
ANTYPE,0
ANTYPE,0
NLGEOM,1
NSUBST,50,1e5,50
OUTRES,ERASE
OUTRES,ALL,ALL  
LNSRCH,1
TIME,load*1.15  
SFSCALE,ALL,1.15*load,1.15*load,
/STATUS,SOLU
SOLVE   
FINISH  
/POST1  
SET,LIST
FINISH  
/POST26
FILE,'asd','rst','.'
/UI,COLL,1  
NUMVAR,200  
SOLU,191,NCMIT  
STORE,MERGE
FILLDATA,191,,,,1,1
REALVAR,191,191
!*  
!*  
!*  
NSOL,2,2,U,Y,   
!*  
XVAR,2  
PLVAR,2,
/AXLAB,X,位移   
/AXLAB,Y,力
/GTHK,AXIS,2
/GRTYP,0
/GROPT,ASCAL,ON
/GROPT,LOGX,OFF
/GROPT,LOGY,OFF
/GROPT,AXDV,1   
/GROPT,AXNM,ON  
/GROPT,AXNSC,1,
/GROPT,DIG1,4,  
/GROPT,DIG2,3,  
/GROPT,XAXO,0,  
/GROPT,YAXO,0,  
/GROPT,DIVX,
/GROPT,DIVY,
/GROPT,REVX,1   
/GROPT,REVY,1   
/GROPT,LTYP,0   
!*  
/XRANGE,DEFAULT
/YRANGE,DEFAULT,,1  
!*  
/AUTO,1
/REP,FAST   
FINISH  
/PREP7  
!*  
TBDE,BISO,1,,,  
TB,BISO,1,1,2,  
TBTEMP,0
TBDATA,,8E+008,0,,,,
FINISH  
! /EXIT,NOSAV   

wang780738

只分析道线性分析也行，分析到非线性分析也行

O2ange

朋友  请问你拱 平面内反对称屈曲  怎么做的
我用abaqus做  直接线性摄动不行
是否要自己导入初始几何缺陷
这个有规定？

南树先生

直接用risks做可以吗？我觉得用屈曲挺好的