adina提取内力.docx

adina中如何查看三维实体单元的截面轴力、弯矩与剪力？1. cut surface方法下面结合一个自由端作用集中荷载的三维悬臂梁实例，讲解如何计算某截面的轴力、弯矩与剪力。1）实例概况一根完全弹性的悬臂梁，截面尺寸为0.10*0.1，长度为1，在自由端作用2个集中力，数值均为1000，需要计算离自由端距离为0.5单位的横截面上的轴力、弯矩与剪力，按照结构力学，该计算截面的轴力为0，剪力为2000，弯矩为1000.下面通过ADINA程序验证上述数值的正确性。2）建模并求解由于模型比较简单，不详细讲解了，需要说明的是，坐标原点位移自由端截面最下边。命令流如下，最终模型如下图：*DATABASE NEW SAVE=NO PROMPT=NOFEPROGRAM ADINACONTROL FILEVERSION=V83*COORDINATES POINT SYSTEM=0CLEAR1 0.00000000000000 0.00000000000000 0.00000000000000 02 0.00000000000000 0.100000000000000 0.00000000000000 03 0.00000000000000 0.100000000000000 0.100000000000000 04 0.00000000000000 0.00000000000000 0.100000000000000 0*SURFACE VERTEX NAME=1 P1=3 P2=4 P3=1 P4=2*VOLUME EXTRUDED NAME=1 SURFACE=1 DX=1.00000000000000,DY=0.00000000000000 DZ=0.00000000000000 SYSTEM=0 PCOINCID=YES,PTOLERAN=1.00000000000000E-05 NDIV=1 OPTION=VECTOR,RATIO=1.00000000000000 PROGRESS=GEOMETRIC CBIAS=NO*FIXBOUNDARY SURFACES FIXITY=ALLCLEAR6 ALL*LOAD FORCE NAME=1 MAGNITUD=1000.00000000000 FX=0.00000000000000,FY=0.00000000000000 FZ=-1.00000000000000*APPLY-LOAD BODY=0CLEAR1 FORCE 1 POINT 3 0 1 0.00000000000000 0 -1 0 0 0 NO,0.00000000000000 0.00000000000000 1 02 FORCE 1 POINT 4 0 1 0.00000000000000 0 -1 0 0 0 NO,0.00000000000000 0.00000000000000 1 0*MATERIAL ELASTIC NAME=1 E=2.00000000000000E+11 NU=0.300000000000000,DENSITY=7800.00000000000 ALPHA=0.00000000000000 MDESCRIP=NONE*EGROUP THREEDSOLID NAME=1 DISPLACE=DEFAULT STRAINS=DEFAULT MATERIAL=1,RSINT=DEFAULT TINT=DEFAULT RESULTS=STRESSES DEGEN=NO FORMULAT=0,STRESSRE=GLOBAL INITIALS=NONE FRACTUR=NO CMASS=DEFAULT,STRAIN-F=0 UL-FORMU=DEFAULT LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA,INCOMPAT=DEFAULT TIME-OFF=0.00000000000000 POROUS=NO,WTMC=1.00000000000000 OPTION=NONE DESCRIPT=NONE PRINT=DEFAULT,SAVE=DEFAULT TBIRTH=0.00000000000000 TDEATH=0.00000000000000*SUBDIVIDE VOLUME NAME=1 MODE=LENGTH SIZE=0.0500000000000000*GVOLUME NODES=27 PATTERN=0 NCOINCID=BOUNDARIES NCFACE=123456 NCEDGE=,123456789ABC NCVERTEX=12345678 NCTOLERA=1.00000000000000E-05,SUBSTRUC=0 GROUP=1 MESHING=MAPPED PREFSHAP=AUTOMATIC,DEGENERA=YES COLLAPSE=NO MIDNODES=CURVED METHOD=DELAUNAY,BOUNDARY=ADVFRONTCLEAR13）进入后处理程序，计算截面轴力、弯矩与剪力（1）首先定义一个切面CUTPLANE，这个切面就是所求内力的截面。点击右上角的工具栏图标CUT SURFACE,选择TYPE为CUTING PLANE,将X-PLANE 的coordinate value 设为0.5，点击save,效果如下图：（2）定义一个模型点积分，此步骤的作用是定义截面积分，表示后续的轴力、弯矩与剪力将由该截面应力积分求得。操作路径：definitions-model point (special)-mesh integration点击add,名称取为INTEGER-M（或其他也可以），三维的 integrate over 选surface，二维的integrate over 选lines，重点是一定Zone Name一定要选择对，保证该切面能顺利切到该Zone。其他参数不需要修改，（3）定义积分的表达式，根据自己需要求的内力种类和方向来写积分表达式。操作路径：definitions-variable-RESULTANT点击add,名称取为STRESS-M（或其他也可以），用来表示弯矩,express为*(-0.05)；*(-2.85)再点击add，名称取为STRESS-N,用来表示轴力，express为；再点击add，名称取为STRESS-Q,用来表示剪力，express为.第一个表达式中的-0.05，中性轴的Z向坐标为0.05。STRESS-M设置对话框如下图：（4）查看操作路径：LIST-VALUE LIST-MODEL POINT ，选择MODEL POINT NAME为INTEGER-M。查看弯矩时，在variales to list 中选择自定义，内力种类选择STRESS-M，点击apply即可。查看轴力时，在variales to list 中选择自定义，内力种类选择STRESS-N，点击apply即可。查看剪力时，在variales to l