连续刚构桥参数优化分析
最近在写论文，做了一个五跨连续刚构的模型，由于是用solid45单元模拟的，不能直接得出内力图，于是用积分的方法计算内力，加自重后求解出两侧对称截面的内力不相等，求高手指教。。。下面给出建模、加载以及内力计算的命令流：
1.建模及网格划分
/PREP7*CREATE,QIANCHULI,macH1=6.4 ！边跨梁高H2=7.0 ！中跨梁高E=2.65 ！两薄壁墩中心距B1=0.7 ！1/2边墩厚度B2=0.7 ！1/2中墩厚度L1=66！边跨长L2=120 ！次中跨长L3=120 ！中跨长K,1,0,0,0,K,2,3,0,0,K,3,3,2.2,0,K,4,3.5,2.7,0,K,5,5.5,2.85,0, K,6,5.5,3,0,K,7,0,3,0,K,8,0,2.7,0,K,9,1.5,2.7,0,K,10,2.5,2.2,0, K,11,2.5,0.8,0, K,12,2,0.4,0, K,13,0,0.4,0, LSTR, 1, 2LSTR, 2, 3LSTR, 3, 4LSTR, 4, 5LSTR, 5, 6LSTR, 6, 7LSTR, 7, 8LSTR, 8, 9LSTR, 10, 9LSTR, 10, 11LSTR, 11, 12LSTR, 12, 13LSTR, 13, 1FLST,2,13,4 FITEM,2,1 FITEM,2,2 FITEM,2,3 FITEM,2,4 FITEM,2,5 FITEM,2,6 FITEM,2,7 FITEM,2,8 FITEM,2,9 FITEM,2,10FITEM,2,11FITEM,2,12FITEM,2,13AL,P51X FLST,3,1,5,ORDE,1 FITEM,3,1 ARSYM,X,P51X, , , ,0,0FLST,2,2,5,ORDE,2 FITEM,2,1 FITEM,2,-2AADD,P51X FLST,2,2,4,ORDE,2 FITEM,2,8 FITEM,2,29LCOMB,P51X, ,0FLST,2,2,4,ORDE,2 FITEM,2,6 FITEM,2,28LCOMB,P51X, ,0FLST,2,2,4,ORDE,2 FITEM,2,12FITEM,2,30LCOMB,P51X, ,0FLST,2,2,4,ORDE,2 FITEM,2,1 FITEM,2,27LCOMB,P51X, ,0KPLOT LPLOT LSTR, 16, 23LSTR, 10, 3LSTR, 24, 15LSTR, 11, 2FLST,2,12,4 FITEM,2,17FITEM,2,16FITEM,2,7 FITEM,2,22FITEM,2,8 FITEM,2,9 FITEM,2,13FITEM,2,3 FITEM,2,4 FITEM,2,5 FITEM,2,6 FITEM,2,18AL,P51X FLST,2,4,4FITEM,2,23FITEM,2,7 FITEM,2,15FITEM,2,14AL,P51X FLST,2,6,4FITEM,2,24FITEM,2,14FITEM,2,1 FITEM,2,19FITEM,2,11FITEM,2,12AL,P51X FLST,2,4,4FITEM,2,19FITEM,2,2 FITEM,2,13FITEM,2,10AL,P51X K,100,3,3-H1,L1-E-B1, K,1000,3,-10,30 , LARC,2,100,1000,370.119 APLOT K,200,5.5,3,L1-E-B1,LSTR, 200, 6LPLOT FLST,2,2,4,ORDE,2 FITEM,2,14FITEM,2,19ADRAG, 19, , , , , , 20 ADRAG, 13, , , , , , 21 LDELE, 20LSTR, 8, 1LSTR, 13, 7FLST,2,4,4FITEM,2,10FITEM,2,29FITEM,2,20FITEM,2,26AL,P51X FLST,2,4,4FITEM,2,2 FITEM,2,30FITEM,2,31FITEM,2,27AL,P51X FLST,2,4,4FITEM,2,28FITEM,2,20FITEM,2,25FITEM,2,31AL,P51X FLST,2,6,5,ORDE,2 FITEM,2,5 FITEM,2,-10 VA,P51X FLST,3,1,6,ORDE,1 FITEM,3,1 VSYMM,X,P51X, , , ,0,0VDRAG, 4, , , , , , 27 VDRAG, 1, , , , , , 21 FLST,3,4,6,ORDE,2 FITEM,3,1 FITEM,3,-4K,300,5.5,3,L1+E+B1 LSTR, 300, 200APLOT FLST,2,4,5,ORDE,4 FITEM,2,10FITEM,2,16FITEM,2,23FITEM,2,36VDRAG,P51X, , , , , , 80 K,2001,3,3-H2,L1+L2-E-B2, K,2002,5.5,3,L1+L2-E-B2,APLOT LSTR, 66, 2002K,2000,3,-10,100LARC,61,2001,2000,370.119 LPLOT ADRAG, 86, , , , , , 134 FLST,2,2,4,ORDE,2 FITEM,2,81FITEM,2,133 ADRAG, 81, , , , , , 133 LSTR, 77, 76LSTR, 78, 75FLST,2,4,4FITEM,2,137 FITEM,2,88FITEM,2,139 FITEM,2,141 AL,P51X FLST,2,4,4FITEM,2,142 FITEM,2,140 FITEM,2,84FITEM,2,136 AL,P51X FLST,2,4,4FITEM,2,138 FITEM,2,141 FITEM,2,135 FITEM,2,142 AL,P51X APLOT FLST,2,6,5,ORDE,3 FITEM,2,41FITEM,2,67FITEM,2,-71 VA,P51X FLST,3,1,6,ORDE,1 FITEM,3,9 VSYMM,X,P51X, , , ,0,0VDRAG, 53, , , , , , 134 VDRAG, 66, , , , , , 133 VPLOT K,2003,5.5,3,L1+L2+E+B2 LSTR, 2002, 2003APLOT FLST,2,4,5,ORDE,4 FITEM,2,71FITEM,2,77FITEM,2,84FITEM,2,97VDRAG,P51X, , , , , , 191 K,2004,3,3-H2,L1+L2+E+B2K,3000,3,-10,200, K,3003,5.5,3,L1+L2+L3-E-B2K,3004,3,3-H2,L1+L2+L3-E-B2 LSTR, 2003, 3003LARC,2004,3004,3000,370.119 LPLOT ADRAG, 197, , , , , , 245 ADRAG, 192, , , , , , 244 LSTR, 134, 133LSTR, 132, 135FLST,2,4,4FITEM,2,252 FITEM,2,250 FITEM,2,199 FITEM,2,248 AL,P51X FLST,2,4,4FITEM,2,253 FITEM,2,251 FITEM,2,195 FITEM,2,247 AL,P51X FLST,2,4,4FITEM,2,252 FITEM,2,249 FITEM,2,253 FITEM,2,246 AL,P51X FLST,2,6,5,ORDE,3 FITEM,2,102 FITEM,2,128 FITEM,2,-132VA,P51X VDRAG, 114, , , , , , 248 VDRAG, 127, , , , , , 244 FLST,3,1,6,ORDE,1 FITEM,3,17VSYMM,X,P51X, , , ,0,0BOPTN,KEEP,0BOPTN,NWARN,0 BOPTN,VERS,RV52 BTOL,0.01,FLST,2,4,6,ORDE,2 FITEM,2,17FITEM,2,-20 VADD,P51X FLST,2,16,6,ORDE,2FITEM,2,1 FITEM,2,-16 VADD,P51X FLST,3,1,6,ORDE,1 FITEM,3,17VSYMM,Z,P51X, , , ,0,0FLST,3,1,6,ORDE,1 FITEM,3,1 VGEN, ,P51X, , , , ,492, , ,1 FLST,2,3,6,ORDE,3 FITEM,2,1 FITEM,2,17FITEM,2,21VADD,P51X BOPTN,KEEP,0BOPTN,NWARN,0 BOPTN,VERS,RV52 BTOL,0.00001, K,4001,3,3-H1,L1-E-B1,K,4002,3,-24,L1-E-B1, K,4003,-3,-24,L1-E-B1,K,4004,3,-24,L1-E+B1, K,4005,-3,-24,L1-E+B1,LSTR, 4001, 4002FLST,2,4,3FITEM,2,4005FITEM,2,4004FITEM,2,4002FITEM,2,4003A,P51XVDRAG, 17, , , , , , 96 FLST,3,1,6,ORDE,1 FITEM,3,1 VGEN,2,P51X, , , , ,5.3, ,0 K,4006,3,3-H2,L1+L2-E-B2, K,4007,3,-24,L1+L2-E-B2,K,4008,-3,-24,L1+L2-E-B2, K,4009,3,-24,L1+L2-E+B2 K,4010,-3,-24,L1+L2-E+B2LSTR, 4006, 4007FLST,2,4,3FITEM,2,4010FITEM,2,4009FITEM,2,4007FITEM,2,4008A,P51XVDRAG, 133, , , , , , 264 FLST,3,1,6,ORDE,1 FITEM,3,4 VGEN,2,P51X, , , , ,5.3, ,0 VPLOT FLST,3,2,6,ORDE,2 FITEM,3,1 FITEM,3,3 FLST,3,4,6,ORDE,3 FITEM,3,1 FITEM,3,3 FITEM,3,-5VSYMM,Z,P51X, , , ,0,0FLST,3,4,6,ORDE,2 FITEM,3,6 FITEM,3,-9VGEN, ,P51X, , , , ,492, , ,1 FLST,2,9,6,ORDE,2 FITEM,2,1 FITEM,2,-9VGLUE,P51XKEYW,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 KEYW,PR_CFD,0 ET,1,SOLID45ET,2,SOLID45MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,3.5E+010 MPDATA,PRXY,1,,0.167MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,2500 MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,3.3E+010 MPDATA,PRXY,2,,0.167MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,2,,2500 TYPE, 2 MAT, 2REAL, ESYS, 0 SECNUM, LESIZE,ALL,1, , , ,1, , ,1, FLST,5,8,6,ORDE,3 FITEM,5,1 FITEM,5,10FITEM,5,-16 CM,_Y,VOLUVSEL, , , ,P51X CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_YMSHAPE,0,3d MSHKEY,1VMESH,_Y1 MSHKEY,0CMDELE,_Y CMDELE,_Y1CMDELE,_Y2TYPE, 1 MAT, 1REAL, ESYS, 0 SECNUM, MSHKEY,0MSHAPE,1,3d CM,_Y,VOLUVSEL, , , , 17 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_YVMESH,_Y1 CMDELE,_Y CMDELE,_Y1CMDELE,_Y2*END
2.2.自重荷载与求解/prep7*CREATE,ZIZHONG,macnsel,s,loc,z,0FINISH/SOLFLST,2,58,1,ORDE,12 FITEM,2,4079FITEM,2,-4093 FITEM,2,4592FITEM,2,-4593 FITEM,2,4720FITEM,2,-4747 FITEM,2,6323FITEM,2,6386FITEM,2,6459FITEM,2,-6460 FITEM,2,7667FITEM,2,-7675 /GO D,P51X, ,0, , , ,UY, , , , ,nsel,s,loc,z,492FLST,2,167,1,ORDE,85FITEM,2,18370 FITEM,2,-18380FITEM,2,18441 FITEM,2,-18443FITEM,2,18504 FITEM,2,-18508FITEM,2,18568 FITEM,2,-18569FITEM,2,18634 FITEM,2,-18638FITEM,2,18698 FITEM,2,18759 FITEM,2,-18760FITEM,2,18880 FITEM,2,-18881FITEM,2,18883 FITEM,2,-18885FITEM,2,18946 FITEM,2,-18948FITEM,2,19010 FITEM,2,-19040FITEM,2,19161 FITEM,2,-19162FITEM,2,19164 FITEM,2,-19166FITEM,2,19228 FITEM,2,-19230FITEM,2,19350 FITEM,2,-19351FITEM,2,19411 FITEM,2,-19421FITEM,2,19490 FITEM,2,19550 FITEM,2,-19556FITEM,2,19973 FITEM,2,-19976FITEM,2,20038 FITEM,2,-20040FITEM,2,20161 FITEM,2,-20162FITEM,2,20164 FITEM,2,-20166FITEM,2,20228 FITEM,2,-20230FITEM,2,20291 FITEM,2,-20292FITEM,2,20354 FITEM,2,-20356FITEM,2,20477 FITEM,2,20537 FITEM,2,-20539FITEM,2,20543 FITEM,2,-20546FITEM,2,20608 FITEM,2,-20610FITEM,2,20612 FITEM,2,-20613FITEM,2,20749 FITEM,2,-20751FITEM,2,20761 FITEM,2,20822 FITEM,2,20888 FITEM,2,21932 FITEM,2,-21940FITEM,2,23328 FITEM,2,-23329FITEM,2,23390 FITEM,2,23392 FITEM,2,-23393FITEM,2,23513 FITEM,2,-23514FITEM,2,23575 FITEM,2,23578 FITEM,2,-23579FITEM,2,23639 FITEM,2,-23641FITEM,2,23644 FITEM,2,-23645FITEM,2,23707 FITEM,2,23779 FITEM,2,23840 FITEM,2,23873 FITEM,2,23911 FITEM,2,28849 FITEM,2,28866 /GO D,P51X, ,0, , , ,UY, , , , ,nsel,s,loc,y,-24FLST,2,168,1,ORDE,84FITEM,2,1 FITEM,2,-7FITEM,2,155 FITEM,2,-156FITEM,2,199 FITEM,2,-204FITEM,2,331 FITEM,2,353 FITEM,2,-357FITEM,2,463 FITEM,2,-469FITEM,2,617 FITEM,2,-618FITEM,2,661 FITEM,2,-666FITEM,2,793 FITEM,2,815 FITEM,2,-819FITEM,2,925 FITEM,2,-931FITEM,2,1072FITEM,2,-1073 FITEM,2,1114FITEM,2,-1119 FITEM,2,1240FITEM,2,1261FITEM,2,-1265 FITEM,2,1366FITEM,2,-1372 FITEM,2,1513FITEM,2,-1514 FITEM,2,1555FITEM,2,-1560 FITEM,2,1681FITEM,2,1702FITEM,2,-1706 FITEM,2,1807FITEM,2,1835FITEM,2,1856FITEM,2,-1860 FITEM,2,1963FITEM,2,1984FITEM,2,2011FITEM,2,2032FITEM,2,-2036 FITEM,2,2138FITEM,2,2159FITEM,2,-2163 FITEM,2,2270FITEM,2,2297FITEM,2,2318FITEM,2,-2322 FITEM,2,2425FITEM,2,2446FITEM,2,2473FITEM,2,2494FITEM,2,-2498 FITEM,2,2600FITEM,2,2621FITEM,2,-2625 FITEM,2,2731FITEM,2,2758FITEM,2,2778FITEM,2,-2782 FITEM,2,2880FITEM,2,2900FITEM,2,2926FITEM,2,2946FITEM,2,-2950 FITEM,2,3047FITEM,2,3067FITEM,2,-3071 FITEM,2,3173FITEM,2,3199FITEM,2,3219FITEM,2,-3223 FITEM,2,3321FITEM,2,3341FITEM,2,3367FITEM,2,3387FITEM,2,-3391 FITEM,2,3488FITEM,2,3508FITEM,2,-3512 /GO D,P51X, ,0, , , ,ALL, , , , , /REPLOT,RESIZE/PLOPTS,MINM,OFF/REPLOT,RESIZE/REPLOT,RESIZE/REPLOT,RESIZEALLSEL,ALL/REPLOT,RESIZE/PLOPTS,MINM,OFF/REPLOT,RESIZE/REPLOT,RESIZE/REPLOT,RESIZEACEL,0,10,0,/STATUS,SOLUSOLVE *END
3.求截面内力宏命令：
*CREATE,Su_NQM,mac
SUCR,su_na,cplane,3
SUSEL,s,su_na
SUMAP,su_sx,s,x
SUMAP,su_sy,s,y
SUMAP,su_sz,s,z
SUMAP,su_sxy,s,xy
SUMAP,su_syz,s,yz
SUMAP,su_sxz,s,xz
SUEVAL,su_Nx,su_sx,intg
SUEVAL,su_Ny,su_sy,intg
SUEVAL,su_Nz,su_sz,intgSUEVAL,su_Qx,su_syz,intg
SUEVAL,su_Qy,su_sxz,intg
SUEVAL,su_Qz,su_sxy,intgSUCALC,sxgcy,su_sx,mult,gcy
SUEVAL,M_XY,sxgcy,intg
SUCALC,sxgcz,su_sx,mult,gcz
SUEVAL,M_XZ,sxgcz,intg
SUCALC,Sxgcz,su_sy,mult,gcz
SUEVAL,M_YZ,sxgcz,intg
SUCALC,sxgcx,su_sy,mult,gcx
SUEVAL,M_YX,sxgcx,intg
SUCALC,sxgcx,su_sz,mult,gcx
SUEVAL,M_ZX,sxgcx,intg
SUCALC,sxgcy,su_sz,mult,gcy
SUEVAL,M_ZY,sxgcy,intg
SUEVAL,A_area,da,sum
SUEVAL,ss_cx,gcx,intg
Su_cx=ss_cx/A_area
SUEVAL,ss_cy,gcy,intg
Su_cy=ss_cy/A_area
SUEVAL,ss_cz,gcz,intg
Su_cz=ss_cz/A_area
M_XY=M_XY-Su_Nx*su_cy
M_XZ=M_XZ-Su_Nx*su_cz
M_YZ=M_YZ-Su_Ny*su_cz
M_YX=M_YX-Su_Ny*su_cx
M_ZX=M_ZX-Su_Nz*su_cx
M_ZY=M_ZY-Su_Nz*su_cy
*END
呵呵 好东西之吃{:4_90:} 应该检查一下模型两侧是不是一致，有时候桥博和midas计算的内力也不一致，可能是截面不完全一样；你这个我只看懂了一半，还没有看完，机器上没装ANSYS没法测试，很遗憾！
页:
[1]