Do you have a verification test of a pressure psd?
In the two inputs below, the 1 sigma result from a pressure psd of 10 psi^2/hz from 50 to 51 hz is compared to a harmonic analysis of a pressure of sqrt(10) psi applied at 50.5 hz. They give the same tip displacement of 0.123. !!! Pressure PSD /PREP7 /title, k8x, ssa, EXAMPLE FOR A PRESSURE PSD LOAD USING PSDUNIT CMD. et,1,187 ex,1,10e6 nuxy,1,.3 dens,1,.00025 block,,20,,1,,2 esize,1 vmesh,1 nsel,s,loc,x, d,all,all nsel,s,loc,y,1 sf,all,pres,1 ! apply a uniform pressure to selected nodes nsel,all finish /solution antype,modal ! run the modal analysis part modopt,lanb,20 mxpand,20,,,yes solve finish /solution antype,spectr spopt,psd,,on dmprat,.02 psdres,disp,abs ! input the press psd table and request output psdres,velo,abs psdres,acel,abs psdunit,1,pres psdfrq,1,,50,51 psdval,1,10,10 lvscale,1 dmprat,.02 ! 2% damping pfact,1,node ! do participation factor calculations psdcom ! new mode combination comm. at rev54 which replaces mcomb,psd to perform psd mode combinations solve finish /post1 set,3,1 uy_max=uy(node(20,1,2)) !!! Equivalent Harmonic with pressure loading /PREP7 /title, k8x, ssa, EXAMPLE FOR A PRESSURE PSD LOAD USING PSDUNIT CMD. et,1,187 ex,1,10e6 nuxy,1,.3 dens,1,.00025 block,,20,,1,,2 esize,1 vmesh,1 nsel,s,loc,x, d,all,all nsel,all finish /solution antyp,harm dmprat,.02 nsel,s,loc,y,1 sf,all,pres,sqrt(10) ! apply a uniform pressure to selected nodes nsel,all harf,,50.5 nsub,1 solve finish /post1 set,1,1 uy_max_real=uy(node(20,1,2)) set,1,1,,1 uy_max_imag=uy(node(20,1,2)) uy_max_srss=sqrt(uy_max_real**2+uy_max_imag**2) |
||
|