找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
查看: 6269|回复: 17

回复逍遥处士的帖子-关于圆管内压

  [复制链接]
发表于 2013-10-10 21:08:18 | 显示全部楼层 |阅读模式
$ C/ d7 M# O/ A2 u2 m/ l
7 h+ m" s4 i7 u/ K1 s, m+ n4 [" Q
逍遥处士有关圆管受压的帖子已经在网上有段时间了。他最近让我从软件角度,看看是否跟数学计算一直。现将分析的结果发出来
. k. m% j  Z7 e  M5 j/ A

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册会员

×

点评

尺寸我回去再修改一下。  发表于 2013-10-11 11:11
非常感谢大侠的回应!另外圆管的内径是φ50,外径是φ100,您这个数据跟那个不符。  发表于 2013-10-11 08:37

评分

参与人数 6威望 +54 收起 理由
憨老马 + 20
老鹰 + 30 热心助人,专业精湛!
逍遥处士 + 1
zh39204128 + 1 热心助人,专业精湛!
の小南灬 + 1 热心助人,专业精湛!
mfka + 1 热心助人,专业精湛!

查看全部评分

本帖被以下淘专辑推荐:

  • · 基础|主题: 799, 订阅: 90
回复

使用道具 举报

发表于 2013-10-10 21:23:58 | 显示全部楼层
发表于 2013-10-10 23:25:47 | 显示全部楼层
楼主,求一个第一和第三主应力随半径变化的曲线,这个可以和拉梅方程的解析解比较一下,看一下误差啦~
发表于 2013-10-10 23:37:54 | 显示全部楼层
我也上一个吧,一个mises,一个第一主应力的。6 _# a4 g) v5 ?4 U
! K9 \% T9 Y" y, p  b: {+ `
$ ^' z" w7 ~) Y+ p4 E( k
! k3 f# M* @; k& T( _; f# v
7 Q( @& ^4 w* t

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册会员

×

点评

这网格!呵呵  发表于 2013-10-11 22:26
这个是abaqus分析的,这个分析是基本功,不难,记得abaqus的教程中有类似的例子吗。abaqus 中help文件肯定可以找到类似的例子,参考下就做出来了。  发表于 2013-10-11 18:54
牛,做机加工艺的,完全看不懂您做的这个东西,就知道是应力分析,感觉很高端。。。  发表于 2013-10-11 15:02
牛,做机加工艺的,完全看不懂您做的这个东西,就知道是应力分析,感觉很高端。。。  发表于 2013-10-11 15:02
弹性模量2E5,泊松比0.28,单元类型C3D20R。  发表于 2013-10-10 23:39

评分

参与人数 1威望 +1 收起 理由
Michael0576 + 1

查看全部评分

发表于 2013-10-11 01:24:57 | 显示全部楼层
茉莉素馨 发表于 2013-10-10 23:37
, Q: ?" q0 w- v. k1 \1 `5 L' d. e$ {& U6 y我也上一个吧,一个mises,一个第一主应力的。
4 F! C: }5 E: J: h
单元类型这个东西有什么讲究?/ A0 r6 x. D  i

点评

ansys的单元大概有一百多种,abaqus的大概四百多吧。需要看手册,这种单元是如何定义的,确定模型如何简化,需要什么样的单元。  发表于 2013-10-11 07:27
发表于 2013-10-11 08:09:06 | 显示全部楼层
本帖最后由 茉莉素馨 于 2013-10-11 08:16 编辑
8 @% J3 d. H$ o, @
7 E- K2 W: v: m2 ?, J: C! W# f( u% @  S9 x举个例子说一下吧~; @' U3 J' e% p/ J, O
我用ABAQUS做了一个圆柱体受扭的模型,很简单,材料力学里面已经有解析解了。8 P- J# o5 L& }( I4 ~
参数定义:圆柱尺寸R25 L100,扭矩100Nm,弹性模量2.11e5,泊松比0.28,单元尺寸:6mm。
. U8 C. B+ P/ q: g2 A8 @
/ r$ w9 F4 u8 B6 G% S1 b5 HC3D4做出来的,基本就是错的,不可用。
6 s4 Q7 |$ l5 O" D4 Y+ N单元类型 R 代表减缩积分,会减低一些精度,但是可以加快求解速度。, e" s, i7 L) G7 U+ g- Q
六面体比四面体网格精度高很多,即便是一阶线性单元。5 G5 _) m9 q9 x& E* d& l; x

+ Z/ Y( @) |$ J5 k( K% D5 i  D+ x: @误差:* G7 s. t) s# [# r' {

% b% h, u. G# W3 T( `& N; p) ~6 F3 ^6 X# g' w
上图:4 k: i  s" c8 o2 J% i2 Y8 N
单元:C3D4
& _& e9 y* z7 l
; r+ x1 X9 b4 M" ?: h6 X) [5 W) P0 Y% Z2 m/ w7 S& ?& R* ~! S
C3D10
0 X0 R' c/ [3 ^) ?( d8 A: m( F9 o+ c6 F+ N4 |1 R, o
% S! o! g2 i' [, j5 Q8 T6 |
C3D8R
& Q4 _& |5 K% T  V- }" V: S% p
5 s7 c" ?( W( O2 p) m. w% }2 t! g
C3D8
: y2 P8 Z$ \3 f- Y. z6 x0 [9 {2 T6 g
9 Q, o4 W5 V/ r2 i4 ~8 m! \
C3D20R
2 f, A! j$ K* Q3 p9 c
' O2 \# A+ [: e& Y/ S' [0 A( b- r+ j3 M: q, e' P: N. k  k
C3D20
" I4 ?! P3 T3 u( @/ }) q4 v- T, q) `! l3 e

9 ~* w5 Y. O4 W" ^0 f& w

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册会员

×

点评

99999大侠说的对,确实每个单元都有自己的作用,有限元博大精深,需要学习的还有很多。  发表于 2013-10-12 09:10
跟受力情况和算法也有关系,比方说扭力,有限元其实跟材料力学里面算扭力变形差不多,所以六面体更合适,你换成拉压力,应该区别不大的  发表于 2013-10-11 22:43
每种单元都有存在的理由,不能一棍子打死,比方说你这个C3D4,如果用H和P方法优化之后结果也不会差,速度一样快  发表于 2013-10-11 22:37
厉害,以前自己鼓捣了几个计算现在看来是不行了。  发表于 2013-10-11 20:43
多谢大侠指点! 进军工企业不太可能了,CAE 在设计的时候,就当做一个参考。如果要用好,确实需要很多知识,现在只是自己随便玩玩就好了。  发表于 2013-10-11 20:29
复杂些的有限元,光划分网格就需要几天,模型的简化则需要专门从事有限元这块的人去弄,分析方法的选择。国内搞有限元最好的地方是军工企业,那边的有限元计算精度与生产实践的误差在1%以内。  发表于 2013-10-11 18:59
有限元分析,看起来难,用心学起来不难,真正学好了又是很难。真正的有限元高手,他力学很牛的,可以看到网上中科院或者其他单位找的都是计算力学的博士搞有限元。这个例子模型比较简单,加载条件也不难。  发表于 2013-10-11 18:57
CATIA和abaqus都是达索公司的,用的核心求解器和abaqus差不多,不过功能可能没有abaqus的多,简单问题应该还是有一些可信度的。ansys学校用的多,线性和简单非线性问题还是可以用的,毕竟ansys是cae,更专业一些。  发表于 2013-10-11 08:49
CATIA的有限元和ansys的差别大吗?听说CATIA比ansys还好点是吗?  发表于 2013-10-11 08:34

评分

参与人数 2威望 +21 收起 理由
Michael0576 + 1
老鹰 + 20

查看全部评分

发表于 2013-10-11 08:36:36 | 显示全部楼层
茉莉素馨 发表于 2013-10-10 23:37 8 I! \. F/ O+ A
我也上一个吧,一个mises,一个第一主应力的。

3 |# P3 Y# X0 R  {6 A: X非常感谢两位的参与,觉得技术要大家一起来研究,你研究一段,我研究一段,然后相互一印证,印证上了,那种感觉是很美妙的。一个人研究是毫无乐趣的。
9 D/ F- N; `6 x& B+ Q' c) [& V& c9 H* J0 [4 B3 g1 I: ^
另外请问圆筒内外径是否是φ50和φ100?还有能否分别显示σθ和σr两种应力?解析解算的是σθ(r=25mm) = 167 MPa,σθ(r=50mm) = 67 MPa,这两个数值如果能对应上,说明解析解的值是可以相信的。1 @; f" T+ q- ^5 Y

' j6 V1 M4 d: v' ^5 X+ r4 g这个圆筒受均布内压,和圆柱受均布外压,其应用是很广的,包括液压缸设计,和过盈配合的设计等等。
7 j# R. W& P2 j6 u# m. }$ L

评分

参与人数 1威望 +20 收起 理由
老鹰 + 20

查看全部评分

发表于 2013-10-11 10:38:49 | 显示全部楼层
逍遥处士 发表于 2013-10-11 08:36 6 Q  ^4 T$ I6 Z; X7 i* ^
非常感谢两位的参与,觉得技术要大家一起来研究,你研究一段,我研究一段,然后相互一印证,印证上了,那 ...

. l/ b' R1 S* |$ U3 I" j( [刚刚用有限元算了一下,按照前辈原帖中的参数。
" m' k) T6 B- b1 \外径100,内径50,长度100,受均布压力100mpa,弹性模量2e5mpa,泊松比0.3。
  w/ |) }! F, j% P# h采用半个模型,对称分析。使用C3D20单元,单元数32500,节点数141627。5 E1 d! v( v: o. X& y8 S
内径单元尺寸1.57,外径单元尺寸3.14,径向单元尺寸2mm。6 t9 b* E8 f3 C: R" B8 ~
. i' R) f. O& q5 w' s/ j# r/ H0 o

4 L" `5 ~0 v& ]# \  i下面上图:
) F! o8 K6 x( c- N7 l第一主应力云图,max 167.2935mpa,min 66.64022mpa。
! k; [! Q$ F( W+ p* N% _- R3 E7 [( s% `6 \& G, m' \* W; J& M( V4 C$ t

0 J& [& W/ `. J. x/ Q3 R# k& Y第三主应力云图,max -99.54505mpa,min 0.07029586mpa。
+ s5 p+ B; }& T3 d; [- ^2 J) y  F6 v1 A; D% B5 P2 Z! Q+ t
' J: g) s: M7 a0 Y5 S9 b+ \- P6 P
这个是沿着直径方向,两个应力的变化曲线,
) G7 X0 a5 A: t- x  Y7 u# z, u( \, l5 Z+ Z. g

- X& h6 v( X; B5 I  Q8 F
* M1 Q! E& U, o

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册会员

×

点评

惭愧。大侠若有兴趣,不妨以汝之有限元,与愚之解析解,来一曲琴箫合奏,何如?  发表于 2013-10-11 12:01
解析解还是最可靠的,数值解总还是有些误差的。前辈的数学和力学功底让人十分佩服!  发表于 2013-10-11 11:56
谢谢!大家都是道友。说明解析解还是比较可靠的。  发表于 2013-10-11 11:50
这个应该可以基本验证前辈解析解的正确性了  发表于 2013-10-11 10:41

评分

参与人数 2威望 +2 收起 理由
逍遥处士 + 1
の小南灬 + 1 问题描述清楚,显得很专业!

查看全部评分

 楼主| 发表于 2013-10-11 11:33:25 | 显示全部楼层
这些都是比较基本的问题。

点评

兄弟,简单的问题应用却很大。  发表于 2013-10-11 11:50
 楼主| 发表于 2013-10-11 11:55:03 | 显示全部楼层
大家可以看看理论分析和计算结果。我虽然尺寸不同但是结果却很一致。逍遥可以继续分析下什么原因。

点评

这个问题,我是知道地,好吧。  发表于 2013-10-11 16:08
这个请前辈参考一下原帖里面的解析解的公式,内外径同时扩大相同的倍数的时候,扩大的倍数被约去了。刚刚我也用有限元算了一下,确实是和您最开始给出的结果是一致的。  发表于 2013-10-11 12:13
您需要登录后才可以回帖 登录 | 注册会员

本版积分规则

Archiver|手机版|小黑屋|机械社区 ( 京ICP备10217105号-1,京ICP证050210号,浙公网安备33038202004372号 )

GMT+8, 2025-9-16 05:32 , Processed in 0.057620 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表