Q) I would like to apply a load for pressure stiffening effects in a PSD analysis with a pressure PSD. However, that pressure loading will be considered in the PSD analysis, so how do I separate the prestress and PSD pressure loads?

A) A prestressed modal analysis can be performed to accomplish this. First, run the static analysis with the prestress pressure and the PSTRES,ON command. Then, remove the preload pressure and apply the PSD pressure in the modal analysis (with PSTRES,ON). This will create the load vector of only the PSD pressure loading. Lastly, the PSD analysis can be run as usual.

Please note that this method includes the stress stiffening effects from the first pressure loading, so the stiffness changes, but the PSD results only include the loading effects from the PSD pressure loading.

An example model is below, where the preload is a membrane pressure loading and the PSD load is an internal pressure loading.
-------------------------------
finish
/clear

MYRADIUS=3

/view,1,1,2,3
/triad,off
/auto

/prep7
et,1,181
keyopt,1,3,2
r,1,0.1
mp,ex ,1,10e6
mp,nuxy,1,.3
mp,dens,1,0.1/386.1

cylind,MYRADIUS,,0,5,0,90

vdele,all
csys,1
asel,s,loc,x,MYRADIUS
asel,inve
adele,all,,,1
csys,0
allsel,all

mshape,0
mshkey,1
esize,MYRADIUS*2*acos(-1)/80
amesh,all

lsel,s,loc,z,0
nsll,s,1
d,all,rotx
d,all,roty
d,all,uz

lsel,s,loc,x,0
nsll,s,1
d,all,ux
d,all,roty
d,all,rotz

lsel,s,loc,y,0
nsll,s,1
d,all,rotx
d,all,uy
d,all,rotz

lsel,s,loc,z,5
nsll,s,1
sf,all,pres,-1e3

allsel,all
eplot
finish

/solu
antype,static
pstres,on
ematwr,yes
solve
finish

/solu
antype,modal
modopt,lanb,10
mxpand,10
pstres,on

sfedele,all,all,pres

sfe,all,1,pres,,1e3

allsel,all
solve
finish

/solu
antype,spectr
spopt,psd

psdunit,1,pres
psdfrq,1,,1,5e3
psdval,1,1e6,1e6
lvscale,1

dmprat,0.05

pfact,1,node
psdres,disp,rel
psdcom
solve
finish

/post26
numvar,20
store,psd,5

nsol,2,node(MYRADIUS,0,5),u,x
nsol,4,node(MYRADIUS,0,5),u,z
rpsd,12,2,,,,UX_NOPRESTRESS
rpsd,14,4


Q) I would like to apply a load for pressure stiffening effects in a PSD analysis with a pressure PSD. However, that pressure loading will be considered in the PSD analysis, so how do I separate the prestress and PSD pressure loads?

A) A prestressed modal analysis can be performed to accomplish this. First, run the static analysis with the prestress pressure and the PSTRES,ON command. Then, remove the preload pressure and apply the PSD pressure in the modal analysis (with PSTRES,ON). This will create the load vector of only the PSD pressure loading. Lastly, the PSD analysis can be run as usual.

Please note that this method includes the stress stiffening effects from the first pressure loading, so the stiffness changes, but the PSD results only include the loading effects from the PSD pressure loading.

An example model is below, where the preload is a membrane pressure loading and the PSD load is an internal pressure loading.
-------------------------------
finish
/clear

MYRADIUS=3

/view,1,1,2,3
/triad,off
/auto

/prep7
et,1,181
keyopt,1,3,2
r,1,0.1
mp,ex ,1,10e6
mp,nuxy,1,.3
mp,dens,1,0.1/386.1

cylind,MYRADIUS,,0,5,0,90

vdele,all
csys,1
asel,s,loc,x,MYRADIUS
asel,inve
adele,all,,,1
csys,0
allsel,all

mshape,0
mshkey,1
esize,MYRADIUS*2*acos(-1)/80
amesh,all

lsel,s,loc,z,0
nsll,s,1
d,all,rotx
d,all,roty
d,all,uz

lsel,s,loc,x,0
nsll,s,1
d,all,ux
d,all,roty
d,all,rotz

lsel,s,loc,y,0
nsll,s,1
d,all,rotx
d,all,uy
d,all,rotz

lsel,s,loc,z,5
nsll,s,1
sf,all,pres,-1e3

allsel,all
eplot
finish

/solu
antype,static
pstres,on
ematwr,yes
solve
finish

/solu
antype,modal
modopt,lanb,10
mxpand,10
pstres,on

sfedele,all,all,pres

sfe,all,1,pres,,1e3

allsel,all
solve
finish

/solu
antype,spectr
spopt,psd

psdunit,1,pres
psdfrq,1,,1,5e3
psdval,1,1e6,1e6
lvscale,1

dmprat,0.05

pfact,1,node
psdres,disp,rel
psdcom
solve
finish

/post26
numvar,20
store,psd,5

nsol,2,node(MYRADIUS,0,5),u,x
nsol,4,node(MYRADIUS,0,5),u,z
rpsd,12,2,,,,UX_NOPRESTRESS
rpsd,14,4,,,,UZ_NOPRESTRESS

/gropt,logx,on
plvar,12,14
-------------------------------





Show Form
No comments yet. Be the first to add a comment!