Swept Sine

Figure 1. Swept Sine event

Figure 2. Swept Sine eventの上面図

初期車両速度、初期ステア値と最大ステア値、初期周波数と最大周波数、周波数変化率、開始時間と終了時間、および時間ステップは、イベントフォームで変更できます。ホイール回転速度は計算された値で、変更できません。

Figure 3.

ステアリングホイール角入力を定義するために使用する簡略式は次のとおりです：(1)

$$\theta ={\theta }_{\max }\ast \sin \left(\left({\omega }_{init}+{\omega }_{rate}\ast time\right)\ast time\right)$$

ステアリングホイール角を定義するために使用する実際の方程式には、角度からラジアンに変換する項、入力周波数を定義済みの最大周波数に制限する項、0のステアリング初期値から正弦関数に移行するための項などが追加されます。この方程式を実現するTemplex式は次のとおりです。

`{ds_swept_sine.max_str.value}*DTOR
*STEP(TIME, {ds_swept_sine.st_tim.value}, 0.0, {ds_swept_sine.st_tim.value}+0.001, 1.0)
*SIN(2*PI
*(min({ds_swept_sine.max_frq.value},{ds_swept_sine.ini_frq.value}+(({ds_swept_sine.frq_rt.value}/2)
*(TIME-{ds_swept_sine.st_tim.value}))))*(TIME-{ds_swept_sine.st_tim.value}))`

この方程式は以下のように評価されます。

45*DTOR*STEP(TIME, 2, 0.0, 2+0.001, 1.0)*SIN(2*PI*(min(3,0.25+((0.275/2)*(TIME-2))))*(TIME-2))

この式で実行される関数を行ごとに分解すると、次のようになります：

このイベントのテンプレートを以下に示します：

<ResOutput
     angle_type          = "YPR"
  />
  <ResOutput
     mrf_file            = "TRUE"
  />
  <ResOutput
     plt_file            = "TRUE"
  />
  <H3DOutput
     switch_on           = "TRUE"
     increment           = "1"
  />
  <ResOutput
     abf_file            = "TRUE"
  />
{if (tire_dataset.opt_omega.ival ==1)}
<!--Initial static analysis -->

<Simulate
	analysis_type = "Static"
	end_time      = "0.0"
/>	
{endif}

 <Deactivate
	element_type = "MOTION"
	element_id = "{mot_frnt_wheel.l.idstring}"
/>

<Deactivate
	element_type = "MOTION"
	element_id = "{mot_frnt_wheel.r.idstring}"
/>

<Deactivate
	element_type = "MOTION"
	element_id = "{mot_rear_wheel.l.idstring}"
/>

<Deactivate
	element_type = "MOTION"
	element_id = "{mot_rear_wheel.r.idstring}"
/>
{if (tire_dataset.opt_omega.ival ==2)}
<!--Initial static analysis -->

<Simulate
	analysis_type = "Static"
	end_time      = "0.0"
/>	
{endif}
<Deactivate
	element_type = "JPRIM"
	element_id = "{j_clamp_1_body.idstring}"
/>

<Deactivate
	element_type = "JPRIM"
	element_id = "{j_clamp_2_body.idstring}"
/>    

{if ds.inp_type.value == "Angle"}
  <Deactivate
	element_type = "FORCE"
	element_id = "{tor_str.idstring}"
  />
  <Motion_Joint
     id                  = "{wh_motion.idstring}"
     expr                = "IF(TIME-.5:{-ds.str_mag.value}D*STEP(TIME,0,0,0.001,1),{-ds.str_mag.value}D,{-ds.str_rate.value}D*(TIME-.5)-{ds.str_mag.value}D)"
  />                                
                                 
 <Simulate
      analysis_type       = "Transient"
      end_time            = "{(ds.max_str.value-ds.str_mag.value)/ds.str_rate.value+.5}"
      num_steps           = "{((ds.max_str.value-ds.str_mag.value)/ds.str_rate.value+.5)*100}"
  />
 
 {else}
 <Deactivate
 	element_type = "MOTION"
 	element_id = "{wh_motion.idstring}"
/>

<Force_Scalar_TwoBody
     id                  = "{tor_str.idstring}"
     type		 = "Torque"
     val_expression      = "IF(TIME-.5:{-ds.str_mag.value},{-ds.str_mag.value},{-ds.str_rate.value}*(TIME-.5)-{ds.str_mag.value})"
  />
  
  <Simulate
        analysis_type       = "Transient"
        end_time            = "{(ds.max_str.value-ds.str_mag.value)/ds.str_rate.value+.5}"
        num_steps      = "{((ds.max_str.value-ds.str_mag.value)/ds.str_rate.value+.5)*100}"
  />
  {endif}
  
 <Stop/>