Merhabalar, ben kod yazmayı ve okumayı bilmiyorum. İnşaat mühendisiyim fakat çalıştığım yerde çözülmesi ve sonuç alınması gereken bir python kodu var. Kod elimde, değerleri girince çalışması lazım fakat şu şekilde bir hata alıyorum.
ValueError: shapes (4,) and (6,6) not aligned: 4 (dim 0) != 6 (dim 0)
Hatanın sebep olduğu satır da şu kısım sanırım:
#
#Strain_Mat=np.zeros((len(elements),4,6))
Strain_Mat={}
#
numnode=4
for el_no,nodes in elements.items():
#coordinates=node_coordinate[nodes-1,1:]
if el_no==230:
print(el_no)
coordinates=[node_coordinate.get(key) for key in nodes]
# Transformation matrix
Tmat=computeLocalCoord(coordinates)
#
R_M=np.zeros((6,6))
R_M[:3,0:3]=Tmat[:,:]; R_M[3:7,3:7]=Tmat[:,:]
#
#
coordinates=np.dot(coordinates,Tmat.T)
#
#
disp_vec=[node_disp.get(key) for key in nodes]
#
disp_vec=np.dot( disp_vec, R_M.T)
#disp_vec = rotation.apply(disp_vec)
disp_vec=np.ravel(disp_vec)
e=flat_shell(4,nodes,coordinates)
p=0
eps_Mat=np.zeros((4,6))
for xi,eta in loc:
wgt=weight[p]
shp, dNdxi, dNdeta=e.calc_shapeFun(xi,eta)
dNdx,dNdy=e.calc_SHP_deriv(dNdxi,dNdeta)
B_m=e.calc_Bmatm(shp,dNdx,dNdy)
B_p=e.calc_Bmatp(shp,dNdx,dNdy)
#
eps_0=np.dot(B_m,disp_vec)
D_eps_p=zcoord*np.dot(B_p,disp_vec)
eps=eps_0+D_eps_p
#
if (abs(abs(Tmat[0,2])-1)<1e-3):
temp=eps[1]
eps[1]=eps[0]
eps[0]=temp
eps_Mat[p][0:3]=eps
#
p=p+1
Strain_Mat[el_no]=eps_Mat
Yardımcı olabilecek birilerini arıyorum.