请教UEL的使用问题:RHS(right-hand-side)的具体涵义
本帖最后由 zsq-w 于 2011-5-7 13:34 编辑RHS(right-hand-side)被解释成剩余力向量,它的具体涵义是什么
它和单元力向量F,即AMATRX和U的乘积,是什么关系? 这个问题归根结底还是对牛顿-拉夫森算法的理解,RHS说白了就是残余力矢量,是个动态的变量,ABAQUS使用这个变量来求解节点位移的,是n步迭代后N+1步残余力矢量,而DDSDDE是N+1步的一致切线刚度矩阵,所以求得的节点位移是N+1步的,所以叫隐式,因为所有的物理量都是在N+1步获得,必须通过迭代计算才能求解,但是计算代价大,这才叫做隐式。
浙江大学 贝壳航母 RHS是一个力。可以读一下theory manual 2.2.1 Nonlinear solution methods in Abaqus/Standard,RHS对应式2.2.1-2的右端项。 这个问题我当时也困扰了很久
rhs可以解释成单元的内力 如下用Fint表示单元内力,Fout表示单元外力
AMATRX*U=Fint-Fout
对于UEL单元用户需要提供给ABAQUS 当前状态下的刚度矩阵 AMATRX 和单元内力 本帖最后由 wzy_hit 于 2009-5-31 10:26 编辑
谢谢您的回复,
在有的用户子程序中,首先将RHS初始化,是一个全是0的列向量,然后再给它附值,等与 -F,这个用意是什么呢 这个问题总算是明白了,其实就是解非线性方程组的一种方法:
In general procedures Abaqus/Standard solves the overall system of equations by Newton's method 谢谢大虾们的指教啊 受教了,谢谢 受益匪浅 看了半天英文文献 头大了 本帖最后由 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用的也是上一个分析步时的切向刚度。
从另一个角度讲,有限元就是将 无限自由度 用 有限自由度 近似处理的,总得有一个近似准则,使我们能接受这种近似结果,就存在一种比较。
不知道能不能这样理解 那为何在程序中 一般都是: rhs(k1,1)=rhs(k1,1)-k*u ,都没有出现外力f 呢? lxj502 发表于 2012-8-2 16:56 static/image/common/back.gif
那为何在程序中 一般都是: rhs(k1,1)=rhs(k1,1)-k*u ,都没有出现外力f 呢?
因为在abaqus中一般集中外力是cae自己算的,而当有分布力的时候就要自己编写程序 hq052244319 发表于 2012-8-3 15:04 static/image/common/back.gif
因为在abaqus中一般集中外力是cae自己算的,而当有分布力的时候就要自己编写程序 ...
原来如此, 谢谢啦~! 多谢分享 多谢分享1!!! 受益了{:soso_e113:} 受教了,谢谢 很好的帖子,一直觉得RHS为什么没有外力项。但是如果考虑多场耦合的时候右侧的RHS该如何考虑呢? 很好的帖子,一直觉得RHS为什么没有外力项 如果要求外力项,以便检验整体平衡,怎么办?
页:
[1]
2