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

[子程序] 一个UEL程序 有些地方看不懂,跪求大虾开点。。

[复制链接]
发表于 2013-8-29 03:01:38 | 显示全部楼层 |阅读模式 来自 英国
程序如下, 问号的地方看不懂,恳请高手用浅显易懂的语言点拨。。
C 2-NODE BEAM C1 ELEMENT
C WITH GENERALIZED SECTION PROPERTIES AND PRESSURE
C
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,PREDF,NPREDF,
3 LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,PERIOD)
C
INCLUDE ’ABA_PARAM.INC’
C
DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),SVARS(*),PROPS(*),
1 ENERGY(7),COORDS(MCRD,NNODE),U(NDOFEL),DU(MLVARX,*),V(NDOFEL),
2 A(NDOFEL),TIME(2),PARAMS(*),JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*),
3 DDLMAG(MDLOAD,*),PREDEF(2,NPREDF,NNODE),LFLAG(4),JPROPS(*)
C
C DEFINE PARAMETER AND VARIABLE
C
DIMENSION H(1,4)
PARAMETER(ZERO=0.D0,ONE=1.D0,TWO=2.D0,THREE=3.D0,FOUR=4.D0,
1 SIX=6.D0,EIGHT=8.D0,TWELVE=12.D0)
C
C CALCULATE LENGTH AND GET PROPERTIES
C
DL=COORDS(1,2)-COORDS(1,1)
DL2=DL**2
E=PROPS(1)
I=PROPS(2)
P=PROPS(3) !ASSIGIN PRESSURE
YMAX=PROPS(4)
C
C INITIALIZE RHS AND AMATRX
C
NSVINT=NSVARS/2
DO K1=1,4
RHS(K1,1)=ZERO
DO K2=1,4
AMATRX(K1,K2)=ZERO
END DO
END DO
DO K1=1,6
SVARS(K1)=ZERO
END DO
C
C AMATRX
Stress analysis in a beam using the UEL 10
C
AMATRX(1,1)=TWELVE
AMATRX(1,2)=SIX*DL
AMATRX(1,3)=-TWELVE
AMATRX(1,4)=SIX*DL
AMATRX(2,1)=SIX*DL
AMATRX(2,2)=FOUR*DL2
AMATRX(2,3)=-SIX*DL
AMATRX(2,4)=TWO*DL2
AMATRX(3,1)=-TWELVE
AMATRX(3,2)=-SIX*DL
AMATRX(3,3)=TWELVE
AMATRX(3,4)=-SIX*DL
AMATRX(4,1)=SIX*DL
AMATRX(4,2)=TWO*DL2
AMATRX(4,3)=-SIX*DL
AMATRX(4,4)=FOUR*DL2
DO K1=1,4
DO K2=1,4
AMATRX(K1,K2)=AMATRX(K1,K2)*E*I/(DL**3)
END DO
END DO
C
C RHS
C
DO K1=1,4
DO K2=1,4
RHS(K1,1)=RHS(K1,1)-AMATRX(K1,K2)*U(K2)  ?U(K2)是啥意思? 整句话是啥意思?
END DO
END DO
RHS(1,1)=RHS(1,1)-(-P*DL/TWO) ?整句话是啥意思?
RHS(2,1)=RHS(2,1)-(-P*DL2/TWELVE)
RHS(3,1)=RHS(3,1)-(-P*DL/TWO)
RHS(4,1)=RHS(4,1)-(P*DL2/TWELVE)
C
C CALCULATE STRESS
C
CALL STRESS(SVARS,PROPS,DL,DL2,AMATRX,U)
C
C OUTPUT
C
IF(ABS(KINC) .LT. ABS(SVARS(5))) THEN
KINC=SVARS(5) ?SVARS(5) 是啥意思?整句话啥意思?
KSTEP=JELEM  ?整句话啥意思?
END
IF
IF(ABS(KINC) .LT. ABS(SVARS(6))) THEN
KINC=SVARS(6)
KSTEP=JELEM+1 ?整句话啥意思?
Stress analysis in a beam using the UEL 11
END IF
OPEN(UNIT=16,FILE=’C:\TEMP\STRESS.DAT’,STATUS=’UNKNOWN’,
$ FORM=’FORMATTED’)
IF(JELEM .EQ. 1) THEN
WRITE(16,*) ’-------------------------NEXT STEP
$ -------------------------’
END IF
WRITE(16,*) ’ELEMENT NUMBER’,JELEM
WRITE(16,*) ’NODE’,JELEM,’FORCE’,SVARS(1)
WRITE(16,*) ’NODE’,JELEM,’MOMENT’,SVARS(2)
WRITE(16,*) ’NODE’,JELEM+1,’FORCE’,SVARS(3)
WRITE(16,*) ’NODE’,JELEM+1,’MOMENT’,SVARS(4)
WRITE(16,*) ’NODE’,JELEM,’NORMAL STRESS’,SVARS(5)
WRITE(16,*) ’NODE’,JELEM+1,’NORMAL STRESS’,SVARS(6)
WRITE(16,*) ’MAX NORMAL STRESS’,KINC,’AT NODE’,KSTEP
WRITE(16,*) ’ ’
C
C END
C
RETURN
END
发表于 2017-4-12 11:00:28 | 显示全部楼层 来自 陕西
Simdroid开发平台
楼主,这个程序我也在做,刚开始学习,同时建立了inp文件,但是始终都运行不了,可能是inp文件了问题,程序应该没什么问题,能帮帮小弟我,谢谢

本帖子中包含更多资源

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

×
回复 不支持

使用道具 举报

 楼主| 发表于 2017-6-9 23:45:44 | 显示全部楼层 来自 英国
旋风小子 发表于 2017-4-12 03:00
楼主,这个程序我也在做,刚开始学习,同时建立了inp文件,但是始终都运行不了,可能是inp文件了问题,程序 ...

看看表头有木有问题。不同版本的Abaqus表头有差异!
回复 不支持

使用道具 举报

发表于 2017-6-10 09:29:45 | 显示全部楼层 来自 陕西
zhuxinyao 发表于 2017-6-9 23:45
看看表头有木有问题。不同版本的Abaqus表头有差异!

我找出问题了,已经可以运行了,
回复 不支持

使用道具 举报

发表于 2017-7-26 11:10:58 | 显示全部楼层 来自 北京
旋风小子 发表于 2017-6-10 09:29
我找出问题了,已经可以运行了,

是什么问题,要改表头么
回复 不支持

使用道具 举报

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

本版积分规则

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

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

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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