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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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

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

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

AaronSpark · 发表于 2006-4-25 06:38

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

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

这个例子很有帮助！

账号		自动登录	找回密码
密码			我要加入

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

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

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