24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1927  |  回复: 9

11hours

铜虫 (小有名气)

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

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

位移



y方向应力
回复此楼
已阅   回复此楼   关注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

铜虫 (小有名气)


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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本9 一志愿西工大085601 324求调剂 +3 wysyjs25 2026-04-10 3/150 2026-04-10 15:19 by sophia_93
[考研] 材料与化工调剂 10+9 下一站上岸@ 2026-04-10 23/1150 2026-04-10 14:56 by Kilig0317
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +17 努力奋斗112 2026-04-08 17/850 2026-04-10 12:17 by pengliang8036
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +19 努力奋斗112 2026-04-04 20/1000 2026-04-10 12:15 by pengliang8036
[考研] 求调剂 +9 翩翩一书生 2026-04-09 9/450 2026-04-10 11:43 by wp06
[考研] 调剂 +15 小张ZA 2026-04-10 16/800 2026-04-10 09:56 by 690616278
[考研] 求调剂 +15 张zic 2026-04-05 16/800 2026-04-10 08:12 by kangsm
[考研] 求调剂 +5 archer.. 2026-04-09 7/350 2026-04-09 22:18 by lbsjt
[考研] 一志愿085502,267分求调剂 +10 再忙也要吃饭啊 2026-04-08 11/550 2026-04-09 19:51 by gong120082
[考研] 337求调剂 +4 Gky09300550, 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
[考研] 一志愿厦大生物学332求调剂 +10 池池池池池池 2026-04-08 10/500 2026-04-09 17:10 by 独醉梦孤城
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 298求调剂 +4 manman511 2026-04-05 4/200 2026-04-08 16:50 by tjzhao
[考研] 298求调剂 +4 残荷新柳 2026-04-07 4/200 2026-04-07 23:02 by lbsjt
[考研] 一志愿南科大生物学297分,求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 081700学硕,323分,一志愿中国海洋大学求调剂学校 +19 披星河 2026-04-04 19/950 2026-04-07 15:00 by 上岸快快
[考研] 求调剂 +11 xzghyuj 2026-04-04 11/550 2026-04-06 11:49 by lijunpoly
[考研] 0817化学工程与技术求调剂,一志愿中海洋319 +14 lv945 2026-04-04 14/700 2026-04-06 10:20 by 蓝云思雨
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
信息提示
请填处理意见