! TITLE: raceo_top_mpb.mac (AAAV) ! !----------------------------------------------------------------------- ! Rev - 04/04/2001 J. Rinaldi !----------------------------------------------------------------------- ! ! This macro computes the equivalent thickness ring for a temperature test ! in an aluminum assembly. ! ! First the actual race deflection is computed with an arbitrary thickness ! aluminum sleeve assembled outside of the race. The race deflection at the ! ball contact points (use the minimum algebraic value) is used to find an ! equivalent uniform ring which compresses the same amount when assembled ! with the aluminum sleeve. ! ! ! CS10: Cartesian central Coordinate System ! CS11: Cylindrical central Coordinate System ! CS12: Spherical central Coordinate System ! !************************** TABLE OF CONTENTS ************************* ! ! Model Assembly Preparation ! ! 1.0 Parameter List ! 2.0 Model Coordinate System ! 3.0 Solid Model ! 4.0 Apply Constraints & Coupling ! 5.0 Solve ! 6.0 Process Thermal Results ! ! !********************* Model Assembly Preparation ********************** ! mod_assemb_prep ! ! /title,MPB Upper Azimuth Outer Race (490A045-101) ! ! /units,bin ! ! !<<<<<<<<<<<<<<<<<<<<<<<< 1. Parameter List >>>>>>>>>>>>>>>>>>>>>>>>>> ! ! raceo_top_mpb.parm ! name of parameter file ! ! !<<<<<<<<<<<<<<<<<<<< 2. Model Coordinate Systems >>>>>>>>>>>>>>>>>>>> ! ! mod_coord_system,CS10,0,0,0,0,0,0 ! ! !<<<<<<<<<<<<<<<<<<<<<<<<<<< 3. Solid Model >>>>>>>>>>>>>>>>>>>>>>>>>> ! /view,1,,,1 /ang,1 ! ! Keypoint Definitions from raceo_top_mpb.parm ! k,,x1,y1,z1 ! km + 1 k,,x2,y2,z2 ! km + 2 k,,x3,y3,z3 ! km + 3 k,,x4,y4,z4 ! km + 4 k,,x5,y5,z5 ! km + 5 k,,x6,y6,z6 ! km + 6 k,,x7,y7,z7 ! km + 7 k,,x8,y8,z8 ! km + 8 k,,x9,y9,z9 ! km + 9 k,,x10,y10,z10 ! km + 10 k,,x11,y11,z11 ! km + 11 (center) k,,x12,y12,z12 ! km + 12 k,,x13,y13,z13 ! km + 13 k,,x14,y14,z14 ! km + 14 (center) ! Save the ball contact KP numbers ktop=km+12 kbot=km+10 ! ! Line definitions from or_mpb.igs ! L,km+1,km+2 ! lm + 1 L,km+2,km+3 ! lm + 2 L,km+4,km+5 ! lm + 3 L,km+6,km+5 ! lm + 4 L,km+3,km+4 ! lm + 5 L,km+7,km+8 ! lm + 6 L,km+1,km+9 ! lm + 7 larc,km+10,km+7,km+11,r1 ! lm + 8 larc,km+12,km+13,km+14,r2 ! lm + 9 L,km+13,km+6 ! lm + 10 larc,km+8,km+12,km+14,r3 ! lm + 11 larc,km+9,km+10,km+11,r4 ! lm + 12 al,all ! ! Add a uniform section of aluminum ! init_entity_num k,,x3,y3,z3 ! Define surface to couple to k,,x2,y2,z2 k,,x2,y2-extra ! Extra mtrl below chamfer k,,r_outer,y2-extra k,,r_outer,y3 k,,x4,y3 create_area ! ! Mesh the race ! ! (1) (2) (3) (4) (5) (6) (7) (8) ! MSH_ALGORITHM, perflsz, perfmsh, mstyp, mshdens, elmszovr, atargov, etpa, etpv, ! ! elmnl, elmxl, mtrl, rl, ecs, volsel ! (9) (10) (11) (12) (13) (14) asel,s,area,,am-1 msh_algorithm,1,1,0,0,,,1,,.02,.03,c440C,CS10 ! ! Mesh the aluminum ring ! asel,s,area,,am msh_algorithm,1,1,0,0,,,1,,.02,.03,alum,CS10 ! !------------------- Apply Boundary conditions ------------------------- ! ! Couple the interface nodes in x direction only ! Contrain OD/ID interface radially nsel,s,loc,x,x2-.001,x2+.001 nsel,r,loc,y,y2,y3 cpintf,ux ! Contrain bearing flange to top of ring nsel,s,loc,y,y3 nsel,r,loc,x,x3,x4 cpintf,uy ! Tie down the aluminum sleeve nsel,s,loc,y,y2-extra d,all,uy alls ! !----------------------- 6 Solve ----------------------- ! fini /solu antype,0 tunif,Ttest tref,Troom solv fini /post1 ksel,s,kp,,ktop nslk,s *get,ntop,node,0,num,max *get,uxctt,node,ntop,u,x ! Total movement of top node ! Compensate for thermal shrinkage of race by itself (no alum ring) uxct=uxctt-nx(ntop)*(Ttest-Troom)*ALPX440C ksel,s,kp,,kbot nslk,s *get,nbot,node,0,num,max *get,uxctb,node,nbot,u,x ! Total movement of bottom node ! Compensate for thermal shrinkage of race by itself (no alum ring) uxcb=uxctb-nx(nbot)*(Ttest-Troom)*ALPX440C ! Worst case deflection *if,uxctt,le,uxctb,then maxdef=uxctt *else maxdef=uxctb *endif /nopr /out,raceo_top_mpb,sum *msg Net Radial deflection of ball contact points with outer race. *msg *msg For use in raceo_top_mpb.tkw to determine equivalent ring thickness: *msg,,uxct,uxcb dri(top) = %g, dri(bottom) = %g *msg *if,uxct,le,uxcb,then *msg Plug the "top" deflection into Tk-Solver *else *msg Plug the "bottom" deflection into Tk-Solver *endif *msg ! Info for calibration below. *msg Total radial deflection of ball contact points: *msg,,uxctt,uxctb dri(top) = %g, dri(bottom) = %g /out /gopr *uilist,raceo_top_mpb,sum ! ! Build a planar model and iterate on it to find the equivalent ! race inner radius dimension ! ! Find: f(r_inner) = maxdef zero,'cal_mpb',1.5,x2-.03,maxdef,1e-7