mario.r
Guest
Hello.
as from title I am performing the design and optimization of a board for launch vehicles reinforced with only double t stringers. Such stringers are designed to have the same size for internal and external flange. the measurement system used is the Anglo-Saxon one. the design of the same happened by parameters (l= panel length, r= panel radius and tskin parameters, tweb,tfe,tfi for various thicknesses, bs for stringer spacing). the resolutive analysis focuses on the buckling and I have to design 3 different panels with lengths fixed to 300,400 and 500 inches. the panel is subject to a compression load of 1000 lb/in. the structural bond is the critical load factor (lambda) which must be greater than 1.2. the code for design and optimization has already been written but I meet various types of problems in the analysis of the first panel (which consequently I have also experienced in the panels from 400 and 500 inches):
1) small mesh leads to a number of high nodes for optimization (about 20000) and in optimization phase, after the first iteration, it stops because the program exceeds 32,000 nodes (limit for student license).
2) based on the values of the parameters, with fairly low values (l and r always remain the same and I cannot change them), ansys does not detect solutions neither static nor for the buckling and does not take the necessary volume for the minimization of the weight, also parameterized (weight=density*volume).
3)when I find an exact combination to obtain solution and volume, the program continues to optimize with solutions for lambda and weight of 99990000, values too high. eventually fails to achieve convergence.
My question is, where am I wrong?
below place the code I am using for the 300 in panel:
as from title I am performing the design and optimization of a board for launch vehicles reinforced with only double t stringers. Such stringers are designed to have the same size for internal and external flange. the measurement system used is the Anglo-Saxon one. the design of the same happened by parameters (l= panel length, r= panel radius and tskin parameters, tweb,tfe,tfi for various thicknesses, bs for stringer spacing). the resolutive analysis focuses on the buckling and I have to design 3 different panels with lengths fixed to 300,400 and 500 inches. the panel is subject to a compression load of 1000 lb/in. the structural bond is the critical load factor (lambda) which must be greater than 1.2. the code for design and optimization has already been written but I meet various types of problems in the analysis of the first panel (which consequently I have also experienced in the panels from 400 and 500 inches):
1) small mesh leads to a number of high nodes for optimization (about 20000) and in optimization phase, after the first iteration, it stops because the program exceeds 32,000 nodes (limit for student license).
2) based on the values of the parameters, with fairly low values (l and r always remain the same and I cannot change them), ansys does not detect solutions neither static nor for the buckling and does not take the necessary volume for the minimization of the weight, also parameterized (weight=density*volume).
3)when I find an exact combination to obtain solution and volume, the program continues to optimize with solutions for lambda and weight of 99990000, values too high. eventually fails to achieve convergence.
My question is, where am I wrong?
below place the code I am using for the 300 in panel:
for optimization instead:/prep7
!introduco parametri geometrici
r=200
l=300
tskin=0.1
tfi=0.1
tfe=0.3
tweb=0.2
hweb=4
fi=0.5
fe=fi
bs=180/10
nstr=nint(180/bs)
!definisco modello (bottom-up)
k,1,0,0,0
k,99,0,0,l
csys,1
k,2,r,90-bs,0
k,3,r,90-bs/2-fi/2,0
k,4,r,90-bs/2,0
k,5,r,90-bs/2+fi/2,0
k,6,r,90,0
k,7,r-hweb,90-bs/2-fi/2,0
k,8,r-hweb,90-bs/2,0
k,9,r-hweb,90-bs/2+fi/2,0
l,1,99
larc,2,3,1,r
larc,3,4,1,r
larc,4,5,1,r
larc,5,6,1,r
lstr,4,8
lstr,8.9
lstr,8,7
adrag,2,,,,,,1
adrag,3,,,,,1
adrag,4,,,,,1
adrag,5,,,,,1
adrag,6,,,,,1
adrag,7,,,,,1
adrag,8,,,,,1
ldele,1
!impostazione pannello completo definendo gruppi di aree
Allsel
asel,s,area,1
asel,a,area,4
agen,nstr,all,,0,-bs
cm,area_skin,area
Allsel
asel,s,area,2
asel,a,area,3
agen,nstr,all,,0,-bs
cm,area_fi,area
Allsel
asel,s,area,5
agen,nstr,all,,0,-bs
cm,area_web,area
Allsel
asel,s,area,6
asel,a,area,7
agen,nstr,all,,0,-bs
cm,area_fe,area
Allsel
nummrg,all
!definizione materiale e costanti reali(ex e gxy in mp)
et,1,shell93
mp,nuxy,1,0.32
mp,ex,1,113800
r,1,tfe
r,2,tweb
r,3,tskin
r,4,tskin+tfi
!genero mesh
esize,20
Allsel
cmsel,s,area_fe
real,1
amesh,all
Allsel
cmsel,s,area_web
real,2
amesh,all
Allsel
cmsel,s,area_skin
real,3
amesh,all
Allsel
cmsel,s,area_fi
real,4
amesh,all
Allsel
nummrg,node
fini
!definizione vincoli e compressione
/solu
csys,1
Allsel
lsel,s,loc,z,0
lsel,r,loc,x,r
nsll,s,1
d,all,ux,0
d,all,uy,0
Allsel
lsel,s,loc,x,r
lsel,r,loc,y,90
nsll,s,1
d,all,ux,0
d,all,uy,0
Allsel
lsel,s,loc,x,r
lsel,r,loc,y,-90
nsll,s,1
d,all,ux,0
d,all,uy,0
Allsel
lsel,s,loc,z,l
lsel,r,loc,x,r
nsll,s,1
d,all,ux,0
d,all,uy,0
Allsel
lsel,s,loc,x,r
lsel,r,loc,y,bs*5
nsll,s,1
d,all,ux,0
d,all,uy,0
Allsel
lsel,s,loc,z,0
lsel,r,loc,x,r
nsll,s,1
sfl,all,pres,1000
Allsel
lsel,s,loc,z,l
lsel,r,loc,x,r
nsll,s,1
sfl,all,pres,1000
Allsel
sftran
!soluzione statica(pstres serve for buckling analysis)
antype,static
pstres,on
solve
fini
/solu
antype,buckle
bucopt,subsp,5
mxpand,5
solve
fini
/post1
!prelevo il fattore di carico critico del primo modo dibuck
*get,lambda,mode,1,freq
!prelevo il volume
Allsel
etable,volu,volu
ssum
*get,volume,ssum,,item,volu
!definisco la funzione peso per l'ottimizzazione
densite=0.15*0.453592
peso=volume*densita
lgwrite,optpan5
Thank you in advance for any answers.*use,optpan5
/opt
opanl,optpan5
opvar,bs,dv,5,20,0.001
opvar,fi,dv,0.01,2,0.001
opvar,hweb,dv,0.1,5,0.001
opvar,tskin,dv,0.01,0.5,0.001
opvar,tfi,dv,0.05,0.5,0.001
opvar,tweb,dv,0.05,0.5,0.001
opvar,tfe,dv,0.05,0.5,0.001
opvar,lambda,sv,1.2,,0.001
opvar,peso,obj,,,0.00001
optype,subp
opsub,50,80
opexe
/out,results,txt
oplist,all
/out,term
fini