真诚请教高手帮忙,谢谢!
命令流:
finish
/clear
/FILNAME,Xxz interphase spacers
/TITLE,Xxz analysis of interphase spacer
/prep7
ET,1,beam188
!定义导线单元类型
l=600
h0=75620
area=2.7596e-4
yingli=0.95*h0/area
q0=9.80665*0.9222
bizai=q0/area
istran=0.95*h0/area/7.3e10
!导线水平应变
ET,2,SOLID185
!定义记忆合金单元类型
ET,3,MASS21
!定义重球单元类型
ET,4,beam188
!间隔棒
!导线材料特性及本构关系
mp,ex,1,7.3e10
!材料属性
mp,prxy,1,0.3
mp,dens,1,3341.8
mp,ex,2,3e6
mp,prxy,2,0.3
mp,ex,4,2.8e11
mp,prxy,4,0.28
mp,dens,4,7850
sectype,1,beam,csolid
!截面
secdata,0.0108
!半径
sectype,4,beam,rect
secdata,0.005,0.03
r,1
r,3
r,3,0,0,1.5
!实常数
*do,i,-40,80,1
a=(i-1)*5+205
k,i+41+101000,10.75,a,-2*yingli/bizai*sinh(bizai*a/(2*yingli))*sinh(bizai*(l-a)/(2*yingli))+sqrt(3)/2*27.075
!节点
*enddo
*do,j,101001,101120,1
l,j,j+1
!连线
latt,1,1,1,,,,1
!赋予属性
lesize,all,,,1
lmesh,all
lsel,u,,,all
*enddo
*do,i,-40,80,1
a=(i-1)*5+205
k,i+41+102000,10.5625,a,-2*yingli/bizai*sinh(bizai*a/(2*yingli))*sinh(bizai*(l-a)/(2*yingli))+sqrt(3)/2*27.45
*enddo
*do,j,102001,102120,1
l,j,j+1
latt,1,1,1,,,,1
lesize,all,,,1
lmesh,all
lsel,u,,,all
*enddo
*do,i,-40,80,1
a=(i-1)*5+205
k,i+41+103000,10.1875,a,-2*yingli/bizai*sinh(bizai*a/(2*yingli))*sinh(bizai*(l-a)/(2*yingli))+sqrt(3)/2*27.45
*enddo
*do,j,103001,103120,1
l,j,j+1
latt,1,1,1,,,,1
lesize,all,,,1
lmesh,all
lsel,u,,,all
*enddo
*do,i,-40,80,1
a=(i-1)*5+205
k,i+41+104000,10,a,-2*yingli/bizai*sinh(bizai*a/(2*yingli))*sinh(bizai*(l-a)/(2*yingli))+sqrt(3)/2*27.075
*enddo
*do,j,104001,104120,1
l,j,j+1
latt,1,1,1,,,,1
lesize,all,,,1
lmesh,all
lsel,u,,,all
*enddo
*do,i,-40,80,1
a=(i-1)*5+205
k,i+41+105000,10.1875,a,-2*yingli/bizai*sinh(bizai*a/(2*yingli))*sinh(bizai*(l-a)/(2*yingli))+sqrt(3)/2*26.7
*enddo
*do,j,105001,105120,1
l,j,j+1
latt,1,1,1,,,,1
lesize,all,,,1
lmesh,all
lsel,u,,,all
*enddo
*do,i,-40,80,1
a=(i-1)*5+205
k,i+41+106000,10.5625,a,-2*yingli/bizai*sinh(bizai*a/(2*yingli))*sinh(bizai*(l-a)/(2*yingli))+sqrt(3)/2*26.7
*enddo
*do,j,106001,106120,1
l,j,j+1
latt,1,1,1,,,,1
lesize,all,,,1
lmesh,all
lsel,u,,,all
*enddo
nummrg,node
nummrg,kp
!第一相导线
l,101061,104061
!间隔棒
l,102061,105061
l,103061,106061
LDIV,ALL,,,2,
LDIV,ALL,,,2,
NUMMRG,NODE
NUMMRG,KP
l,4,5
l,5,6
l,6,7
l,7,8
l,8,9
l,9,4
lsel,s,line,,721,723,1
!赋予属性
latt,4,4,4
lesize,all,,,1
lmesh,721,723,1
lsel,s,line,,730,738,1
latt,4,4,4
lesize,all,,,1
lmesh,730,738,1
lsel,s,line,,724,729,1
latt,2,2,2
lesize,all,,,1
VMESH,724,729,1
nummrg,node
nummrg,kp
type,3
real,3
e,733
d,1,all
!约束
d,122,all
d,243,all
d,364,all
d,485,all
d,606,all
d,121,all
d,242,all
d,363,all
d,484,all
d,605,all
d,726,all
nummrg,all