!btm_blades ,arg1 ,arg2 ,arg3 ,arg4 ,arg5 !btm_blades ,vert load (kg) ,break off above cl(m) ,mass pitch (degs) ,angle of wire (degs) ,bendrad (m) finish /CLEAR,START *abbr,doit,btm_blades /input,start100,ans,'D:\Program Files\Ansys Inc\v100\ANSYS\apdl\',,,,,,,,,,,,,,,,1 /PREP7 arg1=39.5 arg2=-0.00412 arg3=0 arg4=6 !* ! values of parameters !Geometry tlength=0.370 blength=0.365 rootwidth=0.049 tipwidth=0.01 bthick=0.0042 bendrad=0.3959 !0.3915 wire_length=0.5 wire_angle=arg4 !-27 clamp_height=arg2 !.00298 clamp_length= tlength-blength !load weight_carried=arg1 !in kg load_angle=arg4 !-27 !0=vertical, +10=10 degs toward tip, -10=degs toward root !Material Properties maryoung=1.86e11 marpoiss=0.3 mardens=7900 !Maths raddeg=180/3.1415926 halftip=tipwidth/2 halfroot=rootwidth/2 taperl=blength*(1-(halftip/halfroot)) tipload=weight_carried*9.81 hor_tipload=tan(load_angle/raddeg)*tipload lat_wireend_load=tan(arg3/raddeg)*tipload !* thtip=blength/bendrad*raddeg thwaist=taperl/bendrad*raddeg thwire=tlength/bendrad*raddeg !* ET,2,BEAM4 !8 node quadralateral for blade ET,1,SHELL93 !8 node triangle for clamp !Real 1 & 2, shell elements for blade R,1,bthick, , , , , , R,2,bthick*10 !Real 3, beam elements for wire R,3,wire_area,wire_sma,wire_sma,wire_diameter,wire_diameter !material definition !Maraging mp, ex, 1, maryoung mp, dens, 1, mardens mp, prxy, 1, marpoiss /ESHAPE,1 /VIEW, 1, -0.361338990165 , -0.590998239660 , 0.721217869236 /ANG, 1, -68.0920680969 /DIST,1,1.08222638492,1 csys,1 k,31,bendrad,0,-halfroot ,2,bendrad,thwaist,-halftip ,30,bendrad,thtip,-halftip ,4,bendrad,thtip,halftip ,5,bendrad,thwaist,halftip ,6,bendrad,0,halfroot L,31,2 ,2,30 ,30,4 ,4,5 ,5,6 ,6,31 AL,1,2,3,4,5,6 aplot ESIZE,halfroot/8,0 real,1 MAT,1 amesh,1 CLOCAL, 11, 0, bendrad, thtip, 0, 0 csys,11 jang=(-wire_angle/raddeg) jop=wire_length*sin(jang) jadj=wire_length*cos(jang) k, 3, -clamp_height, clamp_length, 0 k, 8, -clamp_height, clamp_length, halftip k, 9, -clamp_height, clamp_length, -halftip k, 17, -jadj-clamp_height, (jop)+clamp_length,0 k, 22, -(-jadj-clamp_height), -((jop)+clamp_length), 0 k, 21, (jadj/10)-clamp_height, (-jop/10)+clamp_length,0 !Blade clamp L,30,9 ,9,3 ,3,8 ,8,4 l,3,17 l,3,21 l,21,22 !blade clamp area AL, 7,8,9,10,3 wire_diameter=0.00110 wire_area=(wire_diameter/2)**2*3.14 wire_sma=3.14*wire_diameter**4/64 R,3,wire_area,wire_sma,wire_sma,wire_diameter,wire_diameter real,3 mat, 1 lmesh,11 ESIZE,halfroot/4 lmesh,12 ESIZE,halfroot/8,0 lmesh,13 real, 2 MAT,1 amesh, 2 csys,0 !*go,:end DL,6,,all,0 FK,21,FX,tipload FK,21,FY,hor_tipload FK,21,Fz,lat_wireend_load FINISH /SOL ANTYPE,0 NLGEOM,1 NSUBST,10,0,0 /STATUS,SOLU SOLVE FINISH /POST1 PLDISP,0 PLDISP,1 nsel,all upcoord,1,ON ksel,s,kp,,30,31 nslk /output,btm_blade,out,,append *Status,argx nlist /output ksel,all nsel,all upcoord,-1,ON !:end