[求救]:获得三唯实体单元某一截面的弯距的最好方法是什么?

beyondboy

我用65号单元做的一个模型,请问如何获得比较精确的截面弯距!```<BR>希望大侠指点~~~~~

AaronSpark

对于3D实体单元的内力计算，可采用面操作（应该是V8以上吧，没有细考！）。例如用悬臂梁采用3D实体单元，其某个截面的的内力计算如下，并具有详细解释。这种方法较原来的更方便，大家不放一用。 <BR>!****************************************************** <BR>!EX4.31 面操作及悬臂梁的内力计算 <BR>FINISH$/CLEAR$/PREP7 <BR>ET,1,SOLID95$MP,EX,1,2E11$MP,PRXY,1,0.3        !定义单元类型、材料特性 <BR>BLC4,2,3,0.2,0.3,4$DA,2,ALL$FK,1,FY,-2E4$FK,3,FY,-2E4  !创建几何模型、加约束和荷载 <BR>FK,3,FX,0.8E4$FK,4,FX,0.8E4$SFA,1,1,PRES,1E6      !施加荷载 <BR>ESIZE,0.05$VMESH,ALL$FINISH$/SOLU$SOLVE      !生成有限元模型并求解 <BR>FINISH$/POST1              !进入后处理层 <BR>WPOFF,,,2$SUCR,SUZ2,CPLANE,3      !移动工作平面、创建面SUZ2 <BR>SUMAP,MYSX,S,X$SUMAP,MYSY,S,Y    !映射X和Y方向应力 <BR>SUMAP,MYSZ,S,Z$SUMAP,MYSXY,S,XY  !映射Z和XY方向应力 <BR>SUMAP,MYSYZ,S,YZ$SUMAP,MYSXZ,S,XZ  !映射YZ和XZ方向应力 <BR>SUPL,SUZ2$SUPL,SUZ2,MYSZ        !显示面本身、面项MYSZ <BR>SUPL,SUZ2,MYSYZ$SUPL,SUZ2,MYS    !显示面项MYSYZ，矢量显示应力 <BR>SUPR,ALL,MYSZ              !列表显示MYSZ面项 <BR>SUEVAL,XFORCE,MYSXZ,INTG  !求截面上FX，理论结果为-16000,误差1% <BR>SUEVAL,YFORCE,MYSYZ,INTG  !求截面上FY，理论结果为40000，误差0.5% <BR>SUEVAL,ZFORCE,MYSZ,INTG    !求截面上FZ，理论结果为-6000，没有误差 <BR>SUEVAL,MYA,DA,SUM      !求截面面积并赋给MYA变量 <BR>SUEVAL,MYYA,GCY,INTG      !求关于X轴的面积矩并赋给变量MYYA <BR>MYYA=MYYA/MYA        !得到面积重心到X轴的距离=面积矩/面积 <BR>SUEVAL,MYXA,GCX,INTG    !求关于Y轴的面积矩并赋给变量MYXA <BR>MYXA=MYXA/MYA        !得到面积重心到Y轴的距离=面积矩/面积 <BR>SUCALC,SZGCY,MYSZ,MULT,GCY  !计算MYSZ×GCY，并赋给面项SZGCY <BR>SUEVAL,MX1,SZGCY,INTG    !对面项SZGCY在面上积分得到MX1 <BR>SUCALC,SZGCX,MYSZ,MULT,GCX  !计算MYSZ×GCZ，并赋给面项SZGCX <BR>SUEVAL,MY1,SZGCX,INTG    !对面项SZGCX在面上积分得到MY1 <BR>!上述弯矩基于总体直角坐标系原点而言的，应对面积重心取矩，将内力简化到面积重心上 <BR>MX1=MX1-ZFORCE*MYYA   !理论结果为80000，误差为0.08% <BR>MY1=MY1-ZFORCE*MYXA  !理论结果-32000，误差为0.2% <BR>

AaronSpark

该方法可以求曲梁和直梁的内力 <BR>-------------------------------------------------------------------------------------------- <BR>文件说明 <BR>计算3根相同的弯梁(梁1用beam44单元，梁2，3用solid45单元)，分别在竖向集中力，纵向集中力， <BR>竖向均布力,重力作用下响应，利用nforce来计算实体单元的内力，并与手算内力，梁单元内力作了比较 <BR>，计算表明用nforce命令计算的实体截面的内力和手算结果相同，计算的beam44单元的内力与手算结果 <BR>基本一致(因为计算中曲梁用直梁代替，直梁的单元的轴向和曲梁单元的轴向有一定夹角，如果把梁单元 <BR>划分得足够小，结果将和手算结果相同) <BR>  <BR>计算内力方法总结 <BR>1:首先选中一个截面上的节点，并选中与之相关的一边的单元，将内力计算点和局部坐标均移到形心上 <BR>使用下面命令,列出内力，内力方向同局部坐标方向一致 <BR>cmsel,s,nodejs  !选择计算内力截面上的点 <BR>NFORCE,ALL  !列出计算内力截面上的点的合力相对局部坐标的合力 <BR>  <BR>内力符合及方向说明,设Z正向和Y正向为梁的上方 <BR>  <BR>对左端截面 <BR>FX "-"表示  及X正向 <BR>FY "-"表示  及Y正向  <BR>FZ "-"表示  ，及Z正向 <BR>MX "-"表示Y→Z  及X正向 <BR>MY "-"表示上端受拉，  及Y正向 <BR>MZ "-"表示上端受压，  及Z负向 <BR>  <BR>右截面和左截面相反 <BR>  <BR>计算结果说明 <BR>端部集中力作用1000N <BR>在梁L/4长度上的内力计算,取截面右端单元计算内力 <BR>        用nforce命令  beam44  手算 <BR>  FX  =  0  0 <BR>  FY  =  0  0 <BR>  FZ  =  -1000.000  -1000  -1000 <BR>  MX  =  1685.300  1739.8  1685.301 <BR>  MY  =  5555.698  5538.9  5555.698 <BR>  MZ  =  0  0  0

sysh320

这个例子很有帮助！