Abaqus UEL 杆单元

刚刚开始接触UEL，先做了一个简单的一维杆拉伸问题。在这过程中遇到一些很有意思的问题，希望有人能答疑解惑^_^

如果我把坐标输入改成只有x坐标输入的话，算出来的结果就不对，而且力的方向会显示向上！！！
如下输入x和y的坐标后，求得的位移没有问题，但是支反力（RF）为零

.inp文件

*HEADINGLINEAR STATIC PERTURBATION OF BAR ELEMENT
*NODE
1   ,   0.0,    0.0
20  ,   1.0,    0.0
*NGEN, NSET=ALLN
1,20
*USER ELEMENT, NODES=2, TYPE=U1001, PROPERTIES=2, COORDINATES=2, VARIABLES=4
1
*ELEMENT, TYPE=U1001, ELSET=ALLEL
1,  1,  2
2,  2,  3
3,  3,  4
4,  4,  5
5,  5,  6
6,  6,  7
7,  7,  8
8,  8,  9
9,  9,  10
10, 10, 11
11, 11, 12
12, 12, 13
13, 13, 14
14, 14, 15
15, 15, 16
16, 16, 17
17, 17, 18
18, 18, 19
19, 19, 20
*UEL PROPERTY, ELSET=ALLEL
3.0E7, 0.1
**
*STEP,PERTURBATION
*STATIC
*BOUNDARY
1, 1,1, 0.0
*CLOAD
20, 1, 1.0
*END STEP

.for文件

      SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS,1     PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,2     KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF,3     NPREDF,LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,4     PERIOD)
CINCLUDE 'ABA_PARAM.INC'
C      DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),1     SVARS(NSVARS),ENERGY(8),PROPS(3),COORDS(MCRD,NNODE),2     U(NDOFEL),DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2),3     PARAMS(*),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*),4     DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAGS(*),5     JPROPS(*)
C
C     DEFINE PARAMETERS
CPARAMETER (ZERO = 0.D0, ONE = 1.0D0)PARAMETER (PI = ACOS(-ONE))
C
C     CALCULATE LENGTH AND GET PROPERTIES
CEMODULE = PROPS(1)              !YOUNG'S MODULEDIA = PROPS(2)                  !DIAMETERAREA = PI*0.25D0*DIA**2         !AREA OF CROSS-SETION OF BARDD = AREA*EMODULE               ALENG = COORDS(1,2) - COORDS(1,1)!LENGTH OF ELEMENTDD = DD/ALENG
C
C     INITIALIZE RHS AND AMATRX
CDO K1 = 1, NDOFELDO KRHS = 1, NRHSRHS(K1, NRHS) = ZEROENDDODO K2 = 1, NDOFELAMATRX(K2, K1) = ZEROENDDOENDDO
C
C     UPDATE AMATRX AND RHS
CAMATRX(1,1) = ONEAMATRX(1,2) = -ONEAMATRX(2,1) = -ONEAMATRX(2,2) = ONEDO K1 = 1, NDOFELDO K2 = 1, NDOFELAMATRX(K1,K2) = AMATRX(K1,K2)*DDENDDOENDDODO K1 = 1, NDOFELDO K2 = 1, NDOFELRHS(K1,1) = RHS(K1, 1) - AMATRX(K1,K2)*U(K2)ENDDOENDDO
C
C     WRITE RHS AND AMATRX OF ELEMENT TO .DAT FILE (OPTIONAL)
CWRITE(6,*) 'RHS:'WRITE(6,100) (RHS(K1,1),K1=1,2)WRITE(6,*) 'AMATRX'DO K1=1,2WRITE(6,100)(AMATRX(K1,K2),K2=1,2)ENDDO
100   FORMAT(9e12.5)RETURNEND

Abaqus UEL 杆单元相关推荐

integer加1_利用Abaqus UEL开发自定义单元1
1. 简介在Abaqus/Standard模块中,用户可以利用子程序UEL来开发自定义单元,用以实现一些通过Abaqus内置单元无法实现的功能.如果编写恰当,用户自定义单元可以正常使用Abaqus/ ...
ABAQUS UEL
这里写自定义目录标题 ABAQUS UEL 用户子程序UEL: Variable that can be updated Variables passed in for information ABA ...
一个简单的UEL梁单元（上）
出于课题需要,学习研究了abaqus uel二次开发,这里通过一个简单的平面线弹性梁单元的UEL实例来复习巩固一下学到的东西,为后续钢筋和混凝土的本构修改做准备. 主要的参考资料为:王涛老师的< ...
UEL手册中变量解释
ABAQUS UEL子程序用于定义单元接口变量定义待定义变量 RHS(右端量) right-hand-side vectors 等号右端的量 AMATRX(刚度矩阵) 一种包含该元素对雅可比矩 ...
软体机器人与拓扑优化
背景介绍软体机器人在军事以及医疗等方面都具有广泛的应用前景,因此,本推文对相关的研究进行简要概括,并针对现存的问题提出一些自己的思考与解决办法. 磁性软体机器在生物医学领域具有广泛的应用,例 ...
【JY】 ABAQUS子程序UEL的有限元原理与应用
不等待即关注 [简述ABAQUS中UEL子程序] ABAQUS作为成熟的商用有限元软件,可为高级用户提供特定的分析需求.ABAQUS常见的二次开发子程序包括:UMAT.VUMAT.UGENS.UEL ...
matlab有限元分析与应用_专栏 | UEL用户子程序开发步骤—有限元理论基础及Abaqus内部实现方式研究系列20...
作者介绍 snowwave02 博士,高级工程师 snowwave02团队:设计仿真领域的软件开发团队,由软件.机械.物理等专业人员组成,10年以上CAD/CAE软件开发经验,精通Abaqus二次开发 ...
Abaqus二次开发——UEL资源大放送！！！
本次给大家分享的是:UEL资源汇总经过之前给大家分享的UEL入门案例,相信大家对于自定义开发单元有一定的理解,最起码会使用,整体的流程大概有一定的了解. 身处科研一线的你们可能会需要一些比较&quo ...
abaqus最大应力准则怎么用_ANSYS与ABAQUS对比，你选择那个？
来自:CAE交流之家 1．在世界范围内的知名度: 由于Ansys产品进入中国市场早于Abaqus,并且在五年前Ansys的界面是当时最好的界面之一,所以在中国,Ansys软件在用户数量和市场推广度方 ...

Abaqus UEL 杆单元

Abaqus UEL 杆单元

Abaqus UEL 杆单元相关推荐

最新文章

热门文章