CN108170984A - 一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法 - Google Patents

一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法 Download PDF

Info

Publication number
CN108170984A
CN108170984A CN201810051101.XA CN201810051101A CN108170984A CN 108170984 A CN108170984 A CN 108170984A CN 201810051101 A CN201810051101 A CN 201810051101A CN 108170984 A CN108170984 A CN 108170984A
Authority
CN
China
Prior art keywords
spring
main spring
main
thickness
order gradient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
CN201810051101.XA
Other languages
English (en)
Inventor
周长城
汪晓
于曰伟
梁勋乾
杨腾飞
杨铖兆
李晓晗
刘灿昌
朱召辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Technology
Original Assignee
Shandong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University of Technology filed Critical Shandong University of Technology
Priority to CN201810051101.XA priority Critical patent/CN108170984A/zh
Publication of CN108170984A publication Critical patent/CN108170984A/zh
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16FSPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
    • F16F1/00Springs
    • F16F1/02Springs made of steel or other material having low internal friction; Wound, torsion, leaf, cup, ring or the like springs, the material of the spring not being relevant
    • F16F1/18Leaf springs
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • Vehicle Body Suspensions (AREA)

Abstract

本发明涉及一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法,属于悬架渐变刚度板簧技术领域。本发明可根据主簧片数和副簧片数,各片板簧的结构参数,开始接触载荷,额定载荷及首片主簧在额定载荷下的许用应力,对一级渐变刚度板簧的各片非等厚主簧的预夹紧应力进行匹配设计。通过样机试验可知,本发明所提供的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法是正确的,为一级渐变刚度板簧的各片非等厚主簧自由切线弧高的设计奠定了可靠的技术基础。利用该方法可确保各片非等厚主簧的预夹紧应力满足设计要求,提高产品的设计水平、可靠性和使用寿命及车辆行驶安全性;同时,降低产品的设计及试验费用,加快产品开发速度。

Description

一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法
技术领域
本发明涉及车辆悬架渐变刚度板簧,特别是一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法。
背景技术
为了满足在不同载荷下的车辆行驶平顺性的设计要求,可采用一级渐变刚度板簧,即在末片主簧和首片主副簧之间设计有一定的主副簧渐变间隙。由于首片主簧受力复杂的要求,通常采用非等厚主簧,即首片主簧的厚度大于其他各片主簧的厚度。为了提高主簧强度和使用寿命,通过主簧各自不同的自由切线弧高,确保装配预夹紧后的首片主簧的初始切线弧高满足设计要求;同时,使首片主簧或前几片主簧受预夹紧压应力,而末片或后几片主簧受预夹紧拉应力,提高板簧可靠性和使用寿命。各片非等厚主簧预夹紧应力的匹配设计,不仅影响各片主簧的预夹紧应力及车辆行驶平顺性和安全性,而且还是各片非等厚主簧自由切线弧高设计的前提。然而,据所查资料可知,由于受各片非等厚主簧等效厚度及各片非等厚主簧之间应力厚度方之差计算的制约,先前一直未曾给出准确可靠的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法,不能满足车辆快速发展及对悬架一级渐变刚度板簧现代化CAD设计的要求。随着车辆行驶速度及其对平顺性要求的不断提高,对一级渐变刚度板簧设计提出了更高要求,因此,必须建立一种精确、可靠的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法,为各片非等厚主簧自由切线弧高的设计奠定可靠的技术基础,满足车辆行业快速发展、车辆行驶平顺性及一级渐变刚度板簧的设计要求,提高产品的设计水平、可靠性和使用寿命及车辆行驶安全性;同时,降低设计及试验费用,加快产品开发速度。
发明内容
针对上述现有技术中存在的缺陷,本发明所要解决的技术问题是提供一种简便、可靠的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法,其匹配设计流程图,如图1所示。一级渐变刚度板簧非等厚主簧和副簧的一半对称夹紧结构示意图,如图2所示,是由主簧1和副簧2构成。一级渐变刚度板簧板簧的宽度为b,弹性模量为E,骑马螺栓夹紧的根部平直段的一半长度L0,各片主簧的厚度不相等,首片主簧的厚度大于其他各片主簧的后,以满足首片主簧受力复杂的要求,即非等厚叠加主簧。主簧1的片数为n,各片主簧的厚度hi,各片主簧的一半作用长度为LiT,各片主簧的自由切线弧高为Hgi0,i=1,2,,…,n。副簧2的片数为m,各片副簧的厚度hAj,各片副簧的一半作用长度为LAjT,j=1,2,…,m。板簧装配预夹紧后的首片主簧的初始切线弧高为HgC1,首片副簧的初始切线弧高为HgAC1,并且在末片主簧和首片副簧之间形成渐变间隙δMA,从而满足板簧渐变刚度的设计要求。通过各片非等厚主簧的各自不同自由切线弧高及自由曲率半径,确保装配预夹紧后的首片主簧初始切线弧高和各片主簧预夹紧应力的设计要求。根据主簧片数和副簧片数,各片主簧和副簧的结构参数,开始接触载荷,额定载荷及首片主簧在额定载荷下的许用应力,对一级渐变刚度板簧各片非等厚主簧的预夹紧应力进行匹配设计,为各片非等厚主簧自由切线弧高设计奠定可靠的技术基础。
为解决上述技术问题,本发明所提供的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法,其特征在于采用以下匹配设计步骤:
(1)一级渐变刚度板簧的非等厚主簧、主副簧的根部重叠部分等效厚度hMe、hMAe计算:
I步骤:主簧根部重叠部分的等效厚度hMe
根据主簧片数n,各片主簧的厚度hi,i=1,2,...,n,对一级渐变刚度板簧的各片非等厚主簧的根部重叠部分的等效厚度hMe进行计算,即
II步骤:主副簧根部重叠部分的等效厚度hMAe
根据副簧片数m,各片副簧的厚度hAj,j=1,2,…m,I步骤中计算得到的hMe,对一级渐变刚度板簧的主副簧根部重叠部分的效厚度hMAe进行计算,即
(2)一级渐变刚度板簧的首片主簧根部最大应力σmax1的计算:
根据板簧的宽度b,根部平直段的一半长度L0,首片主簧的厚度h1,首片主簧的一半作用长度L1T,额度载荷PN,步骤(1)中计算得到的hMe和hMAe,对一级渐变刚度板簧在额定载荷下的首片主簧的根部最大应力σmax1进行计算,即
(3)一级渐变刚度板簧的各片非等厚主簧之间应力厚度方之差的确定:
根据主簧片数n,首片主簧的厚度h1,首片主簧在额度载荷下的许用应力[σ1],步骤(2)中计算得到的σmax1,对一级渐变刚度板簧的各片主簧之间的应力厚度方之差进行确定,即
(4)一级渐变刚度板簧的各片非等厚主簧预夹紧应力σi的匹配设计:
根据主簧片数n,各片主簧的厚度hi,首片主簧在额度载荷下的许用应力[σ1],步骤(2)中计算得到的σmax1,步骤(3)中所确定的对一级渐变刚度板簧各片非等厚主簧的预夹紧应力σi进行匹配设计,i=1,2,…,n,即
本发明比现有技术具有的优点
由于受各片非等厚主簧等效厚度及各片非等厚主簧之间应力厚度方之差计算的制约,先前一直未曾给出准确可靠的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法,不能满足车辆快速发展及对悬架一级渐变刚度板簧现代化CAD设计的要求。本发明可根据主簧片数和副簧片数,各片主簧和副簧的结构参数,开始接触载荷,额定载荷及首片主簧在额定载荷下的许用应力,对一级渐变刚度板簧的各片非等厚主簧的预夹紧应力进行匹配设计。通过样机试验测试可知,本发明所提供的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法是正确的,可得到准确可靠的各片非等厚主簧预夹紧应力的匹配设计值,为一级渐变刚度板簧的各片非等厚主簧自由切线弧高的设计奠定了可靠的技术基础。利用该方法可确保各片非等厚主簧的预夹紧应力满足设计要求,提高产品的设计水平、可靠性和使用寿命及车辆行驶安全性;同时,降低产品的设计及试验费用,加快产品开发速度。
附图说明
为了更好地理解本发明,下面结合附图做进一步的说明。
图1是一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计流程图;
图2是一级渐变刚度板簧非等厚主簧和副簧的一半对称夹紧结构示意图。
具体实施方案
下面通过实施例对本发明作进一步详细说明。
实施例一,某一级渐变刚度板簧的宽度b=60mm,骑马螺栓夹紧的根部平直段的一半长度L0=50mm。主簧片数n=3片,各片非等厚主簧的厚度分别为h1=9mm,h2=8mm,h3=8mm,各片主簧的一半作用长度分别为L1T=525mm,L2T=450mm,L3T=350mm。副簧片数m=2片,各片副簧的厚度hA1=hA2=13mm,各片副簧的一半作用长度分别为LA1T=250mm,LA2T=150mm。开始接触载荷Pk=1900N,额定载荷PN=7500N,首片主簧在额定载荷下的许用应力[σ1]=400MPa。根据主簧片数和副簧片数,各片主簧和副簧的结构参数,开始接触载荷,额定载荷及首片主簧在额定载荷下的许用应力,对该一级渐变刚度板簧的各片非等厚主簧的预夹紧应力进行匹配设计。
本发明实例所提供的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法,其匹配设计流程如图1所示,具体匹配设计步骤如下:
(1)一级渐变刚度板簧的非等厚主簧、主副簧的根部重叠部分等效厚度hMe、hMAe的计算:
I步骤:主簧根部重叠部分的等效厚度hMe
根据主簧片数n=3,各片主簧的厚度h1=9mm,h2=8mm,h3=8mm,对该一级渐变刚度板簧的主簧根部重叠部分的等效厚度hMe进行计算,即
II步骤:主副簧根部重叠部分的等效厚度hMAe
根据副簧片数m=2,各片副簧的厚度hA1=hA2=13mm,I步骤中计算得到的hMe=12.1mm,对该一级渐变刚度板簧的主副簧根部重叠部分的效厚度hMAe进行计算,即
(2)一级渐变刚度板簧的首片主簧根部最大应力σmax1的计算:
根据板簧的宽度b=60mm,根部平直段的一半长度L0=50mm,首片主簧的厚度h1=9mm,首片主簧的一半作用长度L1T=525mm,开始接触载荷Pk=1900N,额度载荷PN=7500N,步骤(1)中计算得到的hMe=12.1mm和hMAe=18.3mm,对该一级渐变刚度板簧在额定载荷下的首片主簧的根部最大应力σmax1进行计算,即
(3)一级渐变刚度板簧的各片非等厚主簧之间应力厚度方之差的确定:
根据主簧片数n=3,首片主簧的厚度h1=9mm,首片主簧在额度载荷下的许用应力[σ1]=400MPa,步骤(2)中计算得到的σmax1=426.4MPa,对该一级渐变刚度板簧的各片非主簧之间的应力厚度方之差进行确定,即
(4)一级渐变刚度板簧的各片非等厚主簧预夹紧应力σi的匹配设计:
根据主簧片数n=3,各片主簧的厚度h1=9mm,h2=8mm,h3=8mm,首片主簧在额度载荷下的许用应力[σ1]=400MPa,步骤(2)中计算得到的σmax1=426.4MPa,步骤(3)中所确定的对该一级渐变刚度板簧的各片非等厚主簧的预夹紧应力σi进行匹配设计,即
通过样机试验测试可知,本发明所提供的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法是正确的,可得到准确可靠的各片非等厚主簧预夹紧应力的匹配设计值。
实施例二,某一级渐变刚度板簧的宽度b=63mm,骑马螺栓夹紧的根部平直段的一半长度L0=55mm。主簧片数n=2片,各片主簧的厚度h1=10mm,h2=9mm,各片主簧的一半作用长度分别为L1T=525mm,L2T=450mm。副簧片数m=2片,各片副簧的厚度hA1=hA2=13mm,各片副簧的一半作用长度分别为LA1T=250mm,LA2T=150mm。开始接触载荷Pk=1800N,额定载荷PN=7000N,首片主簧在额定载荷下的许用应力[σ1]=400MPa。根据主簧片数和副簧片数,各片主簧和副簧的结构参数,开始接触载荷,额定载荷及在额定载荷下首片主簧的许用应力,对该一级渐变刚度板簧的各片非等厚主簧的预夹紧应力进行匹配设计。
采用与实施例一相同的设计方法和步骤,对该一级渐变刚度板簧的各片非等厚主簧的预夹紧应力进行匹配设计,具体匹配设计步骤如下:
(1)一级渐变刚度板簧的非等厚主簧、主副簧的根部重叠部分等效厚度hMe和hMAe的计算:
I步骤:主簧根部重叠部分的等效厚度hMe
根据主簧片数n=2,各片主簧的厚度h1=10mm,h2=9mm,对该一级渐变刚度板簧的各片非等厚主簧根部重叠部分的等效厚度hMe进行计算,即
II步骤:主副簧根部重叠部分的等效厚度hMAe
根据副簧片数m=2,各片副簧的厚度hA1=hA2=13mm,I步骤中计算得到的hMe=12.0mm,对该一级渐变刚度板簧的主副簧根部重叠部分的效厚度hMAe进行计算,即
(2)一级渐变刚度板簧的首片主簧根部最大应力σmax1的计算:
根据板簧的宽度b=63mm,根部平直段的一半长度L0=55mm,首片主簧的厚度h1=9mm,首片主簧的一半作用长度L1T=525mm,开始接触载荷Pk=1900N,额度载荷PN=7500N,步骤(1)中计算得到的hMe=12.0mm和hMAe=18.3mm,对该一级渐变刚度板簧在额定载荷下的首片主簧的根部最大应力σmax1进行计算,即
(3)一级渐变刚度板簧的各片非等厚主簧之间应力厚度方之差的确定:
根据主簧片数n=2,首片主簧的厚度h1=10mm,首片主簧在额度载荷下的许用应力[σ1]=400MPa,步骤(2)中计算得到的σmax1=423.07MPa,对该一级渐变刚度板簧的各片非主簧之间的应力厚度方之差进行确定,即
(4)一级渐变刚度板簧的各片非等厚主簧预夹紧应力σi的匹配设计:
根据主簧片数n=2,各片主簧的厚度h1=10mm,h2=9mm,首片主簧在额度载荷下的许用应力[σ1]=400MPa,步骤(2)中计算得到的σmax1=423.07MPa,步骤(3)中所确定的对该一级渐变刚度板簧的各片非等厚主簧的预夹紧应力σi进行匹配设计,i=1,2,…,n,即
通过样机试验测试可知,本发明所提供的一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法是正确的,可得到准确可靠的各片非等厚主簧预夹紧应力的匹配设计值,为一级渐变刚度板簧的各片非等厚主簧自由切线弧高的设计奠定了可靠的技术基础。利用该方法可确保各片非等厚主簧的预夹紧应力满足设计要求,提高产品的设计水平、可靠性和使用寿命及车辆行驶安全性;同时,降低产品的设计及试验费用,加快产品开发速度。

Claims (1)

1.一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法,其中,一级渐变刚度板簧是由主簧和副簧构成,各片板簧的根部平直段用于骑马螺栓夹紧,装配夹紧后在主簧和副簧之间形成一级渐变间隙,满足板簧一级渐变刚度的要求;各片主簧的厚度不相等,即非等厚叠加主簧;预夹紧应力的匹配设计是各片非等厚主簧自由切线弧高设计的前提;根据主簧片数和副簧片数,各片主簧和副簧的结构参数,开始接触载荷,额定载荷及首片主簧在额定载荷下的许用应力,对一级渐变刚度板簧的各片非等厚主簧的预夹紧应力进行匹配设计,其特征在于采用以下具体匹配设计步骤:
(1)一级渐变刚度板簧的非等厚主簧、主副簧的根部重叠部分等效厚度hMe、hMAe计算:
I步骤:主簧根部重叠部分的等效厚度hMe
根据主簧片数n,各片主簧的厚度hi,i=1,2,...,n,对一级渐变刚度板簧的各片非等厚主簧的根部重叠部分的等效厚度hMe进行计算,即
II步骤:主副簧根部重叠部分的等效厚度hMAe
根据副簧片数m,各片副簧的厚度hAj,j=1,2,…m,I步骤中计算得到的hMe,对一级渐变刚度板簧的主副簧根部重叠部分的效厚度hMAe进行计算,即
(2)一级渐变刚度板簧的首片主簧根部最大应力σmax1的计算:
根据板簧的宽度b,根部平直段的一半长度L0,首片主簧的厚度h1,首片主簧的一半作用长度L1T,额度载荷PN,步骤(1)中计算得到的hMe和hMAe,对一级渐变刚度板簧在额定载荷下的首片主簧的根部最大应力σmax1进行计算,即
(3)一级渐变刚度板簧的各片非等厚主簧之间应力厚度方之差的确定:
根据主簧片数n,首片主簧的厚度h1,首片主簧在额度载荷下的许用应力[σ1],步骤(2)中计算得到的σmax1,对一级渐变刚度板簧的各片主簧之间的应力厚度方之差进行确定,即
(4)一级渐变刚度板簧的各片非等厚主簧预夹紧应力σi的匹配设计:
根据主簧片数n,各片主簧的厚度hi,首片主簧在额度载荷下的许用应力[σ1],步骤(2)中计算得到的σmax1,步骤(3)中所确定的对一级渐变刚度板簧各片非等厚主簧的预夹紧应力σi进行匹配设计,i=1,2,…,n,即
CN201810051101.XA 2018-01-19 2018-01-19 一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法 Withdrawn CN108170984A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810051101.XA CN108170984A (zh) 2018-01-19 2018-01-19 一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810051101.XA CN108170984A (zh) 2018-01-19 2018-01-19 一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法

Publications (1)

Publication Number Publication Date
CN108170984A true CN108170984A (zh) 2018-06-15

Family

ID=62515272

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810051101.XA Withdrawn CN108170984A (zh) 2018-01-19 2018-01-19 一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法

Country Status (1)

Country Link
CN (1) CN108170984A (zh)

Similar Documents

Publication Publication Date Title
CN108253059A (zh) 一级渐变刚度板簧等厚主簧预夹紧应力的匹配设计法
CN108170984A (zh) 一级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法
CN108278297A (zh) 一级渐变刚度板簧非等厚主簧自由切线弧高的设计方法
CN106802996A (zh) 两级副簧式非等偏频型渐变刚度板簧接触载荷的验算方法
CN106812849B (zh) 非等偏频型三级渐变刚度板簧的接触载荷的验算方法
CN108278300A (zh) 两级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法
CN106594139B (zh) 一级渐变刚度板簧的副簧各片厚度及片数的设计方法
CN108278306A (zh) 一级渐变刚度板簧等厚主簧自由切线弧高的设计方法
CN106812848B (zh) 非等偏频型三级渐变刚度板簧刚度特性的计算方法
CN108591328A (zh) 两级渐变刚度板簧等厚主簧预夹紧应力的匹配设计法
CN106777789B (zh) 非等偏频型渐变刚度板簧接触载荷的仿真计算方法
CN108256224A (zh) 三级渐变刚度板簧等厚主簧预夹紧应力的匹配设计法
CN108194549A (zh) 两级渐变刚度板簧非等厚主簧自由切线弧高的设计方法
CN108256222A (zh) 三级渐变刚度板簧非等厚主簧预夹紧应力的匹配设计法
CN106548003B (zh) 非等偏频型三级渐变刚度板簧根部最大应力的仿真计算法
CN106777802B (zh) 非等偏频型三级渐变刚度板簧的各级根部应力的计算方法
CN108278298A (zh) 两级渐变刚度板簧等厚主簧自由切线弧高的设计方法
CN108182332A (zh) 三级渐变刚度板簧等厚主簧初始切线弧高的仿真验算法
CN108006127A (zh) 一级渐变刚度板簧非等厚主簧初始切线弧高的仿真验算法
CN106611091B (zh) 非等偏频一级渐变刚度板簧的初始切线弧高的设计方法
CN108268721A (zh) 两级渐变刚度板簧非等厚主簧初始切线弧高的仿真验算法
CN106650176B (zh) 非等偏频型三级渐变刚度板簧挠度特性的计算方法
CN107967404A (zh) 端部非等构少片变截面板簧预夹紧应力的仿真计算法
CN108268719A (zh) 一级渐变刚度板簧非等厚主簧预夹紧应力的仿真计算法
CN108591327A (zh) 三级渐变刚度板簧等厚主簧自由切线弧高的设计方法

Legal Events

Date Code Title Description
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WW01 Invention patent application withdrawn after publication

Application publication date: 20180615

WW01 Invention patent application withdrawn after publication