Question

我是python的初学者。我有以下代码：

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

基本上，我想要做的是保存矢量＆＃34; newcolumn＆＃34;在Sigmaxy矩阵的新列中。我认为我得到了这个错误，因为＆＃34; newcolumn＆＃34;是象征性的。实际上，E0和K取决于2个变量。

有人可以帮我这个吗？

提前致谢！

Answer 1

如果要以这种方式设置数组的元素，则需要确保形状匹配。

如果你尝试这样的话，你会得到错误：

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

正如错误所解释的那样，您正在尝试将单个元素设置为元素序列，这是不可能的。调试代码，确保.shape和Sigmaxy[i,m]的形状（newcolumn[m]）匹配。

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

ValueError：使用序列符号设置数组元素

1 个答案: