Реферат Моделювання сталих коливань пластин, що контактують з пружним підставою

)))/2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2))/2

k22=(sig1*(ka2**2)*(sig2*sig1*((exp(sig2*h1)+exp(-(sig2*h1)))/2)*((exp(sig1*h1)-exp(-(sig1*h1)))/2)-(alpha**2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2)*((exp(sig1*h1)+exp(-(sig1*h1)))/2)))/(- 4)

del21=(sig1*sig2*((1/4)*(ka2**4)-((gam1+(alpha**2))**2)+((1/4)*(ka2**4)+((gam1+(alpha**2))**2))*((exp(sig1*h1)+exp(-(sig1*h1)))/2)*((exp(sig2*h1)+exp(-(sig2*h1)))/2))-(alpha**2)*((gam1**2)+(sig1**2)*(sig2**2))*((exp(sig1*h1)-exp(-(sig1*h1)))/2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2))/2

del11=((exp (sig2 * h1) + exp (- (sig2 * h1)))/2)

fi=aimag (k22/del21)

function fi

znamr (alpha)! функція

(4), PARAMETER :: i=(0.0,1.0)! уявна одиниця (4) :: znamr, alpha

complex (4) :: eps, eps1, eps2, eps3, ka2, del11, del21, k11, k12, k22, psi, gam1, sig1, sig2! вони комплексні

! real :: mu1=1, mu2=2.6, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=3.8, ro1=1, omega! 1 набір

! real :: mu1=1, mu2=2.6, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=3.8, ro1=1, omega! 2 набір

! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=3.8, ro1=1, omega! 3 набір

! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=3.8, ro1=1, omega! 4 набір

! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=1, ro1=1, omega! 5 набір

real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=1, ro1=1, omega! 6 набір

common/global/omega

eps1=(mu2 * (h ** 3))/(6 * (1-nu2))

eps2=((h ** 3)/12) * ro2 * (omega ** 2)

eps3=h * ro2 * (omega ** 2)

ka2=sqrt ((ro1 * (omega ** 2))/mu1)

eps=(1-2 * nu1)/(2-2 * nu1)

psi=eps1 * (alpha ** 4) -eps2 * (alpha ** 2) -eps3

gam1=(alpha ** 2) - 0.5 * (ka2 ** 2)

sig2=sqrt ((alpha ** 2) - (ka2 ** 2))

sig1=sqrt ((alpha ** 2) -eps * (ka2 ** 2))

k11=(sig2*(ka2**2)*(sig1*sig2*((exp(sig1*h1)+exp(-(sig1*h1)))/2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2)-(alpha**2)*((exp(sig1*h1)-exp(-(sig1*h1)))/2)*((exp(sig2*h1)+exp(-(sig2*h1)))/2)))/(- 4)

k12=(sig1*sig2*(gam1+(alpha**2))*(((exp(sig1*h1)+exp(-(sig1*h1)))/2)*((exp(sig2*h1)+exp(-(sig2*h1)))/2)- 1)-(gam1*(alpha**2)+(sig1**2)*(sig2**2))*((exp(sig1*h1)-exp(-(sig1*h1)))/2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2))/2

k22=(sig1*(ka2**2)*(sig2*sig1*((exp(sig2*h1)+exp(-(sig2*h1)))/2)*((exp(sig1*h1)-exp(-(sig1*h1)))/2)-(alpha**2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2)*((exp(sig1*h1)+exp(-(sig1*h1)))/2)))/(- 4)

del21=(sig1*sig2*((1/4)*(ka2**4)-((gam1+(alpha**2))**2)+((1/4)*(ka2**4)+((gam1+(alpha**2))**2))*((exp(sig1*h1)+exp(-(sig1*h1)))/2)*((exp(sig2*h1)+exp(-(sig2*h1)))/2))-(alpha**2)*((gam1**2)+(sig1**2)*(sig2**2))*((exp(sig1*h1)-exp(-(sig1*h1)))/2)*((exp(sig2*h1)-exp(-(sig2*h1)))/2))/2

del11=((exp (sig2 * h1) + exp (- (sig2 * h1)))/2)

znamr=real ((k22/del21) * psi + 1)

function znamr

znami (alpha)! функція

(4), PARAMETER :: i=(0.0,1.0)! уявна одиниця (4) :: znami, alpha

complex (4) :: eps, eps1, eps2, eps3, ka2, del11, del21, k11, k12, k22, psi, gam1, sig1, sig2! вони комплексні

! real :: mu1=1, mu2=2.6, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=3.8, ro1=1, omega! 1 набір

! real :: mu1=1, mu2=2.6, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=3.8, ro1=1, omega! 2 набір

! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=3.8, ro1=1, omega! 3 набір

! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=3.8, ro1=1, omega! 4 набір

! real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.1, h1=1, ro2=1, ro1=1, omega! 5 набір

real :: mu1=1, mu2=5, nu1=0.3, nu2=0.3, h=0.01, h1=1, ro2=1, ro1=1, omega! 6 набір

common/global/omega

eps1=(mu2 * (h ** 3))/(6 * (1-nu2))

eps2=((h ** 3)/12) * ro2 * (omega ** 2)

eps3=h * ro2 * (omega ** 2)

ka2=sqrt ((ro1 * (omega ** 2))/mu1)

eps=(1-2 * nu1)/(2-2 * nu1)

psi=eps1 * (alpha ** 4) -eps2 * (alpha ** 2) -eps3