找回密码
 注册
Simdroid-非首页
查看: 146|回复: 0

[子程序] 混凝土弹塑性本构 Umat 问题请教!

[复制链接]
发表于 2019-3-8 07:05:46 | 显示全部楼层 |阅读模式 来自 美国
程序来源
1.《ABAQUS在隧道与地下工程中的应用》-陈卫忠 第10章
2.基于修正Mohr-Coulomb准则的弹塑性本构模型及其数值实施 -贾善坡,2010
问题:

  • Line141185291325386,参数m的取值不一样,书中给出的程序可以理解为屈服函数和塑性势函数m取值不同,分别为0.050.02,我的理解是既然屈服函数和塑性势函数一致的话,参数也应该一致,即m取值应相同。或是我理解有错误,只是形状相似,函数并不完全一致。此外,第386行,为什么除以摩擦角的正切?
  • 文中给出的屈服函数F和塑性势函数G并非完全一致,相差一个c*cos(PHI),为什么?
  • 238行,求λ变化率时,屈服函数对粘聚力c的偏导数直接取了1,为什么?
  • 程序中似乎并没有进行奇异点的处理
  • 668行,DDFDDS(K1)/RB(K1),此处为什么用除法?按照文献中的公式应该用乘法。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 16:08 , Processed in 0.027287 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表