忽略法兰部在拉深变形过程中的厚度变化和模具圆角部分的影响。 取法兰部内缘线为初始己知速度边界。设沿模口的质点流动速度等于凸模下行速度，其 方向与法兰内缘线的法向相同。 应力边界条件的选取，模拟计算时一般取法兰外缘为初始应力边界，即沿法兰外缘线·： =0。这里，0。=(24p。)/(L,h)。其中，L,为法兰外缘线长度，p,为压边力，4为板料与 模具接触的摩擦系数。当压边力为待求参数时，取法兰内缘为初始应力边界。 法兰部特征线场的计算，根据边值问题的类型分别按Couchy问题、Riemann)问题和混合 问题求解。 作者编制了不规则直壁件拉深成形的平面应力特征线场模拟计算软件，运用该软件可以 计算拉深变形过程中法兰部的应力分布、速度分布和流动迹线，并能计算最佳毛料形状和合 理的加载条件。 值得指出的是，平面应力特征线场法使用简单，进行一次全程模拟（在IBM-AT机上计 算)只需30~40mn,与有限元法相比计算时间大大缩短。因此，平面应力特征线场法为拉 深过程的模拟提供了一种新的途径，同时使建立在过程数值模拟基础上的工艺参数的优化计 算成为可能。 2.2计算结果及讨论 (1)应力场分析 用本方法计算出的在盒形件拉深变形初期法兰部的应力分布如图1所示。图中的十字线 1.0 35 0.8 0.6 0.4 0.2-1=40mm- 0g.2=250N/mm2 0 3 5 78 Separated point number 图1法兰部应力分布（图彩毛料） 图2沿模口的法向应力分布（圆形毛料） Fig,1 Stress distribution in the flange Fig.2 Normal stress distribution along the region (circular blank shape) dic opening (circular blank shape) 154忽略法兰部在拉 深变形过程中的 厚度变化和模具圆角部分的影响 。 取法兰部内缘线为初始己知速度边界 。 设沿模 口 的 质点流动 速度等于 凸模下行 速度 , 其 方 向与 法兰内缘线的法 向相同 。 应力边界条件的 选取 , 模拟 计算时 一般取 法兰外缘为初始应力边界 , 即沿法兰 外缘线『 、 = 。 : 。 这 里 , a 。 = ( 2产p 。 ) / ( L : h ) 。 其中 , L 二 为 法兰 外缘线长度 , p 。 为压边力 , 声为 板料与 模具接触 的摩擦 系数 。 当压边力为待求参数时 , 取 法兰内缘为初始应力边界 。 法 兰部特征线场 的计算 , 根 据边值 问题的 类型分别按 C o u c h y 向题 、 iR e m a n 问题和混合 问 题求解 。 作 者编制 了不规则直壁件拉深成形的 平面应力特征线 场模拟计算软件 , 运用该软件可以 计算拉深变形过程中法 兰部的应力分布 、 速度分 布和流动迹线 , 并能 计算最佳毛 料形状和合 理的加载 条件 。 算 ) 渲得指 出的是 , 平 面应力特征线场法 使用简单 , 进行 一次全程模拟 (在 BI M 一 A T 机上计 只需 3 0 一 4 0m i n , 与 有限元法相 比 计算时间大大缩短 。 因此 , 平面应力特征线场法 为拉 深过程 的模拟 提供 了一种新的途径 算成为可能 。 , 同时 使建立在过程数值模拟基础上的工艺参 数的优化计 2 。 2 计算结果 及讨论 ( 1) 应力场分析 用 本方 法计算出的 在盒 形件拉深变形初期法兰部的应力分布如图 1所示 。 图中 的十字线 飞 。 0 N 0 . 6 0 . 2 0 . 〔户绪 0 . 吞 。戈` s e P a r a t e d P o i n t n u m b e r 图 1 法兰部应力分布 ( 圆 形毛料 ) 图 2 沿模 口 的法向应力分布 ( 圆 形毛料 ) F 19 . 1 1 5 4 5 t r e s s d i s t r i b u t i o n i n t h e f l a 刀 g e F 19 . 2 N o r m 么 1 s t r e s s d i s t r i b 让 t i o n a l o n g t h e e g i o n ( c i o c u l a r b l a n k s h a p e ) d i e o P e n i n g ( c i r c u l a ` b l a n k S h a P e )