双向密肋楼盖的有限元分析
<P>我现在在做双向密肋楼盖的有限元分析，但在支座处过早出现裂缝。有没有办法解决？</P> 你想怎么做？ <P>FINI<BR>/CLEAR<BR>/UNITS,SI<BR>/CONFIG,NPROC,2<BR>/CONFIG,NRES,5000<BR>/PREP7<BR>PLOAD=-20.436</P><P>ET,1,65<BR>ET,2,8<BR>ET,3,65</P>
<P>!keyopt,1,6,3<BR>!keyopt,1,5,1<BR>!*afun,deg<BR>!R,1,2,0.0026,90,0,2,0.0026,0,0!整体模型，实常数<BR>!定义材料属性<BR>mp,ex,1,3E10!混凝土材料属性<BR>mp,prxy,1,0.2 </P>
<P>TB,MISO,1,1,5, <BR>TBTEMP,0<BR>TBPT,,0.0005,15e6 <BR>TBPT,,0.001,21e6<BR>TBPT,,0.0015,24e6 <BR>TBPT,,0.002,27e6<BR>TBPT,,0.003,24e6<BR> <BR>TB,CONC,1,1,9, <BR>TBTEMP,0<BR>TBDATA,,0.5,0.9,3e6,30e6,,<BR>TBDATA,,,,1,,, <BR>MPTEMP,,,,,,,, <BR>MPTEMP,1,0 </P>
<P>mp,ex,2,2e11!纵向受拉钢筋，钢支座材料属性<BR>mp,prxy,2,0.25<BR>TB,BISO,2,1,2, <BR>TBTEMP,0<BR>TBDATA,,340e6,2e9,,,,</P>
<P>MP,EX,3,3E10<BR>mp,prxy,3,0.2</P>
<P>R,2,3.1415926*0.0035*0.0035,0.25*3.1415926*0.0035*0.0035*0.0035*0.0035</P>
<P>K,1,0.0,0.0,0.024<BR>K,2,0.0,2.710,0.024<BR>K,3,0.0,0.0,0.024<BR>K,4,2.71,0.0,0.024</P>
<P>LSTR,1,2 !L1<BR>LESIZE,1, , ,46<BR>TYPE,2<BR>MAT,2<BR>REAL,2<BR>LMESH,1<BR>*GET,NMAX1,NODE,0,NUM,MAX<BR>*GET,EMAX1,ELEM,0,NUM,MAX<BR>NGEN,18,NMAX1,1,NMAX1,1,0.15<BR>EGEN,18,NMAX1,1,EMAX1,1, , , , , ,0.15<BR>NUMCMP,NODE<BR>NUMCMP,ELEM</P>
<P>LSTR,3,4 !L2<BR>LESIZE,2, , ,46<BR>TYPE,2<BR>MAT,2<BR>REAL,2<BR>*GET,NMAX1,NODE,0,NUM,MAX<BR>*GET,EMAX1,ELEM,0,NUM,MAX<BR>LMESH,2<BR>*GET,NMAX2,NODE,0,NUM,MAX<BR>*GET,EMAX2,ELEM,0,NUM,MAX<BR>NGEN,18,NMAX1,NMAX1+1,NMAX2,1, ,0.15<BR>EGEN,18,NMAX1,EMAX1+1,EMAX2,1, , , , , , ,0.15</P>
<P>BLC4,0.00,0.00,2.730,2.730,0.12 !V1<BR>VSEL, , , ,1<BR>ASLV<BR>LSLA<BR>LSEL,R,LENGTH, ,2.730<BR>LESIZE,ALL, , ,46<BR>VSEL, , , ,1<BR>ASLV<BR>LSLA<BR>LSEL,R,LENGTH, ,0.12<BR>LESIZE,ALL, , ,3</P>
<P>TYPE,1<BR>MAT,1<BR>REAL,1<BR>VMESH,1<BR>ALLS</P>
<P>BLC4,2.55,0,0.18,2.730,-0.06!V2<BR>VSEL, , , ,2<BR>ASLV<BR>LSLA<BR>LSEL,R,LENGTH, ,2.73<BR>LESIZE,ALL, , ,46</P>
<P>VSEL, , , ,2<BR>ASLV<BR>LSLA<BR>LSEL,R,LENGTH, ,0.06<BR>LESIZE,ALL, , ,1</P>
<P>VSEL, , , ,2<BR>ASLV<BR>LSLA<BR>LSEL,R,LENGTH, ,0.18<BR>LESIZE,ALL, , ,3</P>
<P>TYPE,3<BR>VMESH,2</P>
<P>BLC4,0,2.55,2.55,0.18,-0.06 !V3<BR>LSEL,S,LENGTH, ,2.55<BR>LESIZE,ALL, , ,43</P>
<P>VSEL, , , ,3<BR>ASLV</P>
<P>LSLA<BR>LSEL,R,LENGTH, ,0.06<BR>LESIZE,ALL, , ,1<BR>VSEL, , , ,3<BR>ASLV</P>
<P>LSLA<BR>LSEL,R,LENGTH, ,0.18<BR>LESIZE,ALL, , ,3<BR>TYPE,3<BR>VMESH,3</P>
<P>NUMMRG,ALL</P>
<P>/SOLU<BR>NSEL,S,LOC,Z,-0.06<BR>NSEL,R,LOC,Y,2.550,2.730 <BR>NSEL,R,LOC,X,0.00,2.54<BR>!D,ALL, , , , , ,UY,UZ<BR>D,ALL,ALL</P>
<P>NSEL,S,LOC,Z,-0.06 <BR>NSEL,R,LOC,X,2.550,2.730 <BR>NSEL,R,LOC,Y,0.00,2.540<BR>!D,ALL, , , , , ,UX,UZ<BR>D,ALL,ALL</P>
<P>NSEL,S,LOC,Z,-0.06 <BR>NSEL,R,LOC,Y,2.550,2.730 <BR>NSEL,R,LOC,X,2.550,2.730<BR>D,ALL,ALL</P>
<P>ASEL,S,LOC,Y,0.00<BR>ASEL,A,LOC,X,0.00<BR>DA,ALL,SYMM</P>
<P>NSEL,S,LOC,Z,0.12<BR>!SF,ALL,PRES,PLOAD<BR>F,ALL,FZ,PLOAD</P>
<P>NSEL,S,LOC,Z,0.12<BR>NSEL,R,LOC,Y,0<BR>F,ALL,FZ,PLOAD/2.0<BR>NSEL,S,LOC,Z,0.12<BR>NSEL,R,LOC,Y,2.73<BR>F,ALL,FZ,PLOAD/2.0</P>
<P>NSEL,S,LOC,Z,0.12<BR>NSEL,R,LOC,X,0<BR>F,ALL,FZ,PLOAD/2.0<BR>NSEL,S,LOC,Z,0.12<BR>NSEL,R,LOC,X,2.73<BR>F,ALL,FZ,PLOAD/2.0</P>
<P>NSEL,S,LOC,Z,0.12<BR>NSEL,R,LOC,X,0<BR>NSEL,R,LOC,Y,0<BR>F,ALL,FZ,PLOAD/4.0</P>
<P>NSEL,S,LOC,Z,HEIGHT<BR>NSEL,R,LOC,X,2.73<BR>NSEL,R,LOC,Y,0<BR>F,ALL,FZ,PLOAD/4.0</P>
<P>NSEL,S,LOC,Z,0.12<BR>NSEL,R,LOC,X,2.73<BR>NSEL,R,LOC,Y,2.73<BR>F,ALL,FZ,PLOAD/4.0</P>
<P>NSEL,S,LOC,Z,0.12<BR>NSEL,R,LOC,X,0<BR>NSEL,R,LOC,Y,2.73<BR>F,ALL,FZ,PLOAD/4.0</P>
<P>ALLS<BR>!CNVTOL,F, ,0.1,2 !收敛准则<BR>!CNVTOL,U, ,0.002,2 <BR>ANTYPE,0 !分析类型<BR>NLGEOM,1 !大变形开关<BR>NSUBST,200,1000,1 !子步设置</P>
<P>NEQIT,40 !迭代次数<BR>TIME,1 !终止时间<BR>!EQSLV,SPAR !方程求解器<BR>!EQSLV,PCG,1E-4<BR>!EQSLV,ITER,1<BR>!EQSLV,FRONT<BR>LNSRCH,1<BR>!ARCLEN,1,0,0 <BR>!ARCLEN,1,25,0.0001 !弧长法打开<BR>!AUTOTS,-1.0 !关闭自动时间步长<BR>!LNSRCH,-1 !关闭线性搜索</P>
<P>AUTOTS,1 !自动时间步长</P>
<P>SOLVE<BR>SAVE</P>
以上是我的命令流，是块实心板的。当前荷载下支座开裂了，与实验不符。 <P>有没有好心人有相关资料和命令流等，谢谢</P> <P>求救！</P>
页:
[1]