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

[子程序] 请教UEL的使用问题:RHS(right-hand-side)的具体涵义

[复制链接]
发表于 2009-5-21 15:53:02 | 显示全部楼层 |阅读模式 来自 黑龙江哈尔滨
本帖最后由 zsq-w 于 2011-5-7 13:34 编辑

RHS(right-hand-side)被解释成剩余力向量,它的具体涵义是什么
它和单元力向量F,即AMATRX和U的乘积,是什么关系?
发表于 2013-10-1 23:08:44 | 显示全部楼层 来自 浙江杭州
Simdroid开发平台
这个问题归根结底还是对牛顿-拉夫森算法的理解,RHS说白了就是残余力矢量,是个动态的变量,ABAQUS使用这个变量来求解节点位移的,是n步迭代后N+1步残余力矢量,而DDSDDE是N+1步的一致切线刚度矩阵,所以求得的节点位移是N+1步的,所以叫隐式,因为所有的物理量都是在N+1步获得,必须通过迭代计算才能求解,但是计算代价大,这才叫做隐式。

浙江大学 贝壳航母
回复 4 不支持 0

使用道具 举报

发表于 2014-4-10 17:43:48 | 显示全部楼层 来自 挪威
RHS是一个力。可以读一下theory manual 2.2.1 Nonlinear solution methods in Abaqus/Standard,RHS对应式2.2.1-2的右端项。
回复 1 不支持 0

使用道具 举报

发表于 2009-5-26 14:07:43 | 显示全部楼层 来自 浙江杭州
这个问题我当时也困扰了很久
rhs可以解释成单元的内力 如下用Fint表示单元内力,Fout表示单元外力
AMATRX*U=Fint-Fout
对于UEL单元用户需要提供给ABAQUS 当前状态下的刚度矩阵 AMATRX 和单元内力

评分

1

查看全部评分

回复 1 不支持 0

使用道具 举报

 楼主| 发表于 2009-5-27 08:31:56 | 显示全部楼层 来自 黑龙江哈尔滨
本帖最后由 wzy_hit 于 2009-5-31 10:26 编辑

谢谢您的回复,
在有的用户子程序中,首先将RHS初始化,是一个全是0的列向量,然后再给它附值,等与 -F  ,这个用意是什么呢
回复 不支持

使用道具 举报

 楼主| 发表于 2009-5-31 16:16:32 | 显示全部楼层 来自 黑龙江哈尔滨
这个问题总算是明白了,其实就是解非线性方程组的一种方法:
In general procedures Abaqus/Standard solves the overall system of equations by Newton's method

评分

1

查看全部评分

回复 不支持

使用道具 举报

发表于 2010-5-18 15:27:50 | 显示全部楼层 来自 湖北武汉
谢谢大虾们的指教啊
回复 不支持

使用道具 举报

发表于 2010-6-17 16:54:46 | 显示全部楼层 来自 天津
受教了,谢谢
回复 不支持

使用道具 举报

发表于 2011-1-5 19:47:29 | 显示全部楼层 来自 浙江杭州
受益匪浅 看了半天英文文献 头大了
回复 不支持

使用道具 举报

发表于 2011-5-7 12:59:27 | 显示全部楼层 来自 北京
本帖最后由 563999954 于 2011-5-7 13:07 编辑

1# wzy_hit
r=f-I=f-ku
因为首先满足力平衡方程,所以软件内部计算 单元 内力 Ki du  和 施加在单元上的外载 P 之间的偏差(或者叫残差)R——R= P-I ,然后将残差R与收敛准则下的容许值(软件内部定义的)进行比较,这样做就是保证,是每一个加载点位于(F-u曲线上),但愿是符合力平衡方程的,一般的非线性分析是达不到的,需要不断的接近,K用的也是上一个分析步时的切向刚度。
从另一个角度讲,有限元就是将 无限自由度 用 有限自由度 近似处理的,总得有一个近似准则,使我们能接受这种近似结果,就存在一种比较。
不知道能不能这样理解

评分

1

查看全部评分

回复 不支持

使用道具 举报

发表于 2012-8-2 16:56:19 | 显示全部楼层 来自 广东广州
那为何在程序中 一般都是: rhs(k1,1)=rhs(k1,1)-k*u ,  都没有出现外力f 呢?  
回复 不支持

使用道具 举报

发表于 2012-8-3 15:04:26 | 显示全部楼层 来自 武汉大学
lxj502 发表于 2012-8-2 16:56
那为何在程序中 一般都是: rhs(k1,1)=rhs(k1,1)-k*u ,  都没有出现外力f 呢?

因为在abaqus中一般集中外力是cae自己算的,而当有分布力的时候就要自己编写程序
回复 不支持

使用道具 举报

发表于 2012-8-6 09:42:01 | 显示全部楼层 来自 广东广州
hq052244319 发表于 2012-8-3 15:04
因为在abaqus中一般集中外力是cae自己算的,而当有分布力的时候就要自己编写程序 ...

原来如此, 谢谢啦~!
回复 不支持

使用道具 举报

发表于 2013-8-19 20:56:56 | 显示全部楼层 来自 浙江杭州
多谢分享
回复 不支持

使用道具 举报

发表于 2015-10-21 17:54:43 | 显示全部楼层 来自 陕西西安
受益了{:soso_e113:}
回复 不支持

使用道具 举报

发表于 2015-12-7 14:00:54 | 显示全部楼层 来自 上海
受教了,谢谢
回复 不支持

使用道具 举报

发表于 2016-2-29 16:34:10 | 显示全部楼层 来自 欧洲
很好的帖子,一直觉得RHS为什么没有外力项。但是如果考虑多场耦合的时候右侧的RHS该如何考虑呢?
回复 不支持

使用道具 举报

发表于 2016-7-17 11:41:32 | 显示全部楼层 来自 北京
很好的帖子,一直觉得RHS为什么没有外力项
回复 不支持

使用道具 举报

发表于 2016-10-29 23:25:30 | 显示全部楼层 来自 江苏南京
如果要求外力项,以便检验整体平衡,怎么办?
回复 不支持

使用道具 举报

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

本版积分规则

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

GMT+8, 2024-4-28 02:05 , Processed in 0.076382 second(s), 17 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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