有没有大佬帮忙看看这个程序怎么修改,施加载荷的时候位移载荷不会出错,但是压力载荷却总是出错了?

浏览:1199

FINISH

/CLEAR

 

/PREP7                                   !进入前处理器

ET, 1, PLANE182                           !选择单元类型

ET, 2, TARGE169

ET, 3, CONTA171

 

KEYOPT, 1, 3, 0                  !设置单元关键字(1为第一个单元,3为K3,1为所选项(Plane stress=0; Axisymmetric=1; Plane strain=3))

 

MP, EX, 1, 2.06E5                          !输入材料弹性模量(定义材料属性)

MP, PRXY, 1, 0.28                          !输入材料泊松比

TB, BKIN, 1                               !输入材料应力应变特性

TBDATE, 1, 450                            !输入材料屈服强度

TBDATE, 2, 0.8E5                           !输入材料剪切模量

MP, MU, 1, 0.15            !输入摩擦系数(润滑摩擦条件下:钢-钢=0.15,钢-橡胶=0.30)

 

MP, EX, 2, 2.06E5                          !输入材料弹性模量(定义材料属性)

MP, PRXY, 2, 0.28                          !输入材料泊松比

TB, BKIN, 2                               !输入材料应力应变特性

TBDATE, 1, 835                            !输入材料屈服强度

TBDATE, 2, 0.8E5                           !输入材料剪切模量

MP, MU, 2, 0.15            !输入摩擦系数(润滑摩擦条件下:钢-钢=0.15,钢-橡胶=0.30)

 

K, 1, 5, 83, 0                            !创建关键点

K, 2, -80, 83, 0

K, 3, -80, 67.5, 0

K, 4, 5, 67.5, 0

K, 5, 0, 52, 0

K, 6, 0, 12, 0                        

K, 7, -72.26, 12, 0                        

K, 8, -72.26, 28, 0                        

K, 9, -6, 52, 0                        

 

K, 41, -72.26, 65.64, 0                        

K, 42, -72.26, 28, 0                        

K, 43, -6, 52, 0                        

K, 44, -6, 65.64, 0

K, 45, -9.5, 67.64, 0

K, 46, -9.5, 65.64, 0

K, 47, -13, 67.64, 0                        

K, 48, -13, 65.64, 0                        

K, 49, -16.5, 67.64, 0                        

K, 50, -16.5, 65.64, 0                        

K, 51, -20, 67.64, 0                        

K, 52, -20, 65.64, 0                        

K, 53, -23.5, 67.64, 0                        

K, 54, -23.5, 65.64, 0                        

K, 55, -27, 67.64, 0                        

K, 56, -27, 65.64, 0                        

K, 57, -30.5, 67.64, 0                        

K, 58, -30.5, 65.64, 0                        

K, 59, -34, 67.64, 0                        

K, 60, -34, 65.64, 0                        

K, 61, -37.5, 67.64, 0                         

K, 62, -37.5, 65.64, 0                        

K, 63, -41, 67.64, 0                        

K, 64, -41, 65.64, 0                        

K, 65, -44.5, 67.64, 0                        

K, 66, -44.5, 65.64, 0                        

K, 67, -48, 67.64, 0                        

K, 68, -48, 65.64, 0                        

K, 69, -51.5, 67.64, 0                        

K, 70, -51.5, 65.64, 0                        

K, 71, -55, 67.64, 0                        

K, 72, -55, 65.64, 0                         

K, 73, -58.5, 67.64, 0                        

K, 74, -58.5, 65.64, 0                        

K, 75, -62, 67.64, 0                        

K, 76, -62, 65.64, 0                        

K, 77, -65.5, 67.64, 0                        

K, 78, -65.5, 65.64, 0                        

K, 79, -69, 67.64, 0                        

K, 80, -69, 65.64, 0                        

 

* do, i, 41, 79

a = i

b = a + 1

lstr, a, b

* enddo

L, 80, 41

 

* do, j, 4, 39

c = j

d = j + 1

LFILLT, c, d, 0.5                                  !生成倒角线

* enddo

LCOMB, 4, 41

 

* do, k, 5, 39

e = k

f = k + 37

LCOMB, e, f                                      !合并线段

LCOMB, 4, e

* enddo

LCOMB, 4, 40

 

!LSEL, S, , , 1, 40, 1

!AL, ALL                                     !由所有线生成面-锁定滑块

 

L, 1, 2                                        !由关键点生成线段-42

L, 2, 3                                        !43

L, 3, 4                                        !44

L, 4, 1                                        !45

L, 5, 6                                        !46

L, 6, 7                                        !47

L, 7, 8                                        !48

L, 8, 9                                        !49

L, 9, 5                                        !50

 

AL, 1, 2, 3, 4                              !由线生成面

AL, 5, 6, 7, 8  

AL, 9, 10, 11, 12, 13

 

/PNUM, KP, 1                               !显示关键点编号

/PNUM, LINE, 1                             !显示线段编号

/PNUM, AREA, 1                            !显示面编号

 

APLOT                                     !显示面

ALLSEL                                    !选择所有实体

 

 

有没有大佬帮忙看看这个程序怎么修改,施加载荷的时候位移载荷不会出错,但是压力载荷却总是出错了?的图1有没有大佬帮忙看看这个程序怎么修改,施加载荷的时候位移载荷不会出错,但是压力载荷却总是出错了?的图2MAT, 2                                   !锁定滑块

REAL, 3                                  !实常数

TYPE, 1                                  !单元类型

MSHKEY, 0                              !自由网格化分

SMRTSIZE, 2                             !智能网格划分控制

AESIZE, 2                                !对面上线段设置网格划分控制

AMESH, 1                                !面上生成面单元 - 对面进行网格划分

 

MAT, 1                                    !管壁

REAL, 3

TYPE, 1

MSHKEY, 0

SMRTSIZE, 1

AESIZE, 1

AMESH, 2

 

MAT, 2                                    !承压头

REAL, 3

TYPE, 1

MSHKEY, 0

SMRTSIZE, 1

AESIZE, 1

AMESH, 3

 

有没有大佬帮忙看看这个程序怎么修改,施加载荷的时候位移载荷不会出错,但是压力载荷却总是出错了?的图3EPLOT                                          !显示单元

ALLSEL                                         !选择所有实体

 

 

!开始定义接触单元

/COM, CONTACT PAIR CREATION - START

CM,_NODECM,NODE

CM,_ELEMCM,ELEM

CM,_KPCM,KP

CM,_LINECM,LINE

CM,_AREACM,AREA

CM,_VOLUCM,VOLU

/GSAV,cwz,gsav,,temp

MP, MU, 1, 0.15

R, 5

REAL, 5

ET, 6, 169

ET, 7, 172

!KEYOPT, 7, 9, 0

!KEYOPT, 7, 10, 2

KEYOPT, 7, 12, 4

R, 5, , , 0.1, 0.1

! Generate the target surface  

LSEL,S,,, 2

CM,_TARGET,LINE

TYPE,6

MAT, 2

NSLL,S,1

ESLN,S,0

ESURF

CMSEL,S,_ELEMCM

! Generate the contact surface 

LSEL,S,,, 12

CM,_CONTACT,LINE

TYPE,7 

MAT, 2

NSLL,S,1

ESLN,S,0

ESURF  

ALLSEL 

ESEL,ALL

ESEL,S,TYPE,,6 

ESEL,A,TYPE,,7 

ESEL,R,REAL,,5 

/PSYMB,ESYS,1  

/PNUM,TYPE,1

/NUM,1 

EPLOT  

ESEL,ALL

ESEL,S,TYPE,,6 

ESEL,A,TYPE,,7 

ESEL,R,REAL,,5 

CMSEL,A,_NODECM

CMDEL,_NODECM  

CMSEL,A,_ELEMCM

CMDEL,_ELEMCM  

CMSEL,S,_KPCM  

CMDEL,_KPCM

CMSEL,S,_LINECM

CMDEL,_LINECM  

CMSEL,S,_AREACM

CMDEL,_AREACM  

CMSEL,S,_VOLUCM

CMDEL,_VOLUCM  

/GRES,cwz,gsav 

CMDEL,_TARGET  

CMDEL,_CONTACT 

/COM, CONTACT PAIR CREATION - END  

 

/COM, CONTACT PAIR CREATION - START

CM,_NODECM,NODE

CM,_ELEMCM,ELEM

CM,_KPCM,KP

CM,_LINECM,LINE

CM,_AREACM,AREA

CM,_VOLUCM,VOLU

/GSAV,cwz,gsav,,temp

MP,MU, 1, 0.15

R, 6

REAL, 6 

ET, 8, 169

ET, 9, 172

!KEYOPT, 9, 9, 0

!KEYOPT, 9, 10, 2

R,6, , , 0.1, 0.1

RMORE, 

RMORE, , 0

RMORE, 0

! Generate the target surface  

LSEL,S,,,4

CM,_TARGET,LINE

TYPE,8 

MAT, 1

NSLL,S,1

ESLN,S,0

ESURF

CMSEL,S,_ELEMCM

! Generate the contact surface 

LSEL,S,,, 7

CM,_CONTACT, LINE

TYPE,9

MAT, 2

NSLL,S,1

ESLN,S,0

ESURF  

ALLSEL 

ESEL, ALL

ESEL, S, TYPE,, 8 

ESEL, A, TYPE,, 9 

ESEL, R, REAL,, 6 

/PSYMB, ESYS, 1  

/PNUM, TYPE, 1

/NUM,1 

EPLOT  

ESEL,ALL

ESEL,S,TYPE,,8 

ESEL,A,TYPE,,9 

ESEL,R,REAL,,6 

CMSEL,A,_NODECM

CMDEL,_NODECM  

CMSEL,A,_ELEMCM

CMDEL,_ELEMCM  

CMSEL,S,_KPCM  

CMDEL,_KPCM

CMSEL,S,_LINECM

CMDEL,_LINECM  

CMSEL,S,_AREACM

CMDEL,_AREACM  

CMSEL,S,_VOLUCM

CMDEL,_VOLUCM  

/GRES,cwz,gsav 

CMDEL,_TARGET  

CMDEL,_CONTACT 

/COM, CONTACT PAIR CREATION - END  

 

 

/SOLU                                          !进入求解器

ANTYPE, STATIC                                !选择求解类型

NLGEOM, ON

AUTOT, ON

TIME, 1                                        !指定计算终止时间

NSUBST, 25000, 100000, 2500                      !指定迭代子步数

LNSRCH, ON                                   !激活线性搜索

PRED, 1                                        !激活自由度求解预测

 

EPLOT                                          !显示单元

ALLSEL                                         !选择所有实体

 

!LSEL, S, , , 5                                    !选择线段(1,9中间间隔8)

!LSEL, A, , , 10                                    !添加到上面线段选择集

!NSLL, S, 1                                      !选择线段上所有节点

!D, ALL, UY                                     !施加位移约束

 

DL,5,, ALL

DL,10,, ALL

 

!LSEL, S, , , 8

!LSEL, A, , , 9

!NSLL, S, 1

!D, ALL, UX

 

DL,8,, ALL

DL,9,, ALL

DL,1,, UZ

 

 

EPLOT                                          !显示单元

ALLSEL                                         !选择所有实体

 

!LSEL, S, , , 1

!NSLL, S, 1

!D, ALL, UX, 1

 

SFL, 1, PRES, 22.4

 

!LSEL, S, , , 1

!NSLL, S, 1

!SF, ALL, PRES, 22.4

 

ALLSEL

SOLVE                                        !开始求解计算

FINISH

 


邀请回答 我来回答

当前暂无回答

回答可获赠 200金币

没解决?试试专家一对一服务

换一批