2014-10-27 105 views
0

我是初学者,使用python。我有以下代码:ValueError:使用序列符号设置数组元素

import Deformation 
import Qtheta 
import numpy as np 
import sympy as sp 

def champs_de_contraintes(epaisseurs,angles,El,Et,Glt,Nult): 
    E0=np.mat('[0;0;0]') #inizialising E0 
    K=np.mat('[0;0;0]') #inizialising K 
    E0K=np.mat('[0;0;0;0;0;0]') #inizialising K 
    E0K = Deformation.defos(epaisseurs,angles,El,Et,Glt,Nult) #function that return a 6x1 matrix type <class 'numpy.matrixlib.defmatrix.matrix'> 
    E0=E0K[:3] #Slicing E0K into 2 vectors E0 and K 
    K=E0K[3:] 
    nb_composantes = 3 #sigma x, sigma y, taux xy 
    Sigmaxy=np.zeros((nb_composantes,len(epaisseurs))) #sigma x, sigma y, taux xy, Array 

    #===============This bloc calculate the altitude of a ply ================= 

    z=[] 
    e=[] 
    z.append(-(np.sum(epaisseurs))/2) 
    for k in (range(len(epaisseurs))): #initialising e[k] 
     e[len(e):] = [0] #Values to add to get the altitude 
    for i in (range(len(epaisseurs))): 
     e[i]= e[i-1] + epaisseurs[i] #sum e 
     z[len(z):] = [z[0] + e[i]] #Altitude 

    #===========================This bloc calculate a vector=================== 
    for i in (range(len(epaisseurs))): 
     newcolumn= Qtheta.Qtheta(angles[i],El,Et,Glt,Nult)*E0+z[i]*Qtheta.Qtheta(angles[i],El,Et,Glt,Nult)*K #z is the altitude 
     #3x1 = 3x3*3x1 + 1x1*3x3*3x1 
     for m in (range(len(newcolumn))): 
      Sigmaxy[i,m]=newcolumn[m] 
    return Sigmaxy 

这将返回我的错误

Sigmaxy[i,m]=newcolumn[m] 
ValueError: setting an array element with a sequence 

基本上,我想要做的是保存矢量“newcolumn”在Sigmaxy矩阵newcolumn。 我认为我得到这个错误是因为“newcolumn”是符号。实际上,E0和K取决于2个变量。

有人可以帮我这个吗?

在此先感谢!

回答

1

如果你想以这种方式设置数组的元素,你需要确保形状是匹配的。

你被抛出,如果你尝试这样的错误:

my_1d_array = np.array([1,2,3]) 
my_1d_array[0] = [4, 5, 6] 

由于错误解释,你想一个单一的元素设置为元素的序列,这是不可能的。调试您的代码,并确保Sigmaxy[i,m]newcolumn[m]的形状(.shape)匹配。

用正确的大小初始化您的数组,并在开始时为所有元素留出空间。如果您不知道确切的大小,可以通过添加列来创建更大的阵列,请参阅其他问题。

相关问题