vobera
Mechanical
- Feb 15, 2017
- 3
Hello,
I have one problem, I'm trying to connest shell and solid elements with coupling but It do not work. Could zou please help me? I'am attaching my input file.
finish
/clear
*abbr,SEMESTRALKA,'/input,semestr2,inp'
/COM, =======================
/COM, POCATECNI PODMINKY
/COM, =======================
! Rozmery
r_trubka_out=250 !Vnejsi polomer [mm]
t_trubka=10 !Tloušťka trubky [mm]
delka=2*r_trubka_out !Delka trubky [mm]
r_tryska_out=25 !Vnejsi polomer [mm]
t_tryska=5 !Tloušťka trysky [mm]
delka_2=8*r_tryska_out !Delka trysky[mm]
r_inner=3 !Polomer vrubu[mm], musi byt mensi nez t_tryska
r_tryska=5
PI=acos(-1)
! Sit set up
my_esys=0
divider=1/5
divider_thck=1
divider_rad=1.5
divider_vrub=4
lines_of_interest=20
lines_vrub=6
size=1
! Zatizeni
pressure=0.1 !Tlak zatizeni [MPa]
Young=2.1E5 !Younguv modul [MPa]
mu=0.3 !Poisson [-]
/prep7
/COM, =======================
/COM, NASTAVENI GRAFIKY
/COM, =======================
/title, TRYSKA
/rgb,index,100,100,100,0 !Invertuje barvy
/rgb,index,80,80,80,13
/rgb,index,60,60,60,14
/rgb,index,0,0,0,15
/COM, =======================
/COM, PRVEK
/COM, =======================
et,1,solid186 !Vyber objemoveho prvku pro sitovani
et,2,SHELL281 !Vyber shell prvku pro sitovani
keyopt,2,1,1 !Pouze membranova napjatost (ne momenty)
/COM, =======================
/COM, MATERIAL
/COM, =======================
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.1E5
MPDATA,PRXY,1,,0.3
/COM, =======================
/COM, DEFINICE SHELLU
/COM, =======================
! Shellove prvky pro trubku
sect,t_trubka,shell,,t_trubka
secdata,t_trubka,1,0,3
secoffset,BOT
seccontrol,0,0,0,0,1,1,1
! Shellove prvky pro trysku
sect,t_tryska,shell,,t_tryska
secdata,t_tryska,1,0,3
secoffset,BOT
seccontrol,0,0,0,0,1,1,1
/COM, =======================
/COM, TVORBA MODELU
/COM, =======================
! *** Vytvoreni zakladniho modelu -> ***
cylind,r_trubka_out,r_trubka_out-t_trubka,,delka,,180
local,11,cart,r_trubka_out-2*t_trubka,,,,,90,
wpcsys,1,11
cylind,r_tryska_out,r_tryska_out-t_tryska,0,delka_2,90,180
csys,1
asel,s,loc,x,r_trubka_out-t_trubka-0.1,r_trubka_out-t_trubka+0.1
*get,AREA1,area,,num,max
*get,VOLU1,volu,,num,max
vsba,VOLU1,AREA1
*get,VOLU2,volu,,num,max
vdele,VOLU2-1,,,1
local,12,1,,,,,,90
asel,s,loc,x,r_tryska_out-0.1,r_tryska_out+0.1
*get,AREA2,area,,num,max
*get,VOLU3,volu,,num,min
vsba,VOLU3,AREA2
*get,VOLU4,volu,,num,min
vdele,VOLU4,,,1
allsel,all
local,13,cylin,r_trubka_out-2*t_trubka,,,,,90,
csys,13
asel,s,loc,x,r_tryska_out-t_tryska/100,r_tryska_out+t_tryska/100
*get,AREA_FOR_DIV,area,,num,max
*get,AREA_SMALL,area,,num,min
asel,s,,,AREA_SMALL
lsla,s
*get,DIVIDED_LINE,line,,num,min
allsel,all
asbl,AREA_FOR_DIV,DIVIDED_LINE
csys,1
asel,s,loc,x,r_trubka_out-t_trubka/10,r_trubka_out+t_trubka/10
*get,AREA_FILL_1,area,,num,max
csys,13
asel,s,loc,x,r_tryska_out-t_tryska/100,r_tryska_out+t_tryska/100
*get,AREA_FILL_2,area,,num,min
allsel,all
afillt,AREA_FILL_1,AREA_FILL_2,r_tryska
al,1,37,43
al,2,33,42
va,1,4,7,10,22
csys,1
asel,s,loc,x,r_trubka_out-t_trubka,r_trubka_out-t_trubka
*get,AREA_UP,area,,num,max
csys,13
asel,s,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
*get,AREA_DOWN,area,,num,min
allsel,all
afillt,AREA_UP,AREA_DOWN,r_inner
vsba,4,11
*get,VOLU,volu,,num,min
vdele,VOLU+1,,,1
vadd,all
! ***<- Vytvoreni zakladniho modelu ***
! *** Rozdeleni modelu a tvorba ploch -> ***
csys,11
wpcsys,1,11
cylind,1.5*r_tryska_out,,0,delka_2,90,180
csys,13
asel,s,loc,x,1.5*r_tryska_out,1.5*r_tryska_out
*get,AREA_DIV,area,,num,max
vsba,2,AREA_DIV,sepo
vdele,1,,,1
adele,AREA_DIV,,,1
csys,0
wpcsys,1,0
cylind,0.5*r_tryska_out+r_trubka_out,,0,delka_2,0,90
csys,1
asel,s,loc,x,0.5*r_tryska_out+r_trubka_out,0.5*r_tryska_out+r_trubka_out
*get,AREA_DIV_2,area,,num,max
vsba,4,AREA_DIV_2,sepo
vdele,1,,,1
adele,AREA_DIV_2,,,1
allsel,all
local,14,cart,r_trubka_out+r_tryska_out,,,,,90
wpcsys,14
vsbw,all,sepo
local,15,cart,,,2*r_tryska_out
wpcsys,15
vsbw,all,sepo
local,16,cart,,,,(2*r_tryska_out/r_trubka_out)*(180/PI),90,,
wpcsys,16
vsbw,all,sepo
vglue,1,5,9
vdele,4 $ vdele,3 $ vdele,7 $ vdele,8
vsel,s,,,1 $ vsel,a,,,2 $ vsel,a,,,6 $ aslv,s
csys,1
asel,a,loc,x,r_trubka_out-t_trubka,r_trubka_out-t_trubka
csys,13
asel,a,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
asel,inve
adele,all,,,1
vsel,all $ aslv,s $ asel,a,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
asel,inve
aglue,all
allsel,all
! *** <- Rozdeleni modelu a tvorba ploch ***
/COM, =======================
/COM, ZESITOVANI
/COM, =======================
! Sitovani trubky
vsel,all $ aslv,s $ asel,a,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
asel,inve
aatt,1,,2,my_esys,t_trubka
lsel,s,length,,(delka-2*r_tryska_out)
lesize,all,,,divider*(delka-2*r_tryska_out),,1,,,1,
csys,16
lsel,s,length,,(2*r_tryska_out) $ lsel,u,loc,z,1,3*r_trubka_out
lesize,all,,,divider*2*r_tryska_out,,1,,,1,
lsel,s,length,,PI*(r_trubka_out-t_trubka)-(r_trubka_out-t_trubka)*(2*r_tryska_out/r_trubka_out)
lesize,all,,,divider*PI*(r_trubka_out-t_trubka)-(r_trubka_out-t_trubka)*(2*r_tryska_out/r_trubka_out),,1,,,1,
lsel,s,length,,(r_trubka_out-t_trubka)*(2*r_tryska_out/r_trubka_out)
lsel,a,length,,(r_trubka_out)*(2*r_tryska_out/r_trubka_out)
lesize,all,,,divider*(r_trubka_out-t_trubka)*(2*r_tryska_out/r_trubka_out),,1,,,1,
csys,13
vsel,all $ aslv,s $ asel,a,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
asel,inve
mshkey,1 $ amesh,all
! Sitovani trysky
csys,0
lsel,s,length,,PI/2*(r_tryska_out-t_tryska)
lesize,all,,,divider_rad*(r_tryska_out-t_tryska),,1,,,1,
lsel,s,loc,x,r_trubka_out+r_tryska_out+1e-3,r_trubka_out-2*t_trubka+delka_2-1e-3
lesize,all,,,divider*(r_trubka_out-2*t_trubka+delka_2-r_trubka_out+r_tryska_out),,1,,,1,
asel,s,loc,x,r_trubka_out+r_tryska_out+1,r_trubka_out-2*t_trubka+delka_2
aatt,1,,2,my_esys,t_tryska
mshkey,1 $ amesh,all
! Sitovani objemu
vsel,all $ mshkey,0 $ esize,size $ mshape,1 $ vmesh,all
! Coupling
cpintf,all,1e-3
! Constraint Eqn
*get,count_cp,cp,,num
*do,INDEX,1,count_cp-2,3
*get,node_1,cp,INDEX,term,1,node
*get,node_2,cp,INDEX,term,2,node
!ce,INDEX,0,node_1,UX,1,node_2,UX,-1,,UX,1
!ce,INDEX+1,0,node_1,UY,1,node_2,UY,-1,,UY,1
!ce,INDEX+2,0,node_1,UZ,1,node_2,UZ,-1,,UZ,1
cp,count_cp+INDEX,all,node_2,node_1
*enddo
!cpdele,all
cpdele,1,count_cp,1
/COM, =======================
/COM, OKRAJOVE PODMINKY
/COM, =======================
! Symetrie
csys,0
asel,s,loc,y,, $ asel,a,loc,z,, $ da,all,symm
lsel,s,loc,y,, $ lsel,a,loc,z,, $ dl,all,,symm
! Zamezení posuvu
lsel,s,loc,x,r_trubka_out-2*t_trubka+delka_2,r_trubka_out-2*t_trubka+delka_2
dl,all,,ux,
! Definice tlaku
csys,1 $ asel,s,loc,x,r_trubka_out-t_trubka,r_trubka_out-t_trubka
csys,13 $ asel,a,loc,x,r_trryska_out-t_tryska,r_tryska_out-t_tryska
/psymb,adir,1
sfa,all,1,pres,pressure
allsel,all
/COM, =======================
/COM, RESENI
/COM, =======================
/COM, =======================
/COM, DEFINICE CESTY A VYKRESLENI DO GRAFU
/COM, =======================
I have one problem, I'm trying to connest shell and solid elements with coupling but It do not work. Could zou please help me? I'am attaching my input file.
finish
/clear
*abbr,SEMESTRALKA,'/input,semestr2,inp'
/COM, =======================
/COM, POCATECNI PODMINKY
/COM, =======================
! Rozmery
r_trubka_out=250 !Vnejsi polomer [mm]
t_trubka=10 !Tloušťka trubky [mm]
delka=2*r_trubka_out !Delka trubky [mm]
r_tryska_out=25 !Vnejsi polomer [mm]
t_tryska=5 !Tloušťka trysky [mm]
delka_2=8*r_tryska_out !Delka trysky[mm]
r_inner=3 !Polomer vrubu[mm], musi byt mensi nez t_tryska
r_tryska=5
PI=acos(-1)
! Sit set up
my_esys=0
divider=1/5
divider_thck=1
divider_rad=1.5
divider_vrub=4
lines_of_interest=20
lines_vrub=6
size=1
! Zatizeni
pressure=0.1 !Tlak zatizeni [MPa]
Young=2.1E5 !Younguv modul [MPa]
mu=0.3 !Poisson [-]
/prep7
/COM, =======================
/COM, NASTAVENI GRAFIKY
/COM, =======================
/title, TRYSKA
/rgb,index,100,100,100,0 !Invertuje barvy
/rgb,index,80,80,80,13
/rgb,index,60,60,60,14
/rgb,index,0,0,0,15
/COM, =======================
/COM, PRVEK
/COM, =======================
et,1,solid186 !Vyber objemoveho prvku pro sitovani
et,2,SHELL281 !Vyber shell prvku pro sitovani
keyopt,2,1,1 !Pouze membranova napjatost (ne momenty)
/COM, =======================
/COM, MATERIAL
/COM, =======================
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.1E5
MPDATA,PRXY,1,,0.3
/COM, =======================
/COM, DEFINICE SHELLU
/COM, =======================
! Shellove prvky pro trubku
sect,t_trubka,shell,,t_trubka
secdata,t_trubka,1,0,3
secoffset,BOT
seccontrol,0,0,0,0,1,1,1
! Shellove prvky pro trysku
sect,t_tryska,shell,,t_tryska
secdata,t_tryska,1,0,3
secoffset,BOT
seccontrol,0,0,0,0,1,1,1
/COM, =======================
/COM, TVORBA MODELU
/COM, =======================
! *** Vytvoreni zakladniho modelu -> ***
cylind,r_trubka_out,r_trubka_out-t_trubka,,delka,,180
local,11,cart,r_trubka_out-2*t_trubka,,,,,90,
wpcsys,1,11
cylind,r_tryska_out,r_tryska_out-t_tryska,0,delka_2,90,180
csys,1
asel,s,loc,x,r_trubka_out-t_trubka-0.1,r_trubka_out-t_trubka+0.1
*get,AREA1,area,,num,max
*get,VOLU1,volu,,num,max
vsba,VOLU1,AREA1
*get,VOLU2,volu,,num,max
vdele,VOLU2-1,,,1
local,12,1,,,,,,90
asel,s,loc,x,r_tryska_out-0.1,r_tryska_out+0.1
*get,AREA2,area,,num,max
*get,VOLU3,volu,,num,min
vsba,VOLU3,AREA2
*get,VOLU4,volu,,num,min
vdele,VOLU4,,,1
allsel,all
local,13,cylin,r_trubka_out-2*t_trubka,,,,,90,
csys,13
asel,s,loc,x,r_tryska_out-t_tryska/100,r_tryska_out+t_tryska/100
*get,AREA_FOR_DIV,area,,num,max
*get,AREA_SMALL,area,,num,min
asel,s,,,AREA_SMALL
lsla,s
*get,DIVIDED_LINE,line,,num,min
allsel,all
asbl,AREA_FOR_DIV,DIVIDED_LINE
csys,1
asel,s,loc,x,r_trubka_out-t_trubka/10,r_trubka_out+t_trubka/10
*get,AREA_FILL_1,area,,num,max
csys,13
asel,s,loc,x,r_tryska_out-t_tryska/100,r_tryska_out+t_tryska/100
*get,AREA_FILL_2,area,,num,min
allsel,all
afillt,AREA_FILL_1,AREA_FILL_2,r_tryska
al,1,37,43
al,2,33,42
va,1,4,7,10,22
csys,1
asel,s,loc,x,r_trubka_out-t_trubka,r_trubka_out-t_trubka
*get,AREA_UP,area,,num,max
csys,13
asel,s,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
*get,AREA_DOWN,area,,num,min
allsel,all
afillt,AREA_UP,AREA_DOWN,r_inner
vsba,4,11
*get,VOLU,volu,,num,min
vdele,VOLU+1,,,1
vadd,all
! ***<- Vytvoreni zakladniho modelu ***
! *** Rozdeleni modelu a tvorba ploch -> ***
csys,11
wpcsys,1,11
cylind,1.5*r_tryska_out,,0,delka_2,90,180
csys,13
asel,s,loc,x,1.5*r_tryska_out,1.5*r_tryska_out
*get,AREA_DIV,area,,num,max
vsba,2,AREA_DIV,sepo
vdele,1,,,1
adele,AREA_DIV,,,1
csys,0
wpcsys,1,0
cylind,0.5*r_tryska_out+r_trubka_out,,0,delka_2,0,90
csys,1
asel,s,loc,x,0.5*r_tryska_out+r_trubka_out,0.5*r_tryska_out+r_trubka_out
*get,AREA_DIV_2,area,,num,max
vsba,4,AREA_DIV_2,sepo
vdele,1,,,1
adele,AREA_DIV_2,,,1
allsel,all
local,14,cart,r_trubka_out+r_tryska_out,,,,,90
wpcsys,14
vsbw,all,sepo
local,15,cart,,,2*r_tryska_out
wpcsys,15
vsbw,all,sepo
local,16,cart,,,,(2*r_tryska_out/r_trubka_out)*(180/PI),90,,
wpcsys,16
vsbw,all,sepo
vglue,1,5,9
vdele,4 $ vdele,3 $ vdele,7 $ vdele,8
vsel,s,,,1 $ vsel,a,,,2 $ vsel,a,,,6 $ aslv,s
csys,1
asel,a,loc,x,r_trubka_out-t_trubka,r_trubka_out-t_trubka
csys,13
asel,a,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
asel,inve
adele,all,,,1
vsel,all $ aslv,s $ asel,a,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
asel,inve
aglue,all
allsel,all
! *** <- Rozdeleni modelu a tvorba ploch ***
/COM, =======================
/COM, ZESITOVANI
/COM, =======================
! Sitovani trubky
vsel,all $ aslv,s $ asel,a,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
asel,inve
aatt,1,,2,my_esys,t_trubka
lsel,s,length,,(delka-2*r_tryska_out)
lesize,all,,,divider*(delka-2*r_tryska_out),,1,,,1,
csys,16
lsel,s,length,,(2*r_tryska_out) $ lsel,u,loc,z,1,3*r_trubka_out
lesize,all,,,divider*2*r_tryska_out,,1,,,1,
lsel,s,length,,PI*(r_trubka_out-t_trubka)-(r_trubka_out-t_trubka)*(2*r_tryska_out/r_trubka_out)
lesize,all,,,divider*PI*(r_trubka_out-t_trubka)-(r_trubka_out-t_trubka)*(2*r_tryska_out/r_trubka_out),,1,,,1,
lsel,s,length,,(r_trubka_out-t_trubka)*(2*r_tryska_out/r_trubka_out)
lsel,a,length,,(r_trubka_out)*(2*r_tryska_out/r_trubka_out)
lesize,all,,,divider*(r_trubka_out-t_trubka)*(2*r_tryska_out/r_trubka_out),,1,,,1,
csys,13
vsel,all $ aslv,s $ asel,a,loc,x,r_tryska_out-t_tryska,r_tryska_out-t_tryska
asel,inve
mshkey,1 $ amesh,all
! Sitovani trysky
csys,0
lsel,s,length,,PI/2*(r_tryska_out-t_tryska)
lesize,all,,,divider_rad*(r_tryska_out-t_tryska),,1,,,1,
lsel,s,loc,x,r_trubka_out+r_tryska_out+1e-3,r_trubka_out-2*t_trubka+delka_2-1e-3
lesize,all,,,divider*(r_trubka_out-2*t_trubka+delka_2-r_trubka_out+r_tryska_out),,1,,,1,
asel,s,loc,x,r_trubka_out+r_tryska_out+1,r_trubka_out-2*t_trubka+delka_2
aatt,1,,2,my_esys,t_tryska
mshkey,1 $ amesh,all
! Sitovani objemu
vsel,all $ mshkey,0 $ esize,size $ mshape,1 $ vmesh,all
! Coupling
cpintf,all,1e-3
! Constraint Eqn
*get,count_cp,cp,,num
*do,INDEX,1,count_cp-2,3
*get,node_1,cp,INDEX,term,1,node
*get,node_2,cp,INDEX,term,2,node
!ce,INDEX,0,node_1,UX,1,node_2,UX,-1,,UX,1
!ce,INDEX+1,0,node_1,UY,1,node_2,UY,-1,,UY,1
!ce,INDEX+2,0,node_1,UZ,1,node_2,UZ,-1,,UZ,1
cp,count_cp+INDEX,all,node_2,node_1
*enddo
!cpdele,all
cpdele,1,count_cp,1
/COM, =======================
/COM, OKRAJOVE PODMINKY
/COM, =======================
! Symetrie
csys,0
asel,s,loc,y,, $ asel,a,loc,z,, $ da,all,symm
lsel,s,loc,y,, $ lsel,a,loc,z,, $ dl,all,,symm
! Zamezení posuvu
lsel,s,loc,x,r_trubka_out-2*t_trubka+delka_2,r_trubka_out-2*t_trubka+delka_2
dl,all,,ux,
! Definice tlaku
csys,1 $ asel,s,loc,x,r_trubka_out-t_trubka,r_trubka_out-t_trubka
csys,13 $ asel,a,loc,x,r_trryska_out-t_tryska,r_tryska_out-t_tryska
/psymb,adir,1
sfa,all,1,pres,pressure
allsel,all
/COM, =======================
/COM, RESENI
/COM, =======================
/COM, =======================
/COM, DEFINICE CESTY A VYKRESLENI DO GRAFU
/COM, =======================