拉伸结构(国外英文资料).docVIP

拉伸结构(国外英文资料) #包括 stdio. h #包括数学。” #定义PI 3.14159 #定义E1 1e-3 #定义E0 1e-9 #定义M 4 #定义N0 500 / /控制牛顿迭代次数 双F（x [M]，双T2，int n）/角位移代数式组 { 双p [ m ]； p = [ 0 ] = 0.4 *余弦（T2 *π/ 180）+ 1.650 *余弦（x [ 0 ]）+ 1×余弦（x [ 1 ] ] - x [ 3 ]； P [ 1 ] = 0.4 *罪（T2 *π/ 180）+ 1.650 * sin（x [ 0 ]）- 1 *罪（x [ 1 ]）-0.360 p = [ 2 ] = 0.4 *余弦（T2 *π/ 180）+ 1.650 *余弦（x [ 0 ]）+ 0.6×余弦（x [ 1 ]）- 0.800×余弦（x [ 2 ]）- 1 p = [ 3 ] = 0.4 *（T2 *π/ 180）+ 1.650 *（x 0）- 0.6 *（x 1）- 0.800 *（x [ 2 ]） 返回p； } 双（双X [M]，int i，j）/建立雅可比矩阵 { 双p [ m ]； p [ 0 ] [ 0 ]＝1.650 *（x [ 0 ]）；p [ 0 ] 1＝1（* [ 1 ]）；p [ 0 ] 2＝0；p [ 0 ] 3＝1； p 1 [ 0 ]＝1.650×余弦（x [ 0 ] ]；p [ 1 ] 1＝1（1）]；p [ 1 ] 2＝0；p [ 1 ] 3＝0； p [ 2 ] [ 0 ]＝1.650 *（x [ 0 ]）；p [ 2 ] 1＝0.600（* [ 1 ]）；p [ 2 ] 2＝0.800（* [ 2 ]）；p [ 2 ] 3＝0； p 3 [ 0 ]＝1.650×余弦（x [ 0 ] ]；p [ 3 ] 1＝0.600（1）；p [ 3 ] 2＝0.800（2）]；p [ 3 ] 3＝0； 返回p [ i ]； } 双G（双[M] [M]，int n，双B [ M ]）/用高斯消去法求解线性方程组（系数矩阵为方阵） { int i，j，k，I1，J1； 双温； 国际row_maxmod； 双element_maxmod； 对于（j = 0；j；m；j + +） { element_maxmod =一个[ J ] [ J ]； row_maxmod = J； 为（i = j + 1；i 米；i++）/找出第J列最大元素 如果（晶圆厂（element_maxmod）晶圆厂（一个[我] [ J ]）） { element_maxmod =一个[我] [ J ]； row_maxmod =我； } 如果（晶圆厂（element_maxmod）＜E0） { printf（“系数行列式为零，线性方程组无解！“n）； 返回null； } 如果（row_maxmod！= J） { 为（k = J；K＜m；K +）/将主行与第J行交换 { 温度= [ k ]； 一个[ J ] [ k] =一个[ row_maxmod ] [ k ]； 一个[ row_maxmod ] [ k ] =温度； } 温度= b [ j ]； B [ J ] = B [ row_maxmod ]； B【row_maxmod ] =温度； } 为（k = J；K＜m；K +）/高斯消元 一个[ J ] [ k ] / = element_maxmod； B【J】/ = element_maxmod； 对于（I1 = j + 1；I1＜4；I1 + +） { 温度=一个[ 1 ] [ J ]； 对于（J1 = J；J1＜4；J1 + +） 一个[ ] [ ] I1 J1 =一个[ J ] [ 1 ] *温度； B [ 1 ] = B *温度[J].； } } 对（我= m-2；我 = 0；//逆序求线性方程组的解我--） 对于（j = 3；j i；j -） [我]； 返回[ n ]； } 双C（x [M]，双V [M]，双T2，int n）/建立角加速度方程常数项 { 双p [ m ]； p = [ 0 ] = 0.400 *余弦（T2 *π/ 180）+ 1.650×0 [ * ] 0（0）] 1（1）]（1）]； p = [ 1 ] = 0.400 *（T2 *π/ 180）+ 1.650 * v [ 0 ] * [ 0 ] *（x [ 0 ]）- 1×V [ 1 ] * 1 [ * ] *（V [ 1 ]）； p = [ 2 ] = 0.400 *余弦（T2 *π/ 180）+ 1.650×0 [ * ] [ 0 ] * 0.600 [ * ] [ 1 ] * 1（1）] - 0.800（2）]（2）]（2）]； p = [ 3 ] =