24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1935  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 中药学调剂 初试324 +3 洋甘菊、 2026-04-10 5/250 2026-04-10 23:59 by 洋甘菊、
[考研] 342电子信息专硕求调剂 +8 你让我怎么荔枝 2026-04-10 9/450 2026-04-10 21:46 by fxue1114
[考研] 332求调剂 +12 蕉蕉123 2026-04-10 12/600 2026-04-10 19:01 by hmn_wj
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 一志愿西南大学生物学学硕344 求生物学相关调剂/生物与医药 +7 超人不会飞@ 2026-04-08 7/350 2026-04-09 09:35 by gong120082
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +11 emokidd 2026-04-08 12/600 2026-04-09 09:24 by wyf236
[考研] 一志愿北京2,材料与化工308求调剂 +17 熊二想上岸 2026-04-04 18/900 2026-04-09 09:13 by leyan1127
[考研] 计算机408|在校多次国家级竞赛获奖|申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 286求调剂 +20 Faune 2026-04-06 20/1000 2026-04-07 11:33 by 诗与自由
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +11 大火山小火山 2026-04-05 11/550 2026-04-06 22:55 by yunlongyang
[考研] 机械专硕274求调剂,不挑专业学校 +6 泛泛2333 2026-04-05 8/400 2026-04-06 18:06 by 泛泛2333
[考研] 一志愿南航,数一英一学硕317求调剂!! +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by 考研学校招点人
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
信息提示
请填处理意见