! TITLE: stub_shaft_340a168.mac (AAAV) ! !--------------------------------------------------------------------------- ! Rev - 11/08/2000 J. Rinaldi !--------------------------------------------------------------------------- ! ! This macro creates the stub shaft used in the AAAV assembly, config. C. ! ! CSss0: Cartesian Stub-Shaft Coordinate System ! CSss1: Cylindrical Stub-Shaft Coordinate System ! CSss2: Spherical Stub-Shaft Coordinate System ! ! !******************************* TABLE OF CONTENTS ************************** ! ! Model Assembly Preparation ! ! 1. Parameter List ! 2. Model Coordinate System ! 3. Material and Element Type ! 4. Solid Model ! ! !************************** Model Assembly Preparation ********************** ! ! mod_assemb_prep ! prepare database for assembly ! ! /title,Stub Shaft / SK340A168 ! ! /units,bin ! British/inch units ! ! !<<<<<<<<<<<<<<<<<<<<<<<< 1. Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>> ! ! stub_shaft_340a168.parm ! name of parameter file ! ! !<<<<<<<<<<<<<<<<<<<< 2. Model Coordinate Systems >>>>>>>>>>>>>>>>>>>> ! ! mod_coord_system,0,0,0,0,0,0,0 ! ! !<<<<<<<<<<<<<<<<<<<< 3. Material and Element Type >>>>>>>>>>>>>>>>>>> ! ! !---------------------- Material Properties -------------------- ! UIMP,1,EX, , ,28.5E6, ! Material: 416 CRES, cond. T UIMP,1,DENS, , ,.097/386.4, UIMP,1,NUXY, , ,0.28, UIMP,XX,XXXX ! !------------------------- Element Type ------------------------ ! ET,1,shell93 ET,2,solid95 ! ! !<<<<<<<<<<<<<<<<<<<<<<<<<<< 4. Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>> ! csys,CSss0 ! k,,0,0,0 k,,0,0,1 ! init_entity_num ! axp1=km-1 axp2=km ! k,,0,r1,0 k,,0,r2,0 k,,0,r2,l1 k,,0,r3,l1 k,,0,r3,l1+l2 ! k,,0,r4,l1+l2 k,,0,r4,l1+l3 k,,0,r5,l1+l3 ! k,,0,r5,l1+l4-0.74 ! area boundary keypoints k,,0,r5,l1+l4-0.74/2 ! k,,0,r5,l1+l4 k,,0,r6,l1+l4 ! k,,0,r6,l1+l4+l5 k,,0,r7,l1+l4+l5 k,,0,r7,l6 k,,0,r8,l6 k,,0,r8,l6-l7 ! k,,0,r9,l6-l7-((r8-r9)*tan(phi1*rconv)) k,,0,r9,l8+((r10-r9)*tan(phi2*rconv)) k,,0,r10,l8 k,,0,r1,l8 ! create_lines ! lfillt,lm-15,lm-16,rfil2 lfillt,lm-1,lm,rfil1 ! init_entity_num ! lcomb,lm-17,lm-1 ! init_entity_num ! *go,:skip isolate_lines ! al,all ! init_entity_num ! isolate_areas ! revolve_aset,axp1,axp2,360,4 ! nummrg,kp ! init_entity_num ! lplot ! ! ! (1) (2) (3) (4) (5) (6) (7) ! ROTATE_VOLM_MSH, aCS1, stang, mstyp, mshdens, elmszovr, etpa, etpv, ! ! elmnl, elmxl, mtrl, rl, ecs, sprat, radrot, varcang ! (8) (9) (10) (11) (12) (13) (14) (15) ! ! ROTATE_VOLM_MSH,CSss1,0,,,,1,2,0.1,0.2,1,,,,r5,90 ! nummrg,node ! init_entity_num ! csys,CSss0 ! vsel,s,loc,x,0,r2 eslv,s vclear,all vdele,all,,,1 ! init_entity_num ! lplot !*go,:skip csys,CSss0 asel,s,loc,x,0 da,all,symm asel,s,loc,z,0 da,all,all,0 ! /com Apply 100 equivalent pounds to the part in the lateral direction ! Calculate the pressure required to generate F/2 pounds on the model half. ! Apress = (pi/2)*(rL) ! Aproj = rL !asel,s,loc,z,l6 asel,s,area,,am-38,am-37,1 asum,fine *get,Apress,area,0,area Aproj = Apress*(2/pi) press = (F/2)/Aproj /PSF,PRES,NORM,2,0 /PBF,DEFA, ,1 sfa,all,1,PRES,press alls finish !/config,fsplit,250 /solu !****************************** !*** Powerdynamics - PCG Solve !****************************** antype,0 eqslv,pcg,1E-6 solve fini /post26 ! csys,CSss0 nsel,s,loc,x,0 nsel,r,loc,y,r5 nsel,r,loc,z,l1+l4-0.74/2 get_max_entity nsol,2,nm,u,y,del_x :skip