"""Function to Instantiate a JacketSE Assembly: \n

INPUTS \n

All hardwired, so edit the quantities below all the way to the line "#________________ DO NOT MODIFY THE FOLLOWING ________________#" \n

-See JacketOpt_PyOPT for more information. \n

OUTPUTS \n

myjckt -jacket assembly instance \n\n

Optimization parameters: \n\n

f0 -float, target frequency [Hz]

f0epsilon -float, f0*(1+f0epsilon) will not be exceeded \n

jcktDTRmin -Float, minimum jacket member DTR allowed, mostly for 60+waterdepths. \n

mxftprint -Float, max allowed foot print [m]

guesses -Float(n), guesses for all design variables check out DesVar class. \n

bounds -Float(n,2), bounds for all design variables check out DesVar class. \n\n

SAMPLE CALLS: \n

1.OPTIMIZATION: python JacketOpt_ExtCobyla.py C:\RRD\PYTHON\WISDEM\JacketSE\src\jacketse\MyJacketInputs.py \n

2.OPTIMIZATION: python JacketOpt_PyOPT.py C:\RRD\PYTHON\WISDEM\JacketSE\src\jacketse\MyJacketInputs.py True \n

3.BUILD JACKET: python >>> myjacket=C:\RRD\PYTHON\WISDEM\JacketSE\src\jacketse\MyJacketInputs.py \n

"""

#Set inputs ###----ALL USER INPUT----###

Jcktins=JcktGeoInputs()

Jcktins.nlegs =4

Jcktins.nbays =4

Jcktins.batter=30.5495 #=(-43.127+24.614)/(5.939-5.333)

Jcktins.dck_botz =15.651

Jcktins.dck_width= 8.+1.2

Jcktins.weld2D = 0.

Jcktins.VPFlag = True #vertical pile T/F; to enable piles in frame3DD set pileinputs.ndiv>0

Jcktins.clamped= True #whether or not the bottom of the structure is rigidly connected. Use False when equivalent spring constants are being used.

Jcktins.AFflag = False #whether or not to use apparent fixity piles

Jcktins.PreBuildTPLvl = 0 #if >0, the TP is prebuilt according to rules per PreBuildTP

#______________________________________________#

#Soil inputs

Soilinputs=SoilGeoInputs()

Soilinputs.zbots =-np.array([3.,5.,7.,15.,30.,50.])

Soilinputs.gammas =np.array([10000.,10000.,10000.,10000.,10000.,10000.])

Soilinputs.cus =np.array([60000.,60000.,60000.,60000.,60000.,60000.])

Soilinputs.phis =np.array([26.,26.,26.,26.,26.,26])#np.array([36.,33.,26.,37.,35.,37.5])#np.array([36.,33.,26.,37.,35.,37.5])

Soilinputs.delta =25.

Soilinputs.sndflg =True

Soilinputs.PenderSwtch =False #True

Soilinputs.SoilSF =1.

#______________________________________________#

#Water and wind inputs

Waterinputs=WaterInputs()

Waterinputs.wdepth =50.

Waterinputs.wlevel =50. #Distance from bottom of structure to surface THIS, I believe is no longer needed as piles may be negative in z, to check and remove in case

Waterinputs.T=12. #Wave Period

Waterinputs.HW=10. #Wave Height

Waterinputs.Cd=3. #Drag Coefficient, enhanced to account for marine growth and other members not calculated

Waterinputs.Cm=8.#2. #ADded mass Coefficient

Windinputs=WindInputs()

Windinputs.HH=88.15+2.4 #CHECK HOW THIS COMPLIES....

Windinputs.U50HH=30. #assumed gust speed

Windinputs.Cdj=4. #Drag Coefficient for jacket members, enhanced to account for TP drag not calculated otherwise

Windinputs.Cdt=2 #Drag Coefficient for tower, enhanced to account for TP drag not calculated otherwise

#______________________________________________#

#Pile data

Pilematin=MatInputs()

Pilematin.matname=np.array(['RC'])

Pilematin.E=np.array([ 2.1e11])

Pilematin.G=np.array([8.07690e+10])

Pilematin.rho=np.array([3339.12]) *34274.82/36876. #From Test04.txt in SD, to check with official FAST certtest

Pileinputs=PileGeoInputs()

Pileinputs.Pilematins=Pilematin

Pileinputs.ndiv=1 #3 ###----USER INPUT----###

Pileinputs.Dpile=2.082

Pileinputs.tpile=0.491

Pileinputs.Lp=0. #[m] Embedment length

#______________________________________________#

#Legs data

legmatin=MatInputs()

legmatin.matname=(['steel'])

legmatin.rho=np.array([7850.])

Dleg=np.asarray([1.2]).repeat(Jcktins.nbays+1) #e.g., np.array([2.0,1.8,1.6,1.6,1.6])

tleg=np.asarray([0.05,0.05,0.035,0.035,0.035])

leginputs=LegGeoInputs()

leginputs.legZbot = 4.5 ###----USER INPUT----###

leginputs.ndiv=1 ###----USER INPUT----###

leginputs.legmatins=legmatin

leginputs.Dleg=Dleg

leginputs.tleg=tleg

#The following is a passthrough variables

legbot_stmphin =(45.5-43.127) #=2.373 Distance from bottom of leg to second joint along z; must be>0

#______________________________________________#

#Xbrc data

Xbrcmatin=MatInputs()

Xbrcmatin.matname=np.array(['steel']).repeat(Jcktins.nbays)

Xbrcmatin.rho=np.array([7850.])

Dbrc=np.asarray([0.8]).repeat(Jcktins.nbays)

tbrc=np.asarray([0.02]).repeat(Jcktins.nbays)

Xbrcinputs=XBrcGeoInputs()

Xbrcinputs.Dbrc=Dbrc

Xbrcinputs.tbrc=tbrc

Xbrcinputs.ndiv=1 ###----USER INPUT----###

Xbrcinputs.Xbrcmatins=Xbrcmatin

Xbrcinputs.precalc=False #This can be set to true if we want Xbraces to be precalculated in D and t, in which case the above set Dbrc and tbrc would be overwritten

#______________________________________________#

#Mbrc data

Mbrcmatin=MatInputs()

Mbrcmatin.matname=np.array(['steel'])

Mbrcmatin.rho=np.array([7850.])

Mbrcinputs=MudBrcGeoInputs()

Mbrcinputs.Dbrc_mud=0.8 ###----USER INPUT----###

Mbrcinputs.tbrc_mud=0.02

Mbrcinputs.ndiv=2 ###----USER INPUT----###

Mbrcinputs.Mbrcmatins=Mbrcmatin

Mbrcinputs.precalc=False #This can be set to true if we want Mudbrace to be precalculated in D and t, in which case the above set Dbrc_mud and tbrc_mud would be overwritten

#______________________________________________#

#Hbrc data

Hbrcmatin=MatInputs()

Hbrcmatin.matname=np.array(['steel'])

Hbrcmatin.rho=np.array([7850.])

Dbrc_hbrc=1.1 ###----USER INPUT----###

Hbrcinputs=HBrcGeoInputs()

Hbrcinputs.Dbrch=Dbrc_hbrc ###----USER INPUT----###

Hbrcinputs.ndiv=0

Hbrcinputs.Hbrcmatins=Hbrcmatin

Hbrcinputs.precalc=True #This can be set to true if we want Hbrace to be set=Xbrace top D and t, in which case the above set Dbrch and tbrch would be overwritten

#______________________________________________#

#TP data

#Note PrebuildTPLvl is set in JacketIns ###----USER INPUT----###

#TP lumped mass data

TPlumpinputs=TPlumpMass() ###----ALL USER INPUT----###

TPlumpinputs.mass = 666.e3#-98385.33+( 7850.*4*np.pi/4.*(1.2**2-(1.2-2.*0.04)**2)*4) #[kg] TO MODIFY AFTER WE ASSESS OVERALL TP MASS, to be reduced for the steel weight

TPlumpinputs.CMoff= np.array([0.,0.,2.]) #CG of concrete block is 2 m above intersection of diagonal braces

TPlumpinputs.I = 1./12*TPlumpinputs.mass*np.array([8.**2+4.**2,8.**2+4.**2,8.**2+8.**2,0.,0.,0.])

TPstrtmatin=MatInputs()

TPstmpsmatin=MatInputs()

TPgirdmatin=MatInputs()

TPbrcmatin=MatInputs()

TPstemmatin=MatInputs()

TPstmpsmatin.matname=np.array(['steel'])

TPstmpsmatin.rho=np.array([7850.])

TPstrtmatin.matname=np.array(['steel'])

TPstrtmatin.rho=np.array([1350.])

TPbrcmatin.matname=np.array(['steel'])

TPbrcmatin.rho=np.array([1350.])#np.array([7850.]) tpstrucmass=np.pi/4*(1.2**2-(1.2-2*0.04)**2)*7850*4*4 7850.*tpstrucmass/myjckt.TP.TPouts.mass=1462.185

TPgirdmatin.matname=np.array(['steel'])

TPgirdmatin.rho=np.array([1350.])

TPstemmatin.matname=np.array(['steel']).repeat(2) ###----ALL USER INPUT----###

TPstemmatin.rho=np.array([1350.])

TPinputs=TPGeoInputs()

TPinputs.TPstrtmatins=TPstrtmatin

TPinputs.TPbrcmatins=TPbrcmatin

TPinputs.TPstemmatins=TPstemmatin

TPinputs.TPstmpmatins=TPstmpsmatin

TPinputs.TPgirdmatins=TPgirdmatin

#Set TP dimensions as leg and brace dimensions

TPinputs.Dstrut=1.2 ###----ALL USER INPUT----###

TPinputs.tstrut=0.04

TPinputs.Dgir=TPinputs.Dstrut

TPinputs.tgir=TPinputs.tstrut

TPinputs.Dbrc=TPinputs.Dstrut

TPinputs.tbrc=TPinputs.tstrut

###----ALL USER INPUT----###

TPinputs.hstump=0.499#(16.15-15.651)

TPinputs.Dstump=1.2

TPinputs.tstump=0.04

TPinputs.stumpndiv=1

TPinputs.brcndiv=1

TPinputs.girndiv=1

TPinputs.strutndiv=1

TPinputs.stemndiv=1

TPinputs.nstems=3

TPinputs.Dstem=np.array([5.6]).repeat(TPinputs.nstems)

TPinputs.tstem=np.array([0.032,0.032,0.032])

TPinputs.hstem=np.array([(4.)/TPinputs.nstems]).repeat(TPinputs.nstems)

#______________________________________________#

#Tower data

Twrmatin=MatInputs()

Twrmatin.matname=np.array(['steel'])

Twrmatin.rho=np.array([7850.])

Twrinputs=TwrGeoInputs()

Twrinputs.Twrmatins=Twrmatin

#Twrinputs.Htwr=88.15 #Trumped by HH

Twrinputs.Htwr2frac=1./88.15 #fraction of tower height with constant x-section

Twrinputs.ndiv=np.array([1,1]) #ndiv for uniform and tapered section ###----USER INPUT----###

Twrinputs.DeltaZmax= 5. #[m], maximum FE element length allowed in the tower members (i.e. the uniform and the tapered members)

Twrinputs.Db=5.6

Twrinputs.DTRb=Twrinputs.Db/0.032

Twrinputs.Dt=4.

Twrinputs.DTRt=Twrinputs.Dt/0.03

#Set whether or not DTRb and DTRt for the tower are the same. Note if next set to False it will trump DTRt setting above

Twrinputs.DTRsdiff=True ##SET THIS TO TRUE IF YOU WANT DTRs to be different between base and top

#If you use the following 12 lines, The geometry defined above is ignored

ztwr=np.array([20.15,21.15,32.15,42.15,54.15,64.15,74.15,83.15,88.15]) -20.15

Dtwr=np.array([5.6,5.577,5.318,5.082,4.8,4.565,4.329,4.118,4])

ttwr=np.array([0.032,0.032,0.03,0.028,0.024,0.022,0.02,0.03,0.03])

pmtwr=np.array([ np.array([20.15,54.15,88.15])-20.15,[1.9e3,1.4e3,1.0e3]]).T #(3,2) first col z's, second weights

#Interpolate data to refine tower discretization

dz=1. #[m] maximum deltaz allowed in discretization

ztwr2=np.linspace(ztwr[0],ztwr[-1],round((ztwr[-1]-ztwr[0])/dz)) #New discretization

Dtwr_interp=interp1d(ztwr,Dtwr)

ttwr_interp=interp1d(ztwr,ttwr)

Twrinputs.ztwr=ztwr2

Twrinputs.Dtwr=Dtwr_interp(ztwr2)

Twrinputs.ttwr=ttwr_interp(ztwr2)

Twrinputs.TwrlumpedMass=np.zeros([pmtwr.shape[0],11])

Twrinputs.TwrlumpedMass[:,0:2]=pmtwr

TwrRigidTop=False #False=Account for RNA via math rather than a physical rigidmember

#______________________________________________#

#RNA data

RNAins=RNAprops()

RNAins.mass=350.e3 #[kg]

RNAins.I[0]=86.579E+6 #[kg m2]

RNAins.I[1]=53.530E+6 #[kg m2]

RNAins.I[2]=58.112E+6 #[kg m2]

RNAins.CMoff[2]=2.34 #[m]

RNAins.Thoff[2]=2.4 #[m]

RNAins.yawangle=45. #angle with respect to global X, CCW looking from above, wind from left

RNAins.rna_weightM=True

#______________________________________________#

#RNA loads Fx-z, Mxx-zz

RNA_F=np.array([1000.e3,0.,0.,0.,0.,0.]) #unfactored thrust, though accounting for gust and dynamic effects (no IEC PSF though)

#______________________________________________#

# Frame3DD parameters ###----ALL USER INPUT----###

FrameAuxIns=Frame3DDaux()

FrameAuxIns.sh_fg=1 #shear flag-->Timoshenko

FrameAuxIns.deltaz=5.

FrameAuxIns.geo_fg=0

FrameAuxIns.nModes = 6 # number of desired dynamic modes of vibration

FrameAuxIns.Mmethod = 1 # 1: subspace Jacobi 2: Stodola

FrameAuxIns.lump = 0 # 0: consistent mass ... 1: lumped mass matrix

FrameAuxIns.tol = 1e-9 # mode shape tolerance

FrameAuxIns.shift = 0.0 # shift value ... for unrestrained structures

FrameAuxIns.gvector=np.array([0.,0.,-9.8065]) #GRAVITY

#Decide whether or not to consider DLC 6.1 as well

twodlcs=False

#______________________________________________#

#______________________________________________#

# OTHER AUXILIARY CONSTRAINTS AND TARGETS FOR OPTIMIZATION #

#______________________________________________#

#______________________________________________#

#Set Optimization Bounds and guesses for the various variables:

# x= [ batter, Dpile, tpile, Lp, Dleg, tleg, Dbrc, tbrc, Dbrc_mud, tbrc_mud, Dgit, tgir, Db, DTRb Dt, DTRt Htwr2fac dck_widthfact]

MnCnst = np.array([ 8., 1., 1.*0.0254, 20., 1., 1.*0.0254, 1., 1.*0.0254, 1., 1.*0.0254, 1., 1.*0.0254, 5., 120., 3., 120., 0.05, 2.])

MxCnst = np.array([ 15., 2.5, 5.*0.0254, 50., 2.5, 5.*0.0254, 2., 5.*0.0254, 2., 5.*0.0254, 2., 5.*0.0254, 7., 200., 4., 200., 0.25, 3.])

guesses= np.array([ 10., 1.5, 1.5*0.0254, 26., 1.8, 1.5*0.0254, 1.2, 1.5*0.0254, 1.2, 1.5*0.0254, 1.2, 1.5*0.0254, 6., 140., 3.5, 150., 0.2, 2.])

#SET Maximum Footprint [m]

mxftprint =30.

#Set target frequency [Hz] and f0epsilon, i.e. fmax=(1+f0eps)*f0

f0=0.22

f0epsilon=0.1

#Set the minminimum DTR allowed for jacket members; mostly for 60+waterdepths

jcktDTRmin=22.

#_____________________________________________________________#

#________________ DO NOT MODIFY THE FOLLOWING ________________#

#_____________________________________________________________#

bounds=np.vstack((MnCnst,MxCnst))

desvarmeans=np.mean(bounds,1)

# Now Launch the assembly and pass all of the inputs

myjckt=set_as_top(JacketSE(Jcktins.clamped,Jcktins.AFflag,twodlcs=twodlcs))

myjckt.JcktGeoIn=Jcktins

myjckt.Soilinputs=Soilinputs

myjckt.Waterinputs=Waterinputs

myjckt.Windinputs=Windinputs

myjckt.Pileinputs=Pileinputs

myjckt.leginputs=leginputs

myjckt.legbot_stmphin =legbot_stmphin #Distance from bottom of leg to second joint along z; must be>0

myjckt.Xbrcinputs=Xbrcinputs

myjckt.Mbrcinputs=Mbrcinputs

myjckt.Hbrcinputs=Hbrcinputs

myjckt.TPlumpinputs=TPlumpinputs

myjckt.TPinputs=TPinputs

myjckt.Twrinputs=Twrinputs

myjckt.TwrRigidTop=TwrRigidTop #Account for RNA via math rather than a physical rigidmember

myjckt.RNAinputs=RNAins

myjckt.RNA_F=RNA_F

myjckt.FrameAuxIns=FrameAuxIns