最近在写论文，做了一个五跨连续刚构的模型，由于是用solid45单元模拟的，不能直接得出内力图，于是用积分的方法计算内力，加自重后求解出两侧对称截面的内力不相等，求高手指教。。。
下面给出建模、加载以及内力计算的命令流：
1.建模及网格划分
/PREP7 *CREATE,QIANCHULI,mac H1=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, 2 LSTR, 2, 3 LSTR, 3, 4 LSTR, 4, 5 LSTR, 5, 6 LSTR, 6, 7 LSTR, 7, 8 LSTR, 8, 9 LSTR, 10, 9 LSTR, 10, 11 LSTR, 11, 12 LSTR, 12, 13 LSTR, 13, 1 FLST,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,10 FITEM,2,11 FITEM,2,12 FITEM,2,13 AL,P51X FLST,3,1,5,ORDE,1 FITEM,3,1 ARSYM,X,P51X, , , ,0,0 FLST,2,2,5,ORDE,2 FITEM,2,1 FITEM,2,-2 AADD,P51X FLST,2,2,4,ORDE,2 FITEM,2,8 FITEM,2,29 LCOMB,P51X, ,0 FLST,2,2,4,ORDE,2 FITEM,2,6 FITEM,2,28 LCOMB,P51X, ,0 FLST,2,2,4,ORDE,2 FITEM,2,12 FITEM,2,30 LCOMB,P51X, ,0 FLST,2,2,4,ORDE,2 FITEM,2,1 FITEM,2,27 LCOMB,P51X, ,0 KPLOT LPLOT LSTR, 16, 23 LSTR, 10, 3 LSTR, 24, 15 LSTR, 11, 2 FLST,2,12,4 FITEM,2,17 FITEM,2,16 FITEM,2,7 FITEM,2,22 FITEM,2,8 FITEM,2,9 FITEM,2,13 FITEM,2,3 FITEM,2,4 FITEM,2,5 FITEM,2,6 FITEM,2,18 AL,P51X FLST,2,4,4 FITEM,2,23 FITEM,2,7 FITEM,2,15 FITEM,2,14 AL,P51X FLST,2,6,4 FITEM,2,24 FITEM,2,14 FITEM,2,1 FITEM,2,19 FITEM,2,11 FITEM,2,12 AL,P51X FLST,2,4,4 FITEM,2,19 FITEM,2,2 FITEM,2,13 FITEM,2,10 AL,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, 6 LPLOT FLST,2,2,4,ORDE,2 FITEM,2,14 FITEM,2,19 ADRAG, 19, , , , , , 20 ADRAG, 13, , , , , , 21 LDELE, 20 LSTR, 8, 1 LSTR, 13, 7 FLST,2,4,4 FITEM,2,10 FITEM,2,29 FITEM,2,20 FITEM,2,26 AL,P51X FLST,2,4,4 FITEM,2,2 FITEM,2,30 FITEM,2,31 FITEM,2,27 AL,P51X FLST,2,4,4 FITEM,2,28 FITEM,2,20 FITEM,2,25 FITEM,2,31 AL,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,0 VDRAG, 4, , , , , , 27 VDRAG, 1, , , , , , 21 FLST,3,4,6,ORDE,2 FITEM,3,1 FITEM,3,-4 K,300,5.5,3,L1+E+B1 LSTR, 300, 200 APLOT FLST,2,4,5,ORDE,4 FITEM,2,10 FITEM,2,16 FITEM,2,23 FITEM,2,36 VDRAG,P51X, , , , , , 80 K,2001,3,3-H2,L1+L2-E-B2, K,2002,5.5,3,L1+L2-E-B2, APLOT LSTR, 66, 2002 K,2000,3,-10,100 LARC,61,2001,2000,370.119 LPLOT ADRAG, 86, , , , , , 134 FLST,2,2,4,ORDE,2 FITEM,2,81 FITEM,2,133 ADRAG, 81, , , , , , 133 LSTR, 77, 76 LSTR, 78, 75 FLST,2,4,4 FITEM,2,137 FITEM,2,88 FITEM,2,139 FITEM,2,141 AL,P51X FLST,2,4,4 FITEM,2,142 FITEM,2,140 FITEM,2,84 FITEM,2,136 AL,P51X FLST,2,4,4 FITEM,2,138 FITEM,2,141 FITEM,2,135 FITEM,2,142 AL,P51X APLOT FLST,2,6,5,ORDE,3 FITEM,2,41 FITEM,2,67 FITEM,2,-71 VA,P51X FLST,3,1,6,ORDE,1 FITEM,3,9 VSYMM,X,P51X, , , ,0,0 VDRAG, 53, , , , , , 134 VDRAG, 66, , , , , , 133 VPLOT K,2003,5.5,3,L1+L2+E+B2 LSTR, 2002, 2003 APLOT FLST,2,4,5,ORDE,4 FITEM,2,71 FITEM,2,77 FITEM,2,84 FITEM,2,97 VDRAG,P51X, , , , , , 191 K,2004,3,3-H2,L1+L2+E+B2 K,3000,3,-10,200, K,3003,5.5,3,L1+L2+L3-E-B2 K,3004,3,3-H2,L1+L2+L3-E-B2 LSTR, 2003, 3003 LARC,2004,3004,3000,370.119 LPLOT ADRAG, 197, , , , , , 245 ADRAG, 192, , , , , , 244 LSTR, 134, 133 LSTR, 132, 135 FLST,2,4,4 FITEM,2,252 FITEM,2,250 FITEM,2,199 FITEM,2,248 AL,P51X FLST,2,4,4 FITEM,2,253 FITEM,2,251 FITEM,2,195 FITEM,2,247 AL,P51X FLST,2,4,4 FITEM,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,-132 VA,P51X VDRAG, 114, , , , , , 248 VDRAG, 127, , , , , , 244 FLST,3,1,6,ORDE,1 FITEM,3,17 VSYMM,X,P51X, , , ,0,0 BOPTN,KEEP,0 BOPTN,NWARN,0 BOPTN,VERS,RV52 BTOL,0.01, FLST,2,4,6,ORDE,2 FITEM,2,17 FITEM,2,-20 VADD,P51X FLST,2,16,6,ORDE,2 FITEM,2,1 FITEM,2,-16 VADD,P51X FLST,3,1,6,ORDE,1 FITEM,3,17 VSYMM,Z,P51X, , , ,0,0 FLST,3,1,6,ORDE,1 FITEM,3,1 VGEN, ,P51X, , , , ,492, , ,1 FLST,2,3,6,ORDE,3 FITEM,2,1 FITEM,2,17 FITEM,2,21 VADD,P51X BOPTN,KEEP,0 BOPTN,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, 4002 FLST,2,4,3 FITEM,2,4005 FITEM,2,4004 FITEM,2,4002 FITEM,2,4003 A,P51X VDRAG, 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+B2 LSTR, 4006, 4007 FLST,2,4,3 FITEM,2,4010 FITEM,2,4009 FITEM,2,4007 FITEM,2,4008 A,P51X VDRAG, 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,-5 VSYMM,Z,P51X, , , ,0,0 FLST,3,4,6,ORDE,2 FITEM,3,6 FITEM,3,-9 VGEN, ,P51X, , , , ,492, , ,1 FLST,2,9,6,ORDE,2 FITEM,2,1 FITEM,2,-9 VGLUE,P51X 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 KEYW,PR_CFD,0 ET,1,SOLID45 ET,2,SOLID45 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,3.5E+010 MPDATA,PRXY,1,,0.167 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,2500 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,2,,3.3E+010 MPDATA,PRXY,2,,0.167 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,2500 TYPE, 2 MAT, 2 REAL, ESYS, 0 SECNUM, LESIZE,ALL,1, , , ,1, , ,1, FLST,5,8,6,ORDE,3 FITEM,5,1 FITEM,5,10 FITEM,5,-16 CM,_Y,VOLU VSEL, , , ,P51X CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y MSHAPE,0,3d MSHKEY,1 VMESH,_Y1 MSHKEY,0 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 TYPE, 1 MAT, 1 REAL, ESYS, 0 SECNUM, MSHKEY,0 MSHAPE,1,3d CM,_Y,VOLU VSEL, , , , 17 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y VMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 *END
2. 2.自重荷载与求解 /prep7 *CREATE,ZIZHONG,mac nsel,s,loc,z,0 FINISH /SOL FLST,2,58,1,ORDE,12 FITEM,2,4079 FITEM,2,-4093 FITEM,2,4592 FITEM,2,-4593 FITEM,2,4720 FITEM,2,-4747 FITEM,2,6323 FITEM,2,6386 FITEM,2,6459 FITEM,2,-6460 FITEM,2,7667 FITEM,2,-7675 /GO D,P51X, ,0, , , ,UY, , , , , nsel,s,loc,z,492 FLST,2,167,1,ORDE,85 FITEM,2,18370 FITEM,2,-18380 FITEM,2,18441 FITEM,2,-18443 FITEM,2,18504 FITEM,2,-18508 FITEM,2,18568 FITEM,2,-18569 FITEM,2,18634 FITEM,2,-18638 FITEM,2,18698 FITEM,2,18759 FITEM,2,-18760 FITEM,2,18880 FITEM,2,-18881 FITEM,2,18883 FITEM,2,-18885 FITEM,2,18946 FITEM,2,-18948 FITEM,2,19010 FITEM,2,-19040 FITEM,2,19161 FITEM,2,-19162 FITEM,2,19164 FITEM,2,-19166 FITEM,2,19228 FITEM,2,-19230 FITEM,2,19350 FITEM,2,-19351 FITEM,2,19411 FITEM,2,-19421 FITEM,2,19490 FITEM,2,19550 FITEM,2,-19556 FITEM,2,19973 FITEM,2,-19976 FITEM,2,20038 FITEM,2,-20040 FITEM,2,20161 FITEM,2,-20162 FITEM,2,20164 FITEM,2,-20166 FITEM,2,20228 FITEM,2,-20230 FITEM,2,20291 FITEM,2,-20292 FITEM,2,20354 FITEM,2,-20356 FITEM,2,20477 FITEM,2,20537 FITEM,2,-20539 FITEM,2,20543 FITEM,2,-20546 FITEM,2,20608 FITEM,2,-20610 FITEM,2,20612 FITEM,2,-20613 FITEM,2,20749 FITEM,2,-20751 FITEM,2,20761 FITEM,2,20822 FITEM,2,20888 FITEM,2,21932 FITEM,2,-21940 FITEM,2,23328 FITEM,2,-23329 FITEM,2,23390 FITEM,2,23392 FITEM,2,-23393 FITEM,2,23513 FITEM,2,-23514 FITEM,2,23575 FITEM,2,23578 FITEM,2,-23579 FITEM,2,23639 FITEM,2,-23641 FITEM,2,23644 FITEM,2,-23645 FITEM,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,-24 FLST,2,168,1,ORDE,84 FITEM,2,1 FITEM,2,-7 FITEM,2,155 FITEM,2,-156 FITEM,2,199 FITEM,2,-204 FITEM,2,331 FITEM,2,353 FITEM,2,-357 FITEM,2,463 FITEM,2,-469 FITEM,2,617 FITEM,2,-618 FITEM,2,661 FITEM,2,-666 FITEM,2,793 FITEM,2,815 FITEM,2,-819 FITEM,2,925 FITEM,2,-931 FITEM,2,1072 FITEM,2,-1073 FITEM,2,1114 FITEM,2,-1119 FITEM,2,1240 FITEM,2,1261 FITEM,2,-1265 FITEM,2,1366 FITEM,2,-1372 FITEM,2,1513 FITEM,2,-1514 FITEM,2,1555 FITEM,2,-1560 FITEM,2,1681 FITEM,2,1702 FITEM,2,-1706 FITEM,2,1807 FITEM,2,1835 FITEM,2,1856 FITEM,2,-1860 FITEM,2,1963 FITEM,2,1984 FITEM,2,2011 FITEM,2,2032 FITEM,2,-2036 FITEM,2,2138 FITEM,2,2159 FITEM,2,-2163 FITEM,2,2270 FITEM,2,2297 FITEM,2,2318 FITEM,2,-2322 FITEM,2,2425 FITEM,2,2446 FITEM,2,2473 FITEM,2,2494 FITEM,2,-2498 FITEM,2,2600 FITEM,2,2621 FITEM,2,-2625 FITEM,2,2731 FITEM,2,2758 FITEM,2,2778 FITEM,2,-2782 FITEM,2,2880 FITEM,2,2900 FITEM,2,2926 FITEM,2,2946 FITEM,2,-2950 FITEM,2,3047 FITEM,2,3067 FITEM,2,-3071 FITEM,2,3173 FITEM,2,3199 FITEM,2,3219 FITEM,2,-3223 FITEM,2,3321 FITEM,2,3341 FITEM,2,3367 FITEM,2,3387 FITEM,2,-3391 FITEM,2,3488 FITEM,2,3508 FITEM,2,-3512 /GO D,P51X, ,0, , , ,ALL, , , , , /REPLOT,RESIZE /PLOPTS,MINM,OFF /REPLOT,RESIZE /REPLOT,RESIZE /REPLOT,RESIZE ALLSEL,ALL /REPLOT,RESIZE /PLOPTS,MINM,OFF /REPLOT,RESIZE /REPLOT,RESIZE /REPLOT,RESIZE ACEL,0,10,0, /STATUS,SOLU SOLVE *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,intg
SUEVAL,su_Qx,su_syz,intg
SUEVAL,su_Qy,su_sxz,intg
SUEVAL,su_Qz,su_sxy,intg SUCALC,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
|
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
显身卡
