24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 1702  |  回复: 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 的主题更新
信息提示
请填处理意见