24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1932  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工329分求调剂 +8 叶zilin 2026-04-10 8/400 2026-04-10 23:26 by 314126402
[考研] 0703化学 +32 妮妮ninicgb 2026-04-04 37/1850 2026-04-10 22:18 by bljnqdcc
[考研] 吉大计算机技术331分,英语六级,求调剂 +3 峰峰021116 2026-04-09 3/150 2026-04-10 20:01 by chemisry
[硕博家园] 0856材料化工求调剂,一志愿211,初试成绩349 +5 江淮北月 2026-04-05 5/250 2026-04-10 16:26 by 高维春
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 282,电气工程专业,求调剂,不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 调剂 +4 青灯不负 2026-04-09 4/200 2026-04-10 09:55 by LHGeng
[考研] 301求调剂 +4 XYPLR 2026-04-05 6/300 2026-04-10 09:17 by Delta2012
[考研] 087100初试311求调剂 +3 任雅琴 2026-04-09 3/150 2026-04-09 22:42 by lbsjt
[考研] 288求调剂 +15 没有答案_ 2026-04-05 15/750 2026-04-09 10:22 by 5268321
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 22408 调剂材料 +7 我叫ez 2026-04-06 8/400 2026-04-07 17:12 by 蓝云思雨
[考研] 295求调剂 +18 xndjjj 2026-04-04 19/950 2026-04-07 11:02 by wangjy2002
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +11 大火山小火山 2026-04-05 11/550 2026-04-06 22:55 by yunlongyang
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
[考研] 一志愿河北工业大学材料工程,初试344求专硕调剂 +6 15933906766 2026-04-05 6/300 2026-04-06 13:21 by 无际的草原
[考研] 284求调剂 +7 徐同学_001 2026-04-04 13/650 2026-04-05 17:19 by yulian1987
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
信息提示
请填处理意见