beyondboy 发表于 2006-4-24 22:32

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

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

AaronSpark 发表于 2006-4-25 06:37

回复:(beyondboy)[求救]:获得三唯实体单元某一截面...

对于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 发表于 2006-4-25 06:38

回复:(beyondboy)[求救]:获得三唯实体单元某一截面...

该方法可以求曲梁和直梁的内力 <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=00 <BR>FY=00 <BR>FZ=-1000.000-1000-1000 <BR>MX=1685.3001739.81685.301 <BR>MY=5555.6985538.95555.698 <BR>MZ=000

sysh320 发表于 2006-4-28 18:02

这个例子很有帮助!
页: [1]
查看完整版本: [求救]:获得三唯实体单元某一截面的弯距的最好方法是什么?