24小时热门版块排行榜    

查看: 1893  |  回复: 9

11hours

铜虫 (小有名气)

[求助] 自己写的一个小程序,出现问题不知道怎么解决

模型,左边X方向约束,底边Y方向约束。上边均布载荷。左下角有个小洞。
位移出来都好好的。
如图是Y方向位移,
可Y方向应力怎么是这样?
边界上的应力值怎么不对。应力值使用高斯点外推出来的。
是我求应力的过程出错,还是之前分析过程出错?

位移



y方向应力
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

11hours

铜虫 (小有名气)


11hours: 回帖置顶 2012-08-02 06:07:15
xiegangmai: 金币+1, 鼓励讨论交流! 2012-08-02 22:54:05
搞定了,但还是没搞清楚在哪儿出了问题。
原来是想当然的觉得应力磨平要精确,所以没有检查用本构直接解应力的结果如何。
试了各种磨平后应力结果依然不对,但直接解出来的远离孔边的应力就正确了。
孔边的应力直接解应该是差的远的,不过不敢确定应力磨平的结果是否能精确点。怎样评估呢?
10楼2012-08-02 06:07:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

11hours

铜虫 (小有名气)

边界上的应力问题解决了,用单元面积加权平均了一下,边界上的应力和内部差不多了。

可现在的问题是,内部应力比正确值大很多,有两倍多,大部分的应力应该跟上边界的均布力一样才对呀
2楼2012-07-31 16:49:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11hours

铜虫 (小有名气)

位移解跟abaqus相差了10%。是我程序错了还是误差?
3楼2012-07-31 16:51:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

onesupeng

金虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 专家考核, 谢谢应助! 2012-08-02 22:54:19
固体的话差10%的误差确实难以接受

多研究不同网格、不同时间步大小等,可能会容易把程序的问题找出来
长期招收博士生,参见http://fsl-unsw.com
4楼2012-08-01 02:07:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11hours

铜虫 (小有名气)

内容已删除
5楼2012-08-01 03:23:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11hours

铜虫 (小有名气)

放假了,人很少啊
6楼2012-08-01 03:25:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11hours

铜虫 (小有名气)

好像明白点了,刚度矩阵贴错了,应该是
  157106179869.04    ,
  49121509505.371    ,  74301382611.676   
-135867951432.10    , -8460526173.1302    ,  152819724760.43   
  11443319980.716    , -8228229559.4999    , -50360729615.398    ,  75676593946.869   
-156974866704.79    , -49121509505.371    ,  135738532952.75    , -11443319980.716   
  157098907787.68    ,
-49121509505.371    , -74170069447.429    ,  8460526173.1302    ,  8098811080.1506   
  49121509505.371    ,  74294110530.320   
  135736638267.85    ,  8460526173.1302    , -152690306281.08    ,  50360729615.398   
-135862574035.64    , -8460526173.1302    ,  152816242048.87   
-11443319980.716    ,  8096916395.2527    ,  50360729615.398    , -75547175467.520   
  11443319980.716    , -8222852163.0419    , -50360729615.398    ,  75673111235.309
刚好36个,应该是按列排的对称矩阵。如果是这样的话,我计算的刚度矩阵就错了,可是这样的问题刚度矩阵计算应该都一样吧,怎么会错呢?abaqus单元选的是cpe4
7楼2012-08-01 04:08:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11hours

铜虫 (小有名气)

还是这么少人啊。
位移解对了,之前跟abaqus不一样有两个原因
1、形函数少写了一个负号  ——!
欲哭无泪啊。
2、abaqus的刚度矩阵居然代入的是平面应力本构,为什么,为什么?我把本构调成平面应力的之后才和abaqus的单刚差不多。出来的位移解完全一致。可abaqus当中我明明选的是平面应变单元啊,难道还要设置什么东东约束z方向吗?
8楼2012-08-01 16:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

11hours

铜虫 (小有名气)

不能编辑真的有点恼火,我不是刷屏。。。
应力解依然不对,力边界的点应力依然高出正解1.5倍,即正解的2.5倍。怎么回事?
9楼2012-08-01 16:06:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 11hours 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 材料工程(专)一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 278求调剂 +6 烟火先于春 2026-03-17 6/300 2026-03-21 01:57 by JourneyLucky
[考研] 307求调剂 +10 冷笙123 2026-03-17 10/500 2026-03-21 01:54 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 324求调剂 +5 lucky呀呀呀鸭 2026-03-20 5/250 2026-03-20 22:30 by 促天成
[考研] 085600材料与化工 +8 安全上岸! 2026-03-16 8/400 2026-03-20 22:13 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
信息提示
请填处理意见